C语言解决非线性问题.doc

《C语言解决非线性问题.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、#include#include#include#defineN8voidmain(){inti,j;doublex[15][14],y[15][2],sum[20]={0.0},u0=0,b0=0,u1=0,b1=0.0,u2=0,b2=0.0,u3=0,b3=0.0,u4=0,b4=0.0,Q[15],e[15]={0.0},w=0.0,c[15],a,s1=0,s2=0,s3=0,s4=0,d,d0;//FILE*fp;if((fp=fopen("I:\dcsj1.

2、txt","r"))==NULL){printf("Cannotopenfilestrikeanykeyexit!");exit(0);}for(i=0;i<15;i++){for(j=0;j<14;j++){fscanf(fp,"%lf",&x[i][j]);}}fclose(fp);if((fp=fopen("I:\jh.txt","r"))==NULL){printf("Cannotopenfilestrikeanykeyexit!");exit(0);}for(i=0;i<15;i++){for(j=0;j<2;

3、j++){fscanf(fp,"%lf",&y[i][j]);}}sum[0]=15;for(i=0;i

4、um[13]+x[i][2]*x[i][12];sum[14]=sum[14]+x[i][11]*x[i][12];sum[15]=sum[15]+y[i][0];sum[17]=sum[17]+x[i][2]*y[i][0];sum[18]=sum[18]+x[i][11]*y[i][0];sum[19]=sum[19]+x[i][12]*y[i][0];}do{u0=1/sum[0]*(-sum[2]*u2-sum[3]*u3-sum[4]*u4+sum[15]);u2=1/sum[6]*(-sum[2]*u0-sum[1

5、2]*u3-sum[13]*u4+sum[17]);u3=1/sum[7]*(-sum[3]*u0-sum[12]*u2-sum[14]*u4+sum[18]);u4=1/sum[8]*(-sum[4]*u0-sum[13]*u2-sum[14]*u3+sum[19]);w=w+fabs(u0-b0+u2-b2+u3-b3+u4-b4);b0=u0;b2=u2;b3=u3;b4=u4;}while(w>50);for(i=0;i<15-N;i++){Q[i]=b0+b2*x[i][2]+b3*x[i][11]+b4*x[i][

6、12];e[i]=(Q[i]-y[i+N][0])/y[i+N][0];w=w+(Q[i]-y[i+N][0])*(Q[i]-y[i+N][0]);}/**/for(i=0;i<15-N;i++){printf("%lf%lf%lf",Q[i],e[i],w);}if((fp=fopen("I:h.txt","w"))==NULL){printf("Cannotopenfilestrikeanykeyexit!");exit(0);}for(i=0;i<15-N;i++){fprintf(fp,"%lf%lf",

7、c[i+N],y[i+N][0]);}fclose(fp);}