复合形法实验报告..docx

下载文档

46
0
约 6页
2017-01-20 发布于重庆
举报
版权申诉
保障服务

复合形法实验报告..docx

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

复合形法实验报告.

复合形法实验报告课本参考程序#includemath.h#includestdio.h#includestdlib.h#define E1 0.001#define ep 0.00001#define n 2#define k 4double af;int i,j;double X0[n],XX[n],X[k][n],FF[k];double a[n],b[n];double rm=2657863.0;double F(double C[n]){double F;F=pow(C[0]-3,2)+pow(C[1]-4,2);return F;}int cons(double D[n]){if((D[0]=0)(D[1]=0)(D[0]=6)(D[1]=8)((2.5-D[0]+D[1]=0)((5-D[0]-D[1])=0)))return 1;elsereturn 0;}void bou(){a[0]=0;b[0]=6;a[1]=0;b[1]=8;}double r(){double r1,r2,r3,rr;r1=pow(2,35);r2=pow(2,36);r3=pow(2,37);rm=5*rm;if(rm=r3){rm=rm-r3;}if(rm=r2){rm=rm-r2;}if(rm=r1){rm=rm-r1;}rr=rm/r1;return rr;}void produce(double A[n],double B[n]){int jj;double S;s1:for(i=0;in;i++) { S=r(); XX[i]=A[i]+S*(b[i]-A[i]); } if(cons(XX)==0) {goto s1;} for(i=0;in;i++) { X[0][i]=XX[i]; } for(j=1;jk;j++) { for(i=0;in;i++) { S=r(); X[j][i]=A[i]+S*(B[i]-A[i]); } } for(j=1;jk;j++) { for(i=0;in;i++) { X0[i]=0; for(jj=1;jjj+1;jj++) { X0[i]+=X[jj][i]; } X0[i]=(1/j)*(X0[i]); } if(cons(X0)==0) { goto s1; } for(i=0;in;i++) {XX[i]=X[j][i];} while(cons(XX)==0) { for(i=0;in;i++) { X[j][i]=X0[i]+0.5*(X[j][i]-X0[i]); XX[i]=X[j][i]; } } }}main(){double EE,Xc[n],Xh[n],Xg[n],Xl[n],Xr[n],Xs[n],w;int l,lp,lp1;bou();s111:produce(a,b);s222:for(j=0;jk;j++) { for(i=0;in;i++) { XX[i]=X[j][i]; } FF[j]=F(XX); } for(l=0;lk-1;l++) { for(lp=0;lpk-1;lp++) { lp1=lp+1; if(FF[lp]FF[lp1]) { w=FF[lp];FF[lp]=FF[lp1];FF[lp1]=w; for(i=0;in;i++) { XX[i]=X[lp][i];X[lp][i]=X[lp1][i];X[lp1][i]=XX[i]; } } } } for(i=0;in;i++) { Xh[i]=X[0][i];Xg[i]=X[1][i];Xl[i]=X[k-1][i]; } for(i=0;in;i++) { Xs[i]=0; for(j=0;jk;j++) { Xs[i]+=X[j][i]; } Xs[i]=1/(k+0.0)*Xs[i]; } EE=0; for(j=0;jk;j++) { EE+=pow((FF[j]-F(Xs)),2); } EE=pow((1/(k+0.0)*EE),0.5); if(EE=E1) { goto s333; } for(i=0;in;i++) { Xc[i]=0; for(j=1;jk;j++) { Xc[i]+=X[j][i]; } Xc[i]=1/(k-1.0)*Xc[i]; } if(cons(Xc)==1) { af=1.3;ss:for(i=0;in;i++) { Xr[i]=Xc[i]+af*(Xc[i]-Xh[i]); } if(cons(Xr)==1) { if(F(X