《数值分析》计算实习作业。
第二题）算法设计方案：
1、对矩阵a赋值，取计算精度ε=1×10-12；
2、对矩阵a进行拟上三角化，得到a(n-1)，并输出a(n-1)；
对矩阵a的拟上三角化，通过直接调用子函数inftrianglize(a)来实现；拟上三角化得到的矩阵a(n-1)输出至文件中。
3、对a(n-1)进行qr分解并输出q、r及rq矩阵；
qr分解通过直接调用子函数qrdescom(a,q,r, n)实现。
4、运用qr方法求所有的特征值，并输出；
1)初始时令m=n，在m>2的条件下执行；
2)判断如果|amm-1|<ε则得到一个特征值，m=m-1，转(4)；否则转(3)；
3)判断如果|am-1m-2|<ε则得到两个特征值，m=m-2，转(4)；
4)判断如果m≤2，转(6)；否则转(5)；
5)执行相似迭代，转(2)；
6)求出最后的一个或两个特征值；
7)输出全部的特征值至文件中。
5、输出qr分解法迭代结束之后的a(n-1)至文件中；
6、通过反幂法求出所有实特征值的特征向量并输出。
首先令b=(a-λii)，其中λi是实特征值；反幂法通过调用子函数bpowmethod(b,x1)实现，最终λi对应的特征向量就是x1；最后将所有的实特征值的特征向量输出。
北航数值分析第二次大作业
题目 使用带双步位移的qr分解法求矩阵的全部特征值，并对其中的每一个实特征值求相应的特征向量。已知 i,j 1,2,10 一 算法的设计方案 一 总体方案设计 构造矩阵 a，先利用householder矩阵对矩阵a作相似变换，把a化为拟上三角矩阵a n 1 然后进行带双步位移的qr分解求解矩阵的全部...
北航数值分析第二次大作业 QR分解
数值分析a 计算实习题目二。2012年10月。一 算法设计方案。整个程序主要分为四个函数，主函数，拟上三角化函数，qr分解函数以及使用双步位移求解矩阵特征值 特征向量的函数。因为在最后一个函数中也存在qr分解，所以我没有采用参考书上把矩阵m进行的qr分解与矩阵ak的迭代合并的方法，而是在该函数中调用...
数值分析第二次作业
北航2007级研究生。数值分析 计算实习作业。第二题 院系 航空科学与工程学院。学号 sy0705118 姓名 刘德博 日期 2007年12月09日 计算实习说明书。目的 训练运用计算机进行科学与工程计算的能力。要求 1 独立进行算法设计 程序设计和上机运算，并得出正确的结果。2 编制程序时全部采用...