Gauss列主元消去法程序设计

1、精品文档gauss列主元消去法实验报告实验名称：gauss列主元消去法程序设计成绩： 专业班级：数学与应用数学1202班 姓名：王晓阳 学号：2012254010228实验日期：2014年11月10日实验报告日期：2014年11月10日一.实验目的1 .学习gauss消去法的基本思路和迭代步骤.2 .学会运用matlab编写高斯消去法和列主元消去法程序，求解线性方程组.3 .当akg她对值较小时，采用高斯列主元消去法.4 .培养编程与上机调试能力.二、实验内容用消去法解线性方程组的基本思想是用逐次消去未知数的方法把原线性方 程组ax =b化为与其等价的三角形线性方程组，而求解三角形线性方程组可

2、用回 代的方法求解.1.求解一般线性方程组的高斯消去法.(1)消元过程：. 。1欢迎下载精品文档设akkf ),0,第i个方程减去第k个方程的 叫=wk(k)/akk(k第， (i=k+1,|,n),得到 ax=b).w=a/k)-ekaj )(i, j =k +1川,n)h)=nf)-ekbk(k)经过n-1次消元,可把方程组a( x = bc比为上三角方程组an x = b(n).(2)回代过程：xn =bn(n)/ann(n)jnnx = b)-e j xj /a ,i =n-1川,1、lj3)以解如下线性方程组为例测试结果.10x1 -7x2 =7-3x1 2x2 6x3 = 45x1

3、 -x2 5x3 = 62.列主元消去法由高斯消去法可知，在消元过程中可能出现akf)= 0的情况，这是消去法将无法进行，即使主元素akk)r0但很小时，用其作除数，会导致其他元素数量级 的严重增长和舍入误差的扩散，最后也使得计算解不可靠.这时就需要选取主元 素，假定线性方程组的系数矩阵 a是菲奇异的.(1)消元过程：对于k =1,2,|,n-1 ,进行如下步骤：1)按列选主元，记apk =max、k2)交换增广阵a的p,k两行的元素。a(k,j尸a(p,j) (j=k,n+1)3)交换常数项b的p,k两行的元素。b(k尸b(p)4)计算消元卜上冲)=西)-中)(i, j =k+1川,n )

4、ibe=bmkbc(2)回代过程=bn(n)/ann(n) f n 、x = b()_ a xj /aii(i ,i =n-1川,1:j=t+)(3)以解如下线性方程组为例测试结果.0.0000仅 x2 =1j 2x x2 = 2三、实验环境matlab r2014a四、实验步骤1 .高斯列主元消去法流程图：3欢迎下载精品文档2 .程序设计:（一）高斯消去法：a=input（请输入系数阵：）;b=input（请输入常数项：）;n=length（b）;a=a,b；x=zeros（n,1）; % 初始值for k=1:n-1for i=k+1:n%第k次消元m(i,k)=a(i,k)/a(k,k)

5、;for j=k+1:na(i,j)=a(i,j)-a(k,j)*m(i,k);endb(i)=b(i)-m(i,k)*b(k);endendx(n)=b(n)/a(n,n); % 回代for i=n-1:-1:1;s=0;for j=i+1:n;s=s+a(i,j)*x(j);endx(i)=(b(i)-s)/a(i,i)end( 二 ) 高斯列主元消去法：a=input( 请输入系数阵： );b=input( 请输入常数项： );n=length(b);a=a,b;x=zeros(n,1); % 初始值for k=1:n-1if abs(a(k,k)abs(a(k,k)p=i;else p

6、=k;endendif p=k;for j=k:n+1;s=a(k,j);a(k,j)=a(p,j);%交换系数a(p,j)=s;endt=b(k);b(k)=b(p);% 交换常数项b(p)=t;endendfor i=k+1:nm(i,k)=a(i,k)/a(k,k); % 第 k 次消元for j=k+1:na(i,j)=a(i,j)-a(k,j)*m(i,k);endb(i)=b(i)-m(i,k)*b(k);endendx(n)=b(n)/a(n,n); % 回代for i=n-1:-1:1;6 欢迎下载 。精品文档s=0;for j=i+1:n;s=s+a(i,j)*x(j);en

7、dx(i)=(b(i)-s)/a(i,i)end五、实验结果gauss1请输入系数阵：10,-7,0;-3,2,6;5,-1,5请输入常数项：7;4;6x =0-1.00001.0000x =-0.0000-1.00001.0000x=(0， -1 ， 1)gauss2请输入系数阵：10a(-5),1;2,1请输入常数项： 1;2y =1x =0.50001.0000x=(0.5 ， 1)。六、实验讨论、结论本实验通过matlab程序编程实现了高斯消去法及高斯列主元消去法的求解，能加深对高斯消去法基本思路与计算步骤的理解。当主元素特别小时，需要选取主元，否则会影响结果，这时就需要采用高斯列主元消去法。七、参考资料1李庆杨，王能超，易大义.数值分析.