机械优化设计外推 黄金分割法 二次插值法 随机方向法 坐标轮换法 四杆机构.docVIP

《机械优化设计》上机源程序与结果 4 程序一：外推法 #includestdio.h #includemath.h #define B 0.001 double fun(double x) { double min; min=x*x-10*x+36; return min; } void main() { double h0=B,y1,y2,y3,x1,x2,x3,h; x1=0;h=h0;x2=h; y1=fun(x1);y2=fun(x2); if(y2y1) {h=-h; x3=x1; x1=x2; x2=x3; y3=y1; y1=y2; y2=y3; } x3=x2+h;y3=fun(x3); while(y3y2) {h*=1.5; x1=x2; x2=x3; x3=x2+h; y1=y2; y2=y3; y3=fun(x3); } printf(fun(%f)=%f,\nfun(%f)=%f,\nfun(%f)=%f\n,x1,y1,x2,y2,x3,y3); } 运行过程及结果： fun(2.954784)=15.182909, fun(4.432676)=11.321857, fun(6.649513)=13.720895 程序二：黄金分割法 源程序： #includestdio.h #includemath.h #define f(x) x*x