数学建模之规划问答
- 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 标准形式:
min ,
,.,.⋅≤⎧⎪
⋅=⎨⎪≤≤⎩
T s t Aeq beq lb ub f x A x b x x
3、可以转化为线性规划的问题
例:求解下列数学规划问题
1234123412341234min ||2||3||4||,2,..31,123.
2=+++⎧
⎪--+≤-⎪-+-≤-⎨⎪⎪--+≤-⎩
z x x x x x x x x s t x x x x x x x x
解:作変量変换1||||,,1,2,3,4,22
+-=
==i i i i
i x x x x u v i 并把新变量重新排序成一维变量[]1414,,,,,⎡⎤==⎢⎥⎣⎦
L L T
u y u u v v v ,则可把模型转化为线性规划模型
[]min ,
,,..0.⎧⎡⎤-≤⎪⎢⎥⎨⎣⎦⎪≥⎩
T c y u A A b s t v y 其中:[]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、一般整数规划模型
目标函数:
1
max ==
∑ n
j
j j z c
x
约束条件:
1
,1,2,...,,..0,1,2,...,.=⎧==⎪⎨⎪≥=⎩∑n
ij j j j j
a x
b i m s t x j n
例:指派问题的数学模型(0-1型整数规划)
拟分配n 人去做n 项工作,若分配第i 人去做第j 项工作,需花费ij c 单位时间,如何分配工作才能使花费总时间最少? 模型的建立
引入0-1变量1,i j 0第人做第项工作,
,第i 人不做第j 项工作.⎧=⎨⎩ij x
指派问题的数学模型为
1
1
min ,===
∑∑ n
n
ij
ij i j z c
x
111,1,2,...,,..1,1,2,...,,01,,1,2,...,.或==⎧==⎪⎪⎪==⎨⎪⎪= =⎪⎩
∑∑n
ij j n ij i ij x i n s t x j n x i j n 利用匈牙利算法、拍卖算法等求解出最优解。
三、非线性规划
1、简介
目标函数或约束条件中包含非线性函数的规划问题为非线性规划问题。 1.1非线形规划模型的构成
决策变量、目标函数、约束条件。 1.2非线性规划的应用模型
(1)存贮模型 (2)飞行管理问题 (3)森林救火
(4)抽水费用最小问题 (5)钢管下料问题 (6)投资决策问题 (7)供应与选址问题 (8)广告的费用及其效用
2、非线性规划的模型
一般形式:
min (),
()0,1,2,,,..()0,1,2,,.目标函数:约束条件:≥=⎧⎨==⎩L L j
i f h j q s t g i p x x x
其中:[]1,,=L T
n x x x 为模型的决策变量。 Matlab 中非线性规划的数学模型
min (),
,,..()0,()0,.
f x A x b Aeq x beq s t c x ceq x lb x ub ⋅≤⎧⎪⋅=⎪⎪≤⎨⎪=⎪≤≤⎪⎩ 其中:f(x)是标量函数;A,b,Aeq,beq,lb,ub 是相应维数的矩阵和向量;c(x),cex(x)
是非线性向量函数。
3、罚函数法
利用罚函数法可将非线性规划问题的求解转化为求解一系列无约束极值问题。 问题
min (),
()0,1,,,
..()0,1,,,()0,1,,.
i j m f x g x i r s t h x j s k x m t ⎧≤=⎪
≥=⎨⎪
==⎩L L L 取一个充分大的数M>0,构造函数
1
1
1
(,)()max((),0)min((),0)|()|,r s t
i j m i j m p x M f x M g x M h x M k x ====+-+∑∑∑
(或()()(,)()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 ===L L L 可直接利用matlab 中的max 、min 和sum 函数),则增广目标函数P(x,M)为目标函数
的无约束极值问题minP(x,M)的最优解x 即为原问题的最优解。 注意:
1)如果非线性规划问题要求实时算法,则可用罚函数法,但计算精度较低。