現(xiàn)代設(shè)計理論方法黃金分割法算法編程及調(diào)試結(jié)果

《現(xiàn)代設(shè)計理論與方法》實(shí)驗報告實(shí)驗?zāi)康臋C(jī)械優(yōu)化設(shè)計是一門實(shí)踐性較強(qiáng)的課程，學(xué)生通過實(shí)際上機(jī)計算可以達(dá)到以下目的：加深對機(jī)械優(yōu)化設(shè)計方法的基本理論和算法步驟的理解；培養(yǎng)學(xué)生獨(dú)立編制或調(diào)試計算機(jī)程序的能力；掌握常用優(yōu)化方法程序的使用方法；培養(yǎng)學(xué)生靈活運(yùn)用優(yōu)化設(shè)計方法解決工程實(shí)際問題的能力。實(shí)驗項目、學(xué)時分配及對實(shí)驗項目的要求序號實(shí)驗項目學(xué)時實(shí)驗要求1黃金分割法2明確黃金分割法基本原理、計算步驟及程序框圖；編制或調(diào)試黃金分割法應(yīng)用程序；用測試題對所編程序進(jìn)行測試；撰寫實(shí)驗報告。測試題黃金分割法測試題minF(x)=x2–10x+36，取x0=0，h0=0.01，=minF(x)=x4–5x3+4x2–6x+60，取x0=0，h0=0.01，=minF(x)=（x+1），其中x鈮?0，取x0=0，h0=0.01，=10程序及調(diào)試結(jié)果minF(x)=x2–10x+36，取x0=0，h0=0.01，=#include<stdio.h>#include<conio.h>#include<math.h>#definee0.00001#definett0.01floatfunction(floatx){floaty=x*x-10*x+36;return(y);}voidfinding(floata[3],floatf[3]){floatt=tt,a1,f1,ia;inti;a[0]=0;f[0]=function(a[0]);for(i=0;;i++){a[1]=a[0]+t;f[1]=function(a[1]);if(f[1]<f[0])break;if(fabs(f[1]-f[0])>=e){t=-t;a[0]=a[1];f[0]=f[1];}else{if(ia==1)return;t=t/2;ia=1;}}for(i=0;;i++){a[2]=a[1]+t;f[2]=function(a[2]);if(f[2]>f[1])break;t=2*t;a[0]=a[1];f[0]=f[1];a[1]=a[2];f[1]=f[2];}if(a[0]>a[2]){a1=a[0];f1=f[0];a[0]=a[2];f[0]=f[2];a[2]=a1;f[2]=f1;}return;}floatgold(float*ff){floata1[3],f1[3],a[4],f[4];floataa;inti;finding(a1,f1);a[0]=a1[0];f[0]=f1[0];a[3]=a1[2];f[3]=f1[2];a[1]=a[0]+0.382*(a[3]-a[0]);a[2]=a[0]+0.618*(a[3]-a[0]);f[1]=function(a[1]);f[2]=function(a[2]);for(i=0;;i++){if(f[1]>=f[2]){a[0]=a[1];f[0]=f[1];a[1]=a[2];f[1]=f[2];a[2]=a[0]+0.618*(a[3]-a[0]);f[2]=function(a[2]);}else{a[3]=a[2];f[3]=f[2];a[2]=a[1];f[2]=f[1];a[1]=a[0]+0.382*(a[3]-a[0]);f[1]=function(a[1]);}if((a[3]-a[0])<e){aa=(a[1]+a[2])/2;*ff=function(aa);break;}}return(aa);}voidmain(){floatxx,ff;xx=gold(&ff);printf("\nTheOptimalDesignResultIs:\n");printf("\n\tx*=%f\n\tf*=%f",xx,ff);getch();}1）minF(x)=x2–10x+36，取x0=0，h0=0.01，=10-52）minF(x)=x4–5x3+4x2–6x+60，取x0=0，h0=0.01