数值计算A实验报告.docVIP

山东建筑大学 《数值计算A》实验报告 ? ? 二级学院： 理学院 ? 专 业： 信息与计算科学 ? 指导教师： 夏省祥 ? 班级学号： 信计082 2008121273 ? 姓 名： 张建委 ? ? ? ? 实验一 高斯消去法 （1）题目：用选主元素法和高斯消去法求解下列方程组： 已知方程增广矩阵： a = 2 -1 0 0 6 -1 -3 -2 0 1 -1 3 -2 0 0 0 0 -3 5 1 （2）原理： (1)高斯消去法：相对于约当消去法而言，总的来说就是将增广矩阵化为下三角阵。 (2) 顺序选主元素法:相对于高斯消去法的唯一不同是要先按当前要排列元素所在列大小进行排列。 （3）设计思想： (1)高斯消去法：首先让每一行的元素除以该行的主对角线元素。然后利用此行使位于下一行主对角线以前的元素变为0，依次类推。 (2)选主元素法：在高斯消去法的基础上，每次进行化上三角阵之前，重新排列各方程的位置。 （4）对应程序： （a）高斯消去法 function y=gauss1(a,b) %高斯顺序消去法 [m,n]=size(a); if m~=n disp(输入错误，系数矩证阵只能是方阵) end if n~=length(b) disp(输入错误,常数项的个数应与方程的个数相同) end for k=1:n-1 for i=k+1:n a(i,k)=a(i,k)/a(k,k); b(i)=b(i)-a(i,k)*b(k); for j=k+1:n if a(k,k)==0 disp(主元素为零，消去法无法继续) ; break; else a(i,j)=a(i,j)-a(i,k)*a(k,j); end end end end b(n)=b(n)/a(n,n); for i=(n-1):-1:1 w=0; for j=(i+1):n w=w+a(i,j)*b(j); end b(i)=(b(i)-w)/a(i,i); end y=b; （b）高斯列主元消去法 function z=gauss2(a,b,ep) %高斯列主元素消元法 if nargin==2 ep=0.000001 end [m,n]=size(a); if m~=n disp(输入错误，系数矩证阵只能是方阵) end if n~=length(b) disp(输入错误,常数项的个数应与方程的个数相同) end for k=1:n-1 p=a(k,k);I=k; for i=k:n if abs(a(i,k))abs(p) p=a(i,k);I=i; end end if p=ep z=0; end if I~=k for j=k:n w=a(k,j); a(k,j)=a(I,j);a(I,j)=w; end u=b(k);b(k)=b(I); b(I)=u; end for i=k+1:n a(i,k)=a(i,k)/a(k,k); b(i)=b(i)-a(i,k)*b(k); for j=k+1:n a(i,j)=a(i,j)-a(i,k)*a(k,j); end end end b(n)=b(n)/a(n,n); for i=(n-1):-1:1 w=0; for j=(i+1):n w=w+a(i,j)*b(j); end b(i)=(b(i)-w)/a(i,i); end % disp(方程组的解为); z=b; （5）实验结果： a = 2 -1 0 0