数值分析课程设计报告matlab版,数值分析课程设计(DOC毕业设计论文)

数值分析课程设计(DOC毕业设计论文)

[设计名称] 数值分析课程设计[设计时间] 2009 年 12 月 19 日—2009 年 12 月 21 日[机器型号] 东芝 M302[系统环境] Windows XP[报告要求]1.用 Matlab(或 c)语言或你熟悉的其他算法语言编程序，使之尽量具有通用性。2.完成计算后写出计算实验报告，内容包括：计算机型号和所用机时，算法步骤描述，变量说明，程序清单，输出计算结果，结果分析和小结。3．实习时间： 09 年 12 月 19 日----09 年 12 月 21 日。[实习题目一]给定方程组 1235140 xx[实习要求]1 用 G-S 迭代法解该方程组，取初值 X(0)=(0,0,0)T 迭代到‖x(k+1)-x(k)‖∞≤10-4 为止2 讨论迭代法的收敛性3 输出迭代矩阵 BG 迭代公式,迭代次数及计算结果表[算法描述]1. 运用 LU 分解得到 L 矩阵和 U 矩阵，将方程组转化为 LUX=b，令UX=Y，原方程组转化为 LY=b 和 UX=Y，最后解的原方程组得解 X。2. 用打印函数输出系数矩阵 A，b，L，U，X 和行列式 A 的值。3. 根据矩阵条件数估计的定义，通过程序解得 。x[程序清单]function [x,k]=GassSeidelMathod(A,b)k=1;while k>=1x=GSM(A,b,k);y=x(:,k+1)-x(:,k);f=fanshu(y);if fy(i)max=y(i+1);endendfunction Convergence(A)k=0;[m,n]=size(A);for i=1:msum(i)=0;for j=1:nsum(i)=sum(i)+abs(A(i,j));endsum(i)=sum(i)-abs(A(i,i));if abs(A(i,i))>=sum(i)k=k+1;endendif k==3 G-S 迭代法收敛 else G-S 迭代法不收敛 endfunction hanshu1(A,b,v1)D=diag(v1);L=-tril(A-D);U=-triu(A-D);B=inv(D-L)*Uf=inv(D-L)*b[输出结果][实习题目二] ， 即 AX=b1034248516x[实习要求] 1 用平方根法解上述方程组。2 讨论当 b 有微小变化，即 b+ b=(-4.02, 3.5, 10) 时，A 是否为病态矩阵。给T出讨论结果。3 输出分解矩阵 L’及解向量 X 和 detA。[程序清单]function [L,U]=lu(A)U(1,j)=A(1,j);if j>=2L(j,1)=A(j,1)/U(1,1);elseif j>=ifor k=1: i-1sum=0;sum=sum+L(i,k)*U(k,j);endU(i,j)=A(i,j)-sum;elseif i10 A 为病态矩阵 else A 不是病态的 endfunction FS=fanshu(A)for i=1:3sum(i)=0;for j=1:3sum(i)=sum(i)+A(i,j);endendB=[sum(1) sum(2) sum(3)] ;FS=B(1);for k=1:3if B(k)>FSFS=B(k);kFS=k;endend;function L=average(A)for i=1:3sum=0;for j=1:iif i==jfor k=1:j-1sum=sum+L(i,k)^2;endL(i,j)=sqrt(A(i,j)-sum);elseif j