消元法实验报告3.docx

《消元法实验报告3.docx》由会员分享，可在线阅读，更多相关《消元法实验报告3.docx（5页珍藏版）》请在第一文库网上搜索。

1、西京学院数学软件实验任务书课程名称数学软件实验班级*学号*姓名*实验课题线性方程组高斯消去法，高斯列主元消去法，高斯全主元消去法实验目的熟悉线性代数方程组高斯消去法，高斯列主元消去法，高斯全主元消去法实验要求运用Mat1ab/C/C+/Java/Map1e/Mathematica等其中-一种语言完成实验内容线性方程组高斯消去法线性方程组高斯列主元消去法线性方程组高斯全主元消去法成绩教师实验一全主元消去法【实验内容】1 .掌握全主元消去法的基本思路和迭代步骤2 .并能够利用全主元的高斯消去法解任意阶数的线性方程组；3、从课后题中选一题进行验证，得出正确结果，交回实验报告与计算结果。【实验方法与步

2、骤】全主元消去法基本思路(1)消元过程。对于k=1,2,.,nT,进行下列运算：选主元，确定r,t使得若art=O,则系数矩阵为奇异的，停止计算；否则进行下一步。交换A中的r、t两行及t、k两列，并记下交换的号码t、ko对于i=k1,k+2,n；j=k+1,k2,n+1,计算(2)回代过程。对于k=n,n-1,1,计算(3)恢复解的顺序【实验目的】能熟练运用全消元来解线性方程组%【实验程序】functionEx,XA=GuassQuan(A,b)N=size(A);n=N(1);index二O;index_r=0;order=1:n;fori=1:(n-1)me=max(abs(A(1:n,i

3、);fork=i:nforr=i:nif(abs(A(k,r)=me)index_1=k;index_r=r;k=n;break;endendendtemp=A(i,1:n);A(i,1:n)=A(index，1:n);A(index_1,1:n)=temp;bb=b(index_1);b(index)=b(i);b(i)=bb;temp=A(1:n,i);A(1:n,i)=A(1:n,index_r);A(1:n,index_r)=temp;pos=order(i);order(i)=order(index_r);order(index_r)=pos;forj=(1+i):nif(A(i,i)=0)disp(对角元素为0!);return;endI=A(j,i);m=A(i,i);A(j,kn)=A(j,kn)-1*A(i,1:n)/m;b(j)=b(j)-1*b(i)m;endendX=So1veUpTriang1e(A,b);y=zeros(n,1);fori=1:nforj=1:nif(order(j)=i)y(i)=x(j);endendend=y;XA=A;end例：A=138;-529;014运行结果：x二0.3158-1.31580.5789XA二-5.00002.00009.000003.40009.8000001.1176