机械优化设计黄金分割法 外推法.doc

:#include<>#include<>#(doublex){doublem; m=x*x-10*x+36;returnm;}voidmain(){doubleh0=R,y1,y2,y3,x1,x2,x3,h;x1=0;h=h0;x2=h;y1=fun(x1);y2=fun(x2);if(y2>y1){h=-h;x3=x1;y3=y1;x1=x2;y1=y2;x2=x3;y2=y3;}x3=x2+h;y3=fun(x3);while(y3<y2){h*=;x1=x2;y1=y2;x2=x3;y2=y3;x3=x2+h;y3=fun(x3);} printf("fun(%f)=%f,fun(%f)=%f,fun(%f)=%f\n",x1,y1,x2,y2,x3,y3);}运行过程及结果:fun()=,fun()=,fun()=:#include<>#include<>#definef(x)x*x*x*x-5*x*x*x+4*x*x-6*x+60doublehj(double*a,double*b,doublee,int*n){ doublex1,x2,s;if(fabs((*b-*a)/(*b))<=e)s=f((*b+*a)/2);else {x1=*b-*(*b-*a);x2=*a+*(*b-*a);if(f(x1)>f(x2))*a=x1;else*b=x2;*n=*n+1;s=hj(a,b,e,n); }returns;}voidmain(){ doubles,a,b,e,m;intn=0;printf("输入a,b值和精度e值\n");scanf("%lf%lf%lf",&a,&b,&e);s=hj(&a,&b,e,&n);m=(a+b)/2;printf("a=%lf,b=%lf,s=%lf,m=%lf,n=%d\n",a,b,s,m,n);}运行过程及结果:输入a,b值和精度e值-=,b=,s=,m=,n=:#include<>#include<>intmain(void){doublea1,a2,a3,ap,y1,y2,y3,yp,c1,c2,m;doublej[3];inti,h=1;voidfinding(doublea[3]);finding(j);a1=j[0];a2=j[1];a3=j[2];m=;doublef(doublex);y1=f(a1);y2=f(a2);y3=f(a3);for(i=1;1>=1;i++){c1=(y3-y1)/(a3-a1);c2=((y2-y1)/(a2-a1)-c1)/(a2-a3);ap=*(a1+a3-c1/c2);yp=f(ap);if(fabs((y2-yp)/y2)<m)