20100575-黄远洋-《机械优化设计》课程实践报告.docVIP

合肥工业大学 《机械优化设计》课程实践 研究报告 班 级： 机设10-6班 学 号： 姓 名： 黄远洋 授课老师： 王卫荣 日 期： 2013年 5 月 9 日 目录 一维搜索程序作业 【1.1.】例1 【1.1.1】程序文本 【1.1.2】输出结果截图 【1.2】例2 【1.2.1】程序文本 【1.2.2】输出结果截图 单位矩阵程序作业 【2.1】程序文本 【2.2】输出结果截图 连杆机构问题 【3.1】设计变量 【3.2】目标函数 【3.3】约束条件 【3.4】优化方法 【3.5】程序文本 【3.6】数据输入截图 【3.7】输出结果 4、自行选择小型机械设计问题或其他工程优化问题 【4.1】设计变量 【4.2】目标函数 【4.3】约束条件 【4.4】程序文本 【4.5】输入数据截图 5、课程实践心得体会 一维搜索程序作业； 作业：编写用0.618法求函数极小值的程序。 【1.1】例：a=0,b=2,f(x)=cos(x); 【1.1.1】程序文本 #includestdio.h #includemath.h void main() { float a,b,c=0.618,aa[3],y[3],d; scanf(%f,%f,%f,a,b,d); aa[1]=b-c*(b-a); aa[2]=a+c*(b-a); y[1]=cos(aa[1]); y[2]=cos(aa[2]); do{ if(y[1]y[2]) { a=aa[1];aa[1]=aa[2];y[1]=y[2]; aa[2]=a+c*(b-a); y[2]=cos(aa[2]); } else { b=aa[2];aa[2]=aa[1];y[2]=y[1]; aa[1]=b-c*(b-a); y[1]=cos(aa[1]); } }while(fabs((b-a)/b)d); aa[0]=(a+b)/2; y[0]=cos(aa[0]); printf(a*=%f\n,aa[0]); printf(y=%f\n,y[0]); } 【1.1.2】输出结果截图 输入a=0,b=2,精度d=0.0000001,输出极小值点和函数极小值如下： 【1.2】例a=0,b=10,f(x)=(x-2)+3 【1.2.1】程序文本 #includestdio.h #includemath.h void main() { float a,b,c=0.618,aa[3],y[3],d; scanf(%f,%f,%f,a,b,d); aa[1]=b-c*(b-a); aa[2]=a+c*(b-a); y[1]=(aa[1]-2)*(aa[1]-2)+3; y[2]=(aa[2]-2)*(aa[2]-2)+3; do{ if(y[1]y[2]) { a=aa[1];aa[1]=aa[2];y[1]=y[2]; aa[2]=a+c*(b-a); y[2]=(aa[2]-2)*(aa[2]-2)+3; } else { b=aa[2];aa[2]=aa[1];y[2]=y[1]; aa[1]=b-c*(b-a); y[1]=(aa[1]-2)*(aa[1]-2)+3; } }while(fabs((b-a)/b)d); aa[0]=(a+b)/2; y[0]=(aa[0]-2)*(aa[0]-2)+3; printf(a*=%f\n,aa[0]); printf(y=%f\n,y[0]); } 【1.2.2】输出结果截图 输入a=0,b=10,精度d=0.0000001，输入极小值点和函数极小值如下： 2、单位矩阵程序作业； 作业：编写生成单位矩阵的程序。 要求：通用、输出美观、语言少为佳。 【2.1】程序文本 #includestdio.h void main (void) #define m 500 { int i,j,n,a[m][m]; printf(please input a number); scanf(%d,n); for(i=1;i=n;i++) { for(j=1;j=n;j+