机械优化设计惩罚函数内点法

1、#include#include#definem12doublef(doublex,doubler);voidjintuifa(doubleabmm,intn,doublex0,doubleh,intij,doublea,doubleb,doubler0);voidhongjinfa(intn,doublea,doubleb,doubleflag,doublex,doubler0);voidbaoweifa(intn,doublex0,doubleh,doubleflag,doublea,doubleb,doublex,doubler0);doublefahansu(doublett)doub

2、lety;if(tt0)ty=-tt;elsety=0;returnty;doubleyuanhansu(doublex)doubles;/s=x0*x0+x1*x1;s=x0*x0+x1*x1+x2*x2+x3*x3;returns;doublef(doublex,doubler)doubles,t,t2;/t=1-x0;t=1-x0;t2=2-x1;/s=yuanhansu(x)-r*log(fahansu(t);s=yuanhansu(x)-r*log(fahansu(t)-r*log(fahansu(t2);returns;voidjintuifa(doubleabmm,intn,do

3、ublex0,doubleh,intij,doublea,doubleb,doubler0)inti,j,z;doublex1m,x2m,x3m,f1,f2,f3;doublesm;for(i=0;in;i+)si=abiij;for(i=0;i=f1)(h=(-1)*h;for(i=0;in;i+)x3i=x1i;f3=f1;for(i=0;in;i+)x1i=x2i;f1=f2;for(i=0;in;i+)x2i=x3i;f2=f3;)for(i=0;in;i+)x3i=x2i+(double)h*si;f3=f(x3,r0);while(f3f2)(h=2*h;for(i=0;in;i+

4、)x1i=x2i;f1=f2;for(i=0;in;i+)x2i=x3i;f2=f3;for(i=0;in;i+)x3i=x2i+(double)h*si;f3=f(x3,r0);)for(i=0;in;i+)(if(x1ix3i)(ai=x1i;bi=x3i;)else(ai=x3i;bi=x1i;)voidhongjinfa(intn,doublea,doubleb,doubleflag,doublex,doubler0)(inti;doublex1m,x2m,f1,f2;while(1)(for(i=0;in;i+)x1i=bi;f1=f(x1,r0);for(i=0;in;i+)x2i

5、=ai;f2=f(x2,r0);if(fabs(f2-f1)/f2)=flag)(for(i=0;in;i+)xi=(double)(bi+ai)/2.0;break;)else(for(i=0;in;i+)x1i=bi-(double)0.618*(bi-ai);for(i=0;if(x2,r0)for(i=0;in;i+)ai=x1i;elsefor(i=0;in;i+)bi=x2i;)voidbaoweifa(intn,doublex0,doubleh,doubleflag,doublea,doubleb,doublex,doubler0)(inti,j,k,r,sum,p,e;doub

6、lex1m,x2m,abcm,krm,f0,f1,f2,fnm,cm,q,wo;k=0;doubleabmm;for(i=0;in;i+)(for(j=0;jn;j+)(if(i=j)abij=1;elseabij=0;while(1)(wo=f(x0,r0);for(p=0;pn;p+)(q=f(x0,r0);jintuifa(ab,n,x0,h,p,a,b,r0);hongjinfa(n,a,b,flag,x1,r0);for(e=0;en;e+)(x0e=x1e;fnp=(double)(q-f(x1,r0);k=k+1;for(i=0;in;i+)x2i=(double)2*x1i-x

7、0i;for(i=0;in;i+)ci=(double)x1i-x0i;r=0;for(i=0;in;i+)(if(fn0=f0|(double)f0-(double)2*f1+(double)f2)*(double)f0-(double)f1-(double)fn0)*(double)f0-(double)f1-(double)fn0)=0.5*fn0*(double)(f0-f2)*(double)(f0-f2)(if(f1for(i=0;in;i+)x0i=x1i;elsefor(i=0;in;i+)x0i=x2i;)else(intij;ij=n-1;for(i=0;in;i+)kri

8、=x1i;for(i=0;in;i+)(for(j=0;j=r)abij=abij+1;)abin-1=ci;)jintuifa(ab,n,kr,h,ij,a,b,r0);hongjinfa(n,a,b,flag,x1,r0);for(i=0;in;i+)x0i=x1i;)if(fabs(f(x0,r0)-wo)/f(x0,r0)=flag)break;)for(i=0;in;i+)xi=x0i;)intmain()(intv,i,n;doubleh,flag,x0m,am,bm,xm;doublefom,fxo,c,r0;c=0.5;r0=100.0;fom=100;printf(请输入维数：n);scanf(%d,&n);printf(请输入初始点:);for(i=0;iflag)fomufko;r0=c*r0;for(v=0;vn;v+)x0v=xv;