用Jcobi旋转矩阵法法求矩阵A=[2

matlab编程 用乘幂法近似计算矩阵A的主特征值和特征向量,其中A=[2 -1 0;-1 2 0_百度知道
matlab编程 用乘幂法近似计算矩阵A的主特征值和特征向量,其中A=[2 -1 0;-1 2 0
N=20;0 -1 2],急!,要找出向量x中按摸最大分量的位置!,x,其中A=[2 -1 0,误差容限eps=1e-5主程序好像是format longA=[……]x0=ones(3,怎么破;-1 2 0;eps=1e-5?大神们,eps)但里面好像还少了个函数,最大迭代次数为20,N,k]=PowMeth(A,x0;[lemta,1)matlab编程用乘幂法近似计算矩阵A的主特征值和特征向量
提问者采纳
可见基本QR法的迭代精度还是很高的;-2 0 5.1526
0;).3468作50次迭代计算(具体迭代次数可依具体实验矩阵进行)&gt.9877
disp('
A=U&#39.0000由以上结果可得到迭代计算的特征值为.2000对H矩阵作QR分解;
amax=0,自己改构造矩阵&
disp(&#39.4864
0;%最大迭代次数为100m=100;2.5929
0;=Eps & k&lt.9600
U(p;%迭代误差Eps=1E-5.2136
-0,p)-A(q;6 -3 6]A =
6将矩阵A变换为相似的拟上三角矩阵(即为上Hessenberg矩阵)&gt.0000
4,-1,q)=-s,0;H=hess(A)H =
U=eye(n).0000
3;sqrt(2*(1+sqrt(1-w*w))):n
D(i)=A(i:&#39:&);for i=1;r=1,j))&gt.2127R =
U(p;6 -3 6]为例不改了.0000
for i=1;R*Qans =
5;endif k&gt.0000
2;-1;&&A=[1 -1 2;%计算当前迭代误差
%以下为构造正交矩阵U
l=-A(p;%小于迭代次数或迭代误差进入计算while r&gt,q)=c;
amax=abs(A(i.4031
-5,R]=qr(B),q),p)=s:n
for j=1,i):50
B=R*Q.9651
p=i,j));&gt,2,p)=c.8000
s=-w&#47.0332
u=(A(p,R]=qr(H)Q =
if i~=j & abs(A(i;&
DendQR方法以A=[1 -1 2;*A*U%显示每步计算A的计算结果
k=k+1;sqrt(l*l+u*u).4000
4雅克宾方法clc,2]%取矩阵A的维数n=max(size(A)).3162
6;&A矩阵不收敛'k=1,-1;A特征值为;
c=sqrt(1-s*s);
w=sign(u)*l/
q=1;%矩阵AA=[2 ;0,-1;clear all
不管对不对,就拼你打那么长的东西,就选你了*^_^*
提问者评价
太给力了,你的回答完美地解决了我的问题,非常感谢!
其他类似问题
为您推荐:
乘幂的相关知识
等待您来回答
下载知道APP
随时随地咨询
出门在外也不愁【图文】线性代数第二章矩阵及其运算2-3_百度文库
两大类热门资源免费畅读
续费一年阅读会员,立省24元!
评价文档:
线性代数第二章矩阵及其运算2-3
上传于||暂无简介
大小:3.30MB
登录百度文库,专享文档复制特权,财富值每天免费拿!
你可能喜欢算法第2章_图文_百度文库
两大类热门资源免费畅读
续费一年阅读会员,立省24元!
评价文档:
上传于||暂无简介
大小:946.50KB
登录百度文库,专享文档复制特权,财富值每天免费拿!
你可能喜欢求教刘老师设矩阵a=2_百度知道
求教刘老师设矩阵a=2
提问者采纳
☆⌒_⌒☆ 如果问题解决后,只要及时采纳就是对我们最好的回报,您不用添加任何财富,谢谢。若提问人还有任何不懂的地方可随时追问,祝您学业进步,我会尽量解答,补充多些吧 很高兴能回答您的提问题目不完整
来自团队:
其他类似问题
为您推荐:
等待您来回答
下载知道APP
随时随地咨询
出门在外也不愁用分块矩阵的方法求A=[1 -2 6 -3;-2 3 -2 1;0 0 -7 5;0 0 3 -2]的逆A为4x4的矩阵
1 -2 6 -3-2 3 -2 10 0 -7 50 0 3 -2--分块方法:2,2行,2,2列= B C0 D其逆矩阵 =B^-1 -B^-1CD^-10 D^-1B^-1 =-3 -2 -2 -1 D^-1=2 5 3 7 B^-1*C*D^-1=-7 -21 -5 -15 代入得 A^-1=-3 -2 7 21 -2 -1 5 15 0 0 2 5 0 0 3 7
为您推荐:
其他类似问题
扫描下载二维码

我要回帖

更多关于 旋转矩阵法 的文章

 

随机推荐