《现代设计与方法》实验报告.docxVIP

《现代设计与方法》实验报告实验目的机械优化设计是一门实践性较强的课程，学生通过实际上机计算可以达到以下目的：加深对机械优化方法的基本理论和算法步骤的理解培养学生独立编制或调试计算机程序的能力掌握常用优化方法程序的使用方法培养学生灵活运用优化设计方法解决工程实际问题的能力实验要求序号实验项目实验要求1黄金分割法明确黄金分割法基本原理、计算步骤及程序框图；编制或调试黄金分割法应用程序；用测试题对所编程序进行测试；撰写实验报告。2复合形法明确复合形法基本原理、计算步骤及程序框图等；编制或调试复合形法应用程序；用测试题对所编程序进行测试；撰写实验报告。测试题黄金分割法程序测试题min,取min,取min，其中，取复合形法程序测试题1），2）3）程序及运行结果1.黄金分割法流程图（见下页）：程序（见下页）：#includestdio.h#includeconio.h#includemath.h#define e 0.001#define tt 0.01float function(float x){ float y=8*pow(x,3)-2*pow(x,2)-7*x+3; return(y);}void finding(float a[3], float f[3]){ float t=tt, a1,f1,ia; int i; 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;}float gold(float *ff){ float a1[3],f1[3],a[4],f[4]; float aa; int i; 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);}void main(){ float xx, ff; xx=gold(ff); printf(\nThe Optimal Design Result Is:\n); printf(\n\tx*=%f\n\tf*=%f, xx, ff);getch();}运行结果：1)将上述程序中“#define e 0.001”和“float y=8*pow(x,3)-2*pow(x,2)-7*x+3”分别改为“#define e 0.00001”和“float y=pow(x,2)-10*x+36”。结果：2)将上述程序中“#define e 0.001”和“floa