機(jī)械優(yōu)化設(shè)計(jì)黃金分割法_外推法

1、大學(xué)機(jī)械優(yōu)化設(shè)計(jì)部分程序1. 外推法2. 黃金分割法3. 二次插值法4. 坐標(biāo)輪換法5. 隨機(jī)方向法6. 四桿機(jī)構(gòu)優(yōu)化設(shè)計(jì)1.外推法源程序:2.黃金分割法源程序：#in clude#in clude#defi nef(x)x*x*x*x-5*x*x*x+4*x*x-6*x+60double hj(double *a,double *b,double#i nclude#in clude#defi ne R 0.01double fun( double x) double m;m=x*x-10*x+36;return m;void mai n()double h0=R,y1,y2,y3,x1,x2

2、,x3,h;x1=0;h=h0;x2=h;y1=fu n(x1);y2=fu n(x2);if(y2y1)h=-h;x3=x1;y3=y1;x1=x2;y1=y2;x2=x3;y2=y3;x3=x2+h;y3=fu n( x3);while(y3f(x2)*a=x1;else*b=x2;*n=*n+1;s=hj(a,b,e, n);return s;void mai n()double s,a,b,e,m;int n=0;printf(輸入a, b值和精度e值n);sca nf(%lf %lf %lf,&a,&b,&e);prin tf(a=%lf,b=%lf,s=%lf,m=%lf, n=%

3、dn,a,b,s,m, n);運(yùn)行過(guò)程及結(jié)果：輸入a, b值和精度e值-350.0001a=3.279466,b=3.279793,s=22.659008,m=3.279629, n=213.二次插值法源程序：#in clude#in cludeint main(v oid)doublea1,a2,a3,ap,y1,y2,y3,yp,c1,c2,m;double j3;int i,h=1;void findin g(double a3);s=hj(&a,&b,e,&n);a3=ap;y3=yp;else if(y2=yp)a3=a2;y3=y2;a2=ap;y2=yp; elsea1=ap;y

4、1=yp;double x,y;if(y2=1;i+)C1=(y3-y1)/(a3-a1);C2=(y2-y1)/(a2-a1)-c1)/(a2-a3);ap=0.5*(a1+a3-c1/c2);yp=f(ap);if(fabs(y2-yp)/y2)0)if(y2=yp)a仁 a2;y1=y2;a2=ap;y2=yp;elseh=2*h;int h,i;double y3;a0=0;h=1;a1=h;y0=f(a0);y1=f(a1);if(yiyO)h=-h;a2=a0;y2=y0;doa0=a1;a1=a2; y0=y1;y1=y2;a2=a1+h;y2=f(a2);h=2*h;while

5、(y2vy1);elsefor(i=1;i=1;i+)a2=a1+h;y2=f(a2); if(y2=y1) break;a0=a1;y0=y1;a1=a2;y1=y2;return;運(yùn)行過(guò)程及結(jié)果：a*=5.000000y*=11.0000004.坐標(biāo)輪換法源程序：#in clude #in clude #in clude vconi o.hfloat fun 1(float x,float afloat b)float y;y=x+a*b;return y;float fun 2(float x,float y)float 乙z=4*(x-5)*(x-5)+(y-6)*(y-6);retu

6、rn 乙mai n()floatd1003,x1003,xx3,ax1003;floata1,a2,a3,h,t,y1,y2,y3,e,a,b,l,fi;int i,k;printf(輸入初始點(diǎn)坐標(biāo)n);sea nf(%f%f, &x01, &x02);e=0.000001;1=0.618;x21=x01;x22=x02;k=0;k-;dox01=x21;x02=x22;k+;for(i=1;iyi)h=-h;a3=a1;y3=y1;a1=a2;a2=a3;y1=y2;y2=y3;a3=a2+h;xi1=fu n1(xi-11,di1,a3);xi2=fu n1(xi-12,di2,a3);

7、y3=fu n2(xi1,xi2);doa仁 a2;y1=y2;a2=a3;y2=y3;a3=a2+h;xi1=fu n1(xi-11,di1,a3);xi2=fu n1(xi-12,di2,a3); y3=fu n2(xi1,xi2);while(y3a3;)t=a3;a3=a1;a1=t;y3=y1;y1=t;a=a1;b=a3;a仁 b-l*(b-a);a2=a+l*(b-a);xi1=fu n1(xi-11,di1,a1);xi2=fu n1(xi-12,di2,a1);y1=fu n2(xi1,xi2);xi1=fu n1(xi-11,di1,a2);xi2=fu n1(xi-12,

8、di2,a2);y2=fu n2(xi1,xi2);if(be;)if(y1=y2)a=a1;a仁 a2;y仁 y2;a2=a+l*(b-a);xi1=fu n1(xi-11,di1,a2);t=y1;yi=y2；xi2=fu n1(xi-12,di2,a2);y2=fu n2(xi1,xi2);elseb=a2;a2=a1;y2=y1;a仁 b-l*(b-a);xi1=fu n1(xi-11,di1,a1);xi2=fu n1(xi-12,di2,a1);y1=fu n2(xi1,xi2);elsefor(;fabs(b-a)/b)=e|fabs(y2-y1)/y2)=e;)if(y1=y2

9、)a=a1;a2=a+l*(b-a);xi1=fu n1(xi-11,di1,a2);xi2=fu n1(xi-12,di2,a2);y2=fu n2(xi1,xi2);elseb=a2;a2=a1;y2=y1;a仁 b-l*(b-a);xi1=fu n1(xi-11,di1,a1);xi2=fu n1(xi-12,di2,a1);y1=fu n2(xi1,xi2);axki=0.5*(a+b);xi1=fu n1(xi-11,di1,axki);a1=a2;xi2=fu n1(xi-12,di2,axki)；5.隨機(jī)方向法源程序：#in clude#in cludewhile(sqrt(po

10、w(x21-x01),2)+pow(x22-x02),2)=1e-6);xx1=x21;xx2=x22;fi=fu n2(xx1,xx2);printf(最優(yōu)解為 nx1*=%fnx2*=%fn f*=%fnk=%dn,xx1,xx2,fi,k);運(yùn)行過(guò)程及結(jié)果：輸入初始點(diǎn)坐標(biāo)89最優(yōu)解為x1*=5.000000x2*=6.000000#in cludefloat f(float x,float y)float 乙z=(x-2)*(x-2)+(y-1)*(y-1);return 乙float g1(float x,float y)float 乙z=x*x-y;return 乙float g2(

11、float x,float y)float 乙 z=x+y-2; return 乙k=2void main()floatk=8,c=0.000001,a0=-3,b0=3,a 1=-3,b 1=3;float x10,x010,xl10,e10,r10,d10, h,fl,f0,fx;While(g1(x00,x01)0|g2(x00,x01)0)x00=a0+(ra nd()/32767.00)*(b0-a0);x01=a1+(ra nd()/32767.00)*(b1-a1);fl=f(x00,x01);f0=f(x00,x01);while(1)h=0.01;e0=r0/sqrt(r0*

12、r0+r1*r1);e1=r1/sqrt(r0*r0+r1*r1);x0=x00+h*e0;x1=x01+h*e1;if(g1(x0,x1)=0&g2(x0,x1)=0) fx=f(x0,x1);if(fxfl)fl=fx;for(i=0;i2;i+)di=ei;xli=xi;while(j=k)j+;r0=-1+(ra nd()/32767.00)*(1-(-1);r0=-1+(ra nd()/32767.00)*(1-(-1);e0=r0/sqrt(r0*r0+r1*r1);e1=r1/sqrt(r0*r0+r1*r1); x0=x00+h*e0; x1=x01+h*e1;if(g1(x0

13、,x1)v=0&g2(x0,x1)v=0) fx=f(x0,x1);if(fxfl)fl=fx;for(i=0;i0|g2(x0,x1)0)break;fx=f(x0,x1);if(fxfl) fl=fx;else break;dox0=x0-h*d0;x1=x1-h*d1;h=0.7*h;if(h0|g2(x0,x1)0) con ti nue;fx=f(x0,x1);while(fx=fl);if(fabs(f0-fx)/f0)=c)xOO=xO;xO1=x1;fl=fx;fO=fx;elsebreak;printf(輸出最優(yōu)解為 nx1*=%f,x2*=%f, y*=%fn,xO,x1,

14、fx);運(yùn)行過(guò)程及結(jié)果：輸出最優(yōu)解為x1*=0.995421,x2*=1.004521,y*=1.0092006.四桿機(jī)構(gòu)優(yōu)化設(shè)計(jì)源程序：#in clude#i nclude#i nclude#define Pai 3.1415926&( -I2=O)&( 6-l1-l2=0)&(1-l2-4=0)&( I2-I1-4=O)&(I1*I1+I2*I2-1.414*I1*I2-16=O)&(36-I1*I1-I2*I2-1.414*I1*I2v=0)return (1);elsereturn (0);float fun (float x02)floatf,a31,b31,r31,p31,q31,

15、w31,x12;int i;p0=acos(1+x00)*(1+x00)-x0 1*x01+25)/(10+10*x00);q0=acos(1+x00)*(1+x00)-x01*x01-25)/(10*x01);f=0;for(i=1;i=30;i+)Pi=P0+(Pai/60)*i;if(-I1=0)ri=sqrt(26-10*cos(pi);ai=acos(ri*ri+x01*x01-x00*x00)/(2*ri*x01);while(g(z0,z1)=0);for(i=0;i=99;i+)bi=acos(ri*ri+24)/(10*ri)）；d0i=-1+2*(ra nd()/32767

16、.00);qi=Pai-ai-bi;for(j=0;jv=99;j+)wi=q0+(2*(pi-p0)*(pi-p0)/(3*Pai);d1j=-1+2*(ra nd()/32767.00);f=f+(Pai/60)*(qi-wi)*(qi-wijf0=fu n(z);)*(pi-pi-1);fl=fun( z);ss:return f;a=0.01;for(i=0,j=0;i=99&j=99;i+,j+)void mai n()n=1/sqrt(d0i)*(d0i)+d1j*d1floatj);a,q,f,fl,f0,l2,z2,d0100,d1100.d0i=n*d0i;x2,xi2,fx,m0,m1,e;d1j=n*d1j;int i,j, n,k;x0=z0+a*d0i;printf(輸入精度);x1=z1+a*d1j;scan f(%f,&e);if(g(x0,x1)=1)doz0=0+5*(ra nd()/32767.00);f=fun (x);doif(ffl)x0=x0-a*m0;x1=x1-a*m1;fl=f;a=0.7*a;m0=d0i;if(a=fl);if(fabs(f0-f)/f0)e)x0=l0;