非线性规划及MATLAB解法及其应用.docVIP

题 目 非线性规划的MATLAB解法及其用线性规划的MATLAB解法 题二： 某厂生产一种产品有甲、乙两个牌号，讨论在产销平衡的情况下如何确定各自的产量，使总利润最大. 所谓产销平衡指工厂的产量等于市场上的销量.符号说明：z(x1,x2)表示总利润；p1，q1，x1分别表示甲的价格、成本、销量； p2，q2，x2分别表示乙的价格、成本、销量； aij，bi，λi,ci（i，j =1，2）是待定系数. 题三：设有400万元资金, 要求4年内使用完, 若在一年内使用资金x万元, 则可得效益万元(效益不能再使用),当年不用的资金可存入银行, 年利率为10%. 试制定出资金的使用计划, 以使4年效益之和为最大. 数据结构 题一:设剪去的正方形的边长为x，则水槽的容积为：;建立无约束优化模型为：min y=-， 0x1.5 题二：总利润为： z(x1,x2)=(p1-q1)x1+(p2-q2)x2 若根据大量的统计数据,求出系数b1=100,a11=1,a12=0.1,b2=280, a21=0.2,a22=2,r1=30,λ1=0.015,c1=20, r2=100,λ2=0.02,c2=30,则 问题转化为无约束优化问题：求甲,乙两个牌号的产量x1，x2，使 总利润z最大. 为简化模型,先忽略成本,并令a12=0,a21=0,问题转化为求: z1 = ( b1 - a11x1 ) x1 + ( b2 - a22x2 ) x2 的极值. 显然其解为x1 = b1/2a11 = 50, x2 = b2/2a22 = 70, 我们把它作为原问题的初始值. 题三：设变量表示第i年所使用的资金数,则有 源程序 题一： 编写M文件fun0.m: function f=fun0(x) f=-(3-2*x).^2*x; 主程序为wliti2.m: [x,fval]=fminbnd(fun0,0,1.5); xmax=x fmax=-fval 题二： 建立M-文件fun.m: function f = fun(x) y1=((100-x(1)- 0.1*x(2))-(30*exp(-0.015*x(1))+20))*x(1); y2=((280-0.2*x(1)- 2*x(2))-(100*exp(-0.02*x(2))+30))*x(2); f=-y1-y2; 输入命令: x0=[50,70]; x=fminunc(‘fun’,x0), z=fun(x) 题三： 建立M文件 fun44.m,定义目标函数: function f=fun44(x) f=-(sqrt(x(1))+sqrt(x(2))+sqrt(x(3))+sqrt(x(4))); 建立M文件mycon1.m定义非线性约束： function [g,ceq]=mycon1(x) g(1)=x(1)-400; g(2)=1.1*x(1)+x(2)-440; g(3)=1.21*x(1)+1.1*x(2)+x(3)-484; g(4)=1.331*x(1)+1.21*x(2)+1.1*x(3)+x(4)-532.4; ceq=0 主程序youh4.m为: x0=[1;1;1;1];vlb=[0;0;0;0];vub=[];A=[];b=[];Aeq=[];beq=[]; [x,fval]=fmincon(fun44,x0,A,b,Aeq,beq,vlb,vub,mycon1) 运行结果 题一: 运算结果为: xmax = 0.5000,fmax =2.0000.即剪掉的正方形的边长为0.5米时水槽的容积最大,最大容积为2立方米. 题二： 运行结果为:x=23.9025, 62.4977, z=6.4135e+003即甲的产量为23.9025,乙的产量为62.4977,最大利润为6413.5. 题三： 运行结果为：x1=86.2;x2=104.2;x3=126.2;x4=152.8;z=43.1 相关知识 用Matlab解无约束优化问题 一元函数无约束优化问题 常用格式如下： （1）x= fminbnd (fun,x1,x2) （2）x= fminbnd (fun,x1,x2 ，options) （3）[x，fval]= fminbnd（...） （4）[x，fval，exitflag]= fminbnd（...） （5）[x，fval，exitflag，output]= fminbnd（...） 其中（3）、（4）、（5）的等式右边可选用（1）或（2）的等式右边。 函数fminbnd的算法基于黄金分割法和二次插值法，它要求目标函数必须是连续