下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、数值分析实验报告专业 信息与计算科学 班级 信计101 姓名 学号 协作队员 实验日期 2013 年 1 月 5 日 星期 六 成绩评定 教师签名 批改日期 题目一、 问题提出给定矩阵,(1)用Matlab函数“eig”求矩阵全部特征值。(2)用幂法求A的主特征值及对应的特征向量。(3)用基本QR算法求全部特征值(可用Matlab函数“qr”实现矩阵的QR分解)。二、 模型建立用幂法求A的主特征值及对应的特征向量的模型:选取v0=(1,1,1,1,1)T,按照下列公式构造向量序列ukvkv0=u00vk=Auk-1k=maxvkuk=vk/k则有limkuk=x1maxx1 limkk=1循环
2、足够多次后,可以近似得出,1=k x1=uk三、 求解方法(1) 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;a=eig(A)(2)pmethod.mfunction l,v,s=pmethod(A,x0,eps)if nargin=2 eps = 1.0e-6;endv = x0; %v为主特征向量M = 5000; %迭代步数限制m = 0; l = 0;for(k=1:M) y = A*v; m = max(y); %m为按模最大的分量 v = y/m; if(abs(m - l)eps) l = m; %到所需精度,退出,l为
3、主特征值 s = k; %s为迭代步数 return; else if(k=M) disp(迭代步数太多,收敛速度太慢!); l = m; s = M; else l = m; end endend(3)function l = rqrtz(A,M)%QR算法求矩阵全部特征值%已知矩阵:A%迭代步数:M%求得的矩阵特征值:lA = hess(A);for i=1:MN = size(A);n = N(1,1);u = A(n,n);q,r=qr(A-u*eye(n,n);A = r*q+u*eye(n,n);l = diag(A);end四、 输出结果(1) a = 13.1724 6.5519 1.5957 -0.3908 -0.9291(2)主特征值x = 13.1724对应的特征向量y = 0.7250 1.0000 0.7930 0.3533 0.0268(3)所有特征值x = 13.1724 6.551
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 记忆的价值阅读理解
- 馆舍租赁合同4篇
- 装修合同书完整版六篇
- 销售合同汇编九篇
- 施工员:设备安装工程知识学习(三)
- 危险化学品储存管理的注意事项
- 药店周年庆促销活动方案
- 2024年新高考高中新教材文言句子翻译技巧(附答案解析)
- 物业项目规划设计方案(2篇)
- 中班舞蹈教学计划方案(2篇)
- 2024年05月内蒙古包头师范学院招考聘用16人笔试历年高频备考试题库含答案带难、易错考点黑钻版解析
- 2024年注册安全工程师考试题库附完整答案【名师系列】
- 2024年山东高速集团限公司上半年校园招聘(1168人)公开引进高层次人才和急需紧缺人才笔试参考题库(共500题)答案详解版
- MOOC 学术英语写作与演讲-上海交通大学 中国大学慕课答案
- 2022版义务教育信息科技新课程标准试题(含答案)
- 中国旅游地理智慧树知到期末考试答案章节答案2024年平凉职业技术学院
- 中华民族发展史智慧树知到期末考试答案2024年
- 江苏省无锡市锡山区2024年中考化学全真模拟试卷含解析
- 硕博研究生英语综合教程完整版电子课件
- 9.1 科学立法 课件-高中政治统编版必修三政治与法治
- 乌拉俄语АБВ智慧树知到期末考试答案2024年
评论
0/150
提交评论