解线性方程组--数值分析.docVIP

下载本文档

0
0
约4.01千字
约 5页
2017-01-12 发布于江苏
举报
版权申诉

解线性方程组--数值分析.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

解线性方程组--数值分析

# includeiostream # includecmath # define MAXN 4 using namespace std; void M(double e[MAXN]); void Jacobi(double a[MAXN][MAXN], double b[MAXN], double TOL); void Gauss( double a[MAXN][MAXN], double b[MAXN], double TOL); void SOR( double a[MAXN][MAXN], double b[MAXN], double w, double TOL); int main(){ int i,j,choice; double a[MAXN][MAXN], b[MAXN],TOL,w; TOL=pow(10,-6); cout\t**************************************************endl; cout\t\t本程序运用三种方法解简单的线性方程组!endl; cout\t\t 线性方程组的阶数暂定为4阶!endl; cout\t\t 源程序中可以更改方程阶数!endl; cout\t\t 温馨提示：更改MAXN!endl; cout\t**************************************************endlendl; cout请先输入有关该线性方程组的数据:endlendl; cout线性方程组的系数矩阵:endl; for(i=0;iMAXN;i++) for(j=0;jMAXN;j++) cina[i][j]; cout线性方程的右端项:endl; for(i=0;iMAXN;i++) cinb[i]; coutendl; cout以下有三种方法可供选择:endl; cout1--Jacobi 2--Gauss seidl 3--SOR endl; cout请选择要运用的方法:endl; cinchoice; switch(choice){ case 1: Jacobi(a,b,TOL); break; case 2: Gauss(a,b,TOL); break; case 3: cout请输入w的值：; cinw; SOR(a,b,w,TOL); break; default: cout该线性方程组无法运用以上的三种方法求的解!endl; break; } return 0; } void M(double e[MAXN]){ for(int i=1;iMAXN;i++) if(e[i]e[0]) e[0]=e[i]; } void Jacobi(double a[MAXN][MAXN],double b[MAXN], double TOL){ int i,j; //声明一个中间变量 double s=0; //第一步：对线性方程组的系数矩阵进行赋值 //右端项以及初始向量的初始化 double x0[MAXN]={0,0,0,0}; double x[MAXN]; double e[MAXN]; //第二步：利用while循环进行求解 do{ for( i=0;iMAXN;i++){ //中间插入的 for循环该循环用来模仿求和公式 for(j=0;jMAXN j!=i;j++) s+=a[i][j]*x0[j]; //外层循环用来求新的解 x[i]=(b[i]-s)/a[i][i]; //求误差向量 e[i]=fabs(x[i]-x0[i]); } //调用子函数取得该误差向量中的最大分量 M(e); //进行判断将符合条件的解向量输出 if(e[0]TOL){ for( i=0;iMAXN;i++) coutx[i]=x[i]endl; } //如果不满足条件则更新初始向量 for(i=0;iMAXN;i++) x0[i]=x[i]; //清空中介值 s s=0; }while(e[0]TOL); } void Ga