gauss列主元消去法程序设计.docxVIP

Gauss列主元消去法 序设计 《Gauss列主元消去法》实验报告 实验名称：Gauss列主元消去法程序设计 成绩： 专业班级：数学与应用数学1202班 姓名：王晓阳学号： 实验日期：2014年11月10 口 实验报告日期：2014年11月10日 实验目的 学习Gauss消去法的基本思路和迭代步骤. 学会运用matlab编写高斯消去法和列主元消去法程序，求解线性方程组. 当他⑷绝对值较小时，釆用高斯列主元消去法. 培养编程与上机调试能力. 二、 实验内容 用消去法解线性方程组的基本思想是用逐次消去未知数的方法把原线性方程 组化为与其等价的三角形线性方程组，而求解三角形线性方程组可用回代 的方法求解. 求解一般线性方程组的高斯消去法. (1)消元过程: 设他⑷工0,第i个方程减去第k个方程的叫=畋叽叫咅, (i = R +1,…,n),得到 A(k^]x =沪州. 经过n-l次消元，可把方程组A叭⑴化为上三角方程组A叫= 0). (2)回代过程： 以解如下线性方程组为例测试结果. 列主元消去法 由高斯消去法可知，在消元过程中可能出现伦⑷=0的情况，这是消去法将 无法进行，即使主元素他的H0但很小时，用其作除数，会导致其他元素数量级 的严重增长和舍入误差的扩散，最后也使得计算解不可靠.这时就需要选取主元 素，假定线性方程组的系数矩阵A是菲奇异的. (1)消元过程： 对于k = l,2,?.?,n-l,进行如下步骤： 按列选主元，记 交换增广阵A的p,k两行的元素。 A(k, j)=A(p, j)(丿％ …，卅 1) 交换常数项b的p, k两行的元素。 b(k) =b (p) 计算消元 (2)回代过程 ⑶以解如下线性方程组为例测试结果. 三、 实验环境 MATLAB R2014a 四、 实验步骤 高斯列主元消去法流程图： 程序设计: (一)高斯消去法： a=input(请输入系数阵：); b=input(，请输入常数项：)； n二length (b); A=[a, b]; x二zeros (n, 1) ; %初始值 for k=l:n-l for i =k+1:n %第k次消元 m(i, k) =A(i, k)/A(k, k); for j=k+l：n A(i, j)=A(i, j)-A(k, j)*m(i, k); end b(i)=b(i)-m(i, k) *b (k); end end x(n)=b (n)/A(n, n) ; %回代 for i=n~l:~1:1; s=0; for j=i+l:n; end x(i) = (b(i)-s)/A(i, i) end (二)高斯列主元消去法： a二input (请输入系数阵：)； b=input(，请输入常数项：)； n=length (b); A=[a, b]; x=zeros (n, 1) ; %初始值 for k=l:n-l if abs(A(k, k)) 107-4) ;%判断是否选主元 y=l else y=0; end if y; %选主元 for i=k+l:n; if abs (A (i, k)) abs (A (k, k)) p二 i; else p=k; end end if pk; for j二k:n+l; s=A(k,j); A(k, j)=A(p, j);% 交换系数 A(p, j)=s; end t=b (k); b(k)=b(p);%交换常数项 b(p)=t; end end for i二k+l：n m(i,k)=A(i,k)/A(k, k); %第1＜次消元 for j二k+l:n A(i, j)=A(i, j)-A(k, j)*m(i, k); end b(i)=b(i)-m(i, k) *b (k); end end x(n)=b(n)/A(n, n) ; %回代 for i二nT:T:l; s 二0; s=s+A(i, j)*x(j); end x(i) = (b(i)-s)/A(i,i) end 五、实验结果 Gauss1 请输入系数阵:[10, -7,0;-3, 2,6; 5, -1, 5] 请输入常数项：[7;4;6] X 二 0 X二(0, -1, 1) Gauss2 请输入系数阵:[10-(-5),1；2,1] 请输入常数项：[1;2] X= (, 1) 六、实验讨论、结论 本实验通过matlab程序编程实现了高斯消去法及高斯列主元消去法的求 解，能加深对高斯消去法基本思路与计算步骤的理解。当主元素特别小时，需要 选取主元，否则会影响结果，这时就需要采用高斯列主元消去法。 七、参考资料 [1]李庆杨，王能超，易大义.数值分析?清华大学出版社，2008, P142 [2]刘卫国.MATLAB程序设计与应用.高等教育出版社，2006