双参数法模型
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
分布参数法建模
前面建立的模型都用了考察对象在系统中的均匀分布假设。这种方法建模被称为集中参数法。
考虑个体差异(或分布差异)的建模方法被称为分布参数法。分布参数法用于连续变量的问题时,得到的通常都是偏微分方程,无论建模还是求解都比较困难。仅举两个简单例子,来说明这种方法的应用。
例8 人口问题的偏微分方程模型
人有年龄、性别等区别,本例中考虑到这些因素,用分布参数法来建立人口问题的数学模型。
令p (t ,x )为t 时刻年龄为x 的人口密度,则t 时人口总数为:
A 为人的最大寿命。
设t 时刻年龄为x 的人的死亡率为d (t ,x ),则有:
=dt ,由上式可导出:
3.38)初始条件: P (0,x )=P 0(x )(3.39)边界条件: 3.40)
k (t ,x )女性性别比b (t ,x )女性生育率[x 1,x 2]妇女生育期0()(,)A
P t p t x dx =⎰(,)(,)(,)(,)p t dt x dx p t x dt dx d t x dt p t x dxdt +--=--(,)(,)p p d t x p t x t x
∂∂+=-∂∂2
1(,0)(,)(,)(,)x x P t b t x k t x p t x dx =⎰
对(3.38)式关于x 从0到A 积分,得:
(t )、D (t )分别为t 时刻的生育率和死亡率。则有:B (t )、D (t )与t 无关,则可得:Malthus 模型
2100(,0)(,)(,)(,)(,)(,)(,)(,)A x A
x dP P t d t x p t x dx dt
b t x k t x p t x dx d t x p t x dx =-=-⎰⎰⎰00(,)(,)(,)()()
(,)(,)()()
A A b t x k t x p t x dx
B t P t d t x p t x dx D t P t ==⎰⎰(()())()dP B t D t P t dt =-0
()()(0)dP B D P t dt P P ⎧=-⎪⎨⎪=⎩
例9 交通流问题
问题的两个角度:
司机或旅客安全、快速地到达目的地
交通管理部门尽可能多的人安全地通过集中参数法:
假设车流量是均匀分布
目标使车流密度保持在安全的范围之内,让司机尽可能开得快些即可,必要时司机自己会刹车。
现实生活中可能吗?
车流密度和车速不可能是常数
分布参数法:
x 轴表示公路,x 轴正向表示车流方向。
如果采用连续模型,设u (t ,x )为时刻t 时车辆按x 方向分布的密度,再设q (t ,x )为车辆通过x 点的流通率。车辆数守恒,有:
(,)(,)0u q t x t x t x
∂∂+=∂∂(3.41)由于安全上的原因,q 是u 的函数,该函数关系称为基本方程或结构方程。
(,)(,)(,)(,)u t dt x dx u t x dx q t x dt q t x dx dt
+-=-+
利用经验公式导出基本方程。q
0u
u m u j 图3-28图3-28是根据美国公路上的车辆情况而统计出来的曲线,其中u 的单位是车辆数/每英里,q 的单位为车辆数/每小时。图中可以看出:
(2)u 增大到一定程度(达到u m )时,q 达到最大;u 继续增大时,车辆流q 将减小,这表示车辆密度太大反而会影响车辆率,使之下降,(出现堵塞)。
(1)当u 的值较小时,公路利用率较低,q 较小(u =0时公路是空置的,车辆率q 为零);随着u 的增大,公路利用率逐渐提高,q 逐渐增大。
根据美国公路实际统计:
当u ≈75辆/每英里可达到最大车辆流
当u ≈225辆/英里时,q ≈0,即堵塞。
根据图3-28中曲线的特征,可用多种函数来拟合q =q (u )。u f 为自由速度,u j 为出现完全堵塞时的车流密度。
Greenshields 用二次函数来拟合。
他令:(1/)f j q u u u u =-0≤u ≤u j
u m =u j /2,q m =u f u m
/2有:将Greenshields 的基本方程代入(3.41),利用复合函数求导法则并注意到u f 、u j 均为常数,可得:2(,)()(,)0f f j u u u u t x u t x t u x
∂∂+-=∂∂令,方程可简化为:2f
f j u h u u u =-(,)(,)0h h t x h t x t x ∂∂+=∂∂初值条件:02(0,)()f
f j u h x u u x u =-