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

《机械优化设计惩罚函数内点法.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、#include#include#definem12doublef(doublex[],doubler);voidjintuifa(doubleab[m][m],intn,doublex0[],doubleh,intij,doublea[],doubleb[],doubler0);voidhongjinfa(intn,doublea[],doubleb[],doubleflag,doublex[],doubler0);voidbaoweifa(intn,doublex0[],doubleh,doubleflag,double

2、a[],doubleb[],doublex[],doubler0);doublefahansu(doublett){doublety;if(tt<0)ty=-tt;elsety=0;returnty;}doubleyuanhansu(doublex[]){doubles;//s=x[0]*x[0]+x[1]*x[1];s=x[0]*x[0]+x[1]*x[1]+x[2]*x[2]+x[3]*x[3];returns;}doublef(doublex[],doubler){doubles,t,t2;//t=1-x[0];t=1-x[0];t2=2-x[1];/

3、/s=yuanhansu(x)-r*log(fahansu(t));s=yuanhansu(x)-r*log(fahansu(t))-r*log(fahansu(t2));returns;}voidjintuifa(doubleab[m][m],intn,doublex0[],doubleh,intij,doublea[],doubleb[],doubler0){inti,j,z;doublex1[m],x2[m],x3[m],f1,f2,f3;doubles[m];for(i=0;i

4、++){x1[i]=x0[i];x2[i]=x0[i]+(double)h*s[i];}f1=f(x1,r0);f2=f(x2,r0);if(f2>=f1){h=(-1)*h;for(i=0;i

5、1[i]=x2[i];f1=f2;for(i=0;i

6、[m],f1,f2;while(1){for(i=0;i

7、[i]);if(f(x1,r0)>f(x2,r0))for(i=0;i

8、++){for(j=0;j