文档介绍：(课程实践报告封面模版)合肥工业大学《机械优化设计》课程实践研究报告班级:机设10-04学号:姓名:李健授课老师:王卫荣日期:2012年月日目录一主要内容1、一维搜索程序作业 =....................................................................22、单位矩阵程序作业............................................43、其他工程优化问题..................................94连杆机构问题.....................................12二实践心得体会...............................15一::λ=(y2>y1)证明:、α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],根据插入点的对称性,也能推得同样的λ的值。::(1)a=0,b=2π,f(x)=cox(x)(2)a=0,b=10,f(x)=(x-2)2+3(1)#include<>#include<>voidmain(void){inti;floata1,a2,aa,y1,y2,ymin,e;floata=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<>voidmain(void){inti;floata1,a2,aa,y1,y2,ymin,e;floata=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;}}}运行结果:#include<>voidmain(void){inta[100][100]; intN,i,j; printf("请输入所要输出矩阵的阶数(最多100阶):");scanf("%d",&N); printf("输出的矩阵阶数为%d\n",N); printf("N");/*****制作表头*****/ for(i=0;i<N;i++)/*****横行序号*****/ printf("%3d",i+1); printf("\n"); for(i=0;i<N+1;i++) printf("---");/*****分割线*****/ printf("\n"); for(i=0;i<100;i++)/*****数组赋值*****/ for(j=0;j<100;j++) { if(i==j) a[i][j]=1; else a[i][j]=0; } for(i=0;i<N;i++)/*****输出所需数组*****/ { prin