数值分析实验报告_5.docVIP

应用数学学院12应数《数值分析》课程设计 2014-2015学年第一学期 PAGE 系 别： 应用数学学院 专业名称： 数学与数学应用 指导教师： 陈玉成 姓 名： 汪芹 学 号： 1210052138 二〇一五年一月 课程设计评分表 评 分 内 容 满分 实际得分 1 出勤情况，旷课一节课扣2分，迟到一节课扣1分。 20 2 Word论文层次清晰，按要求的几个方面叙述，语言准确，行文流畅，书写格式规范。 设计能力强，能较好地独立设计出程序，代码清晰，有适当的注释，运行无误。 50 3 答辩时对答流畅，对主要问题叙述清楚准确。 30 总分 100 等级 优（90-100），良（80-89），中（70-79），及格（60-69），不及格（60以下）。 五级制 评语： 指导教师（签名） 年 月 日 PAGE 4 设计题目1 比较Jacobi方法和G-S方法用于方程 基本原理 使用迭代法求解线性方程组的一般做法是先给定解的初始近似值,然后产生一系列列向量.迭代法首先把方程组转化成具有形式的等价方程组.然后,对于给定的初值,由可以得到一系列近似解向量,如果存在,称迭代收敛,否则迭代发散. 程序代码 Jacobi迭代法程序： function x=jacobi(A,b,x0,e) %x0为初始近似值，e为题目要求的误差 D=diag(diag(A));%由A的对角线元素构成对角矩阵D L=-tril(A,-1);%由A的严格下三角元素的相反数构成 %矩阵L U=-triu(A,1); %由A的严格上三角元素的相反数构成矩阵U B=(inv(D))*(L+U);%inv表示矩阵的逆 a=max(abs(eig(B))) f=(inv(D))*b; if a1 x=B*x0+f; n=1; while max(abs(x-x0))=e %相邻两向量差的无穷范数大于e时循环继续 x0=x; x=B*x0+f; n=n+1; end else disp(迭代矩阵谱半径大于1，jacobi迭代不收敛); end G-S迭代法程序： function x=gauss_seidel(A,b,x0,e) %x0为初始近似值，e为题目要求的误差 D=diag(diag(A));%由A的对角线元素构成对角矩阵D L=-tril(A,-1); %由A的严格下三角元素的相反数构成矩阵L U=-triu(A,1); %由A的严格上三角元素的相反数构成矩阵U B=(inv(D-L))*U; a=max(abs(eig(B))) f=(inv(D-L))*b; if a1 x=B*x0+f; n=1; while max(abs(x-x0))=e %相邻两向量差的无穷范数大于e时循环继续 x0=x; x=B*x0+f; n=n+1; end else disp(迭代矩阵谱半径大于1，jacobi迭代不收敛); end 其中，矩阵A,b由以下程序得到： a1=5*ones(1,20); a2=-2*ones(1,19); a3=ones(1,18); A=diag(a1)+diag(a2,1)+diag(a3,2)+… diag(a2,-1)+diag(a3,-2) b=[1;zeros(19,1)] 运算结果 x0=ones(20,1); jacobi(A,b,x0,0.001) a =1.174334860339060 迭代矩阵谱半径大于1，jacobi迭代不收敛 gauss_seidel(A,b,x0,0.001) a =0.450143116282673 x1 x2 x3 x4 x5 0.242 0.094 -0.022 -0.031 -0.007 x6 x7 x8 x9 x10 0.005 0.003 0 -0.001 -0.001 x11 x12 x13 x14 x15 0 0.001 0 0.001 0 x16 x17 x18 x19 x20 0 0 0 0 0 结果分析 根据运行