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

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