- 1、本文档共28页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
第五章矩阵的特征值与特征向量的计算
第五章矩阵的特征值与特征向量的计算
5.2幂法及其MATLAB程序
幂法的MATLAB程序
用幂法计算矩阵的主特征值和对应的特征向量的MATLAB主程序
function[k,lambda,Vk,Wc]=mifa(A,V0,jd,max1)
lambda=0;k=1;Wc=1;,jd=jd*0.1;state=1;V=V0;
while((k=max1)(state==1))
Vk=A*V;[mj]=max(abs(Vk));mk=m;
tzw=abs(lambda-mk);Vk=(1/mk)*Vk;
Txw=norm(V-Vk);Wc=max(Txw,tzw);V=Vk;lambda=mk;state=0;
if(Wcjd)
state=1;
end
k=k+1;Wc=Wc;
end
if(Wc=jd)
disp(请注意:迭代次数k,主特征值的近似值lambda,主特征向量的近似向量Vk,相邻两次迭代的误差Wc如下:)
else
disp(请注意:迭代次数k已经到达最大迭代次数max1,主特征值的迭代值lambda,主特征向量的迭代向量Vk,相邻两次迭代的误差Wc如下:)
end
Vk=V;k=k-1;Wc;
例用幂法计算以下矩阵的主特征值和对应的特征向量的近似向量,精度.并把〔1〕和〔2〕输出的结果与例5.1.1中的结果进行比拟.
(1);(2);(3);〔4〕.
解〔1〕输入MATLAB程序
A=[1-1;24];V0=[1,1];
[k,lambda,Vk,Wc]=mifa(A,V0,0.00001,100),
[V,D]=eig(A),Dzd=max(diag(D)),wuD=abs(Dzd-lambda),wuV=V(:,2)./Vk,
运行后屏幕显示结果
请注意:迭代次数k,主特征值的近似值lambda,主特征向量的近似向量Vk,相邻两次迭代的误差Wc如下:
k=lambda=Wc=
Vk=V=wuV=
Dzd=wuD=
31.738368038406435e-006
由输出结果可看出,迭代33次,相邻两次迭代的误差Wc8.6919e-007,矩阵的主特征值的近似值lambda3.00000和对应的特征向量的近似向量Vk〔-0.50000,1.00000,lambda与例中的最大特征值近似相等,绝对误差约为1.73837e-006,Vk与特征向量的第1个分量的绝对误差约等于0,第2个分量的绝对值相同.由wuV可以看出,的特征向量V(:,2)与Vk的对应分量的比值近似相等.因此,用程序mifa.m计算的结果到达预先给定的精度.
(2)输入MATLAB程序
B=[123;213;336];V0=[1,1,1];
[k,lambda,Vk,Wc]=mifa(B,V0,0.00001,100),[V,D]=eig(B),
Dzd=max(diag(D)),wuD=abs(Dzd-lambda),wuV=V(:,3)./Vk,
运行后屏幕显示结果
请注意:迭代次数k,主特征值的近似值lambda,主特征向量的近似向量Vk,相邻两次迭代的误差Wc如下:
k=lambda=Wc=Dzd=wuD=
39090
Vk=wuV=
V=
(3)输入MATLAB程序
C=[122;1-11;4-121];V0=[1,1,1];
[k,lambda,Vk,Wc]=mifa(C,V0,0.00001,100),[V,D]=eig(C),
Dzd=max(diag(D)),wuD=abs(Dzd-lambda),
Vzd=V(:,1),wuV=V(:,1)./Vk,
运行后屏幕显示
请注意:迭代次数k已经到达最大迭代次数max1,主特征值的迭代值lambda,主特征向量的迭代向量Vk,相邻两次迭代的误差Wc如下:
k=lambda=Wc=
1000.
您可能关注的文档
- 竖脊肌协助稳定骶髂关节(运动康复).doc
- 第11章三角形单元测试卷.doc
- 农村社会问题及其解决路径.ppt
- 入大学后感言.ppt
- 简笔画涂色大全.doc
- 做诚信中学生.ppt
- 糕点厂生产工艺文件.doc
- 第4课文艺复兴与宗教改革教案.doc
- 公共管理的范畴与要素.ppt
- 空间几何中的平行与垂直关系.doc
- 2025-2026年度国内抗酸药及治疗消化性溃疡和胃肠胀气用药市场发展规划及投资前景咨询报告.doc
- 2025-2026年度第一学期园务计划.doc
- 2025-2026年度第一学期学校工作计划.docx
- 2025-2026年度口服轮状病毒活疫苗市场深度分析及产业链投资价值研究咨询报告.doc
- 高考历史一轮复习 世界史 第03讲 走向整体的世界(原卷版).docx
- 高考历史一轮复习 世界现代史单元检测(学生版).docx
- 高考历史一轮复习 世界史 第03讲 走向整体的世界(解析版).docx
- 中考数学总复习第七模块图形的变化练习题整理.docx
- 《幼儿园管理条例》专业解读课件.pptx
- 人教版二年级数学上册第一二单元综合素质达标测试题课件.pptx
文档评论(0)