matlab分析雅克比矩阵,科学网—数值分析---雅克比求特征值matlab程序

%%%程序编写者西北工业大学自动化学院 Email： yincwxa2013@mail.nwpu.edu.cn

%雅克比求特征值

clear

clc

a=input('输入一个n阶对称方阵A=')

n=length(a(1,:))

q=2;

Q=eye(n);

while (q>0.0001)

i=1;

j=2;

for t=1:n-1

for s=i+1:n

if abs(a(t,s))>abs(a(i,j))

i=t

j=s

end

if a(i,i)==a(j,j) cosk=cos(pi/4);sink=sign(a(i,i))*cosk;

else

d=(a(i,i)-a(j,j))/(-2*a(i,j));

t=sign(d)/(abs(d)+sqrt(d^2+1));

cosk=1/sqrt(t^2+1);

sink=t*cosk;

end

for s=1:n

p(s,s)=1;

end

p(i,i)=cosk;

p(j,j)=cosk;

p(i,j)=sink;

p(j,i)=-sink;

a1=p'*a*p

Q=Q*p

q=abs(sum(sum(a.^2))-sum(diag(a.^2)))

a=a1;

end

tezhengzhi=diag(a1)

%% 输入一个n阶对称方阵A=[3.5 -6 5;-6 8.5 -9;5 -9 8.5]

转载本文请联系原作者获取授权，同时请注明本文来自殷春武科学网博客。

链接地址：http://blog.sciencenet.cn/blog-748007-751156.html

上一篇：数值分析----三角分解matlab程序

下一篇：LMIs矩阵内部信息定义法

matlab分析雅克比矩阵,科学网—数值分析---雅克比求特征值matlab程序 - 殷春武的博文...相关推荐

雅克比矩阵（上）-----雅克比推导
1.前言回顾前面几期的内容,在第一期中介绍了机器人的正/逆运动学建模,正运动学解决的问题是如何从关节空间的关节变量描述操作空间的位姿,反之则是逆运动学的内容.将操作空间和关节的空间的关系用以下关系式 ...
matlab调和级数求和,科学网—疯狂的绝技------级数加速收敛的艺术 - 张江敏的博文...
很多时候,我们需要计算一个无穷级数之和.比如,历史上著名的Basel问题是要计算级数之和.这个问题之所以叫巴塞尔问题,是因为来自巴塞尔的约翰-伯努利和雅克比-伯努利为之苦恼了很久,尔后解决之的数学家 ...
matlab比Excel好在哪,科学网—再讨论浮点运算精度问题（IDL与Excel和Matlab相比） - 董彦卿的博文...
论坛里面一个哥们的问题,原帖地址:http://bbs.esrichina-bj.cn/ESRI/viewthread.php?tid=90950&page=1&extra=#pid8 ...
非线性动力学 matlab,科学网—非线性脉冲耦合动力学以及Matlab 源程序 - 王又法的博文...
在这篇文章中,我们给出了非线性脉冲和孤子耦合的理论,它包含了2阶3阶色散,模间色散,损耗增益,色散不匹配,高阶非线性等等.文章主要说明了几个非线性耦合的基本问题: 非线性脉冲耦合行为主要依赖于参量Lc ...
matlab根据乐谱播放音乐,科学网-[原][Matlab][04] Midi音乐键盘-王楠的博文
以前在网上读到一个用matlab播canon音乐的源码感觉很有意思,但声音还不够好. 于是想做个自己的音乐播放器,有键盘,有对应的钢琴声音,读取乐谱就能播.如下. (1)按键的音高与频率从1到下一个 ...
matlab call lapack,科学网—fortran lapack的dgesvd和matlab SVD结果不一致，请大神指点！ - 周锋的博文...
对于A矩阵,matlab SVD分解结果如下: A = 2.2700 -1.5400 1.1500 -1.9400 0.2800 -1.6700 0.9400 -0.780 ...
matlab批量生成灰度图像_科学网—matlab彩色图像的批处理转换为灰度、二值和主成分图图像 - 金秀良的博文...
这段代码主要用来进行图像的批处理转换为灰度.二值和主成分图图像,希望给大家借鉴.其中 RGB to bw可以直接实现,但是效果不好,所以先用RGB to gray,之后再gray to bw.RGB ...

matlab分析雅克比矩阵,科学网—数值分析---雅克比求特征值matlab程序 - 殷春武的博文...

matlab分析雅克比矩阵,科学网—数值分析---雅克比求特征值matlab程序 - 殷春武的博文...相关推荐

最新文章

热门文章