|
微分方程在数学建模中的应用举例 数学建模是数学方法解决各种实际问题的桥梁,随着计算机技术的快速发展,数学的应用日益广泛,数学建模的作用越来越重要,而且已经应用到各个领域。用微分方程解决实际问题的关键是建立实际问题的数学模型——微分方程。这首先要根据实际问题所提供的条件,选择确定模型的变量,再根据有关学科,如物理、化学、生物、经济等学科理论,找到这些变量遵循的规律,用微分方程的形式将其表示出来。9 W5 z2 d- X% R7 c0 ?- g+ \& z
一、交通红绿灯模型
2 q0 J! H! q- t 在十字路口的交通管理中,亮红灯之前,要亮一段时间的黄灯,这是为了让那些正行驶在十字路口的人注意,告诉他们红灯即将亮起,假如你能够停住,应当马上刹车,以免冲红灯违反交通规则。这里我们不妨想一下:黄灯应当亮多久才比较合适?
$ r4 b: `2 D5 M C# d2 }. ~4 a: W1 E 停车线的确定,要确定停车线位置应当考虑到两点:一是驾驶员看到黄灯并决定停车需要一段反应时间 ,在这段时间里,驾驶员尚未刹车。二是驾驶员刹车后,车还需要继续行驶一段距离,我们把这段距离称为刹车距离。驾驶员的反应时间(实际为平均反应时间) 较易得到,可以根据经验或者统计数据求出,交通部门对驾驶员也有一个统一的要求(在考驾照时都必须经过测试)。例如,不失一般性,我们可以假设它为1秒,(反应时间的长短并不影响到计算方法)。停车时,驾驶员踩动刹车踏板产生一种摩擦力,该摩擦力使汽车减速并最终停下。设汽车质量为m,刹车摩擦系数为f,x(t)为刹车后在t时刻内行驶的距离,更久刹车规律,可假设刹车制动力为fmg(g为重力加速度)。由牛顿第二定律,刹车过程中车辆应满足下列运动方程:
$ l" s0 @- J. b8 F2 Q0 m md2xdt2=-fmg
# D1 V+ O* X9 Z2 d8 H x(0)=0, dxdtt=0=v0
7 k( g. k7 ^( g6 I (1), t: S; l+ z) @+ s3 T, ]$ |
在方程(1)两边同除以 并积分一次,并注意到当t=0时dxdt=V0,得到
/ T$ s5 L/ x$ c% `1 o, R% ~ dxdt=-fgt+v0! L; J5 n9 y( x" c. i6 W/ E8 c/ ]
(2)( v+ g% E* R- D9 X h1 h; ?% h* ?
刹车时间t2可这样求得,当t=t2时,dxdt=0,故1 i g! h1 V0 R2 j7 V+ J1 L2 U: d& j
t2=v0fg) `/ u: L9 [) t& H/ o
将(2)再积分一次,得
$ R/ t x. ?6 Q1 D x(t)=-12fgt2+v0t7 b/ c6 d, m8 E- H
将t2=v0fg代入,即可求得停车距离为6 s0 F7 I* Z# k. ]3 `% A* D' u. B
x(t2)=1v202fg& b) c) \6 x( E( K& V, _6 a* }
据此可知,停车线到路口的距离应为:; j: ~/ `# S, K9 x) [0 `- X2 f
L=v0t1+12v20fg* f% n$ X ^; }
等式右边的第一项为反应时间里驶过的路程,第二项为刹车距离。
) W0 L5 @2 ~ c; l: Z& n' o 黄灯时间的计算,现在我们可以来确定黄灯究竟应当亮多久了。在黄灯转为红灯的这段时间里,应当能保证已经过线的车辆顺利地通过街口,记街道的宽度为D(D很容易测得),平均车身长度为 ,这些车辆应通过的路程最长可达到L+D+l,因而,为保证过线的车辆全部顺利通过,黄灯持续时间至少应当为:
2 B; q4 b5 ~8 Q0 Y& C T=L+D+lv0
# R8 e" c1 w! |/ } 二、市场价格调整模型6 s- S$ n5 N% J
对于纯粹的市场经济来说,商品市场价格取决于市场供需之间的关系,市场价格能促使商品的供给与需求相等这样的价格称为(静态)均衡价格。也就是说,如果不考虑商品价格形成的动态过程,那么商品的市场价格应能保证市场的供需平衡,但是,实际的市场价格不会恰好等于均衡价格,而且价格也不会是静态的,应是随时间不断变化的动态过程。 X& v! V. R7 P0 a/ k( x7 m5 a# Y
如果设某商品在时刻t的售价为P,社会对该商品的需求量和供给量分别是P的函数D(P),S(P),则在时刻t的价格p(t)对于时间t的变化率可认为与该商品在同时刻的超额需求量D(P)-S(P)成正比,即有微分方程
7 E1 }: I" [0 M# d2 Z& C dPdt=k[D(P)-S(P)] (k>0)
& Q; F, S6 [) R' [: ^ (3)) f' T, U- Z8 K: w$ s6 t( A: L! j
在D(P)和S(P)确定情况下,可解出本文由论文联盟http://wWw.LWlm.cOm收集整理价格与t的函数关系,这就是商品的价格调整模型。7 n: I, R& G! h/ @9 G
某种商品的价格变化主要服从市场供求关系。一般情况下,商品供给量 是价格 的单调递增函数,商品需求量Q是价格P的单调递减函数,为简单起见,分别设该商品的供给函数与需求函数分别为& c# V% z8 k3 @9 |
S(P)=a+bP,Q(p)=α-βP' }$ T5 x {4 `$ o- C! `8 Q
(4)/ r r1 E. E& t! E# N
其中a,d,α,β均为常数,且b>0,β>0。4 g0 W# y9 `6 c% \+ v3 m$ J7 F
当供给量与需求量相等时, 由(4)可得供求平衡时的价格
4 B& E. m5 q9 M! ~2 G Pe=α-aβ+b5 c2 H+ D& X, c0 c
并称Pe为均衡价格。% N; ~4 x- }/ V
一般地说,当某种商品供不应求,即S<Q时,该商品价格要涨,当供大于求,即S>Q时,该商品价格要落。因此,假设t时刻的价格P(t)的变化率与超额需求量Q-S成正比,于是有方程 dPdt=k[Q(P)-S(P)]
( V! j0 o2 Y7 U6 f& x3 ^4 g$ h 其中k>0,用来反映价格的调整速度。! J4 ^1 W4 g& G7 b a7 z5 P8 Y
将(4)代入方程,可得
# q$ i& y& I4 p dPdt=λ(pe-P)
' {: | m+ D, \+ e% g (5)7 n4 T* y. Z* y3 h/ H4 F5 E
其中常数λ=(b+β)k>0,方程(5)的通解为2 X8 d; |0 L0 N. w* F* `# b
P(t)=Pe+Ce-λt t' ~( k+ p' {( `$ |
假设初始价格P(0)=P0,代入上式,得C=P0-Pe,于是上述价格调整模型的解为
: [1 G7 S4 a4 t3 A6 Y7 F2 e$ \* f P(t)=Pe+(P0-Pe)eλt
: Z2 Q* \4 D! V8 a' H. W 由于λ>0知,t→+∞时,P(t)→Pe。
+ B1 j$ r: U5 o+ l: b! n 说明随着时间不断推延,实际价格P(t)将逐渐趋近均衡价格Pe。这符合我们实际生活中具体事实。
, F, f9 ]) c: H% u8 b$ S. C g1 `3 z/ B3 n. m& @* G
|