数值分析计算方法实验报告.docVIP

课程 名称 计算方法 实验项目 名 称 线性方程组的数值解法 实验项目类型 验证 演示 综合 设计 其他 指导 教师 胡小兵 成 绩 √ 1．实验目的： （1）高斯列主元消去法求解线性方程组的过程 （2）熟悉用迭代法求解线性方程组的过程 （3）设计出相应的算法，编制相应的函数子程序 2．实验内容 分别用高斯列主元消去法 ,Jacobi迭代法,Gauss--Saidel迭代法,超松弛迭代法求解线性方程组 3．实验环境及实验文件存档名 写出实验环境及实验文件存档名 4.实验结果及分析 输出计算结果，结果分析和小结等。format long; A=[2,10,0,-3;-3,-4,-12,13;1,2,3,-4;4,14,9,-13] b=[10,5,-2,7] [m,n]=size(A); if m~=n error(矩阵A的行数和列数必须相同); return; end if m~=size(b) error(b的大小必须和A的行数或A的列数相同); return; end if rank(A)~=rank([A,b]) error(A矩阵的秩和增广矩阵的秩不相同,方程不存在唯一解); return; end c=n+1; A(:,c)=b; for k=1:n-1 [r,m]=max(abs(A(k:n,k))); m=m+k-1; if(A(m,k)~=0) if(m~=k) A([k m],:)=A([m k],:); end A(k+1:n, k:c)=A(k+1:n, k:c)-(A(k+1:n,k)/ A(k,k))*A(k, k:c); end end x=zeros(length(b),1); x(n)=A(n,c)/A(n,n); for k=n-1:-1:1 x(k)=(A(k,c)-A(k,k+1:n)*x(k+1:n))/A(k,k); end disp(X=); disp(x); format short; 输出结果： 2．Jacobi迭代法： %Jacobi迭代法求解实验1 % A为方程组的增广矩阵 clc; A=[2 10 0 -3 10;-3 -4 -12 13 5;1 2 3 -4 -2;4 14 9 -13 7] MAXTIME=50; eps=1e-5; [n,m]=size(A); x=zeros(n,1); y=zeros(n,1); k=0; disp(迭代过程X的值情况如下:) disp(X=); while 1 disp(x); for i=1:1:n s=0.0; for j=1:1:n if j~=i s=s+A(i,j)*x(j); end y(i)=(A(i,n+1)-s)/A(i,i); end end for i=1:1:n maxeps=max(0,abs(x(i)-y(i))); end if maxeps=eps for i=1:1:n x(i)=y(i); end return; end for i=1:1:n x(i)=y(i); y(i)=0.0; end k=k+1; if kMAXTIME error(超过最大迭代次数,退出); return; end end 输出结果： 由于不收敛，故出现上述情况。 3.Gauss--Saidel迭代法 %Gauss_Seidel迭代法求解实验1 % A为方程组的增广矩阵 clc; format long; A=[2 10 0 -3 10;-3 -4 -12 13 5;1 2 3 -4 -2;4 14 9 -13 7] [n,m]=size(A); Maxtime=50; Eps=10E-5; x=zeros(1,n); disp(x=); for k=1:Maxtime disp(x); for i=1:n s=0.0; for j=1:n if i~=j s=s+A(i,j)*x(j); end