实验2_求解线性方程组直接法.pdfVIP

数值分析实验报告二 求解线性方程组的直接方法（ 2 学时） 班级专业 统计一班 姓名 何斌 学号 06 日期 2011-3-30 一 实验目的 1．掌握求解线性方程组的高斯消元法及列主元素法； 2. 掌握求解线性方程组的克劳特法； 3. 掌握求解线性方程组的平方根法。 二 实验内容 1．用高斯消元法求解方程组（精度要求为 10 6 ）： 3 x1 x2 2 x3 7 x1 2 x2 2 x3 1 2x1 2 x2 4 x3 0 2．用克劳特法求解上述方程组（精度要求为 10 6 ）。 3. 用平方根法求解上述方程组（精度要求为 10 6 ）。 4. 用列主元素法求解方程组（精度要求为 10 6 ）： 3 x x 4 x 3 1 2 3 x1 2 x2 2x 3 2 2 x1 3 x2 2x 3 5 三 实验步骤（算法）与结果 1. 高斯消元法求解： 设所求方程组的增广矩阵为 A0 ，系数矩阵为 A01 ，（择定值） L= （l11,l22,l33 ）， 方程组右边的值用有序数组 B 表示 . clear A0=[3 -1 2 7;-1 2 -2 -1;2 -2 4 0];L=[1,1,1]; A01=[3 -1 2 ;-1 2 -2 ;2 -2 4 ];B=[7 -1 0]; %a=zeros(size(A0)); a=A0; %l=zeros(1,3); l=L; 1 %消元过程的计算公式 u11=a(1,1)/l(1,1);u12=a(1,2)/l(1,1);u13=a(1,3)/l(1,1);z1=a(1,4)/l(1,1); a(2,1)=a(2,1)/u11;a(2,2)=(a(2,2)-a(2,1)*u12)/l(1,2); a(2,3)=(a(2,3)-a(2,1)*u13)/l(1,2);z2=(a(2,4)-a(2,1)*z1)/l(1,2); a(3,1)=a(3,1)/u11;a(3,2)=(a(3,2)-a(3,1)*u12)/a(2,2); a(3,3)=(a(3,3)-a(3,1)*u13-a(3,2)*a(2,3))/l(1,3); z3=(a(3,4)-a(3,1)*z1-a(3,2)*z2)/l(1,3); A1=[u11 u12 u13 ;0 a(2,2) a(2,3) ; 0 0 a(3,3)]; % 得到上三角矩阵 Z=[z1;z2;z3]; %迭代求解 x3=z3/a(3,3); x2=(z2-a(2,3)*x3)/a(2,2); x1=(z1-u13*x3-u12*x2)/u11; x=[x1 x2 x3]; for k=1:3 fprintf( x(%d) 的结果为 %.6f\n\n ,k,x(k)), end %控制精度 %残差法检验答案 C=A01*X-B 程序运行结果： x(1) 的结果为 3.500000 x(2) 的结果为 -1.000000 x(3) 的结果为 -2.250000 残差法检验结果： C = 1.0e-01