机械优化设计黄金分割法程序.doc

《机械优化设计黄金分割法程序.doc》由会员分享，可在线阅读，更多相关《机械优化设计黄金分割法程序.doc（3页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、机械优化设计黄金分割法程序机械优化设计实验报告黄金分割法：已知：F(x）=（x+1）（x2）2， include stdio.h#include conio。h#include #define e 0。00001#define tt 0.1float f（double x）float y=pow（x,3)-3pow（x,2)+4;return（y）;void finding（float a3,float f3)float h=tt，a1，f1,ia;a0=0；f0=function(a0）；for（int i=0；i+)a1=a0+t; f1=function(a1);if(f1f0) bre

2、ak；if（fabs(f1-f0)=e） t=-t; a0=a1; f0=f1；else if(ia=1) return；t=t/2; ia=1;for(i=0；i+）a2=a1+t; f2=function(a2);if（f2f1） break；t=2t;a0=a1；f0=f1；a1=a2；f1=f2;if（a0a2)a1=a0;f1=f0;a0=a2;f0=f2;a2=a1；f2=f1；return；float gold(float *ff)float a13，f13,a4，f4；float aa;finding（a1，f1）;a0=a10; f0=f10；a3=a12; f3=f12;a

3、1=a0+0.382(a3a0);a2=a0+0。618（a3a0）;f1=function（a1）； f2=function(a2）；for(int i=0;；i+）if（f1=f2)a0=a1； f0=f1;a1=a2; f1=f2；a2=a0+0。618*(a3a0)； f2=function(a2）；else a3=a2;f3=f2；a2=a1；f2=f1；a1=a0+0。382*（a3-a0)； f1=function(a1);if(（a3-a0)e）aa=(a1+a2)/2; *ff=function(aa);break；return(aa); void main（）float xx，ff; xx=gold(ff）;printf(“n The Optimal Design Result Is：n”）; printf（“n tx* =fntf =f”，xx，ff）；getch（）；