- 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
- DB23_T 3866-2024冰上龙舟赛事组织服务规范.docx
- DB32 2163-2012 棉纱单位可比综合电耗限额及计算方法.docx
- DB32∕T 2429-2013 棉花田间生长发育观察记载规范.docx
- DB23_T 3820-2024 工业互联网综合平台数据质量管理规范.docx
- DB23_T 3474-2023非煤智慧矿山信息系统技术规范.docx
- DB3201_T 1106-2022 工程地质层划分技术规范.docx
- DB33T 2515-2022公共机构“零碳”管理与评价规范.docx
- DB3301_T 0415-2023 装修垃圾收运处置管理规范.docx
- DB3311/T 106―2019“丽水山居”民宿服务要求与评价规范.docx
- DB3201T 1127-2022 慈善捐赠受赠管理规范.docx
最近下载
- 亿赛通电子文档安全管理系统V5.0--客户端使用手册V1113379.pdf VIP
- [标准规范]煤矿窄轨铁道维修质量标准及检查评级办法(精编).doc VIP
- 人教统编版小升初语文总复习专题七:句式变换课件(共28张PPT).ppt VIP
- 2025重庆永川区板桥镇村(社区) 后备干部选拔笔试备考试题及答案解析.docx VIP
- 深度解析《GBT 44026-2024预制舱式锂离子电池储能系统技术规范》.pptx
- 品质部各岗位职责说明书.pdf VIP
- 2022.01 深圳南山红花岭工业区城市更新项目一期建筑方案设计 华艺.pdf VIP
- 电冰箱安全技术规程培训.pptx VIP
- 07FK02 防空地下室 通风设备安装.pdf VIP
- 解决我国农民负担问题的路线图设想-中国发展.pdf
文档评论(0)