北航研究生数值分析编程大作业.docx

优质参考文档 数值分析大作业 一、算法设计方案 1、 矩阵初始化 ■ s 15011，在矩这样所需要的存储矩阵A= ? 501 501的下半带宽r=2,上半带宽s=2,设置矩阵C 阵C中检索矩阵A中的带内元素aj的方法是：a ■ s 15011，在矩 这样所需要的存储 单元数大大减少，从而极大提高了运算效率。 2、 利用幕法求出1，501 幕法迭代格式： 非零向量 Uo R n k」 yk」 其中 其中 B 二 A- I (I ，比较 后来存放中间向 Uk 二 Ayk」 当歸-円/Bk兰1012时，迭代终止。 首先对于矩阵A利用幕法迭代求出一个■，然后求出矩阵B， 为单位矩阵)，对矩阵B进行幕法迭代，求出■，之后令 与 ”的大小，大者为-501，小者为-1 。 3、利用反幕法求出-s, ik 反幕法迭代格式： 非零向量Uo ￡ Rn ^k』=Juk』uk v yk」=Uk 4/^4 A? = yk 4 忆=yJk 当\ 1 / \ 10^2时，迭代终止，’s =1 \ o 每迭代一次都要求解一次线性方程组 Au^yk4，求解过程为: 作分解A二LU 对于k =1,2,..., n执行 k 4 ck -j s 1,j ?— ck 4:S 1, j ■— ck -t s 1,t ct -j s 1, j t :max(1,k -r, j -s) 〔j = k, k 1,..., mi