机械优化设计程序.docVIP

序号13—18的同学 1、 参考《MATLAB实用数值分析》 张德丰 清华大学出版社 P361 二次规划 数学模型 min1/2xTHx+cx Ax=b s.t. Aeqx=beq Lb=x=ub 将目标函数化为标准形式： 代码： clear all; clc; H=[4,-4;-4,6]; f=[3,-4]; A=[2,1;-1,2]; b=[4;4]; lb=zeros(2,1); [x,fval,exitflag]=quadprog(H,f,A,b,[],[],lb) 运行结果： x = 0 0.6667 fval = -1.3333 exitflag = 1 2、 3、 这是个多元约束优化问题 首先编写目标函数的M文件，代码如下： function f=li9_12fun(x) f=-x(1)*x(2)*x(3); 由于约束条件为非线性不等式约束，因此，编写一个约束条件的M文件，代码如下： function[c,ceq]=li9_12non(x) c=[]; ceq=2*x(1)*x(2)+2*x(1)*x(3)+2*x(2)*x(3)-150; 实现代码： clear all; clc; x0=[1;1;1]; [x,fval,exitflag]=fmincon(@fun,x0,[],[],[],[],[],[],@mycon) 运行结果： x = 5.0000 5.0000 5.0000 fval = -125.0000 exitflag = 5 表示方向导数的级小于给定的容忍度 4、求函数的极小值。 参考《MATLAB实用数值分析》 张德丰 清华大学出版社 P350 这是一个无约束优化问题；实现代码如下： clear all; clc; f=@(x)2*x(1)^2+x(2)^2+2*x(1)*x(2)+x(1)-x(2); x0=[0;0]; [x,fval,exitflag]=fminsearch(f,x0) 结果如下： x = -1.0000 1.5000 fval = -1.2500 exitflag = 1 5、 参考《MATLAB实用数值分析》 张德丰 清华大学出版社 P337 这是一个线性规划问题；实现代码如下： clear all; clc; F=[1,1,1]; A=[1,2,3;2,1,5;1,-4,2]; b=[15;20;10]; lb=zeros(2,1); [x,fval,exitflag]=linprog(F,A,b,[],[],lb) 结果如下： x = 1.0e+018 * 0.0000 0.0000 -2.2904 fval = -2.2904e+018 exitflag = -3 表示所求问题是无界的 6、有一块边长为6m的正方形铝板，四角截去相等的边长为x的方块并折转，造一个无盖的箱子，问如何截法（x取何值）才能获得最大容器的箱子。 7、任务分配问题：某车间有甲、乙两台机床，可用于加工三种工件。假定这两台车床的可用台时数分别为800和900，三种工件的数量分别为400、600和500，且已知用三种不同车床加工单位数量不同工件所需的台时数和加工费用如下表。问怎样分配车床的加工任务，才能既满足加工工件的要求，又使加工费用最低？ 车床类 型 单位工件所需加工台时数 单位工件的加工费用 可用台时数 工件1 工件2 工件3 工件1 工件2 工件3 甲 0.4 1.1 1.0 13 9 10 800 乙 0.5 1.2 1.3 11 12 8 900 8、已知：轴的一端作用载荷 P=1000N，扭矩 M=100N·m；轴长不得小于8cm；材料的许用弯曲应力 [σw]=120MPa，许用扭剪应力 [τ]= 80MPa，许用挠度 [f] = 0.01cm；密度[ρ] = 7.8t /m，弹性模量E=2×105MPa。 要求：设计销轴，在满足上述条件的同时，轴的质量应为最轻。