C++-数值分析-三次样条插值-自动选取步长梯形法-Romberg求积法-列主元高斯消去法-列主元LU分解法-Jacobi迭_第1页
C++-数值分析-三次样条插值-自动选取步长梯形法-Romberg求积法-列主元高斯消去法-列主元LU分解法-Jacobi迭_第2页
C++-数值分析-三次样条插值-自动选取步长梯形法-Romberg求积法-列主元高斯消去法-列主元LU分解法-Jacobi迭_第3页
C++-数值分析-三次样条插值-自动选取步长梯形法-Romberg求积法-列主元高斯消去法-列主元LU分解法-Jacobi迭_第4页
C++-数值分析-三次样条插值-自动选取步长梯形法-Romberg求积法-列主元高斯消去法-列主元LU分解法-Jacobi迭_第5页
已阅读5页，还剩10页未读， 继续免费阅读

下载本文档

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

文档简介

1、 测试矩阵为：2-10-12-10-12测试矩阵为：2-10-12-10-1200-11010源代码：cpp#include#include9松弛法（SOR迭代）取松弛系数w=1.46#defineN4staticstaticstaticstaticstaticdoubledoubledoubledoubledouble/矩阵的维数，可按需更改ANN=2厂1,0,0,-1,2厂1,0,0,-1,2,-1,0,0,-1,2;系数矩阵BN二YN;XN;GN1;（1,0,1,0;右端项输出比较项输出项/X#defineN4staticstaticstaticstaticstaticdoubledou

2、bledoubledoubledouble/矩阵的维数，可按需更改ANN=2厂1,0,0,-1,2厂1,0,0,-1,2,-1,0,0,-1,2;系数矩阵BN二YN;XN;GN1;（1,0,1,0;右端项输出比较项输出项/X二BX+G的G矩阵/计数器inti,j,k;doubleeps=0.001;/最小容许误差intM二100;最大容许比较次数doublew=1.46:/松地因子booldistance0/求两输出项的差的范数是否满足精度要求doubleten?=0;for(i=0;ieps)returnfalse;elsereturntrue;)voidmain()满足精度要求则结束程序/

3、形成迭代矩阵B,存放到A中for(i=0;iN;i+)if(fabs(Aii)eps)(cout打印失败z，endl;wturn;)doubleT=Aii;for(j=0;jN;j+)(Aij=-w*Aij/T;Aii=1-w;Gi=w*Bi/T;intcounter=0:vdiile(counterM)(迭代for(i二0;iN;i+)doubleten)=0;for(j=0;jN;j+)temp=ten?+Aij*Xj;Xi=Gi+ten)；if(distance0二二true)break;else(交换X,Y向量；for(i=0;iN:i+)Yi二Xil;)counter+;)cout迭代次数为:，zcounter，z次。该线性方程组的解(XI,X2,X3Xn)为:zzendl;for(i二0;i*2,X3Xn)为:1

人人文库> 全部分类> 行业资料 > 信息产业

温馨提示

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

最新文档

评论

0/150

提交评论