实验1 线性方程组的数值解法.doc

重 庆 大 学 学 生 实 验 报 告 实验课程名称 计算方法 开课实验室 数统学院实验室 学 院 机械工程 年级 2011 专业班 机自一班 学 生 姓 名 林凯 学 号 开 课 时 间 2012 至 2013 学年第 二 学期 总 成 绩 教师签名 课程 名称 计算方法 实验项目 名 称 线性方程组的数值解法 实验项目类型 验证 演示 综合 设计 其他 指导 教师 胡小兵 成 绩 √ 1．实验目的： （1）高斯列主元消去法求解线性方程组的过程 （2）熟悉用迭代法求解线性方程组的过程 （3）设计出相应的算法，编制相应的函数子程序 2．实验内容 分别用高斯列主元消去法 ,Jacobi迭代法,Gauss--Saidel迭代法,超松弛迭代法求解线性方程组 3． 计算实例、数据。 1）%Gauss列主元 A=[2 10 0 -3;-3 -4 -12 13;1 2 3 -4;4 14 9 -13]; B=[10;5;-2;7]; A1=A; [n,n]=size(A); A(:,n+1)=B; for k=1:n [l,m]=max(abs(A(k:n,k))); if m~=1 temp=A(k,:); A(k,:)=A(k+m-1,:); A(k+m-1,:)=temp; end for i=1:n if i~=k A(i,:)=A(i,:)-A(k,:)*A(i,k)/A(k,k); end end end for i=n:(-1):1 A(i,:)=A(i,:)/A(i,i); end disp(A(:,n+1)); 2）%Jacobi 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); x=zeros(n,1); y=zeros(n,1); eps=1e-5; k=0; 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 k500 error(e wrong); return; end end 3）Gauss—Saidel 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); x=zeros(n,1); y=zeros(n,1); z=zeros(n,1); eps=1e-5; k=0; 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)*z(j); end y(i)=(A(i,n+1)-s)/A(i,i); z(i)=y(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