【数学建模】第4讲_线性规划模型VIP免费

线性规划模型数学建模基本内容2.用数学软件包MATLAB求解线性规划问题.4.实验作业.1.两个引例.3.建模案例：投资的收益与风险.问题一：任务分配问题：某车间有甲、乙两台机床，可用于加工三种工件.假定这两台车床的可用台时数分别为800和900，三种工件的数量分别为400、600和500，且已知用三种不同车床加工单位数量不同工件所需的台时数和加工费用如下表.问怎样分配车床的加工任务，才能既满足加工工件的要求，又使加工费用最低？单位工件所需加工台时数单位工件的加工费用车床类型工件1工件2工件3工件1工件2工件3可用台时数甲0.41.11.013910800乙0.51.21.311128900两个引例解设在甲车床上加工工件1、2、3的数量分别为x1、x2、x3，在乙车床上加工工件1、2、3的数量分别为x4、x5、x6,可建立以下线性规划模型：6543218121110913minxxxxxxz142536123456x400600500s.t.0.41.18000.51.21.39000,1,2,,6ixxxxxxxxxxxxi解答问题二：某厂每日8小时的产量不低于1800件.为了进行质量控制，计划聘请两种不同水平的检验员.一级检验员的标准为：速度25件/小时，正确率98%，计时工资4元/小时；二级检验员的标准为：速度15件/小时，正确率95%，计时工资3元/小时.检验员每错检一次，工厂要损失2元.为使总检验费用最省，该工厂应聘一级、二级检验员各几名？解设需要一级和二级检验员的人数分别为x1、x2人,则应付检验员的工资为：212124323848xxxx因检验员错检而造成的损失为：21211282)%5158%2258(xxxx故目标函数为：2121213640)128()2432(minxxxxxxz约束条件为：0,0180015818002581800158258212121xxxxxx线性规划模型：213640minxxz12121253459s.t.150,0xxxxxx解答返回线性规划模型的一般形式11min,1,2,...,.s.t.0,1,2,...,.niiinikkikiucxaxbimxin目标函数和所有的约束条件都是设计变量的线性函数.min.s.tucxAxbvlbxvub矩阵形式：实际问题中的优化模型T1min(max)(),(,,)s.t.()0,1,2,,nizfxxxxgxim或x是决策变量f(x)是目标函数gi(x)0是约束条件数学规划线性规划(LP)二次规划(QP)非线性规划(NLP)纯整数规划(PIP)混合整数规划(MIP)整数规划(IP)0-1整数规划一般整数规划连续规划优化模型的分类用MATLAB优化工具箱解线性规划minz=cXs.t.AXb1.模型：命令：x=linprog（c,A,b）2.模型：minz=cXs.t.AXbbeqXAeq命令：x=linprog（c，A，b，Aeq,beq）注意：若没有不等式：存在，则令A=[]，b=[].bAX3.模型：minz=cXs.t.AXbbeqXAeqVLB≤X≤VUB命令：[1]x=linprog（c，A，b，Aeq,beq,VLB，VUB）[2]x=linprog（c，A，b，Aeq,beq,VLB，VUB,X0）注意：[1]若没有等式约束:,则令Aeq=[],beq=[].[2]其中X0表示初始点beqXAeq4.命令：[x,fval]=linprog(…)返回最优解ｘ及ｘ处的目标函数值fval.解编写M文件xxgh1.m如下：c=[-0.4-0.28-0.32-0.72-0.64-0.6];A=[0.010.010.010.030.030.03;0.02000.0500;00.02000.050;000.03000.08];b=[850;700;100;900];Aeq=[];beq=[];vlb=[0;0;0;0;0;0];vub=[];[x,fval]=linprog(c,A,b,Aeq,beq,vlb,vub)例1max6543216.064.072.032.028.04.0xxxxxxz123456s.t.0.010.010.010.030.030.03850xxxxxx70005.002.041xx10005.002.052xx90008.003.063xx01,2,,6jxj例2321436minxxxz123s.t.120xxx301x5002x203x解:编写M文件xxgh2.m如下：c=[634];A=[010];b=[50];Aeq=[111];beq=[120];vlb=[30,0,20];vub=[];[x,fval]=linprog(c,A,b,Aeq,beq,vlb,vub)123min(634)xzxx32120030xxx123111120s.t.01050xxxs.t.Xz8121110913min9008003.12.15.000000011.14.0X50060040010010001001...