機(jī)械優(yōu)化設(shè)計課程設(shè)計-閱讀頁
2025-01-27 05:45本頁面
【正文】 for(i=0。i++) x[i]=(double *)malloc(n*sizeof(double))。 for(i=0。i++) *(x[0]+i)=x0[i]。 for(i=0。i++) *(x[1]+i)=*(x[0]+i)+h*s[i]。 if(f2=f1) {h=h0。in。 f3=f1。in。 *(x[1]+i)=*(x[2]+i)。 f2=f3。) {h=2*h。in。 f3=objf(x[2])。 else { for(i=0。i++) {*(x[0]+i)=*(x[1]+i)。 } f1=f2。 } } if(h0) for(i=0。i++) {a[i]=*(x[2]+i)。 } else for(i=0。i++) {a[i]=*(x[0]+i)。 } for(i=0。i++) free(x[i])。 double f1,f2,*x[2],ff,q,w。i2。 for(i=0。i++) {*(x[0]+i)=a[i]+*(b[i]a[i])。 } f1=objf(x[0])。 do {if(f1f2) {for(i=0。i++) {b[i]=*(x[0]+i)。 } f1=f2。in。 f2=objf(x[1])。in。 *(x[1]+i)=*(x[0]+i)。 for(i=0。i++) *(x[0]+i)=a[i]+*(b[i]a[i])。 } q=0。in。 w=sqrt(q)。 for(i=0。i++) xx[i]=*(a[i]+b[i])。 for(i=0。i++) free(x[i])。 } double oneoptim(double x0[],double s[],double h0,double epsg,int n,double x[]) {double *a,*b,ff。 b=(double *)malloc(n*sizeof(double))。 ff=gold(a,b,epsg,n,x)。 free(b)。 } double powell(double p[],double h0,double eps,double epsg,int n,double x[]) {int i,j,m。 double f,f0,f1,f2,f3,fx,dlt,df,sdx,q,d。 s=(double *)malloc(n*sizeof(double))。in。j=n。 *(ss+i*(n+1)+i)=1。i4。 for(i=0。i++) *(xx[0]+i)=p[i]。) {for(i=0。i++) {*(xx[1]+i)=*(xx[0]+i)。 } f0=f1=objf(x)。 for(j=0。j++) {for(i=0。i++) {*(xx[0]+i)=x[i]。 } f=oneoptim(xx[0],s,h0,epsg,n,x)。 if(dfdlt) {dlt=df。 } } sdx=0。in。 if(sdxeps) {free(ss)。 for(i=0。i++) free(xx[i])。 } for(i=0。i++) *(xx[2]+i)=x[i]。 for(i=0。i++) {*(xx[3]+i)=2*(*(xx[2]+i)(*(xx[1]+i)))。 } fx=objf(x)。 q=(f12*f2+f3)*(f1f2dlt)*(f1f2dlt)。 if((f3f1)||(qd)) {if(f2=f3) for(i=0。i++) *(xx[0]+i)=*(xx[2]+i)。in。 } else {for(i=0。i++) {*(ss+(i+1)*(n+1))=x[i](*(xx[1]+i))。 } f=oneoptim(xx[0],s,h0,epsg,n,x)。in。 for(j=m+1。j++) for(i=0。i++) *(ss+i*(n+1)+j1)=*(ss+i*(n+1)+j)。 double ff,x[2]。 printf(x[0]=%f,x[1]=%f,ff=%f\n,x[0],x[1],ff)。 getchar()。相信在以后的工作當(dāng)中，會有所
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1