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

 PAGE \* MERGEFORMAT 4 《现代设计理论与方法》实验报告 实验目的 机械优化设计是一门实践性较强的课程，学生通过实际上机计算可以达到以下目的： 加深对机械优化设计方法的基本理论和算法步骤的理解； 培养学生独立编制或调试计算机程序的能力； 掌握常用优化方法程序的使用方法； 培养学生灵活运用优化设计方法解决工程实际问题的能力。 实验项目、学时分配及对实验项目的要求 序号实验项目学时实验要求 1 黄金分割法 2明确黄金分割法基本原理、计算步骤及程序框图； 编制或调试黄金分割法应用程序； 用测试题对所编程序进行测试； 撰写实验报告。测试题 黄金分割法测试题 minF (x) = x2 – 10x + 36，取  QUOTE x0 x0=0，  QUOTE h0 h0=0.01， QUOTE = QUOTE 10-5 10-5 minF (x) = x4 – 5 x3 +4 x2 – 6x + 60，取  QUOTE x0 x0=0，  QUOTE h0 h0=0.01， QUOTE = QUOTE 10-5 10-5 minF (x) =（x+1） QUOTE ，其中 QUOTE x鈮? x鈮?0， 取  QUOTE x0 x0=0，  QUOTE h0 h0=0.01， QUOTE = QUOTE 10-5 10-5 程序及调试结果 minF (x) = x2 – 10x + 36，取  QUOTE x0 x0=0，  QUOTE h0 h0=0.01， QUOTE = QUOTE 10-5 10-5 #include stdio.h #include conio.h #include math.h #define e 0.00001 #define tt 0.01 float function(float x) { float y=x*x-10*x+36; 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