機械優(yōu)化設(shè)計上機實踐報告(編輯修改稿)
2025-06-10 01:02 本頁面
【文章內(nèi)容簡介】 ]=*(x[2]+i)。 }for(i=0。i3。i++) free(x[i])。}double gold(double a[],double b[],double eps,int n,double xx[]){ int i。 double f1,f2,*x[2],ff,q,w。 for(i=0。i2。i++) x[i]=(double*)malloc (n*sizeof(double))。 for(i=0。in。i++) {*(x[0]+i)=a[i]+*(b[i]a[i])。 *(x[1]+i)=a[i]+*(b[i]a[i])。 } f1=objf(x[0])。 f2=objf(x[1])。 do {if(f1f2) {for(i=0。in。i++) {b[i]=*(x[0]+i)。 *(x[0]+i)=*(x[1]+i)。 } f1=f2。 for(i=0。in。i++) *(x[1]+i)=a[i]+*(b[i]a[i])。 f2=objf(x[1])。 } else {for(i=0。in。i++) {a[i]=*(x[1]+i)。 *(x[1]+i)=*(x[0]+i)。 } f2=f1。 for(i=0。in。i++) *(x[0]+i)=a[i]+*(b[i]a[i])。 f1=objf(x[0])。 } q=0。 for(i=0。in。i++) q=q+(b[i]a[i])*(b[i]a[i])。 w=sqrt(q)。 }while(weps)。 for(i=0。in。i++) xx[i]=*(a[i]+b[i])。 ff=objf(xx)。 for(i=0。i2。i++) free(x[i])。 return(ff)。}double oneoptim(double x0[],double s[],double 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)。elsefor (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[]={1,1}。double ff,x[2],x1,x2,f。ff=powell(p,2,x)。printf(shuchuzuiyoujie：\n)。x1=x[1]。x2=x[2]。f=ff。printf(x1=%f,x2=%f,f=%f\n,x1,x2,f)。getchar()。}（2）運行結(jié)果為：3《約束優(yōu)化搜索方法》上機實踐報告寫出所選擇的約束優(yōu)化搜索算法的基本過程、原理（可附流程圖說明）。程序計算結(jié)果分析：中間各步驟的結(jié)果分析及與理論計算結(jié)果分析對比。算例1:minf=（x(1)2)^2+(x(2)1)^2。 g1(x)=x(1)^2x(2)=0 g2(x)=x(1)+x(2)2 =0 初始點：x0=[3。3],f(x0)=5 最優(yōu)解：x*=[1。1] f(x*)=1 (1)MATLAB程序清單：function f=myfun_constrain(x)f=（x(1)2)^2+(x(2)1)^2。function [c,ceq]=mycon(x)c=[x(1)^2x(2)。x(1)+x(2)2 ]ceq=[] [x,fval]=fmincon(@myfun_constrain,x0,A,b,[],[],[],[],@mycon)（2） 運行結(jié)果： A=[1,0。0,1]b=[0。0]x0=[3。3]A = 1 0 0 1b = 0 0x
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1