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

--WORD格式 -- 可编辑 《数值分析 B 》大作业一 SY1103120 朱舜杰 一． 算法设计方案： 1.矩阵 A 的存储与检索 将带状线性矩阵 A[501][501] 转存为一个矩阵 MatrixC[5][501] . 由于 C 语言中数组角标都是从 0 开始的，所以在数组 MatrixC[5][501] 中检索 A 的带内元素 a ij 的方法是： ij i-j+2,j A 的带内元素 a =C 中的元素 c 1 501 s 2.求解 λ，λ ，λ ① 首先分别使用幂法和反幂法迭代求出矩阵按摸最大和最小的特征 max min min s 值 λ 和 λ 。λ 即为 λ； max 501 max max max 如果 λ 0, 则 λ = λ ；如果 λ 0, 则 λ1= λ 。 ②使用带原点平移的幂法（ mifa （）函数），令平移量 p= λ ，求出 max ， max 对应的按摸最大的特征值 λ ， ， ， max 1 max max 501 max 如果 λ 0, 则 λ= λ +p ；如果 λ 0, 则 λ = λ +p 。 k 1 501 1 ik 3.求解 A 的与数 μ = λ+k （λ - λ ）/40 的最接近的特征值 λ （k=1,2 ，?， 39 ）。 k k 使用带原点平移的反幂法，令平移量 p= μ，即可求出与 μ最 ik 接近的特征值 λ 。 2 4. 求解 A 的（谱范数）条件数 cond （A ） 和行列式 detA 。 2 1 n 1 n ① cond （A ） =| λ / λ | ，其中 λ和 λ分别是矩阵 A 的模最大和最小特 征值。 -- 专业资料分享 -- --