常用方法MATLAB求解.pptVIP

* [x,fval]=fmincon(‘fun’,X0,A,b,Aeq,beq,VLB,VUB,’nonlcon’) 说明：fan为目标函数的表达式字符串；如果没有不等式：AX≤b存在，则令A=[],b=[]；如果没有等式约束： AeqX=beq存在，则令Aeq=[],beq=[]；如果没有不等式：vlb≤x≤vub约束，则令vlb=[],ulb=[]；nonlcon的作用是通过接收的向量x来计算非线性不等约束G(x)≤0和等式约束Ceq(x)=0分别在x处的估计G和Ceq, 通过指定函数柄来使用；函数可能会出现局部最优解，这与初值x0的选取有关。 * 例４ 1）先建立M文件 fun4.m,定义目标函数: function f=fun4(x); f=exp(x(1))*(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1); 2）再建立M文件mycon.m定义非线性约束： function [g,ceq]=mycon(x) g=[ 1.5+x(1)*x(2)-x(1)-x(2);-x(1)*x(2)-10]; ceq=[]; * 3）主程序youh3.m为: x0=[-1;1]; A=[];b=[]; Aeq=[1 1];beq=[0]; vlb=[];vub=[]; [x,fval]=fmincon(fun4,x0,A,b,Aeq,beq,vlb,vub,mycon) 3. 运算结果为： x = -1.2250 1.2250 fval = 1.8951 * （2）作出三次插值海底地形图 输入命令 ： x0=[129 140 103.5 88 185.5 195 105.5 157.5 107.5 77 81 162 162 117.5]; y0=[7.5 141.5 23 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.5 84 -33.5]; z0=[-4 -8 -6 -8 -6 -8 -8 -9 -9 -8 -8 -9 -4 -9]; y=-70:0.5:150; z=griddata(x0,y0,z0,x,y,cubic); meshz(x,y,z); xlabel(x轴); ylabel(y轴); zlabel(z轴) title(三次插值海底地形图) x=75:0.5:200; * 结果见图 * （3）作出水深低于5英尺的海域范围 输入命令 ： contour(x,y,z,[-5,-5],’-k’); grid, xlabel(x轴); ylabel(y轴); title(船不能进入区域) 结果见图 * 从图可以看出，船应避免进入危险区域 （-10，110）。 （110，190） * 四、MATLAB软件在优化中的应用 在研究与解决具体问题中，经常遇到有关优化问题，下面介绍用MATLAB软件求解一些优化问题，包括求解线性规划和非线性规划等问题。 线性规划是运筹学的一个重要分支，它起源于工业生产组织管理的决策问题。在数学上它用来确定多变量线性函数在变量满足线性约束条件下的最优值；随着计算机的发展，出现了如单纯形法等有效算法，它在工农业、军事、交通运输、决策管理与规划等领域中有广泛的应用。 * 问题一 ： 任务分配问题：某车间有甲、乙两台机床，可用于加工三种工件。假定这两台车床的可工作时间分别为800和900，三种工件的数量分别为400、600和500，且已知车床甲加工单位数量三种工件所需的时间和加工费分别为0.4、1.1、1和13、9、10，车床乙加工单位数量三种工件所需的时间和加工费分别为0.5、1.2、1.3和11、12、8。问怎样分配车床的加工任务，才能既满足加工工件的要求，又使加工费用最低？ * 解 设在甲车床上加工工件1、2、3的数量分别为x1、x2、x3，在乙车床上加工工件1、2、3的数量分别为x4、x5、x6。 可建立以下线性规划模型： * 问题二： 某厂每日8小时的产量不低于1800件。为了进行质量控制，计划聘请两种不同水平的检验员。一级检验员的标准为：速度25件/小时，正确率98%，计时工资4元/小时；二级检验员的标准为：速度15件/小时，正确率95%，计时工资3元/小时。检验员每错检一次，工厂要损失2元。为使总检验费用最省，该工厂应聘一级、二级检验员各几名？ 解 设需要一级和二级检验员的人数分别为x1、x2人, 则应付检验员的工资为： 因检验员错检而造成的损失为： 故目标函数为： * 约束条件为： 线性规划模型： * 丁的蛙泳成绩退步到1’15”2；戊的自由泳成绩进步到57”5,