现代设计理论方法黄金分割法算法编程及调试结果

《现代设计理论与方法》实验报告实验目的机械优化设计是一门实践性较强的课程，学生通过实际上机计算可以达到以下目的：加深对机械优化设计方法的基本理论和算法步骤的理解；培养学生独立编制或调试计算机程序的能力；掌握常用优化方法程序的使用方法；培养学生灵活运用优化设计方法解决工程实际问题的能力。实验项目、学时分配及对实验项目的要求序号实验项目学时实验要求1黄金分割法2明确黄金分割法基本原理、计算步骤及程序框图；编制或调试黄金分割法应用程序；用测试题对所编程序进行测试；撰写实验报告。测试题黄金分割法测试题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程序及调试结果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