北航数值分析大作业一精品.pdfVIP

《数值分析》大作业（1） 一、算法的设计方案 1. 总体算法设计 题中给出的矩阵 A 为上半带宽和下半带宽都等于 2 的 501 阶的带状稀 疏矩阵。故可以将矩阵 A 压缩存储到一个 5 ×501 矩阵 B 中，对 501×501 矩阵的运算就可以减少成对 5 ×501矩阵的运算，减少了计算量并节省了存 储空间。 总体的算法设计方案如下： （1）计算最小和最大特征值 、 ：  501 1 ① 先通过幂法计算出原矩阵 A 的按摸最大特征值m 1 ，则m 1 为1 或 者501 ； ② 将原矩阵 A 作原点平移，平移量为m 1 ，即令 C=A ﹣m 1 I，再由幂 法计算矩阵 C 的按摸最大特征值m 2 ； ③ 计算  ，则 min{ ,}， max{ ,} 。   m 1 m 2 1 m1 501 m1 （2 ）计算s ，因为s 为按摸最小特征值，对原矩阵 A 使用反幂法即 可算得到。    （3 ）计算矩阵 A 的与数 k 1 k 501 1 最接近的特征值 ik 40 （k=1,2,…,39 ），可以通过带原点平移的反幂法求得该值，平移量分别为    k 1 k 501 1 （k=1,2,…,39 ）。 40 m 1 （4 ）计算 A 的（谱范数）条件数cond (A ) ，cond A ，其中 、 2  2 m 1 s s 分别为（1）、（2 ）中得到的原矩阵A 的按摸最大特征值和按模最小特征 值。 （5 ）计算 A 的行列式的值，对A 进行 Doolittle 分解，A=LU ，L 是单 位下三角阵，U 是上三角阵，则det A A LU U ，由于U 为上三角阵， 故 A 的行列式的值为U 的对角线元素的乘积。 2. 各算法的实现 本次作业中使用的向量范数均为向量的 2-范数，算法的具体实现如下： （1）幂法的实现 n ① 取任意非零向量u0 R ； ②  uT u ； k 1 (k 1) k 1 ③ y k 1 uk 1 k 1 ； ④ uk Ay k 1 ； ⑤ k y T1uk (k 1,2, ) 。 k 当  = 1012 时，迭代终止，以当前的 作为的近似值。   