数学模型姜启源.ppt

《数学模型姜启源.ppt》由会员分享，可在线阅读，更多相关《数学模型姜启源.ppt（259页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、数数 学学 模模 型型 数学模型数学模型数学模型数学模型 姜启源姜启源 主编主编 课 程 名 称学 时 36数学模型与数学建模Mathematical Modeling学 分3课程类别专业选修课先 修 课 程微积分、线性代数、概率论与数理统计课 程 简 介本本课课程是程是计计算机及管理算机及管理专业专业的一的一门专业选门专业选修修课课。也是本科生参加数学建。也是本科生参加数学建模模竞赛竞赛的的辅导课辅导课程。数学模型是架于数学理程。数学模型是架于数学理论论和和实际问题实际问题之之间间的的桥桥梁。梁。数学建模是数学建模是应应用数学解决用数学解决实际问题实际问题的重要手段和途径。本的重要手段和途径

2、。本书书介介绍绍数学建数学建模中常用的一些基本概念、理模中常用的一些基本概念、理论论和典型的数学模型，包括：数据和典型的数学模型，包括：数据拟拟合，合，网网络络模型，模型，优优化模型，离散模型、随机模型，化模型，离散模型、随机模型，时间时间序列序列预报预报模型，回模型，回归归分析及其分析及其试验设计试验设计。通。通过过数学模型和数学建模有关数学模型和数学建模有关问题问题的的论论述和模型述和模型实实例的介例的介绍绍，使学生，使学生应应用数学解决用数学解决实际问题实际问题的能力有所提高。的能力有所提高。教 材 及 参 考 书 目数学模型，姜启源主编，高等教育出版社 课课 程程 简简 介介 数学模型

3、数学模型数学模型数学模型 姜启源姜启源 主编主编 数学模型数学模型数学模型数学模型 姜启源姜启源 主编主编第一章第一章 建立数学模型建立数学模型第二章第二章 初等模型初等模型第三章第三章 简单的优化模型简单的优化模型第四章第四章 数学规划模型数学规划模型 第五章第五章 微分方程模型微分方程模型第六章第六章 稳定性模型稳定性模型第七章第七章 差分方程模型差分方程模型第八章第八章 离散模型离散模型第九章第九章 概率模型概率模型第十章第十章 统计回归模型统计回归模型附录附录:数学建模实验数学建模实验周次节次教学内容课时作业执行情况1五 56 1.1-1.5数学模型的介绍 1.6数学模型的基本方法步骤

4、、特点和分类22五 562.1公平的席位分配(讨论课)2.2录像机计数器的用途2.3双层玻璃的功效23五 562.7实物交换3.2生猪的出售时机24五 563.3森林救火(讨论课)3.4最优价格25五 563.6消费者的选择4.3汽车生产与原油采购26五 564.5饮料厂的生产与检修5.1传染病模型(讨论课)27五 56 5.2经济增长模型 5.6人口的预测和控制 28五 56 6.1捕鱼业的持续收获 6.2军备竞赛(讨论课)2 教教 学学 进进 度度 数学模型数学模型数学模型数学模型 姜启源姜启源 主编主编9五 56 6.4种群的相互依存 7.1市场经济中的蛛网模型2 10五 56 7.2减

5、肥计划-节食与运动 8.3层次分析模型212五 56 8.4效益的合理分配 9.2报童的诀窍(讨论课)213五 56 9.5随机人口模型 9.6航空公司的预定票策略214五五56 10.1牙膏的牙膏的销销售量售量2评估周评估周15五 56 Mtlab,Mathematcia数学软件学习(上机)216五 56 数学建模实验(上机)217五 56 数学建模实验(上机)218 考试 数学模型数学模型数学模型数学模型 姜启源姜启源 主编主编 数学模型数学模型数学模型数学模型 姜启源姜启源 主编主编第第一一章章 建立数学模型建立数学模型1.1从现实对象到数学模型从现实对象到数学模型1.2数学建模的重要意

6、义数学建模的重要意义1.3数学建模示例数学建模示例1.4数学建模的方法和步骤数学建模的方法和步骤1.5数学模型的特点和分类数学模型的特点和分类1.6怎样学习数学建模怎样学习数学建模玩具、照片、飞机、火箭模型玩具、照片、飞机、火箭模型 实物模型实物模型水箱中的舰艇、风洞中的飞机水箱中的舰艇、风洞中的飞机 物理模型物理模型地图、电路图、分子结构图地图、电路图、分子结构图符号模型符号模型模型模型是为了一定目的，对客观事物的一部分是为了一定目的，对客观事物的一部分进行简缩、抽象、提炼出来的进行简缩、抽象、提炼出来的原型原型的替代物的替代物模型模型集中反映了集中反映了原型原型中人们需要的那一部分特征中人

7、们需要的那一部分特征1.1 从现实对象到数学模型从现实对象到数学模型我们常见的模型我们常见的模型 第一章第一章 建立数学模型建立数学模型数学模型数学模型 姜启源姜启源 主编主编你碰到过的数学模型你碰到过的数学模型“航行问题航行问题”用用 x 表示船速，表示船速，y 表示水速，列出方程：表示水速，列出方程：答：船速每小时答：船速每小时20千米千米/小时小时.甲乙两地相距甲乙两地相距750千米，船从甲到乙顺水航行需千米，船从甲到乙顺水航行需30小时，小时，从乙到甲逆水航行需从乙到甲逆水航行需50小时，问船的速度是多少小时，问船的速度是多少?x=20y=5求解求解 第一章第一章 建立数学模型建立数学

8、模型数学模型数学模型 姜启源姜启源 主编主编航行问题航行问题建立数学模型的基本步骤建立数学模型的基本步骤作出简化假设（船速、水速为常数）；作出简化假设（船速、水速为常数）；用符号表示有关量（用符号表示有关量（x,y表示船速和水速）；表示船速和水速）；用物理定律（匀速运动的距离等于速度乘以用物理定律（匀速运动的距离等于速度乘以时间）列出数学式子（二元一次方程）；时间）列出数学式子（二元一次方程）；求解得到数学解答（求解得到数学解答（x=20,y=5）；）；回答原问题（船速每小时回答原问题（船速每小时20千米千米/小时）。小时）。第一章第一章 建立数学模型建立数学模型数学模型数学模型 姜启源姜启源

9、 主编主编数学模型数学模型(MathematicalModel)和和数学建模（数学建模（MathematicalModeling)对于一个对于一个现实对象现实对象，为了一个，为了一个特定目的特定目的，根据其根据其内在规律内在规律，作出必要的，作出必要的简化假设简化假设，运用适当的运用适当的数学工具数学工具，得到的一个，得到的一个数学结构数学结构。建立数学模型的全过程建立数学模型的全过程（包括表述、求解、解释、检验等）（包括表述、求解、解释、检验等）数学模型数学模型数学数学建模建模 第一章第一章 建立数学模型建立数学模型数学模型数学模型 姜启源姜启源 主编主编1.2 数学建模的重要意义数学建模的

10、重要意义电子计算机的出现及飞速发展；电子计算机的出现及飞速发展；数学以空前的广度和深度向一切领域渗透。数学以空前的广度和深度向一切领域渗透。数学建模作为用数学方法解决实际问题的第一步，数学建模作为用数学方法解决实际问题的第一步，越来越受到人们的重视。越来越受到人们的重视。在一般工程技术领域数学建模仍然大有用武之地；在一般工程技术领域数学建模仍然大有用武之地；在高新技术领域数学建模几乎是必不可少的工具；在高新技术领域数学建模几乎是必不可少的工具；数学进入一些新领域，为数学建模开辟了许多处女地。数学进入一些新领域，为数学建模开辟了许多处女地。第一章第一章 建立数学模型建立数学模型数学模型数学模型

11、姜启源姜启源 主编主编数学建模的具体应用数学建模的具体应用分析与设计分析与设计预报与决策预报与决策控制与优化控制与优化规划与管理规划与管理数学建模计算机技术知识经济知识经济如虎添翼如虎添翼 第一章第一章 建立数学模型建立数学模型数学模型数学模型 姜启源姜启源 主编主编1.3 数学建模示例数学建模示例1.3.1 椅子能在不平的地面上放稳吗椅子能在不平的地面上放稳吗问题分析问题分析模模型型假假设设通常通常三只脚着地三只脚着地放稳放稳四只脚着地四只脚着地四条腿一样长，椅脚与地面点接触，四脚四条腿一样长，椅脚与地面点接触，四脚连线呈正方形连线呈正方形;地面高度连续变化，可视为数学上的连续地面高度连续变

12、化，可视为数学上的连续曲面曲面;地面相对平坦，使椅子在任意位置至少三地面相对平坦，使椅子在任意位置至少三只脚同时着地。只脚同时着地。第一章第一章 建立数学模型建立数学模型数学模型数学模型 姜启源姜启源 主编主编模型构成模型构成用数学语言把椅子位置和四只脚着地的关系表示出来用数学语言把椅子位置和四只脚着地的关系表示出来椅子位置椅子位置利用正方形利用正方形(椅脚连线椅脚连线)的对称性的对称性xBADCODCBA用用(对角线与对角线与x轴的夹角轴的夹角)表示椅子位置表示椅子位置四只脚着地四只脚着地距离是距离是 的函数的函数四个距离四个距离(四只脚四只脚)A,C两脚与地面距离之和两脚与地面距离之和f(

13、)B,D两脚与地面距离之和两脚与地面距离之和g()两个距离两个距离 椅脚与地面距离为零椅脚与地面距离为零正方形正方形ABCD绕绕O点旋转点旋转正方形正方形对称性对称性 第一章第一章 建立数学模型建立数学模型数学模型数学模型 姜启源姜启源 主编主编用数学语言把椅子位置和四只脚着地的关系表示出来用数学语言把椅子位置和四只脚着地的关系表示出来f(),g()是是连续函连续函数数对任意对任意,f(),g()至少一个为至少一个为0数学数学问题问题已知：已知：f(),g()是是连续函数连续函数;对任意对任意，f()g()=0;且且g(0)=0，f(0)0.证明：存在证明：存在 0，使，使f(0)=g(0)=

14、0.模型构成模型构成地面为连续曲面地面为连续曲面椅子在任意位置椅子在任意位置至少三只脚着地至少三只脚着地 第一章第一章 建立数学模型建立数学模型数学模型数学模型 姜启源姜启源 主编主编模型求解模型求解给出一种简单、粗糙的证明方法给出一种简单、粗糙的证明方法将椅子将椅子旋转旋转900，对角线，对角线AC和和BD互换。互换。由由g(0)=0，f(0)0，知，知f(/2)=0,g(/2)0.令令h()=f()g(),则则h(0)0和和h(/2)p2/n2，对对不公平不公平A p1/n1p2/n2=5 第二章第二章 初等模型初等模型数学模型数学模型 姜启源姜启源 主编主编公平分配方案应公平分配方案应使

15、使 rA,rB 尽量小尽量小设设A,B已分别有已分别有n1,n2席，若增加席，若增加1席，问应分给席，问应分给A,还是还是B不妨设分配开始时不妨设分配开始时p1/n1p2/n2，即对即对A不公平不公平对对A的相对不公平度的相对不公平度将绝对度量改为相对度量将绝对度量改为相对度量类似地定义类似地定义rB(n1,n2)将一次性的席位分配转化为动态的席位分配将一次性的席位分配转化为动态的席位分配,即即“公平公平”分配方法分配方法若若p1/n1p2/n2，定义定义 第二章第二章 初等模型初等模型数学模型数学模型 姜启源姜启源 主编主编1）若）若p1/(n1+1)p2/n2，则这席应给则这席应给A2）若

16、）若p1/(n1+1)p2/(n2+1)，应计算应计算rB(n1+1,n2)应计算应计算rA(n1,n2+1)若若rB(n1+1,n2)p2/n2问：问：p1/n1rA(n1,n2+1),则这席应给则这席应给B 第二章第二章 初等模型初等模型数学模型数学模型 姜启源姜启源 主编主编当当rB(n1+1,n2)640g=0.1 第三章第三章 简单的优化模型简单的优化模型数学模型数学模型 姜启源姜启源 主编主编敏感性分析敏感性分析研究研究r,g变化时对模型结果的影响变化时对模型结果的影响估计估计r=2，g=0.1设设g=0.1不变不变t 对对r 的（相对）敏感度的（相对）敏感度生猪每天体重增加量生猪

17、每天体重增加量r 增加增加1%，出售时间推迟，出售时间推迟3%。rt 第三章第三章 简单的优化模型简单的优化模型数学模型数学模型 姜启源姜启源 主编主编敏感性分析敏感性分析估计估计r=2，g=0.1研究研究r,g变化时对模型结果的影响变化时对模型结果的影响设设r=2不变不变t 对对g的（相对）敏感度的（相对）敏感度生猪价格每天的降低量生猪价格每天的降低量g增加增加1%，出售时间提前，出售时间提前3%。gt 第三章第三章 简单的优化模型简单的优化模型数学模型数学模型 姜启源姜启源 主编主编强健性分析强健性分析保留生猪直到利润的增值等于每天的费用时出售保留生猪直到利润的增值等于每天的费用时出售由由

18、S(t,r)=3建议过一周后建议过一周后(t=7)重新估计重新估计,再作计算。再作计算。研究研究r,g不是常数时对模型结果的影响不是常数时对模型结果的影响w=80+rt w=w(t)p=8-gt p=p(t)若若(10%),则则（30%）每天利润的增值每天利润的增值每天投入的资金每天投入的资金 第三章第三章 简单的优化模型简单的优化模型数学模型数学模型 姜启源姜启源 主编主编3.3 森林救火森林救火森林失火后，要确定派出消防队员的数量。森林失火后，要确定派出消防队员的数量。队员多，森林损失小，救援费用大；队员多，森林损失小，救援费用大；队员少，森林损失大，救援费用小。队员少，森林损失大，救援费

19、用小。综合考虑损失费和救援费，确定队员数量。综合考虑损失费和救援费，确定队员数量。问题问题分析分析问题问题记队员人数记队员人数x,失火时刻失火时刻t=0,开始救火时刻开始救火时刻t1,灭火时刻灭火时刻t2,时刻时刻t森林烧毁面积森林烧毁面积B(t).损失费损失费f1(x)是是x的减函数的减函数,由烧毁面积由烧毁面积B(t2)决定决定.救援费救援费f2(x)是是x的增函数的增函数,由队员人数和救火时间由队员人数和救火时间决定决定.存在恰当的存在恰当的x，使，使f1(x),f2(x)之和最小之和最小 第三章第三章 简单的优化模型简单的优化模型数学模型数学模型 姜启源姜启源 主编主编关键是对关键是对

20、B(t)作出作出合理的简化假设合理的简化假设.问题问题分析分析失火时刻失火时刻t=0,开始救火时刻开始救火时刻t1,灭火时刻灭火时刻t2,画出时刻画出时刻t 森林烧毁面积森林烧毁面积B(t)的大致图形的大致图形t1t20tBB(t2)分析分析B(t)比较困难比较困难,转而讨论森林烧毁转而讨论森林烧毁速度速度dB/dt.第三章第三章 简单的优化模型简单的优化模型数学模型数学模型 姜启源姜启源 主编主编模型假设模型假设3）f1(x)与与B(t2)成正比，系数成正比，系数c1(烧毁单位面积损失费）烧毁单位面积损失费）1）0 t t1,dB/dt与与t成正比，系数成正比，系数(火势蔓延速度）火势蔓延速

21、度）2）t1 t t2,降为降为-x(为队员的平均灭火为队员的平均灭火速度）速度）4）每个）每个队员的单位时间灭火费用队员的单位时间灭火费用c2,一次性费用一次性费用c3假设假设1）的解释的解释 rB火势以失火点为中心，火势以失火点为中心，均匀向四周呈圆形蔓延，均匀向四周呈圆形蔓延，半径半径r与与t 成正比成正比面积面积B与与t2成正比，成正比，dB/dt与与t成正比成正比.第三章第三章 简单的优化模型简单的优化模型数学模型数学模型 姜启源姜启源 主编主编模型建立模型建立b0t1tt2假设假设1）目标函数目标函数总费用总费用假设假设3）4）假设假设2）第三章第三章 简单的优化模型简单的优化模型

22、数学模型数学模型 姜启源姜启源 主编主编模型建立模型建立目标函数目标函数总费用总费用模型求解模型求解求求x使使C(x)最小最小结果解释结果解释 /是火势不继续蔓延的最少队员数是火势不继续蔓延的最少队员数b0t1t2t其中其中c1,c2,c3,t1,为已知参数为已知参数 第三章第三章 简单的优化模型简单的优化模型数学模型数学模型 姜启源姜启源 主编主编模型模型应用应用c1,c2,c3已知已知,t1可估计可估计,c2 x c1,t1,x c3,x 结果结果解释解释c1烧毁单位面积损失费烧毁单位面积损失费,c2每个每个队员单位时间灭火费队员单位时间灭火费,c3每个每个队员一次性费用队员一次性费用,t

23、1开始救火时刻开始救火时刻,火火势蔓延速度势蔓延速度,每个每个队员平均灭火队员平均灭火速度速度.为什么为什么?,可可设置一系列数设置一系列数值值由模型决定队员数量由模型决定队员数量x 第三章第三章 简单的优化模型简单的优化模型数学模型数学模型 姜启源姜启源 主编主编3.4 最优价格最优价格问题问题根据产品成本和市场需求，在产销平根据产品成本和市场需求，在产销平衡条件下确定商品价格，使利润最大衡条件下确定商品价格，使利润最大假设假设1）产量等于销量，记作）产量等于销量，记作x2）收入与销量）收入与销量x 成正比，系数成正比，系数p 即价格即价格3）支出与产量）支出与产量x 成正比，系数成正比，系

24、数q 即成本即成本4）销量）销量x 依赖于价格依赖于价格p,x(p)是减函数是减函数建模建模与求解与求解收入收入支出支出利润利润进一步设进一步设求求p使使U(p)最大最大 第三章第三章 简单的优化模型简单的优化模型数学模型数学模型 姜启源姜启源 主编主编使利润使利润U(p)最大的最优价格最大的最优价格p*满足满足最大利润在边际收入等于边际支出时达到最大利润在边际收入等于边际支出时达到建模建模与求解与求解边际收入边际收入边际支出边际支出 第三章第三章 简单的优化模型简单的优化模型数学模型数学模型 姜启源姜启源 主编主编结果结果解释解释q/2成本的一半成本的一半b 价格上升价格上升1单位时销量的下

25、降单位时销量的下降幅度（需求对价格的敏感度）幅度（需求对价格的敏感度）a 绝对需求绝对需求(p很小时的需求很小时的需求)b p*a p*思考：如何得到参数思考：如何得到参数a,b?第三章第三章 简单的优化模型简单的优化模型数学模型数学模型 姜启源姜启源 主编主编q2U(q1,q2)=cq103.6 消费者均衡消费者均衡问题问题消费者对甲乙两种商品的偏爱程度用无差别消费者对甲乙两种商品的偏爱程度用无差别曲线族表示，问他如何分配一定数量的钱，曲线族表示，问他如何分配一定数量的钱，购买这两种商品，以达到最大的满意度。购买这两种商品，以达到最大的满意度。设甲乙数量为设甲乙数量为q1,q2,消消费者的无

26、差别曲线族费者的无差别曲线族(单调减、下凸、不相单调减、下凸、不相交），记作交），记作U(q1,q2)=cU(q1,q2)效用函数效用函数已知甲乙价格已知甲乙价格p1,p2,有钱有钱s，试分配，试分配s,购买甲乙数量购买甲乙数量q1,q2,使使U(q1,q2)最大最大.第三章第三章 简单的优化模型简单的优化模型数学模型数学模型 姜启源姜启源 主编主编s/p2s/p1q2U(q1,q2)=cq10模型模型及及求解求解已知价格已知价格p1,p2,钱钱s,求求q1,q2,或或p1q1/p2q2,使使U(q1,q2)最大最大几几何何解解释释直线直线MN:最优解最优解Q:MN与与l2切点切点斜率斜率MQ

27、N第三章第三章 简单的优化模型简单的优化模型数学模型数学模型 姜启源姜启源 主编主编结果结果解释解释边际效用边际效用消费者均衡状态在两种商品消费者均衡状态在两种商品的边际效用之比恰等于它们的边际效用之比恰等于它们价格之比时达到。价格之比时达到。效用函数效用函数U(q1,q2)应满足的条件应满足的条件A.U(q1,q2)=c所确定的函数所确定的函数q2=q2(q1)单调减、下凸单调减、下凸解释解释B的实际意义的实际意义第三章第三章 简单的优化模型简单的优化模型数学模型数学模型 姜启源姜启源 主编主编效用函数效用函数U(q1,q2)几种常用几种常用的形式的形式消费者均衡状态下购买两种商品费用之比消

28、费者均衡状态下购买两种商品费用之比与二者价格之比的平方根成正比。与二者价格之比的平方根成正比。U(q1,q2)中参数中参数,分别表示消费者对甲乙分别表示消费者对甲乙两种商品的偏爱程度。两种商品的偏爱程度。第三章第三章 简单的优化模型简单的优化模型数学模型数学模型 姜启源姜启源 主编主编购买两种商品费用之比与二者价格无关。购买两种商品费用之比与二者价格无关。U(q1,q2)中参数中参数,分别表示对甲乙分别表示对甲乙的偏爱程度。的偏爱程度。思考：如何推广到思考：如何推广到m(2)种商品的情况种商品的情况效用函数效用函数U(q1,q2)几种常用几种常用的形式的形式第三章第三章 简单的优化模型简单的优

29、化模型数学模型数学模型 姜启源姜启源 主编主编第四章第四章 数学规划模型数学规划模型 4.3 汽车生产与原油采购汽车生产与原油采购4.5 饮料厂的生产与检修饮料厂的生产与检修 数学模型数学模型数学模型数学模型 姜启源姜启源 主编主编数学规划模型数学规划模型 实际问题中实际问题中的优化模型的优化模型x决策变量决策变量f(x)目标函数目标函数gi(x)0约束条约束条件件多元函数多元函数条件极值条件极值决策变量个数决策变量个数n和和约束条件个数约束条件个数m较大较大最优解在可行域最优解在可行域的边界上取得的边界上取得数数学学规规划划线性规划线性规划非线性规划非线性规划整数规划整数规划重点在模型的建立

30、和结果的分析重点在模型的建立和结果的分析第四章第四章 数学规划模型数学规划模型数学模型数学模型 姜启源姜启源 主编主编 如果生产某一类型汽车，则至少要生产如果生产某一类型汽车，则至少要生产8080辆，辆，那么最优的生产计划应作何改变？那么最优的生产计划应作何改变？例例1 汽车厂生产计划汽车厂生产计划 汽车厂生产三种类型的汽车，已知各类型每辆车对钢汽车厂生产三种类型的汽车，已知各类型每辆车对钢材、劳动时间的需求，利润及工厂每月的现有量。材、劳动时间的需求，利润及工厂每月的现有量。小型小型中型中型大型大型现有量现有量钢材（吨）钢材（吨）1.535600劳动时间（小时）劳动时间（小时）2802504

31、0060000利润（万元）利润（万元）234制订月生产计划，使工厂的利润最大。制订月生产计划，使工厂的利润最大。4.3 汽车生产与原油采购汽车生产与原油采购第四章第四章 数学规划模型数学规划模型数学模型数学模型 姜启源姜启源 主编主编设每月生产小、中、大型设每月生产小、中、大型汽车的数量分别为汽车的数量分别为x1,x2,x3汽车厂生产计划汽车厂生产计划 模型建立模型建立 小型小型中型中型大型大型现有量现有量钢材钢材1.535600时间时间28025040060000利润利润234线性线性规划规划模型模型(LP)第四章第四章 数学规划模型数学规划模型数学模型数学模型 姜启源姜启源 主编主编模型模

32、型求解求解 3）模型中增加条件：模型中增加条件：x1,x2,x3均为整数，重新求解。均为整数，重新求解。OBJECTIVEFUNCTIONVALUE1)632.2581VARIABLEVALUEREDUCEDCOSTX164.5161290.000000X2167.7419280.000000X30.0000000.946237ROWSLACKORSURPLUSDUALPRICES2)0.0000000.7311833)0.0000000.003226结果为小数，结果为小数，怎么办？怎么办？1）舍去小数：取）舍去小数：取x1=64，x2=167，算出目标函数值，算出目标函数值z=629，与，与

33、LP最优值最优值632.2581相差不大。相差不大。2）试试探探：如如取取x1=65，x2=167；x1=64，x2=168等等，计计算算函函数数值值z，通过比较可能得到更优的解。，通过比较可能得到更优的解。但必须检验它们是否满足约束条件。为什么？但必须检验它们是否满足约束条件。为什么？第四章第四章 数学规划模型数学规划模型数学模型数学模型 姜启源姜启源 主编主编IP可用可用LINDO直接求解直接求解整数规划整数规划(IntegerProgramming,简记简记IP)“gin3”表示表示“前前3个变量个变量为整数为整数”，等价于：，等价于：ginx1ginx2ginx3IP的最优解的最优解x

34、1=64，x2=168，x3=0，最优值，最优值z=632max2x1+3x2+4x3st1.5x1+3x2+5x3600280 x1+250 x2+400 x360000endgin3OBJECTIVEFUNCTIONVALUE1)632.0000VARIABLEVALUEREDUCEDCOSTX164.000000-2.000000X2168.000000-3.000000X30.000000-4.000000模型求解模型求解 IP结果输出结果输出第四章第四章 数学规划模型数学规划模型数学模型数学模型 姜启源姜启源 主编主编其中其中3个个子模型应子模型应去掉，然后去掉，然后逐一求解，比较目

35、标函数值，逐一求解，比较目标函数值，再加上整数约束，得最优解：再加上整数约束，得最优解：方法方法1：分解为：分解为8个个LP子模型子模型汽车厂生产计划汽车厂生产计划 若生产某类汽车，则至少生产若生产某类汽车，则至少生产8080辆，求生产计划。辆，求生产计划。x1,x2,x3=0或或 80 x1=80，x2=150，x3=0，最优值，最优值z=610第四章第四章 数学规划模型数学规划模型数学模型数学模型 姜启源姜启源 主编主编LINDO中中 对对 0-1变量的限定：变量的限定：inty1inty2inty3方法方法2：引入引入0-1变量，化为整数规划变量，化为整数规划M为大的正数，为大的正数，可

36、取可取1000OBJECTIVEFUNCTIONVALUE1)610.0000VARIABLEVALUEREDUCEDCOSTX180.000000-2.000000X2150.000000-3.000000X30.000000-4.000000Y11.0000000.000000Y21.0000000.000000Y30.0000000.000000 若生产某类汽车，则至少生产若生产某类汽车，则至少生产8080辆，求生产计划。辆，求生产计划。x1=0 或 80 x2=0 或 80 x3=0 或 80最优解同前最优解同前第四章第四章 数学规划模型数学规划模型数学模型数学模型 姜启源姜启源 主编

37、主编NLP虽虽然然可可用用现现成成的的数数学学软软件件求求解解(如如LINGO,MATLAB)，但是其结果常依赖于初值的选择。，但是其结果常依赖于初值的选择。方法方法3：化为非线性规划化为非线性规划非线性规划（非线性规划（Non-LinearProgramming，简记，简记NLP）实实践践表表明明，本本例例仅仅当当初初值值非非常常接接近近上上面面方方法法算算出出的最优解时，才能得到正确的结果。的最优解时，才能得到正确的结果。若生产某类汽车，则至少生产若生产某类汽车，则至少生产8080辆，求生产计划。辆，求生产计划。x1=0 或 80 x2=0 或 80 x3=0 或 80第四章第四章 数学规

38、划模型数学规划模型数学模型数学模型 姜启源姜启源 主编主编应如何安排原油的采购和加工应如何安排原油的采购和加工？例例2 原油采购与加工原油采购与加工 市场上可买到不超过市场上可买到不超过1500吨的原油吨的原油A：购买量不超过购买量不超过500吨时的单价为吨时的单价为10000元元/吨；吨；购买量超过购买量超过500吨但不超过吨但不超过1000吨时，超过吨时，超过500吨的吨的 部分部分8000元元/吨；吨；购买量超过购买量超过1000吨时，超过吨时，超过1000吨的部分吨的部分6000元元/吨。吨。售价售价4800元元/吨吨售价售价5600元元/吨吨库存库存500吨吨库存库存1000吨吨汽油

39、甲汽油甲(A 50%)原油原油A原油原油B汽油乙汽油乙(A 60%)第四章第四章 数学规划模型数学规划模型数学模型数学模型 姜启源姜启源 主编主编决策决策变量变量 目标目标函数函数问题问题分析分析 利润：销售汽油的收入利润：销售汽油的收入 -购买原油购买原油A的支出的支出 难点：原油难点：原油A的购价与购买量的关系较复杂的购价与购买量的关系较复杂甲甲(A 50%)AB乙乙(A 60%)购买购买xx11x12x21x224.8千元千元/吨吨5.6千元千元/吨吨原油原油A的购买量的购买量,原油原油A,B生产生产汽油汽油甲甲,乙的数量乙的数量c(x)购买原油购买原油A的支出的支出利润利润(千元千元)

40、c(x)如何表述？如何表述？第四章第四章 数学规划模型数学规划模型数学模型数学模型 姜启源姜启源 主编主编原油供应原油供应 约束约束条件条件 x 500吨单价为吨单价为10千千元元/吨；吨；500吨吨 x 1000吨，超过吨，超过500吨的吨的8千千元元/吨；吨；1000吨吨 x 1500吨，超过吨，超过1000吨的吨的6千千元元/吨。吨。目标目标函数函数购买购买x ABx11x12x21x22库存库存500吨吨库存库存1000吨吨第四章第四章 数学规划模型数学规划模型数学模型数学模型 姜启源姜启源 主编主编目标函数中目标函数中c(x)不是线性函数，是非线性规划；不是线性函数，是非线性规划；对

41、于用分段函数定义的对于用分段函数定义的c(x)，一般的非线性规划软一般的非线性规划软件也难以输入和求解；件也难以输入和求解；想办法将模型化简，用现成的软件求解。想办法将模型化简，用现成的软件求解。汽油含原油汽油含原油A的比例限制的比例限制约束约束条件条件甲甲(A 50%)AB乙乙(A 60%)x11x12x21x22第四章第四章 数学规划模型数学规划模型数学模型数学模型 姜启源姜启源 主编主编x1,x2,x3以价格以价格10,8,6(千元千元/吨吨)采购采购A的吨数的吨数目标目标函数函数 只有当以只有当以10千元千元/吨的价格购买吨的价格购买x1=500(吨吨)时，才能以时，才能以8千元千元/

42、吨的价格购买吨的价格购买x2方法方法1 非线性规划模型非线性规划模型，可以用，可以用LINGO求解求解模型求解模型求解x=x1+x2+x3,c(x)=10 x1+8x2+6x3 500吨吨 x 1000吨，超过吨，超过500吨的吨的8千千元元/吨吨增加约束增加约束x=x1+x2+x3,c(x)=10 x1+8x2+6x3 第四章第四章 数学规划模型数学规划模型数学模型数学模型 姜启源姜启源 主编主编方法方法1：LINGO求解求解Model:Max=4.8*x11+4.8*x21+5.6*x12+5.6*x22-10*x1-8*x2-6*x3;x11+x12x+500;x21+x220;2*x1

43、2-3*x220;x=x1+x2+x3;(x1-500)*x2=0;(x2-500)*x3=0;x1500;x2500;x30;x110;x120;x210;x220;x10;x20;x30;endObjectivevalue:4800.000VariableValueReducedCostX11500.00000.0000000E+00X21500.00000.0000000E+00X120.0000000E+000.0000000E+00X220.0000000E+000.0000000E+00X10.1021405E-1310.00000X20.0000000E+008.000000X3

44、0.0000000E+006.000000X0.0000000E+000.0000000E+00LINGO得到的是局部最优解，还得到的是局部最优解，还能得到更好的解吗？能得到更好的解吗？用库存的用库存的500吨原油吨原油A、500吨原油吨原油B生产汽油甲，不购买新的原油生产汽油甲，不购买新的原油A，利润为利润为4,800千千元。元。第四章第四章 数学规划模型数学规划模型数学模型数学模型 姜启源姜启源 主编主编y1,y2,y3=1以价格以价格10,8,6(千元千元/吨吨)采购采购A增增加加约约束束方法方法2 0-1线性规划模型线性规划模型，可用，可用LINDO求解求解y1,y2,y3=0或或1O

45、BJECTIVEFUNCTIONVALUE1)5000.000VARIABLE VALUE REDUCEDCOSTY11.0000000.000000Y21.0000002200.000000Y31.0000001200.000000X110.0000000.800000X210.0000000.800000X121500.0000000.000000X221000.0000000.000000X1500.0000000.000000X2500.0000000.000000X30.0000000.400000X1000.0000000.000000购买购买1000吨原油吨原油A，与，与库存的库

46、存的500吨原油吨原油A和和1000吨原油吨原油B一起，生一起，生产汽油乙，利润为产汽油乙，利润为5,000千元千元。x1,x2,x3以价格以价格10,8,6(千元千元/吨吨)采购采购A的吨数的吨数y=0 x=0 x0y=1优于方法优于方法1的结果的结果第四章第四章 数学规划模型数学规划模型数学模型数学模型 姜启源姜启源 主编主编b1b2b3b4方法方法3 b1 x b2，x=z1b1+z2b2，z1+z2=1，z1,z2 0,c(x)=z1c(b1)+z2c(b2).c(x)x1200090005000050010001500b2 x b3，x=z2b2+z3b3，z2+z3=1，z2,z3

47、 0,c(x)=z2c(b2)+z3c(b3).b3 x b4，x=z3b3+z4b4，z3+z4=1，z3,z4 0,c(x)=z3c(b3)+z4c(b4).直接处理处理分段线性函数直接处理处理分段线性函数c(x)第四章第四章 数学规划模型数学规划模型数学模型数学模型 姜启源姜启源 主编主编IP模型，模型，LINDO求求解，得到的结果与解，得到的结果与方法方法2相同相同.处理分段线性函数，方法处理分段线性函数，方法3更具一般性更具一般性bk x bk+1yk=1,否则否则,yk=0方法方法3 bk x bk+1,x=zkbk+z k+1bk+1zk+zk+1=1，zk,zk+1 0,c(x

48、)=zkc(bk)+zk+1c(bk+1).c(x)x1200090005000050010001500b1b2b3b4对于对于k=1,2,3第四章第四章 数学规划模型数学规划模型数学模型数学模型 姜启源姜启源 主编主编4.5 饮料厂的生产与检修饮料厂的生产与检修单阶段生产计划单阶段生产计划多阶段生产计划多阶段生产计划生产批量问题生产批量问题企业生产计划企业生产计划考虑与产量无关的固定费用考虑与产量无关的固定费用给优化模型求解带来新的困难给优化模型求解带来新的困难外部需求和内部外部需求和内部资源随时间变化资源随时间变化第四章第四章 数学规划模型数学规划模型数学模型数学模型 姜启源姜启源 主编主

49、编 安排生产计划安排生产计划,满足每周的需求满足每周的需求,使使4周总费用最小。周总费用最小。存贮费存贮费:每周每千箱饮料每周每千箱饮料 0.2千元。千元。例例1 饮料厂的生产与检修计划饮料厂的生产与检修计划 在在4周内安排一次设备检修，占用当周周内安排一次设备检修，占用当周15千箱生产能千箱生产能力，能使检修后每周增产力，能使检修后每周增产5千箱，检修应排在哪一周千箱，检修应排在哪一周?周次周次需求量需求量(千箱千箱)生产能力生产能力(千箱千箱)成本成本(千元千元/千箱千箱)115305.0225405.1335455.4425205.5合计合计100135某种饮料某种饮料4周的需求量、生产

50、能力和成本周的需求量、生产能力和成本第四章第四章 数学规划模型数学规划模型数学模型数学模型 姜启源姜启源 主编主编问题分析问题分析除第除第4周外每周的生产周外每周的生产能力超过每周的需求；能力超过每周的需求；生产成本逐周上升；生产成本逐周上升；前几周应多生产一些。前几周应多生产一些。周次周次需求需求能力能力11530225403354542520合计合计100135成本成本5.05.15.45.5 饮料厂在第饮料厂在第1周开始时没有库存；周开始时没有库存；从费用最小考虑从费用最小考虑,第第4周末不能有库存；周末不能有库存；周末有库存时需支出一周的存贮费；周末有库存时需支出一周的存贮费；每周末的