中北大学数字分析5-8实验报告报告.doc

下载文档 降价啦

6
0
约7.64千字
约 26页
2017-01-14 发布于湖北
举报
版权申诉
保障服务

中北大学数字分析5-8实验报告报告.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

中北大学数字分析5-8实验报告报告

实验类别：数值分析专业：信息与计算科学班级学号：1308024120 姓名：杨燕中北大学理学院实验五线性方程组的迭代解法【】Jacobi迭代法，Gauss-Seidol迭代法和SOR方法计算下列方程组的解。【实验方法与步骤】用C语言对Jacobi迭代法编程如下： #includestdio.h #includemath.h #define Max_N 10 #define Max_M 100 #define epsilon 1e-6 main() {int i,j,k=0,n; double a[Max_N][Max_N+1],x[Max_N],y[Max_N+1]; double s,temp; printf(\nPlease input n value(dim of Ax=b):); do { scanf(%d,n); if(nMax_N) printf(\nPlease re-input n value:); } while(nMax_N||n=0); /*输入Ax=b的[A|b]矩阵*/ printf(Input the A(i,j):\n); for(i=0;in;i++) for(j=0;jn+1;j++) scanf(%lf,a[i][j]); for(i=0;in;i++) x[i]=0; while(1) { temp=0.0;k++; for(i=0;in;i++) { s=0; for(j=0;jn;j++) { if(j==i) continue; s+=a[i][j]*x[j]; } y[i]=(a[i][n]-s)/a[i][i]; temp+=fabs(y[i]-x[i]); } if(tempepsilon) { printf(Diedaicishuis:%d\n,k); printf(Slove is:); for(i=0;in;i++) printf(x[%d]=%f\n,i,x[i]); break; } if(k=Max_M) { printf(The method is disconvergent!); return; } for(i=0;in;i++) x[i]=y[i]; } } 用C语言对Gauss-Seidol迭代法编程如下： #includestdio.h #includemath.h #define Max_N 10 #define Max_M 100 #define epsilon 1e-6 main() { int i,j,k=0,n; double a[Max_N][Max_N+1],x[Max_N],y[Max_N+1]; double s,temp; printf(\nPlease input n value(dim of Ax=b):); do { scanf(%d,n); if(nMax_N) printf(\nPlease re-input n value:); } while(nMax_N||n=0); /*输入Ax=b的[A|b]矩阵*/ printf(Input the A(i,j):\n); for(i=0;in;i++) for(j=0;jn+1;j++) scanf(%lf,a[i][j]); for(i=0;in;i++) x[i]=0; do { for(i=0;in;i++) { s=0; for(j=0;jn;j++) { s+=a[i][j]*x[j]; } temp=(a[i][n]-s)/a[i][i]; x[i]+=temp; } if(k=Max_M) { printf(The method is disconvergent!); return; } k++; } while(fabs(temp)e