数学建模中规划方案问题的MATLAB求解培训.pptVIP

第二章 规划问题的MATLAB求解 2.1 线性规划 2.2 非线性规划 2.3 整数规划 数学规划模型 序： 规划问题是常见的数学建模问题，离散系统的优化问题一般都可以通过规划模型来求解。因此快速求解规划问题是数学建模的基本素质。利用MATLAB提供的规划模型求解命令，可以快速得到想要的结果。 2.1 线性规划 人们在生产实践中经常会遇到如何利用现有资源安排生产，取得最大经济效益的问题。此类问题构成了运筹学的一个重要分支——数学规划。而线性规划是数学规划中的重要组成部分。 1947年美国数学家G.B. Dantzig提出线性规划的一般数学模型和求解线性规划的通用方法——单纯形法，为学科奠定基础。 1947年美国数学家J.von诺依曼提出对偶理论，开创了线性规划的许多新领域，扩大的应用范围和解题能力。 1951年美国经济学家T.C.库普曼斯把线性规划应用到经济领域，并与康拓罗维奇一起获得1975年诺贝尔经济学奖。 2.1.1 线性规划的实例与定义 【例2-1】某电视台播放两套宣传片，宣传甲片播映时间为3.5min，广告时间为30s，观众60万；宣传乙片播映时间为1min，广告时间为1min，收视观众为20万。广告公司规定每周至少有3.5min广告，电视台只能为该栏目宣传不小于16min的节目。电视台每周至少要播两套宣传片多少次，才能使人使收视观众最多？ 类别 播甲片 播放乙片 节目要求 片集时间（min） 3.5 1 =16 广告时间（min） 0.5 1 =3.5 收视观众/万 60 20 2.1.2线性规划的MATLAB标准形式 c=[2 3 1]; a=[1 4 2;3 2 0]; b=[8 6]; [x,y]=linprog(c,-a,-b,[],[],zeros(3,1)) 2.1.4 求解线性规划的MATLAB解法 [x,fval]=linprog(c,A,b,Aeq,beq,LB,UB,X0,OPTIONS) 其中的fval为返回目标函数的值；A，b对应不等式约束；Aeq，beq对应等式约束；LB，UB为变量x的下界和上界。X0为x的初始值；OPTIONS是控制参数。 f=[-5 -4 -6]; A=[1 -1 1;3 2 4;3 2 0]; b=[20 42 30]; lb=zeros(3,1); [x,fval,exitflag,output,lambda]=linprog(f,A,b,[],[],lb) 2.2 非线性规划 如果规划问题中目标函数或约束条件中包含非线性函数，这种规划问题称为非线性规划问题。 2.2.1 非线性规划 特点：非线性规划问题比线性规划问题的求解要复杂得多，没有统一的解法。各种方法有自己特定的适用范围。 2.2.2 非线性规划的MATLAB求解 MATLAB命令： [x,y] = fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon,options) function f=fun1(x) f=x(1)^2+x(2)^2+8; function [g,h]=fun2(x) g=-x(1)^2+x(2); h=-x(1)-x(2)^2+2; options=optimset; [x,y]=fmincon(fun1,rand(2,1),[],[],[],[],zeros(2,1),[],... fun2,options) function f=fun44(x) f=-(sqrt(x(1))+sqrt(x(2))+sqrt(x(3))+sqrt(x(4))); x0=[1;1;1;1];lb=[0 0 0 0];ub=[];A=[];b=[];Aeq=[];beq=[]; [x,fval]=fmincon(fun44,x0,A,b,Aeq,beq,lb,ub,mycon1) 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; 结果：x = 86.2890 103.7586 127.1018 152.1895 fval = -43.0858 2.2.3 二次规划 若非线性规划的目标函数为自变量x的二次函数，约束条件是线性的，则称这种规划为二次规划。 求解命令一般格式 x = quadprog(H,f,A,b,Aeq,beq,