北航数值分析大业一.docVIP

《数值分析》大作业一 1、算法设计方案： ⑴矩阵A的存储与检索： 将矩阵A[501][501],转存为新矩阵A[5][501]。在存入时，按每一行依次存入，存入时原矩阵与新矩阵的列号保持不变，其中原矩阵的主对角线上的元素存为新矩阵的第三行，最后所有元素存入新矩阵中。 ⑵求λ1，λ501，λs λs：λs表示矩阵的按模最小特征值，为求得λs直接对待求矩阵A应用反幂法即可。 λ1、λ501：已知矩阵A的特征值满足关系λ1＜…＜…λ501，要求λ1、及λ501时，可按如下方法求解： 对矩阵A用幂法，求得按模最大的特征值λn1。 按平移量λn1对矩阵A进行原点平移得矩阵，对矩阵B=A+λn1I用反幂法求得B的按模最小特征值λn2。 λn3=λn2一λn1 则：λ1=min（λn1，λn3），λ501=max（λn1，λn3） ⑶求A中与数μk=λ1+k（λ501-λ1）/40最接近的特征值λik（k=1，…39）： 先求矩阵 B=A-I 对应的按模最小特征值，则+即为矩阵A与最接近的特征值。 重复以上过程39次即可求得（k=0，1，…39）的值。 ⑷求A的（谱范数）条件数和行列式： 在（1）中用反幂法求矩阵A的按模最小特征值时，要用到Doolittle分解方法，在Doolittle分解完成后得到的两个矩阵分别为L和U，detA等于U所有对角线上元素的乘积。 ，λmax和λmin分别为模最大