数学建模之规划问题
数学建模目标规划方法
30
x1
2x1
12x2 x2
d1 d2
d1 d2
2500 140
x1
d
3
d3
60
a x (,)b
ij j
i
j 1
(i 1,2, , m)
绝对约束
x 0 ( j 1,2, , n) j
d , d 0 (l 1,2, , L) ll
非负约束
K
L
min Z
pk
(kl
d
l
kl
dl
)
k 1
l 1
n
c(l) x d d g ( l 1,2, , L)
三 目标规划方法
通过前面的介绍和讨论，我们知道，目标规划方法 是解决多目标规划问题的重要技术之一。
这一方法是美国学者查恩斯（A.Charnes）和库 伯（W.W.Cooper）于1961年在线性规划的基础上提 出来的。后来，查斯基莱恩（U.Jaashelainen）和李 （Sang.Lee）等人，进一步给出了求解目标规划问题 的一般性方法——单纯形方法。
34
4
所以目标规划模型为：
min Z p d p (7d 12d ) p (d d )
11
2
2
3
34
4
70x 120x d d 50000
1
2
1
1
x 1
d d 200
2
2
x d d 250
生产甲、乙两种产品，
数学建模中规划问题的MATLAB求解
min P ( x, M )
的最优解 x 也是原问题的最优解。
【例 2-10】求非线性规划问题的解
min f ( x) x1 x2 8
2 x1 x2 0 2 s.t. x1 x2 2 0 x , x 0 1 2
2 2
function g=test(x) M=50000; f=x(1)^2+x(2)^2+8; g=f-M*min(x(1),0)-M*min(x(2),0)-M*min(x(1)^2-x(2),0)... +M*abs(-x(1)-x(2)^2+2);
命令窗输入[x,y]=fminunc('test',rand(2,1))
第二章
规划问题的MATLAB求解
2.1 线性规划 2.2 非线性规划
2.3 整数规划
数学规划模型
序：
规划问题是常见的数学建模问题，离散系统的优化问题 一般都可以通过规划模型来求解。因此快速求解规划问题是 数学建模的基本素质。利用MATLAB提供的规划模型求解命 令，可以快速得到想要的结果。
2.1 线性规划
以下介绍外罚函数法，内罚函数法参考PPT“内罚函数 法”。
设规划模型为 min f ( x)
gi ( x ) 0, i 1 ~ r s.t. hi ( x ) 0, i 1 ~ s k ( x ) 0, i 1 ~ t i
取充分大的数 M 0 ，构造函数 P ( x , M ) f ( x ) M max( g i ( x ), 0) M min( hi ( x ), 0)
数学建模测试题-线性规划部分
313数学教育1、2班，510数学教育1、2、3班数学建模上机测试题，需要把运行结果写出来。
模型包括目标函数、约束条件，编写的程序和程序运行结果四部分内容。
写在作业本上。
按学号顺序做，如35号同学做习题35习题1：某厂计划生产甲、乙、丙三种零件，有机器、人工工时和原材料的限制，有关数据1、2、若原材料为2元/公斤，试建立获得最大利润生产计划的线性规划模型。
习题2：一塑料厂利用四种化工原料合成一种塑料产品。
这四种原料含A、B、C的成分见下表，这种塑料产品要求含A为25%，含B、C都不得少于30%。
问各种原料投放比例为习题3：建立以下线性规划模型1）某家具厂生产桌椅，每张桌子耗用木材0.28立方米、2小时人工，售价288元；每把椅子耗用木材0.13立方米、0.8小时人工，售价147元。
且1张桌子必须配4把椅子。
已知木材本月供应量不得超过52立方米，且每立方米成本价为500元。
本月人工工时上限为288小时，且每小时成本为20元。
（1）写出最大月收益线性规划模型；（2）写出月收益不低于8000元而动用木材最省的线性规划模型（其余条件不变）。
习题4 某工厂要用三种原料1、2、3混合调配出三种不同规格的产品甲、乙、丙，数据如右表。
问：该厂应如何安排生产，使利润收入为最大？习题5、某部门现有资金200万元，今后五年内考虑给以下的项目投资。
已知：项目A ：从第一年到第五年每年年初都可投资，当年末能收回本利110%；项目B ：从第一年到第四年每年年初都可投资，次年末能收回本利125%，但规定每年最大投资额不超过30万元；项目C ：需在第三年年初投资，第五年末能收回本利140%，但规定最大投资额不能超过80万元；项目D ：需在第二年年初投资，第五年末能收回本利155%，但规定最大投资额不能超过100万元；问：a.应如何确定这些项目的每年投资额，使得第五年年末拥有资金的本利金额为最大？ b.应如何确定这些项目的每年投资额，使得第五年年末拥有资金的本利在330万元的基础上使得其投资总的风险系数为最小？习题6 某公司计划在三年的计划期内，有四个建设项目可以投资：项目Ⅰ从第一年到第三年年初都可以投资。
数学建模lingo作业-习题讲解
基础题：1．目标规划问题最近，某节能灯具厂接到了订购16000套A 型和B 型节能灯具的订货合同，合同中没有对这两种灯具的各自数量做要求，但合同要求工厂在一周内完成生产任务并交货。
根据该厂的生产能力，一周内可以利用的生产时间为20000min ，可利用的包装时间为36000min 。
生产完成和包装一套A 型节能灯具各需要2min ；生产完成和包装完成一套B 型节能灯具各需要1min 和3min 。
每套A 型节能灯成本为7元，销售价为15元，即利润为8元；每套B 型节能灯成本为14元，销售价为20元，即利润为6元。
厂长首先要求必须按合同完成订货任务，并且即不要有足量，也不要有超量。
其次要求满意销售额达到或者尽量接近275000元。
最后要求在生产总时间和包装总时间上可以有所增加，但过量尽量地小。
同时注意到增加生产时间要比包装时间困难得多。
试为该节能灯具厂制定生产计划。
解：将题中数据列表如下：根据问题的实际情况，首先分析确定问题的目标级优先级。
第一优先级目标：恰好完成生产和包装完成节能灯具16000套，赋予优先因子p1；第二优先级目标：完成或者尽量接近销售额为275000元，赋予优先因子p2； 第三优先级目标：生产和包装时间的增加量尽量地小，赋予优先因子p3； 然后建立相应的目标约束。
在此，假设决策变量12,x x 分别表示A 型，B 型节能灯具的数量。
（1） 关于生产数量的目标约束。
用1d -和1d +分别表示未达到和超额完成订货指标16000套的偏差量，因此目标约束为1111211min ,..16000z d d s t x x d d -+-+=+++-=要求恰好达到目标值，即正、负偏差变量都要尽可能地小（2） 关于销售额的目标约束。
用2d -和2d +分别表示未达到和超额完成满意销售指标275000元的偏差值。
因此目标约束为221222min ,..1520-275000.z d s t x x d d --+=++=要求超过目标值，即超过量不限，但必须是负偏差变量要尽可能地小，(另外：d +要求不超过目标值，即允许达不到目标值，就是正偏差变量要尽可能地小) （3） 关于生产和包装时间的目标约束。
数学建模公交线路规划问题
3． 我校教职员工、学生的出行特点：上班、上课我校师生往返两校区的首要需求，结合我校教职 员工、学生的居住分布特点，因此我校教职员工、学生的出行特点十分明显，表现为时间空间 上的集中，具体特征如下： （1） 时间特点：上下课、上下班时间段（沙河校区—清水河校区：7:20、9:10、13:20、 15:10 、 18:20 ；清水河校区 — 沙河校区： 10:30 、 12:20 、 16:30 、 18:20 、 22 ： 20）出行人数骤增，其他时间段出行人数较少，甚至没有。 （2） 路线特点：起点、终点绝大多数为清水河校区、沙河校区两站。 本着 “保障教学科研工作开展， 满足师生往返两校” 的原则， 利用快速公交系统 （Bus Rapid Transit ——BRT）的便利因素、技术特点，结合我校师生出行特点，统筹便利性、社会效益、经济效益， 兼顾公交公司利益，进行方案制定。 2.1 线路选择 本线路以服务科大师生往返新老校区为初衷，所以在选择线路时，要使往返新老校区的时间最 短。由于交管部门数据不足，本文忽略由路况产生的拥塞、限速等情况，即认为路径最短时间最短。 2.2 站点设置 对于选择好的公交线路，在普通时段，与普通公交相同，按既定站点运行。在我校师生集中出 行时段，采用线路组合，即线路组合这种调度方式。首先我们对线路调度进行说明。 2.2.1 线路组合 此调度方式从普通线路按既定站点运行，站站停靠的方式派生出来。线路组合分标准线路、大 站快线、直达线路 ，并根据客流情况选择不同的方式（标准线路、大站快线、直达线路） 。它适用 于客流量大且集中，同时适用于开发分散的市郊区域。 其次对标准线路、大站快线、直达线路三种调度方式进行说明。 （1）标准线路：与普通公交线路相同，每站都停。
摘要
为配合我校和成都市公交规划部门，开设往返新老校区的快速公交线路。以高效便捷地保障广 大师生往返两校的交通需求。 本文解决了该公交线路的路线走向、站点设置、运行时长，发车间隔等设计问题，分析了拟定 的方案对学校的校车运行方案的影响，并作为向公交公司提供的策划论证的技术材料。本设计运用 Dijskra 算法，寻找到最快捷的路线走向。引入站点选择向量，发车间隔两个变量，结合客流量 OD 矩阵和站点距离矩阵，从出行时间成本和线路运营成本两个方面建立目标函数，运用遗传算法，求 解使目标函数最小的站点选择向量和发车间隔。 设计方案为：路线走向，沙河校区，一环路、蜀汉路、蜀西路、土龙路、金辉路、西源大道至 清水河校区。设置站点：电子科技大学沙河校区、苏宁电器建设路店、萤门口立交桥、蜀西路、土 龙路、金辉路、电子科大清水河校区。运行时间：7:30 首发车，21：30 末班车，共 14 小时。发车 间隔：11.43 分钟。
数学建模——规划问题
3.5 习题P 54.2.某工厂用21.A A 两台机床加工B B B 321,,三种不同零件。
已知在一个生产周期内A 1只能工作80机时；A 2只能工作100机时。
一个生产周期内计划加工B 1为70件、B 2为50件、B 3为20件。
两台机床加工每个零件的时间和加工每个零件的成本，分别如下列各表所示：加工每个零件时间表（单位：机时/个）加工每个零件成本表（单位：元/个）问：怎么样安排两台机床一个周期的加工任务，才能使加工成本最低？解：设在A 1机床加工零件B B B 321,,的数量分别为x x x 321,,，在A 2机床加工零件B B B 321,,的数量分别为x x x 654,,，建立如下线性规划模型：x x x x x x Z 654321633532min +++++=s.t.6,5,4,3,2,1,020507010038032635241654321=≥=+=+=+≤++≤++i x x x x x x x x x x x x x i改写成：[]⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=654321633532min x x x x x x Zs.t. 020507010010001001000100110080311000000321654321654321654321≥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡≤⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡x x x x x x x x x x x x x x x x x x结果： 解得7,0,40,0,25,30654321======x x x x x x297760340305253302min =⨯+⨯+⨯+⨯+⨯+⨯=Z 即机床A 1在一个周期加工零件B B B 321,,的数量分别为30件,25件,0件；A 2机床加工零件B B B 321,,的数量分别为40件，0件，7件；加工成本最低为297元。
数学建模作业数学规划模型----供应与选址的问题
再编写主程序liaochang2.m为：
clear
x0=[3 5 4 7 1 0 0 0 0 0 5 11 5 4 7 7];
A=[1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 1 1 1 1 1 1 0 0 0 0];
B=[20;20];
Aeq=[1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0
0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0
0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0
0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0
0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0
0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0];
使用临时料场的情形：
使用两个临时料场A(5,1)，B(2,7)．求从料场j向工地 的运送量 .在各工地用量必须满足和各料场运送量不超过日储量的条件下，使总的吨千米数最小，这是线性规划问题。线性规划模型为：
其中 ，i=1,2,…,6,j=1,2,为常数
设X11=X1,X21=X 2,,X31=X 3,X41=X 4,X51=X 5,,X61=X 6
程序截图如下：
程序的运行结果为：
xx =
3.0000
5.0000
0.0000
7.0000
0.0000
1.0000
0.0000
0.0000
4.0000
0.0000
6.0000
10.0000
fval =
136.2275
运行结果截图如下：
即由料场A、B向6个工地运料方案为：
数学建模-整数规划
整数规划
Integer Programming
数信学院 任俊峰
2012-4-15
数学建模之整数规划
整数规划模型（IP）
如果一个数学规划的某些决策变量或全部决策 变量要求必须取整数，则称这样的问题为整数规 划问题，其模型称为整数规划模型。 如果整数规划的目标函数和约束条件都是线性 的，则称此问题为整数线性规划问题.
松弛问题最优解满足整数要求，则该最优解为整数 规划最优解；
数学建模之整数规划
整数线性规划的求解方法
从数学模型上看整数规划似乎是线性规划的 一种特殊形式，求解只需在线性规划的基础上，通 过舍入取整，寻求满足整数要求的解即可。 但实际上两者却有很大的不同，通过舍入得到
的解（整数）也不一定就是最优解，有时甚至不能
1 xj 0
选中第j个项目投资 不 选中第j个项目投资
max Z 160 x 1 210 x 2 60 x 3 80 x 4 180 x 5 210 x 1 300 x 2 150 x 3 130 x 4 260 x 5 600 x x2 x3 1 1 x3 x4 1 x x 1 5 x1 , x 2 , x 3 , x 4 , x 5 0 或 1
1 2
14 x1 9 x 2 51 6 x1 3 x 2 1 x1 , x 2 0
数学建模之整数规划
用图解法求出最优解 x1＝3/2, x2 = 10/3 且有 z = 29/6 现求整数解（最优解）： 如用“舍入取整法”可得到4 个点即(1，3) (2，3) (1，4) (2，4)。显然，它们都不可能 是整数规划的最优解。
数学建模之整数规划
例5 固定费用问题
数学建模论文-投资规划问题
数学建模一周论文课程设计题目：投资规划问题摘要目前，证券在我国得到了迅速健康的发展，并且为我国的经济发展作出了很大贡献。
本文针对目前流行的各种不同的证券发行方案，建立线性规划模型，得出最佳的证券组合投资方案。
问题一中假设该经理有1000万资金可以进行投资支配，在满足题目给出的各限制范围内，以最大收益为目标函数，建立三个线性规划模型，分别为冒险模型、保守模型和一个折中模型，但是前两个不符合题目给出的约束条件，综合考虑，应选用折中模型，用Lingo求解得出了最大收益为29.83636万元，各种证券的投资方案见表二。
问题二中假设能以2.75%的利率借到不超过100万元资金，在相同的约束条件下，仍然建立线性规划模型，采用Lingo求解，得出最大收益为32.82000万元，投资方案见表五。
问题三中在1000万元资金情况下，若证券A的税前收益增加为4.5%，仍然建立线性规划模型，通过Lingo解得最大收益相对问题一中增加了，为30.27273万元，投资方案见表六；若证券C的税前收益减少为4.8%，用同样的方法求出最大收益相对问题一中减少了，为29.42400万元，投资方案见表七。
关键字：证券投资、线性规划、Lingo求解软件、投资风险某银行经理计划用一笔资金进行有价证券的投资，可供购进的证券及其信用等级、到期年限、收益如下表所示。
按照规定，市政证券的收益可以免税，其他证券的收益需按50%的税率纳税。
此外还有以下限制：●政府及代办机构的证券总共至少要购进400万元●所购证券的平均信用等级不超过1.4（信用等级数字越小，信用程度越高）●所购证券的平均到期年限不超过5年（2）如果能够以2.75%的利率借到不超过100万元资金，该经理应如何操作？（3）在1000万元资金情况下，若证券A的税前收益增加为4.5%，投资应否改变？若证券C的税前收益减少为4.8%，投资应否改变？模型假设1.假设在有价证券到期前，该经理不会中断投资。
数学建模案例之线性规划
资，每天最多购买多少桶牛奶？ 2）若可以聘用临时工人以增加劳动时间，付给临时工人的
工资最多是每小时多少元？ 3）由于市场需求变化，每公斤A1的获利增加到30元，应否
问题分析
企业内部的生产计划有各种不同的情况。
空间层次
数学建模案例之线性规划 奶制品的生产与销售
引言
优化问题及其一般模型：
优化问题是人们在工程技术、经济管理和科学研究等领域中 最常遇到的问题之一。例如： 设计师要在满足强度要求等条件下选择材料的尺寸， 使
结构总重量最轻； 公司经理要根据生产成本和市场需求确定产品价格，使所获
利润最高； 调度人员要在满足物质需求和装载条件下安排从各供应点
Step 2. 确定决策变量 第一来源：Step 1的结果，用变量固定需要回答的决策 第二来源：由决策导出的变量（具有派生结构） 其它来源：辅助变量（联合完成更清楚的回答）
Step 3. 确定优化目标 用决策变量表示的利润、成本等。
Step 4. 寻找约束条件 决策变量之间、决策变量与常量之间的联系。 第一来源：需求； 第二来源：供给； 其它来源：辅助以及常识。
引言
建立数学规划模型的步骤:
当你打算用数学建模的方法来处理一个优化问题的时候，首先要确定寻 求的决策是什么，优化的目标是什么，决策受到那些条件的限制（如果 有限制的话），然后用数学工具（变量、常数、函数等）表示它们，最 后用合适的方法求解它们并对结果作出一些定性、定量的分析和必要的 检验。
引言
Step 1. 寻求决策，即回答什么？必须清楚，无歧义。 阅读完题目的第一步不是寻找答案或者解法，而是……
例1：加工奶制品的生产计划
Lingo软件编程求解规划问题数学建模
例2——军队供给问题
用Lingo软件求解
结果
Lingo软件
Lingo 是一个可以简洁地阐述、解决和分析复杂问题的简便工具。
其特点是程序执行速度很快，易于输入、修改、求解和分析一个数 学规划问题。
LINGO模型的构成：4个段
• 集合段（SETS ENDSETS） • 目标与约束段 • 数据段（DATA ENDDATA） • 初始段（INIT ENDINIT）
用Lingo软件求解
求 z 72 x1 64 x2的最大值？
x1 x2 50
12
x1 8x2 480 3x1 100
x1 0
x2 0
Lingo结果中其它项的意思！ 更多信息
例 1 求 f = 2ex sin x 在 0<x<8 中的最小值与最大值
max=@sin(x)*2*@exp(-x); x<8; x>0;
获利24元/公斤 获利16元/公斤
每天： 50桶牛奶 时间480小时 至多加工100公斤A1
例1——加工奶制品的生产计划
1桶牛 奶或
12小时 8小时
3公斤A1 4公斤A2
每天：50桶牛奶 时间480小时
获利24元/公斤 获利16元/公斤 至多加工100公斤A1
决策变量 目标函数
约束条件
x1桶牛奶生产A1 获利 24×3x1 每天获利 原料供应 劳动时间 加工能力 非负约束
x2桶牛奶生产A2 获利 16×4 x2
Max z 72 x1 64 x2
x1 x2 50
12 x1 8x2 480 3x1 100
x1, x2 0
线性规 划模型 (LP)
例1——加工奶制品的生产计划
x1 x2 50
数学建模——规划模型
假设：料 场和工地 之间有直 线道路
1)现有 2 料场，位于 A (5, 1), B (2, 7),记为 (xj,yj),j=1,2, 日储量 ej 各有 20 吨。
i 1 i
n
i
a ik x k bi , i 1, 2 ,..., n. s.t . k 1 x 0 , i 1, 2 ,..., n. i
（3）二次规划问题
目标函数为二次函数，约束条件为线性约束
1 n min u f ( x ) ci xi bij xi x j 2 i , j 1 i 1 n a ij x j bi , i 1, 2,..., n. s.t . j 1 x 0 .i 1, 2,..., n. i
改写为： S.t.
min z 13 9 10 11 12 8X
0 0 800 0.4 1.1 1 0 X 0 0 0 0 . 5 1 . 2 1 . 3 900
x1 x2 x 3 ，X 0 x4 x 5 x 6
编写M文件xxgh4.m如下： c = [40 36]; A=[-5 -3]; b=[-45]; Aeq=[]; beq=[]; vlb = zeros(2,1); vub=[9;15]; %调用linprog函数： [x,fval] = linprog(c,A,b,Aeq,beq,vlb,vub)
（一）规划模型的数学描述
u f ( x)
数学建模动态规划问题
个阶段的决策过程有 个状态变量， 表示 演变的结果。在例1中 取 ，或定义为 ，即 。
根据过程演变的具体情况，状态变量可以是离散的或连续的。为了计算的方便有时将连续变量离散化；为了分析的方便有时又将离散变量视为连续的。
状态变量简称为状态。
2.1.3决策
当一个阶段的状态确定后，可以作出各种选择从而演变到下一阶段的某个状态，这种选择手段称为决策（decision），在最优控制问题中也称为控制（control）。
.
决策 Байду номын сангаас允许集合为
.
状态转移方程和阶段指标应对 的每个取值 和 的每个取值 计算，即 ， 。最优值函数应对 的每个取值 计算。基本方程可以表为
（4）
按照（3），（4）逆向计算出 ，为全过程的最优值。记状态 的最优决策为 ，由 和 按照状态转移方程计算出最优状态，记作 。并得到相应的最优决策，记作 。于是最优策略为 。
描述决策的变量称决策变量（decision variable），变量允许取值的范围称允许决策集合（set of admissible decisions）。用 表示第 阶段处于状态 时的决策变量，它是 的函数，用 表示 的允许决策集合。在例1中 可取 或 ，可记作 ，而 。
决策变量简称决策。
2.1.4策略
（ ）写出基本方程即最优值函数满足的递归方程，以及端点条件。
数学建模之农场规划问题
农场规划问题问题重述：某农户拥有100亩土地和15000元可供投资，每年冬季（9月中旬至来年5月中旬），该家庭的成员可以贡献3500小时的劳动时间，而夏季为4000小时。
如果这些劳动时间有富裕，该家庭中的年轻成员将去附近的农场打工，冬季每小时6.8元，夏季每小时7.0元。
现金收入来源于三中农作物（大豆、玉米和燕麦）以及奶牛和母鸡。
农作物不需要付出投资，但每头奶牛需要400元的初始投资，可产奶3年，每只母鸡需要3元的吃食投资，只饲养1年。
每头奶牛需要1.5亩的土地，并且冬季需要付出100小时劳动时间，夏季付出50小时劳动时间，每年产生的净现金收入为1350元；每只母鸡的对应数字为：不占用土地，冬季0.6小时，夏季0.3小时，年净现金收入10.5元。
养鸡厂房最多容纳3000只母鸡，栅栏的大小限制了最多能饲养32头奶牛。
根据统计，三种农作物每种植一亩所需要的劳动时间和收入数据分别为：大豆：冬季20小时，夏季30小时，年净收入360.0元；玉米：冬季35小时，夏季75小时，年净收入600.0元；燕麦：冬季10小时，夏季40小时，年净收入400.0元。
基本假设：1、假设该农户每年都能及时获得现金收入，即本年度所获得的利润可及时用于下一年的投资；2、第五年的投资也考虑到计算中。
问题分析:这个问题的目标是使得5年内净现金收入最大，要做的决策是生产规划，即确定每种农作物应该种植多少亩，奶牛和鸡各应蓄养多少只，决策受到6个变量的限制，即土地总面积、投资资金、劳动力时间（夏季和冬季）以及奶牛和鸡的总饲养量。
模型建立:决策变量：设用i=0,1,2,3,4,5表示年数，用j=1,2,3,4,5分别表示三种农作物（大豆、玉米、燕麦）及奶牛和母鸡。
可表示第i年种植三种农作物的亩数或者蓄养奶牛和母鸡的个数，表示第i 年的总现金收入。
目标函数：设第i年的总获利为元，因农作物不用投资，则第i年种植大豆为亩，每亩收入360元，获利360元；第i年种植玉米亩，每亩收入600元，获利600；第i年种植燕麦亩，每亩收入400元，获利400元；第i年买奶牛头，每头收入1350元，获利1350（++）元；第i年鸡购买只，每只收入10.5元，获利10.5元；若劳动力有剩余，则第i年夏季劳动力收入［4000-(3075)］元，冬季劳动力收入［3500-(2035)］元。
数学建模第1章线性规划
数学
建模
例 1.6
min{max
xi
yi
|
ei
|}，其中e i
=
xi -
yi 。
取v
=
max yi
|
e
i
|，这样，上面的问题就变换成
min v，
s.t.
ìïïíïïî
x1 y1
-
y1 ? x1 ?
v,L , xn v,L , yn
yn ? v, n ? v.
25/39
基础部数学教研室
数学 建模
2x1 - 5x2 + x3 ? 10， x1 + 3x2 + x3 ? 12， x1, x2 , x3 ³ 0.
11/39
基础部数学教研室
数学 建模
解 （1）化成 Matlab 标准型
min w = - 2x1 - 3x2 + 5x3，
s.t.
轾 犏- 2 犏 臌1
5 3
-1 1
轾 犏x1 犏 犏x2 犏 臌x3
a=1 -1 -1 1 1 -1 1 -3 1 -1 -2 3;
enddata
min=@sum(col:c*@abs(x));
@for(row(i):@sum(col(j):a(i,j)*x(j))<b(i));
@for(col:@free(x)); !x的分量可正可负;
end
24/39
基础部数学教研室
@for(row(i):@sum(col(j):a(i,j)*x(j))<b(i));
@sum(col:x)=7;
14/39
end
基础部数学教研室
数学 建模
例 1.2 求解下列线性规划问题 max z = 2x1 + 3x2 - 5x3， s.t. x1 + x2 + x3 = 7， 2x1 - 5x2 + x3 ? 10， x1 + 3x2 + x3 ? 12， x1, x2 , x3 ³ 0.
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
一、线性规划
1.简介
1.1适用情况
用现有资源来安排生产，以取得最大经济效益的问题。
如： （1）资源的合理利用
（2）投资的风险与利用问题 （3）合理下料问题 （4）合理配料问题 （5）运 输 问 题 （6）作物布局问题
（7）多周期生产平滑模型 （8）公交车调度安排 1.2建立线性规划的条件
（1）要求解问题的目标函数能用数值指标来反映，且为线性函数；
（2）要求达到的目标是在一定条件下实现的，这些约束可用线性等式或不等式描述。
1.3线性规划模型的构成
决策变量、目标函数、约束条件。
2、一般线性规划问题 数学标准形式：
目标函数：
1
max ==
∑ n
j
j
j z c
x
约束条件：1
,1,2,...,,..0,1,2,...,.=⎧==⎪⎨⎪≥=⎩∑n
ij j i j j
a x
b i m s t x j n
matlab 标准形式：
3、可以转化为线性规划的问题 例：求解下列数学规划问题
解：作変量変换1||||,,1,2,3,4,22
+-===i i i i
i x x x x u v i 并把新变量重新排序成一维变量
[]1414,,,,
,⎡⎤
==⎢⎥⎣⎦
T
u y u u v v v ，则可把模型转化为线性规划模型
其中：[]1,2,3,4,1,2,3,4;=T c 12,1,;2⎡⎤=---⎢⎥⎣
⎦T
b 111111131 - - ⎡⎤
⎢⎥= - -⎢⎥⎢⎥ -1 -1 3⎣⎦
A 。
利用matlab 计算得最优解：12342,0,=-===x x x x 最优值z=2。
程序如下： 略
二、整数规划 1.简介
数学规划中的变量(部分或全部)限制为整数时称为整数规划。
目前流行求解整数规划的方法一般适用于整数线性规划。
1.1整数规划特点
1）原线性规划有最优解，当自变量限制为整数后，出现的情况有
①原线性规划最优解全是整数，则整数规划最优解与线性规划最优解一致。
②整数规划无可行解。
③有可行解（存在最优解），但最优解值变差。
2）整数规划最优解不能按照实数最优解简单取整获得。
1.2求解方法分类
（1）分枝定界法—可求纯或混合整数线性规划。
（2）隔平面法—可求纯或混合整数线性规划。
（3）隐枚举法—可求“0-1”整数规划。
（4）匈牙利法—解决指派问题。
（5）蒙特卡洛法—求解各种类型规划. 1.3整数规划的应用模型 （1）固定费用的问题。
（2）指派问题。
（3）合理下料问题。
（4）流动推销员问题。
（5）生产与销售计划问题。
2、一般整数规划模型
目标函数： 约束条件：
例：指派问题的数学模型(0-1型整数规划)
拟分配n 人去做n 项工作，若分配第i 人去做第j 项工作，需花费ij c 单位时间，如何分配工作才能使花费总时间最少？ 模型的建立
引入0-1变量1,i j 0第人做第项工作，
，第i 人不做第j 项工作.
⎧=⎨⎩ij x
指派问题的数学模型为
利用匈牙利算法、拍卖算法等求解出最优解。
三、非线性规划 1、简介
目标函数或约束条件中包含非线性函数的规划问题为非线性规划问题。
1.1非线形规划模型的构成
决策变量、目标函数、约束条件。
1.2非线性规划的应用模型
（1）存贮模型
（2）飞行管理问题 （3）森林救火
（4）抽水费用最小问题 （5）钢管下料问题 （6）投资决策问题 （7）供应与选址问题 （8）广告的费用及其效用 2、非线性规划的模型 一般形式： 其中：[]1,
,=T
n x x x 为模型的决策变量。
Matlab 中非线性规划的数学模型
其中：f(x)是标量函数；A,b,Aeq,beq,lb,ub 是相应维数的矩阵和向量；c(x),cex(x)是非线性向量函数。
3、罚函数法
利用罚函数法可将非线性规划问题的求解转化为求解一系列无约束极值问题。
问题
取一个充分大的数M>0，构造函数
（或()()(,)()max min ||()||,00m G x H x p x M f x Msum Msum M k x ⎛⎫⎛⎫
⎛⎫⎛⎫=+-+ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭⎝⎭这里
[][][]111()(),
,(),()(),
,(),()(),
,(),r s t G x g x g x H x h x h x K x k x k x ===可直接利用matlab
中的max 、min 和sum 函数），则增广目标函数P(x,M)为目标函数的无约束极值问题
minP(x,M)的最优解x 即为原问题的最优解。
注意：
1）如果非线性规划问题要求实时算法，则可用罚函数法，但计算精度较低。
2）如果非线性规划问题不要求实时算法，但要求精度高，可使用Lingo 软件编程求解或使用Matlab 的fmincon 命令求解。
四、目标规划 1、简介
1.1求解目标规划的思路
（1）加权系数法
为每一目标赋一个权系数，把多目标模型转化成单一目标的模型。
但困难是要确定合理的权系数，以反映不同目标之间的重要程度。
（2）优先等级法
将各目标按其重要程度不同的优先等级，转化为单目标模型。
（3）有效解法
寻求能够照顾到各个目标，并使决策者感到满意的解。
由决策者来确定选取哪一个解，即得到一个满意解。
但有效解的数目太多而难以将其一一求出。
1.2建立目标规划的条件
（1）正、负偏差变量。
（2）绝对（刚性）约束和目标约束。
（3）优先因子（优先等级）与权系数。
1.3 目标规划的目标函数
目标规划的目标函数基本三种形式为
（1）第i 个目标要求恰好达到目标值，即正、负偏差变量都要尽可能地小，这时
i i i i min w d w d --++ +.
（2）第i 个目标要求不超过目标值，即允许达不到目标值，就是正偏差变量要尽可能地小，这时
（3）第i 个目标要求超过目标值，即超过量不限，但必须是负偏差变量要尽可能地小，这时
1.4目标规划的模型应用
（1）求多目标下产品利润最优的决策方案。
（2）求多目标下总运费最小的运输调度方案。
2、目标规划的一般数学模型
设(1,2,
,)j x j n =是目标规划的决策变量，共有m 个约束是刚性约束，可能
是等式约束，也可能是不等式约束。
设有l 个柔性目标约束，其目标规划约束的偏差为,(1,2,
,)i i d d i l +-=。
设有q 个优先级别，分别为12,,,q p p p ⋯。
在同一
个优先级k p 中，有不同的权重，分别记为,(1,2,,)ki
ki w w i l +-=。
目标规划模型的
一般数学表达式如下
可用序贯算法求解目标规划。
3、数据包网络分析（DEA ） 3.1适用范围
DEA 特别适用于具有多输入多输出的复杂系统，如技术进步、技术创新、资源配置、金融投资等领域，特别对非单纯利益公共部门，如学校、医院、某些文化设施的评价方面。
3.2特点
1）DEA 以决策单位各输入/输出的权重为变量， 1）数据包络分析的C2R 模型 设有n 个DMU ，每个DMU 都有m 种投入和s 种产出，设1,
,;1,,ij x i m j n ==（）表示第
j 个DMU 的第i 种投入量， 1,,1,,rj y r s j n ==（；）表示第 j 个DMU 的第r 种产出量，(1,
,)i v i m =表示第i 种投入的权值，1,
,r u r s =（）表示第r 种产出的权值。
向量,(1,
,)j j X Y j n =分别表示决策单元 j 的输入和输出向量，v 和u 分别表示输入输
出权值向量，则12(,,
,)T j j j mj X x x x =，12(,,,),T j j j sj Y x x x =12(,,,)T m u u u u =，
12(,,
,)T S v v v v =。
定义决策单元j 的效率评价指数为
评价决策单元0j 效率的数学模型为
00
max
,
1,1,2,,,..0,0,0,0.
T j T
j T j
T
j
u Y v X u Y j n s t v X u v u v ⎧≤=⎪⎨⎪
≥≥≠≠⎩ （1） 通过Charnes ?Cooper 变换：01
,,()
T j tv tu t v X ωμ===可以将模型（1）转化为等
价的线性规划问题
对于C2R 模型，有如下定义：
（1）若线性规划问题的最优目标01j v =,则称决策单元0j 是弱DEA 有效的。
（2）若线性规划问题存在最优解**0,0,μω>>并且其最优目标值01j V =，则称决策单元0j 是EDA 有效的。