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

合肥工业大学《机械优化设计》课程实践研究报告班级:学号:姓名:授课教师:日期:=、+自行选择小型机械设计问题或其他工程优化问题(1)分析优化对象,根据设计问题的要求,选择设计变量,确立约束条件,建立目标函数,建立优化设计的数学模型并编制问题程序;(2)选择适当的优化方法,简述方法原理,进行优化计算;(3)进行结果分析,并加以说明。=、:aα1α2b1-λλ1aα3α1α2λ(1-λ)λ2λ黄金分割法要求插入点α1,α2的位置相对于区间[a,b]两端点具有对称性,即α1=b-λ(b-a)α2=b+λ(b-a)其中λ为待定常数。除了对称要求外,黄金分割法还要求在保留下来的区间内再插入一点,所形成的新三段与原来区间的三段具有相同的比例分布,故有1-λ=λ2取方程正数解,得λ=5-12≈:(以正弦函数y=sinx为例)#include<>#include<>intmain(){doublea,b,c=,x[3],y[3],d;printf("请输入区间[a,b]的值以及精度:\n");scanf("%lf,%lf,%lf",&a,&b,&d);x[1]=b-c*(b-a);x[2]=a+c*(b-a);y[1]=sin(x[1]);y[2]=sin(x[2]);do{if(y[1]>y[2]){a=x[1];x[1]=x[2];y[1]=y[2];x[2]=a+c*(b-a);y[2]=sin(x[2]);}else{b=x[2];x[2]=x[1];y[2]=y[1];x[1]=b-c*(b-a);y[1]=sin(x[1]);}} while(fabs((b-a)/b)>d);x[0]=(a+b)/2;y[0]=sin(x[0]);printf("极小点x*=%lf\n",x[0]);printf("极小值y=%lf\n",y[0]);}C语言程序运行结果::#include<>intmain(){ inti,j,N,y; printf("请输入单位矩阵阶数N:"); scanf("%d",&N); for(i=1;i<=N;i++) for(j=1;j<=N;j++) { if(i==j) y=1; else y=0; printf("%d",y); if(j==N) printf("\n"); }}C语言程序运行结果::clc;clearall;N=input('请输入单位矩阵阶数N:');eye(N)matlab程序运行结果:(FORTRAN语言)C连杆机构问题函数子程序C 目标函数============== SUBROUTINEFFX(N,X,FX);计算X点的目标函数值的子程序C ====================== DIMENSIONX(N);MON/ONE/I1,I2,I3,I4,NFX,I6;公用语句,程序段传输数据 NFX=NFX+1;目标函数子程序调用次数计数器 P0=ACOS(((+X(1))**2-X(2)**2+)/(*(+X(1))))Q0=ACOS(((+X(1))*+*2-X(2)**2-)/(*X(2))) T=*/(*);将考察范围输入角30等分后每一份的值 FX=;连加先置0 DO10K=0,30;循环程序入口,循环30次 PI=P0+K*T;计算每一次循环后的输入角 QE=Q0+*(PI-P0)**2/(*) D=SQRT(-*COS(PI)) AL=ACOS((D*D+X(2)*X(2)-X(1)*X(1))/(*D*X(2))) BT=ACOS((D*D+)/(*D))IF()THENQI=-AL-BT;输入角在0到pi的实际输出角ELSEQI=-AL+BT;输入角在大于pi的实际输出角ENDIFIF()THENFX=FX+(QI-QE)**2*T;数值积分的非首点,非末点ELSEFX=FX+(QI-QE)**2*T/;数值积分的首端点,末端点ENDIF