简单迭代法求解线性方程组_第1页
简单迭代法求解线性方程组_第2页
简单迭代法求解线性方程组_第3页
简单迭代法求解线性方程组_第4页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

简单迭代法求解线性方程组1原理：将原线性方程组Ax=b中系数矩阵的主对角线移到一边并将其系数化为一，然后通过迭代的方法求解线性方程组的值。 (1)计算迭代矩阵：将系数矩阵的所有值分别处以各自所在行的主对角线值，然后将主对角线赋值为0。(2)输入迭代初值，进行迭代将迭代初值存入yn矩阵，然后利用迭代式nn=nn+xij*yj; yi=nn+bi;经过有限次迭代得到误差要求以内的值2源程序如下：#include #include #include using namespace std; #define kk 50 /定义最大方程元数 int n,i,c,j,hh,gg,mm; double Akkkk,xkkkk,bkk,ykk,akk,zkk,m,nn,d,e=1,w,fff ; void main() cout输入的方程元数n; cout请输入方程系数矩阵:endl; for(i=0;in;i+)for(j=0;jAij; cout请输入右边向量:endl; for(i=0;ibi; cout输入你想要的迭代精度(建议1e-5以上)！fff; cout输入最大迭代次数（建议300次以上）！mm; /计算出迭代矩阵for(i=0;in;i+) bi=bi/Aii; for(j=0;jn;j+) if(i=j) xii=0; else xij=-Aij/Aii; /输出迭代矩阵 cout计算出迭代矩阵为：endl; for(i=0;in;i+) for(j=0;jn;j+) coutxij ; coutbi ; coutendl; /赋迭代初值cout输入迭代初值endl; for(i=0;iyi; int f=1; /简单迭代法cout ; for(i=1;in+1;i+) coutXi ; cout精度; coutendl; cout迭代初值为： ; coutsetiosflags(ios:fixed); for(i=0;in;i+) coutyi ; coutfff) for(i=0;in;i+) zi=yi; nn=0; for(j=0;je) e=fabs(zi-yi); if(i=0) coutsetiosflags(ios:fixed); cout第setw(3)setprecision(3)f+次迭代 ; coutsetiosflags(ios:fixed); coutsetw(8)setprecision(8)yi ; coute; coutmm) cout迭代次数大于mm次endl; cout认为方程发散，迭代不收敛endl; exit(1); coutendl; coutendl; cout方程迭代了f-1次,达到你所要求的精度fffendl; cout最后结果为：endl; coutendl; for(i=0;in;i+) coutXi
人人文库> 全部分类> 行业资料 > 管理策划

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论