1
简单的优化模型
3.1 存贮模型
3.2 生猪的出售时机
3.3 森林救火
3.4 最优价格
3.5 血管分支
3.6 消费者均衡
3.7 冰山运输
2
现实世界中普遍存在着优化问题
静态优化问题指最优解是 数 (不是函数 )
建立静态优化模型的 关键之一 是根据建模目的确定恰当的目标函数
求解静态优化模型一般用 微分法静 态 优 化 模 型
3
1 存贮模型问 题配件厂为装配线生产若干种产品,轮换产品时因更换设备要付生产准备费,产量大于需求时要付贮存费。该厂生产能力非常大,即所需数量可在很短时间内产出。
已知某产品日需求量 100件,生产准备费 5000元,贮存费每日每件 1元。试安排该产品的生产计划,即多少天生产一次(生产周期),每次产量多少,使总费用最小。
要求不只是回答问题,而且要建立生产周期、产量与需求量、准备费、贮存费之间的关系。
4
问题分析与思考
每天生产一次,每次 100件,无贮存费,准备费 5000元。
日需求 100件,准备费 5000元,贮存费每日每件 1元。
10天生产一次,每次 1000件,贮存费 900+800+…+100 =4500
元,准备费 5000元,总计 9500元。
50天生产一次,每次 5000件,贮存费 4900+4800+…+100
=122500元,准备费 5000元,总计 127500元。
平均每天费用 950元平均每天费用 2550元
10天生产一次平均每天费用最小吗?
每天费用 5000元
5
这是一个优化问题,关键在建立目标函数。
显然不能用一个周期的总费用作为目标函数目标函数 —— 每天总费用的平均值
周期短,产量小
周期长,产量大问题分析与思考贮存费少,准备费多准备费少,贮存费多存在最佳的周期和产量,使总费用 (二者之和 )最小
6
模 型 假 设
1,产品每天的需求量为常数 r;
2,每次生产准备费为 c1,每天每件产品贮存费为 c2;
3,T天生产一次(周期),每次生产 Q件,当贮存量为零时,Q件产品立即到来(生产时间不计);
建 模 目 的设 r,c1,c2 已知,求 T,Q 使每天总费用的平均值最小。
4,为方便起见,时间和产量都作为连续量处理。
7
模 型 建 立
0 t
q
贮存量表示为时间的函数 q(t)
T
Q
rt=0生产 Q件,q(0)=Q,q(t)以需求速率 r递减,q(T)=0.
一周期总费用 TQccC 2~ 21
每天总费用平均值(目标函数) 2~)( 21 rTcTcTCTC
离散问题连续化
Acdttqc T 202 )(
一周期贮存费为
A=QT/2
2
2
21
rTcc
rTQ?
8
模型求解
M in2)( 21 rTcTcTC
求 T 使
0?dTdC
2
12
c
rcrTQ
2
12
rc
cT?
模型分析
QTc,1 QTc,2 QTr,
模型应用 c1=5000,c2=1,r=100
T=10(天 ),Q=1000(件 ),C=1000(元 )? 回答问题
9
经济批量订货公式 ( EOQ公式 )
2
12
rc
cT?
2
12
c
rc
rTQ
每天需求量 r,每次订货费 c1,每天每件贮存费 c2,
用于订货、供应、存贮情形不允许缺货的存贮模型
问:为什么不考虑生产费用?在什么条件下才不考虑?
T天订货一次 (周期 ),每次订货 Q件,当贮存量降到零时,Q件立即到货。
10
允许缺货的存贮模型
A
B0
q
Q
r
T1 t
当贮存量降到零时仍有需求 r,
出现缺货,造成损失原模型假设:贮存量降到零时 Q件立即生产出来 (或立即到货 )
现假设:允许缺货,每天每件缺货损失费 c3,缺货需补足
T
1rTQ?
Acdttqc T 202 1 )(
一周期贮存费
Bcdttqc TT 33
1
)(
一周期缺货费周期 T,t=T1贮存量降到零
2
)(
2
2
1
3
1
21
TTrcQTccC
一周期总费用
11
rT
QrTc
rT
Qc
T
c
T
CQTC
2
)(
2),(
2
3
2
21
0,0 QCTC
每天总费用平均值
(目标函数)
2
13121 )(2
1
2
1 TTrcQTccC一周期总费用
M in),(?QTC求 T,Q 使
3
32
2
12
c
cc
rc
cT
32
3
2
12
cc
c
c
rcQ

为与 不允许缺货的存贮模型相比,T记作 T ’,Q记作 Q’
12
2
12
rc
cT?
2
12
c
rcrTQ
不允许缺货模型
QQTT,
3
32
c
cc记
1 QQTT ','
13c QQTT,
3
32
2
12'
c
cc
rc
cT
32
3
2
12'
cc
c
c
rcQ
允许缺货模型不允许缺货
3c
13
3
32
2
12
c
cc
rc
cT
32
3
2
12
cc
c
c
rcQ

允许缺货模型
0
q
Q?
r
T1 tT注意:缺货需补足
Q?~每周期初的存贮量
R
每周期的生产量
R (或订货量)
3
32
2
12
c
cc
c
rcTrR
Q~不允许缺货时的产量 (或订货量 ) QQR
14
2 生猪的出售时机饲养场每天投入 4元资金,用于饲料、人力、设备,估计 可使 80千克重的生猪体重增加 2公斤。
问题市场价格目前为每千克 8元,但是 预测 每天会降低 0.1元,问生猪应何时出售。
如果 估计 和 预测 有误差,对结果有何影响。
分析投入资金使生猪体重随时间增加,出售单价随时间减少,故存在最佳出售时机,使利润最大
15
trtgttQ 4)80)(8()(
求 t 使 Q(t)最大
rg
grt 2404
10天后出售,可多得利润 20元建模及求解生猪体重 w=80+rt
出售价格 p=8-gt
销售收入 R=pw
资金投入 C=4t
利润 Q=R-C=pw -C
估计 r=2,
若当前出售,利润为 80× 8=640(元)
t 天出售
=10
Q(10)=660 > 640
g=0.1
16
敏感性分析研究 r,g变化时对模型结果的影响 估计 r=2,g=0.1
rg
grt 2404
设 g=0.1不变
5.1,6040 rrrt
t 对 r 的(相对)敏感度
rr
ttrtS

/Δ),(?
t
r
dr
dt?
36040 60),( rrtS
生猪每天体重增加量 r 增加 1%,出售时间推迟 3%。
1,5 2 2,5 3
0
5
10
15
20
r
t
17
敏感性分析估计 r=2,g=0.1
rg
grt 2404
研究 r,g变化时对模型结果的影响
设 r=2不变
15.00,203 gg gt
t 对 g的(相对)敏感度
t
g
dg
dt
gg
ttgtS

/Δ),(
3203 3),( ggtS
生猪价格每天的降低量 g增加 1%,出售时间提前 3%。
0,0 6 0,0 8 0,1 0,1 2 0,1 4 0,1 6
0
10
20
30
g
t
18
强健性分析保留生猪直到利润的增值等于每天的费用时出售由 S(t,r)=3
建议过一周后 (t=7)重新估计,再作计算 。wwpp,,,
研究 r,g不是常数时对模型结果的影响
w=80+rt?w = w(t)
4)()()()( twtptwtp
p=8-gt? p =p(t)
若 (10%),则 ( 30%)2.28.1 w 137 t
0)( tQ
每天利润的增值 每天投入的资金
ttwtptQ 4)()()(
25
4 最优价格问题根据产品成本和市场需求,在产销平衡条件下确定商品价格,使利润最大假设 1)产量等于销量,记作 x
2)收入与销量 x 成正比,系数 p 即价格
3)支出与产量 x 成正比,系数 q 即成本
4)销量 x 依赖于价格 p,x(p)是减函数建模与求解
pxpI?)(收入 qxpC?)(支出
)()()( pCpIpU利润进一步设 0,,)( babpapx
求 p使 U(p)最大
26
0
*
ppdp
dU
使利润 U(p)最大的最优价格 p*满足
** pppp dp
dC
dp
dI

最大利润在边际收入等于边际支出时达到
pxpI?)(
qxpC?)(
bpapx)(
))(( bpaqp
)()()( pCpIpU
b
aqp
22
*
建模与求解边际收入 边际支出
27
结果解释
b
aqp
22
* 0,,)( babpapx
q / 2 ~ 成本的一半
b ~ 价格上升 1单位时销量的下降幅度(需求对价格的敏感度)
a ~ 绝对需求 ( p很小时的需求 )
b p*?
a p*?
思考:如何得到参数 a,b?
28
3.5 血 管 分 支背景机体提供能量维持血液在血管中的流动给血管壁以营养 克服血液流动的阻力消耗能量取决于血管的几何形状在长期进化中动物血管的几何形状已经达到能量最小原则研究在能量最小原则下,血管分支处粗细血管半径比例和分岔角度问题
29
模型假设一条粗血管和两条细血管在分支点对称地处于同一平面血液流动近似于粘性流体在刚性管道中的运动血液给血管壁的能量随管壁的内表面积和体积的增加而增加,管壁厚度近似与血管半径成正比 q q
1
q1
A
B
B′
C
H
L
l
l1
r
r1
q=2q1 r/r1,
考察血管 AC与 CB,CB′
30
粘性流体在刚性管道中运动
l
prq
8
4?
p~A,C压力差,
~粘性系数克服阻力消耗能量
4
2
1
8
d
lqpqE

提供营养消耗能量
21,2 lbrE
管壁内表面积 2?rl
管壁体积?(d2+2rd)l,
管壁厚度 d与 r成正比模型假设
q q1
q1
A
B
B′
C
H
L
l
l1
r
r1
31
模型建立
q q1
q1
A
B
B′
C H
L
l
l1
r
r1
4
2
1
8
d
lqpqE

克服阻力消耗能量
21,2 lbrE
提供营养消耗能量
1141214221 2)/()/( lbrrkqlbrrkqEEE
s in/,/ 1 HLltgHLl

s i n/2)/(
)t a n/)(/(),,(
1
4
1
2
1
42
1
Hbrrkq
HLbrrkqrrE

机体为血流提供能量
32
模型求解
q q1
q1
A
B
B′
C H
L
l
l1
r
r1

0,0
1
rErE
0/4
0/4
5
1
21
1
521


rkqrb
rkqrb
4
1
1
4
r
r
0E
4
1
2c o s



r
r?
4
4
2c o s?

21 001 4937,32.1/26.1rr

s in/2)/(
)t a n/)(/(),,(
1
4
1
2
1
42
1
Hbrrkq
HLbrrkqrrE

33
模型解释生物学家:结果与观察大致吻合大动脉半径 rmax,毛细血管半径 rmin
大动脉到毛细血管有 n次分岔
4
1
1
4
r
r
4
m in
m a x 4
n
r
r
5m i nm a x 41000/rr
21
00
1
4937
32.1/26.1


rr
观察:狗的血管
)4(5n 30~25?n
血管总条数 973025 10~1032~22n
推论 n=?
34
q2
U(q1,q2) = c
q10 1l
2l
3l
3.6 消费者均衡问题 消费者对甲乙两种商品的偏爱程度用无差别曲线族表示,问他如何分配一定数量的钱,购买这两种商品,以达到最大的满意度。
设甲乙数量为 q1,q2,消费者的无差别曲线族
(单调减、下凸、不相交),记作 U(q1,q2)=c
U(q1,q2) ~ 效用函数已知甲乙价格 p1,p2,有钱 s,试分配 s,
购买甲乙数量 q1,q2,使 U(q1,q2)最大,
35
s/p2
s/p1
q2
U(q1,q2) = c
q10 1l
2l
3l
模型及求解已知价格 p1,p2,钱 s,
求 q1,q2,或 p1q1 / p2q2,
使 U(q1,q2)最大 sqpqpts
qqUZ

2211
21
..
),(m a x
),( 2211 qpqpUL )2,1(0
i
q
L
i
2
1
2
1
p
p
q
U
q
U

1
2
2 dq
dqK
l
几何解释
sqpqp 2211直线 MN,
最优解 Q,MN与 l2切点
21 / ppK MN
斜率 ·M
Q

·
21
/
q
U
q
U

36
0,0,0,0,0.B
21
2
2
2
2
2
1
2
21





qq
U
q
U
q
U
q
U
q
U
2
1
2
1
p
p
q
U
q
U
结果解释 21,q
U
q
U
—— 边际效用消费者均衡状态在两种商品的边际效用之比恰等于它们价格之比时达到。
效用函数 U(q1,q2) 应满足的条件
A,U(q1,q2) =c 所确定的函数 q2=q2(q1)单调减、下凸
解释 B的实际意义AB?
37
0,,)(.1 1
21

qq
U
效用函数 U(q1,q2) 几种常用 的形式
2
1
2
1
p
p
q
U
q
U
2
1
22
11
p
p
qp
qp

消费者均衡状态下购买两种商品费用之比与二者价格之比的平方根成正比。
U(q1,q2)中参数?,? 分别表示消费者对甲乙两种商品的偏爱程度。
38
1,0,.2 21 qqU
0,,)(.3 221 baqbqaU
2
1
2
1
p
p
q
U
q
U

22
11
qp
qp
购买两种商品费用之比与二者价格无关。
U(q1,q2)中参数?,? 分别表示对甲乙 的偏爱程度。
思考:如何推广到 m ( > 2) 种商品的情况效用函数 U(q1,q2) 几种常用 的形式
39
3.7 冰山运输背景? 波斯湾地区水资源贫乏,淡化海水的成本为每立方米 0.1英镑。
专家建议从 9600千米远的南极用拖船运送冰山,取代淡化海水
从经济角度研究冰山运输的可行性。
建模准备 1,日租金和最大运量船 型 小 中 大日租金(英镑)
最大运量(米 3)
4.0 6.2 8.0
5?105 106 107
40
2,燃料消耗(英镑 /千米)
3,融化速率(米 /天)
与南极距离 (千米 )
船速 (千米 /小时 ) 0 1000 >4000
1
3
5
0 0.1 0.3
0 0.15 0.45
0 0.2 0.6
冰山体积 (米 3)
船速 (千米 /小时 ) 10
5 106 107
1
3
5
8.4 10.5 12.6
10.8 13.5 16.2
13.2 16.5 19.8
建模准备
41
建模目的选择船型和船速,使冰山到达目的地后每立米水的费用最低,并与淡化海水的费用比较模型假设
航行过程中船速不变,总距离 9600千米
冰山呈球形,球面各点融化速率相同
到达目的地后,每立方米冰可融化 0.85立方米水建模分析目的地水体积 运输过程融化规律总费用目的地冰体积初始冰山体积燃料消耗租金船型,船速船型船型,船速船型
42


4000),1(
40000),1(
2
1
dbua
dbuda
r
4.0,2.0,105.6 251 baa
模型建立
1,冰山融化规律船速 u (千米 /小时 )
与南极距离 d(千米 )
融化速率 r(米 /天)
r是 u 的线性函数;
d<4000时 u与 d成正比
d>4000时 u与 d无关,
utd 24?
航行 t 天



u
tu
u
ttuu
r
t
6
1 0 0 0
),4.01(2.0
6
1 0 0 0
0,)4.01(1056.1
3
第 t天融化速率
0 1000 >4000
1
3
5
0 0.1 0.3
0 0.15 0.45
0 0.2 0.6
u r
d
43
1,冰山融化规律

t
k
kt rRR
1
0
冰山初始半径 R0,航行 t天时半径冰山初始体积
3
00 3
4 RV
3
3
4
tt RV
t天时体积总航行天数
3
1
3 0
0 4
3
3
4),,(



t
k
kr
VtVuV
选定 u,V0,航行
t天时冰山体积
3
1
3 0
0 4
3
3
4),(



T
t
tr
VVuV
到达目的地时冰山体积
uuT
4 0 0
24
9 6 0 0
44
1,6,3.0 321 ccc


1
4
3
3
4
lo g)6(2.7
]),,() [ lo g(24),,(
3
1
3 0
10
3010210
t
k
kr
V
uu
ctVuVcucutVuq
),) ( l o g( 310211 cVcucq
2,燃料消耗
105 106 107
1
3
5
8.4 10.5 12.6
10.8 13.5 16.2
13.2 16.5 19.8
V
u q1燃料消耗 q1(英镑 /千米 )
q1对 u线性,对 log10V线性选定 u,V0,航行第 t天燃料消耗 q (英镑 /天 )
燃料消耗总费用
T
t
tVuqVuQ
1
00 ),,(),(
45
V0 5?105 106 107
f(V0) 4.0 6.2 8.0
3,运送每立方米水费用 冰山初始体积 V0的日租金 f(V0)(英镑)
uT
400?航行天数总燃料消耗费用拖船租金费用
uVfVuR
400)(),(
00
冰山运输总费用 ),(),(),(
000 VuQVuRVuS




1
4
3
3
4
l o g)6(2.7),(
3
1
3 0
10
1
0
t
k
k
T
t
r
V
uuVuQ
46
冰山到达目的地后得到的水体积 ),(85.0),( 00 VuVVuW?
3,运送每立方米水费用冰山运输总费用运送每立方米水费用 ),(
),(),(
0
0
0 VuW
VuSVuY?
3
1
3 0
0 4
3
3
4),(



T
t
tr
VVuV
到达目的地时冰山体积
),(),(),( 000 VuQVuRVuS
47
模型求解选择船型和船速,使冰山到达目的地后每立方米水的费用最低求 u,V0使
Y(u,V0)最小
u=4~5(千米 /小时 ),V0= 107 (米 3),Y(u,V0)最小
V0只能取离散值经验公式很粗糙
3 3.5 4 4.5 5
107 0.0723 0.0683 0.0649 0.0663 0.0658
0.2251 0.2013 0.1834 0.1842 0.1790
106 78.9032 9.8220 6.2138 5.4647 4.5102
V0 u
5?106
取几组( V0,u)用 枚举法 计算
48
结果分析由于未考虑影响航行的种种不利因素,冰山到达目的地后实际体积会显著小于 V(u,V0)。
有关部门认为,只有当计算出的 Y(u,V0)显著低于淡化海水的成本时,才考虑其可行性。
大型拖船 V0= 107 (米 3),船速 u=4~5(千米 /小时 ),冰山到达目的地后每立米水的费用 Y(u,V0)约 0.065(英镑 )
虽然 0.065英镑略低于淡化海水的成本 0.1英镑,
但是模型假设和构造非常简化与粗糙。