版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、精选优质文档-倾情为你奉上一、实验名称:用QR算法求矩阵的特征值二、实验目的:1、通过实验进一步熟悉掌握求矩阵特征值的QR方法及原理。 2、理解QR方法的计算流程。 3、能够编程实现QR方法。三、实验内容:给定矩阵 , ,采用QR方法计算A和H矩阵的全部特征值。四、实验要求:(1) 根据QR算法原理编写程序求矩阵及矩阵的全部特征值(要求误差10)。(2) 直接用MATLAB的内部函数eig求矩阵及矩阵的全部特征值,并与(1)的结果比较。五、QR方法计算矩阵特征值的程序:专心-专注-专业function namda,time,data_na=qr_tz(A,tol)if nargin=1;tol
2、=1e-5;endwucha=1;time=0;while (wucha>tol)&(time<500)q,r=qr(A);A1=r*q;tz0=diag(A1);tz1=diag(A);wucha=norm(tz0-tz1);A=A1;time=time+1;data_na(time,:)=tz1;endnamda=tz1;disp(特征值为)namdadisp(第一个特征在值)timen1=length(data_na);n2=(1:n1);temp1=n2,data_na;subplot(2,2,1:2)plot(date_na(:,1)title(迭代次数为)gri
3、dsubplot(2,2,3)plot(data-na(:,2)title(第二个特征值)gridsubplot(2,2,4)plot(data-na(:,3)title(第三个特征值)grid六、实验结果:>> A=6,2,1;2,3,1;1,1,1;namda,time,data_na=qr_tz(A,1e-5);特征值为namda = 7.5196 2.5311 0.9494迭代次数为time = 6图 1>> A=6,2,1;2,3,1;1,1,1;V,D=eig(A,'nobalance'),V = -0.4791 -0.1127 -0.475
4、5 -0.2482 0.1505 -0.2585 0.7125 0.6152 -0.6346D = 0.1052 0 0 0 2.8525 0 0 0 7.0417>> A=2,3,4,5,6;4,4,5,6,7;0,3,6,7,8;0,0,2,8,9;0,0,0,1,0;namda,time,data_na=qr_tz(A,1e-5);特征值为namda = 13.8078 6.8945 1.8729 -0.7164 -0.8589迭代次数为time =22图 2>> A=2,3,4,5,6;4,4,5,6,7;0,3,6,7,8;0,0,2,8,9;0,0,0,1,
5、0;V,D=eig(A,'nobalance'),V = 0.1240 0.5961 0.8703 0.3744 -0.6599 1.0000 1.0000 0.5544 -1.0000 0.5090 0.8331 0.8813 -1.0000 0.8985 -0.4224 0.9037 -0.4478 0.4509 0.9684 -0.2297 0.8908 -0.3626 0.1483 -0.3073 1.0000D = 13.3192 0 0 0 6.5665 0 0 0 1.9939 0 0 0 0 0 0 0 0 0 0 -0.6488 0 0 -0.2297 表1 用两种方法求得矩阵A的全部特征值特征值QR法直接用eig法17.51967.0417 22.53112.8525 30.94940.1052表2 用两种方法求得矩阵H的全部特征值特征值QR法直接用eig法113.807813.319226.89456.5665 31.87291.9939 4-0.7164-0.6488 5-0.8589-0.2297七、实验结果分析:从图1和图2中可以看出在迭代前几次可能会有一些波动,但逐渐趋于平稳,并且收敛速度快,算法稳定。从表1和表2可以看出直接用MATLAB的内部函数eig求矩阵及矩阵的全部特征值与QR方法
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 英语演讲稿题目
- 请给自己设计一份英语寒假一天的时间计划表
- 红楼梦人物介绍大全
- 高二文科班政治哲学周考卷(五)
- 1 例胸壁完全植入式输液港针头滑脱致长春瑞滨药物外渗患者
- 部编版道德与法治二年级上册第6课《班级生活有规则》教学课件
- 江西省南昌市重点中学2021-2022学年高考考前提分物理仿真卷含解析
- 苏教版四年级下册数学分课时练习作业
- 江西省吉安市新干中学2021-2022学年高三六校第一次联考物理试卷含解析
- 新苏教版五年级数学下册教案(二)
- 实体书店与文化产业市场前景分析
- 慢性肾脏病高血压指南解读
- 民族宗教法律知识讲座
- 2024年重庆联合产权交易所集团招聘笔试参考题库含答案解析
- 关于胃肠道健康知识讲座
- 人文关怀与老年护理
- 文明出行拍摄方案
- 市立医院院级患者安全专项行动方案根据患者安全专项行动方案2023-2025年制订
- 小学项目化学习案例:数学《用我的“尺”解决问题》迭代案例
- Python网络爬虫与数据采集答案1
- 检验科培训课件
评论
0/150
提交评论