简单迭代法求解线性方程组.doc

《简单迭代法求解线性方程组.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、简单迭代法求解线性方程组1．原理：将原线性方程组Ax=b中系数矩阵的主对角线移到一边并将其系数化为一，然后通过迭代的方法求解线性方程组的值。(1)计算迭代矩阵：将系数矩阵的所有值分别处以各自所在行的主对角线值，然后将主对角线赋值为0。(2)输入迭代初值，进行迭代将迭代初值存入y[n]矩阵，然后利用迭代式nn=nn+x[i][j]*y[j];y[i]=nn+b[i];经过有限次迭代得到误差要求以内的值2．源程序如下：#include#include#includeusingnamespacestd;#definekk50//定义最大方

2、程元数intn,i,c,j,hh,gg,mm;doubleA[kk][kk],x[kk][kk],b[kk],y[kk],a[kk],z[kk],m,nn,d,e=1,w,fff;voidmain(){cout<<"输入的方程元数"<>n;cout<<"请输入方程系数矩阵:"<>A[i][j];cout<<"请输入右边向量:"<>b[i];cout<<"输入你想要的迭代精度(建议1e-5以上)！"<

3、>fff;cout<<"输入最大迭代次数（建议300次以上）！"<>mm;//计算出迭代矩阵for(i=0;i

4、dl;for(i=0;i>y[i];intf=1;//简单迭代法cout<<"";for(i=1;ifff){for(i=0;i

5、}e=fabs(z[0]-y[0]);if(fabs(z[i]-y[i])>e)e=fabs(z[i]-y[i]);if(i==0){cout<mm){cout<<"迭代次数大于"<

6、不收敛"<