平衡原理及机理模型.doc

§3.3 平衡原理与机理模型 一. 平衡原理 自然界任何物质在其运动变化过程中一定受到某种平衡关系的支配。 二. 机理模型 在一定的假设下，根据主要因素相互作用的机理，对它们之间的平衡关系的数学描述。 三. 连续模型 连续模型组建的微元法 在自变量的微小的区间内以简单的形式描述有关变量之间的平衡关系, 再利用微分学的思想进一步处理它,得到以微分方程的形式描述的数学模型。 例1. 人口的自然增长. 建模描述一个地区内人口的自然增殖的过程。即考虑由于人口的生育和死亡所引起的人群数量变化的过程。 假设1. 人群个体同质。令N(t)表示t时刻的人口数。 假设2. 群体规模大。 N(t) 连续可微. 假设3. 群体封闭，只考虑生育和死亡对人口的影响。 平衡关系：人口数在区间[t，t+ (t ]内的改变量等于这段时间内出生的个体数与死亡的个体数之差。 令B(t, (t, N), D(t, (t , N) 分别表示在时间区间[t，t+ (t ]内生育数和死亡数, 则有 N(t+(t)-N(t)=B(t, ( t,N)-D(t, ( t,N) 假设4. 从大群体的平均效应考虑生育和死亡对人口的影响。（生育率和死亡率） 生育率b(t, (t, N) = B(t, (t, N)/N, 死亡率d(t, (t, N) = D(t, (t, N)/N 记增长率为 R(t, ( t,N)= b(t, ( t,N)-d(t, ( t,N) 则有 N(t+(t)-N(t)=R(t, ( t,N)N 将R(t, (t,N)关于(t展开. 由于R(t, h, N)|h=0=0，所以 N(t+ (t ) - N(t)=r(t,N)N (t +o((t). 两边除以(t, 并令(t →0, 得到 dN/dt=r(t, N)N 假设5. 群体增长恒定。（r与 t 无关） dN/dt=r(N) N 假设6. 个体增长独立。（r 与 N 无关） dN/dt=r N 给定初值 N(0)=N0，可得人口增长的指数模型(Maithus 模型) N(t)=N0ert 在离散时间点k=0, 1, 2, …, 上有 N(k+1) = er N(k ) Maithus: “若我的两个假设是成立的,那么,我认为人口繁殖的能量是无限地大于自然界为人类提供资料的能量的。人口如果不受控制，它会以几何比率增长。而生活资料只能以算术比率增长。只要稍微看一下数字，就将明确第一种能量比之第二种能量是无比巨大的。” 《论人口原理》 总结对人口指数增长模型的假设, 1. 人群个体同质。 2. 群体规模大。 3. 群体封闭，只考虑生育和死亡对人口的 影响。 4. 从大群体的平均效应考虑生育和死亡对人口的影响。（生育率和死亡率） 5. 群体增长恒定。 6. 个体增长独立。 由这些假设可分析这个模型的作用. 例2 池水含盐 池中有一定体积的盐水，从池的上部向池中注入一定浓度的盐水。混合后的盐水将从池的下部流出。建模描述池中盐水浓度的动态。 假设： 盐水注入池中后迅速混合, 使得盐水浓度均匀。 变量、参量： 池中盐水体积 V(t), 池中盐水浓度 p(t); 池中原有盐水体积 V0, 原有盐水浓度p0; 流入盐水速度 rI(t), 流入盐水浓度 pI(t); 流出盐水速度rO(t), 流出盐水浓度 p(t). 平衡关系 在时间段[t+ ( t]内, 池中(纯)盐的改变量=这段时间内流入的(纯)盐的量与流出的(纯)盐的量之差。 池中盐的改变量： p(t+(t)V(t+ ( t)-p(t)V(t) 流入盐量： 流出盐量： 利用积分中值定理可得 在时间段[t+ ( t]内, 池中盐水体积的改变量=这段时间内流入盐水的体积与流出盐水体积之差； 令 ( t (0 得 模型 进一步问题：池中有水 2000 m3，含盐 2 kg，以 6m3 / 分 的速率向池中注入浓度为 0.5 kg / m3 的盐水，又以 4 m3 / 分的速率从池中流出混合后的盐水。问欲使池中盐水浓度达到 0.2 kg / m3，需要多长时间？ 此时 V(t)=2000+2*t. dp/dt=3/V(t)-6*p(t)/V(t), p(0)=0.001. 用MATLAB求 p(t) 求表达式（符号运算） S=dsolve(‘Dx=(3-6*x)/(2000+2*t)’); 求数值解 建立M文件 fun . M, function y=fun(t,x) y=(3-6*x)