机械优化设计实验报告 - 图文 

}

{ a1=a[0]; f1=f[0]; a[0]=a[2]; f[0]=f[2]; a[2]=a1; f[2]=f1; } return;

double lagrange(double *xk,double *ft,double *s) {

int i;

double a[3],f[3]; double b,c,d,aa; finding(a,f,xk,s); for(i=0;;i++) {

if(ia==1)

{ }

aa=a[1]; *ft=f[1]; break;

d=(pow(a[0],2)-pow(a[2],2))*(a[0]-a[1])-(pow(a[0],2)-pow(a[1],2))*(a[0]-a[2]); if(fabs(d)==0) break;

c=((f[0]-f[2])*(a[0]-a[1])-(f[0]-f[1])*(a[0]-a[2]))/d; if(fabs(c)==0) break;

b=((f[0]-f[1])-c*(pow(a[0],2)-pow(a[1],2)))/(a[0]-a[1]); aa=-b/(2*c); *ft=func(xk,aa,s); if(fabs(aa-a[1])<=ad)

{ }

if(*ft>f[1]) { }

aa=a[1]; break;

if(aa>a[1]) {

if(*ft>f[1]) { }

a[2]=aa; f[2]=*ft;

else if(*ft

{

}

a[0]=a[1]; a[1]=aa; f[0]=f[1]; f[1]=*ft;

else if(*ft==f[1]) {

a[2]=aa;a[0]=a[1]; f[2]=*ft;f[0]=f[1]; a[1]=(a[0]+a[2])/2; f[1]=func(xk,a[1],s);

} } else {

if(*ft>f[1])

{ }

a[0]=aa; f[0]=*ft;

else if(*ft

{

a[2]=a[1]; a[1]=aa;

}

}

f[2]=f[1]; f[1]=*ft;

else if(*ft==f[1]) {

a[0]=aa;a[2]=a[1];

f[0]=*ft;f[2]=f[1]; a[1]=(a[0]+a[2])/2; f[1]=func(xk,a[1],s);

} }

if(*ft>f[1]) { } return aa; }

double *gradient(double *xk) {

double *g,f1,f2,q; int i;

g=(double*)malloc(n*sizeof(double)); f1=fny(xk); *ft=f[1];aa=a[1];

}

for(i=0;i

q=ff;

xk[i]=xk[i]+q; f2=fny(xk); g[i]=(f2-f1)/q; xk[i]=xk[i]-q; } return g;

double * bfgs(double *xk) {

double u[n],v[n],h[n][n],dx[n],dg[n],s[n]; double aa,ib;

double *ft,*xk1,*g1,*g2,*xx,*x0=xk; double fi; int i,j,k;

ft=(double *)malloc(sizeof(double)); xk1=(double *)malloc(n*sizeof(double)); for(i=0;i