機械優(yōu)化設(shè)計課程設(shè)計-資料下載頁
2025-01-12 05:45本頁面
【正文】 ouble h0,double epsg,int n,double x[]) {double *a,*b,ff。 a=(double *)malloc(n*sizeof(double))。 b=(double *)malloc(n*sizeof(double))。 jtf(x0,h0,s,n,a,b)。 ff=gold(a,b,epsg,n,x)。 free(a)。 free(b)。 return (ff)。 } double powell(double p[],double h0,double eps,double epsg,int n,double x[]) {int i,j,m。 double *xx[4],*ss,*s。 double f,f0,f1,f2,f3,fx,dlt,df,sdx,q,d。 ss=(double *)malloc(n*(n+1)*sizeof(double))。 s=(double *)malloc(n*sizeof(double))。 for(i=0。in。i++) {for(j=0。j=n。j++) *(ss+i*(n+1)+j)=0。 *(ss+i*(n+1)+i)=1。 } for(i=0。i4。i++) xx[i]=(double *)malloc(n*sizeof(double))。 for(i=0。in。i++) *(xx[0]+i)=p[i]。 for(。) {for(i=0。in。i++) {*(xx[1]+i)=*(xx[0]+i)。 x[i]=*(xx[1]+i)。 } f0=f1=objf(x)。 dlt=1。 for(j=0。jn。j++) {for(i=0。in。i++) {*(xx[0]+i)=x[i]。 *(s+i)=*(ss+i*(n+1)+j)。 } f=oneoptim(xx[0],s,h0,epsg,n,x)。 df=f0f。 if(dfdlt) {dlt=df。 m=j。 } } sdx=0。 for(i=0。in。i++) sdx=sdx+fabs(x[i](*(xx[1]+i)))。 if(sdxeps) {free(ss)。 free(s)。 for(i=0。i4。i++) free(xx[i])。 return(f)。 } for(i=0。in。i++) *(xx[2]+i)=x[i]。 f2=f。 for(i=0。in。i++) {*(xx[3]+i)=2*(*(xx[2]+i)(*(xx[1]+i)))。 x[i]=*(xx[3]+i)。 } fx=objf(x)。 f3=fx。 q=(f12*f2+f3)*(f1f2dlt)*(f1f2dlt)。 d=*dlt*(f1f3)*(f1f3)。 if((f3f1)||(qd)) {if(f2=f3) for(i=0。in。i++) *(xx[0]+i)=*(xx[2]+i)。 else for(i=0。in。i++) *(xx[0]+i)=*(xx[3]+i)。 } else {for(i=0。in。i++) {*(ss+(i+1)*(n+1))=x[i](*(xx[1]+i))。 *(s+i)=*(ss+(i+1)*(n+1))。 } f=oneoptim(xx[0],s,h0,epsg,n,x)。 for(i=0。in。i++) *(xx[0]+i)=x[i]。 for(j=m+1。j=n。j++) for(i=0。in。i++) *(ss+i*(n+1)+j1)=*(ss+i*(n+1)+j)。 } } } void main() {double p[]={8,9}。 double ff,x[2]。 ff=powell(p,2,x)。 printf(x[0]=%f,x[1]=%f,ff=%f\n,x[0],x[1],ff)。 printf(學(xué)號： 202212701041 姓名：朱振 )。 getchar()。 } 4）分析結(jié)果 四、課程總結(jié) 優(yōu)化設(shè)計課程的學(xué)習(xí)，初步掌握了一些對機械工程設(shè)計優(yōu)化的方法，能夠應(yīng)用于機械的設(shè)計中。相信在以后的工作當(dāng)中，會有所應(yīng)用。 參考文獻 《機械優(yōu)化設(shè)計》第四版 孫靖民 機械工業(yè)出版社
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1