数值方法源代码 C 插值积分微分迭代.doc

下载文档 降价啦

9
0
约6.32千字
约 10页
2018-06-21 发布于河南
举报
版权申诉
保障服务

数值方法源代码 C 插值积分微分迭代.doc

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

数值方法源代码 C 插值积分微分迭代

高斯算法和追赶法求解线性方程组高斯： #includestdio.h #includestdlib.h #includemath.h void guass(double a[],int n,double b[],double x[]) { int i,j,m; double s; for(i=1,j=0;in;i++,j++) { for(m=i;mn;m++){ s=a[j*n+j]/a[m*n+j]; //printf(%f ,s); //printf(\n); for(int k=j;kn;k++) { a[m*n+k]=a[m*n+k]*s-a[(i-1)*n+k]; } b[m]=b[m]*s-b[i-1]; } } for(i=0;in;i++) for(int j=0;jn;j++){ printf(%f ,a[n*i+j]); if(j==n-1){ printf(%f ,b[i]); printf(\n);} } for(i=n-1;i-1;i--) { if(i==n-1) x[i]=b[i]/a[i*n+i]; else { s=b[i]; for(j=1;jn-i;j++) { s=s-a[i*n+n-j]*x[n-j]; } x[i]=s/a[i*n+i]; } } } void main() { int n=3; double x[3]; double a[3][3]={{2.5,2.3,-5.1}, {5.3,9.6,1.5}, {8.1,1.7,-4.3}, }; double b[4]={3.7,3.8,5.5}; //double x[3]; /*for(int i=0;in;i++) for(int j=0;jn;j++){ printf(%f ,a[i][j]); if(j==n-1) printf(\n); }*/ guass(a[0][0],n,b[0],x[0]); for(int i=0;in;i++) printf(x[%d]=%f ,i,x[i]); } 追赶法： #includestdio.h #includestdlib.h #includemath.h void chase(float a[],int n,float l[],float u[]) { int i; for(i=1;i=n-1;i++) { l[i*n+i-1]=a[i*n+i-1]; if(i==1) l[0]=a[0]; else l[(i-1)*n+i-1]=a[(i-1)*n+i-1]-l[(i-1)*n+i-2]*u[(i-2)*n+i-1]; u[(i-1)*n+i-1]=1; u[(i-1)*n+i]=a[(i-1)*n+i]/l[(i-1)*n+i-1]; } l[(i-1)*n+i-1]=a[(i-1)*n+i-1]-l[(i-1)*n+i-2]*u[(i-2)*n+i-1]; u[(i-1)*n+i-1]=1; } void solve(int n,float l[],float u[],float f[]) { float y[5],x[5]; for(int i=1;i=n;i++) { if(i==1) y[i-1]=f[0]/l[0]; else y[i-1]=(f[i-1]-l[(i-1)*n+i-2]*y[i-2])/l[(i-1)*n+i-1]; } for(i=n;i0;i--) { if(i==n) x[i-1]=y[i-1]; else x[i-1]=y[i-1]-x[i]*u[(i-1)*n+i]; } printf(\n); for(i=0;in;i++) printf(%f ,x[i]); } void main() { int n=5; //double x[4]; float a[5][5]= { {-2,0,0,0,0}, {1,-2,0,0,0}, {0,1,-2,0,0}, {0,0,1,-2,0}, {0,0,0,1,-2} }; float f[]={-10,0,0,0,0}; float l[5][5],u[5][5]; //double b[4]={1,0,1,0}; //double x[3]; /*for(int i=0;in;i++)