《机械优化设计》课程实践报告.docVIP

（课程实践报告封面模版） 合肥工业大学 《机械优化设计》课程实践 研究报告 班 级： 机设10 -04 学 号： 姓 名： 李健 授课老师： 王卫荣 日 期： 2012年 月 日  目录 一主要内容 1、一维搜索程序作业 A.λ = 0.618的证明..........................................1 B.编写用0.618法求函数极小值的程序..........................2 2、单位矩阵程序作业............................................4 3、其他工程优化问题..................................9 4连杆机构问题.....................................12 二实践心得体会...............................15 一: 主要内容 1. 一维搜索程序作业： λ = 0.618的证明 (y2 y1) 证明：0.618法要求插入点α1、α2 的位置相对于区间 [a，b] 两端点具有对称性，即 已知 a1=a2 , 要求α11=α22 由于α1=b-λ(b-a) α2=a+λ(b-a) 若使α11=α22 则有：b1-λ(b1-a1)=a2+λ(b2-a2)= a1+λ2(b1-a1) 因此: b1- a1=(λ2+λ)( b1- a1) ( b1- a1)( λ2+λ-1)=0 因为: b1= a1 所以: λ2+λ-1=0 则有: 取方程正数解得 若保留下来的区间为 [α1，b] ，根据插入点的对称性，也能推得同样的λ的值。 其0.618法的程序框图如下: B.编写用0.618法求函数极小值的程序 例：（1）a=0 ，b=2π，f（x）=cox（x） （2）a=0 ，b=10， f（x）=（x-2）2+3 （1） #includestdio.h #includemath.h void main(void) { int i; float a1,a2,aa,y1,y2,ymin,e; float a=0,b=2*3.14159,n=0.618; a1=b-n*(b-a); a2=a+n*(b-a); print(“输入精度：”)； scanf（“%f”,e）; for(i=0;i=10000;i=i++) { y1=cos(a1); y2=cos(a2); if(y1y2) { a=a1; a1=a2; a2=a+n*(b-a); } If(y1y2) b=a2; a2=a1; a1=b-n*(b-a); } if(fabs(b-a)/be)(fabs(y1-y2)/y2e)) {aa=(a+b)/2; ymin=cos(aa); printf(“x=%7.4f\tf(x)=%7.4f\n”),aa,ymin); break; } } } 运行结果： (2) #includestdio.h #includemath.h void main(void) { int i; float a1,a2,aa,y1,y2,ymin,e; float a=0,b=10,n=0.618; a1=b-n*(b-a); a2=a+n*(b-a); print(“输入精度：”)； scanf（“%f”,e）; for(i=0;i=10000;i=i++) { y1=（a1-2）*（a1-2）+3; y2=（a2-2）*（a2-2）+3; if(y1=y2) { a=a1; a1=a2; a2=a+n*(b-a); } If(y1y2) b=a2; a2=a1; a1=b-n*(b-a); } if(fabs(b-a)/be)(fabs(y1-y2)/y2e)) {aa=(a+b)/2; ymin=（aa-2）*（aa-2）+3; printf(“x=%6.3f\tf(x)=%6.3f\n”),aa,ymin); break; } } } 运行结果： 2.单位矩阵程序作业 编写生成单位矩阵的程序 程序文本 #include stdio