数值分析 全主元消去法实验 程序代码

《数值分析 全主元消去法实验 程序代码》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、一、实验目的本次试验研究的对象是n阶线性代数方程组Ax=b，其中A=(aij)是方程组的系数aij构成的n×n阶矩阵，叫做系数矩阵。B=（ai（n+1）），x为所求的解。主元素消去法是为控制舍入误差而提出来的一种算法，在Gauss消去法的消元过程中，若出现akk(k)=0,则消元无法进行，即使akk(k)≠0，但很小，把它作为除数，就会导致其他元素量级的巨大增长和舍入误差的扩散，最后使计算结果不可靠，而全主元消去法正式解决这种问题的算法。二、算法描述抑制舍入误差的增长，通常有两个途径，一是增加参加计算的数字位数，从而使最后结果中积累起来的误差随之减小。但这样做会使计算的时间增加，

2、我们这里要讲的是另一种途径，在做除法运算时，分母的绝对值越小，舍入误差影响就越大，因此在做除法运算时，要选取绝对值比较大的做分母，这就是主元素消去法的基本思想。三，源代码#include"iostream.h"doublejueduizhi(doublex){if(x<0)x=x*(-1);returnx;8}intmain(){inti,j,k,l,n,t=1,s=1,r=1,v,w;doublea[50][51],y,x[50],temp;//cout<<"请输入数组的维数n然后依照下面的格式输入数据";//cout<<"x[]x[]x[]x[]";//cout<

3、<"a11a12....a1nb1a21a22....a2nb2.....an1an2....annbn";cin>>n;for(i=0;i>a[i][j];cout<<"原矩阵"<

4、i(a[j][k]);if(y

5、w]<<"";cout<1;i--){y=((a[i-1][t])/(a[t][t]));for(j=1;j

6、

7、ut<