1、,3.3 平衡原理 与机理模型,3.3 平衡原理与机理模型,一. 平衡原理 自然界任何物质在其运动变化过程中 一定受到某种平衡关系的支配。二. 机理模型 在一定的假设下， 根据主要因素相互作用的机理， 对它们之间的平衡关系的数学描述,关于假设 是对实际问题的抽象、化简和规范 是组建数学模型的基础和前提。 假设较强，模型简单易分析和操作 但与实际差距大。 假设较弱，模型比较接近实际。 但模型复杂不易操作。,关于平衡关系1.平衡关系是数学模型的核心，建模的关键。2.有些平衡关系是明显的。 有些平衡关系隐藏在问题的背后 需要在化简之后逐渐明确出来。3.有些平衡关系本身就直接构成了模型。 有些平衡关系
2、还需要经过数学上 的加工整理才能得到理想的模型。,三. 模型举例例1 买房贷款：银行可以向购房人提供个人住房贷款的业务。偿还贷款时要求借款人在借款期间内每月以相等的月均还款额偿还银行贷款本金和利息。试组建计算月均还款额的数学模型。 假设: 1. 每月月底还款； 2. 每月还款金额相等； 3. 按月计算利息； 4. 到期欠款全部还清。,参量、变量 贷款额：A(万元)， 贷款期限：N年(n=12N月) , 月利率：r， 月均还款额：x。,模型平衡关系： 相邻两月欠款余额之的关系本月月底还款后的欠款余额等于上月欠款余额的本利和扣除月还款后的金额。 令Ck表示第k月月底还款后的欠款余额，则有 Ck=
3、(1+r)Ck-1-x,令 C*= (1+r)C* - x, 则有 C*= x/r. 可得 Ck-C* = (1+r) (Ck-1-C*) = (1+r)k(C0-C*)利用 C0=A, Cn=0 和 C*=x/r 可得 (1+r)n(A-x/r)-x/r=0 (1+r)nA-(1+r)n-1/rx = 0,例 : 一位教师筹措月利息为0.3675%的15年公积金贷款20万元买房。1. 每月等额偿还贷款。计算它的月均还款额。2. 他已经还款20个月，想知道他还欠多少贷款？,月均还款额,k月后的欠款余额,练习：一项买房贷款期限15年，每月支付0.51%利息的20万元商品房贷款，每月应还款额是多少
4、？1 试组建计算等额本息（等额月均）还款的数学模型。2 如果每月偿还当月利息后，再以每月相等的额度偿还本金（等额本金还款）。 试组建等额本金还款的数学模型。 比较 2种还款方式的优劣。,问题P88 15，14，13,例2. 兔子的繁殖 I由一对兔子开始，一年可以繁殖成多少对兔子？假设兔子的生殖力是这样的：一对兔子每一个月可以生一对兔子，兔子在出生两个月以后就具有繁殖后代的能力。假设： 1. 每对兔子每一个月定生一对兔子。 2. 兔子出生两个月后都具有繁殖能力。 3. 兔子每经过一个月底就增加一个月令。 4. 兔子不离开群体（不考虑死亡）。,变量、参量 月份：n，幼兔：a0(n)，成兔: a1(
5、n) 平衡关系： 本月初(一月龄)的幼兔是上月成兔繁殖的后代。 本月的成兔是上月的成兔和上个月(一月龄)的幼兔发育结果的总和。 模型 I a0(n) = a1(n-1) a1(n) = a0(n-1) + a1(n-1) a0(1)=1 a1(n) = 0 目标：求到第12个月底兔子总数 a(12)= a0(12) + a1(12),相邻两月兔群数量的变化,1. 模拟. a0(1)=1, a1(1)=0 n 1 2 3 4 5 6 7 8 9 10 11 12 a0(n) 1 0 1 1 2 3 5 8 13 21 34 55 a1(n) 0 1 1 2 3 5 8 13 21 34 55 8
6、9 a (n) 1 1 2 3 5 8 13 21 34 55 89 144 a(12)=144 斐波那契数列(黄金数) a(n+1) = a(n) + a(n-1),2. 证明斐波那契数列(黄金数) 的性质 a(n+1) = a(n) + a(n-1) 因为 a0(n+1)=a1(n) a1(n+1)=a0 (n)+a1 (n)=a(n) 所以 a (n+1) = a (n)+a1 (n) = a(n) +a(n-1) 于是 a(n) /a(n+1)= 1/a(n-1)/a(n)+1 记 xn= a(n) /a(n+1), 则 xn=1/1+xn-1 可以证明xn 收敛，记其极限为x0， 则
7、由 x0=1/1+x0 得到 x0=0.618,3. 模型的作用机理 令 a(n) = (a0(n), a1(n), 则模型为 a(n) = A a(n-1) a11幼兔的繁殖能力， a12成兔的繁殖能力， a21幼兔的发育为成兔的比例， a22成兔存活的比例。,4. 群体的渐近性质 A 有主特征值 =1.618 相应的右特征向量 L=(0.382 0.618)， 于是，当n 时, a(n)/a(n) L群体的渐进增长速率为1.618幼兔将占群体总数的 0.382成兔将占群体总数的 0.618称向量 L为种群的稳定分布。问题：如果每对兔子每月可生两对兔子，求这个种群的稳定分布。,兔子的繁殖II
8、 一对兔子每月可生一对幼兔，幼兔出生二个月后就具有繁殖能力，三个月后就离开群体。问一对幼兔一年后繁殖的群体多大？求这个种群的稳定分布。,假设 1. 2. 3 同上 假设4. 兔子在三个月 后生完一对幼兔就离开群体。参量、变量 月份: n, 幼兔: a0(n), 成兔: a1(n), 老兔: a2(n)平衡关系 本月初的幼兔是上月成兔老兔繁殖的后代。 本月初的成兔是上月幼兔发育的结果。 本月初的老兔是上月成兔发育的结果。,模型 a0(n+1)=a1(n)+a2(n) a1(n+1)=a0(n) a2(n+1)=a1(n)令 a(n) = (a0(n), a1(n), a2(n), 则 a(n)
9、= A a(n-1)其中,分析 1. 模拟. a0(1)=1, a1(1)=0, a2(1)=0 n 1 2 3 4 5 6 7 8 9 10 11 12a0(n) 1 0 1 1 1 2 2 3 4 5 7 9a1(n) 0 1 0 1 1 1 2 2 3 4 5 7a2(n) 0 0 1 0 1 1 1 2 2 3 4 5a(n) 1 1 2 2 3 4 5 7 9 12 16 21 2. 证明 a(n+1) = a(n-1) + a(n-2) a(n+1) = a(n) + a(n-4) Padovan 数列(塑料数),矩阵A的特征方程为 31=0有特征值（渐进增长率）=1.3247特征
10、向量（稳定的年龄结构） （0.4302，0.3247，0.2451）,离散模型（差分方程）的组建利用平衡原理，找出每一步对前一步或前几步的依赖关系，得到以差分方程的形式描述的数学模型。,例3. 人口的自然增长. 建模描述一个地区内人口的自然增殖的过程。即只考虑由于人口的生育和死亡所引起的人群数量变化的过程。 令N(t)表示t时刻的人口数。 假设1. 人群个体同质。 N(t) 连续可微. 假设2. 群体规模大。,平衡关系：人口数在区间t，t+t内的改变量等于这段时间内出生的个体数与死亡的个体数之差。假设3. 群体封闭，只考虑生育和死亡对人口的影响。令B(t, t, N), D(t, t, N)
11、分别表示生育数和死亡数, 则有,假设4. 从大群体的平均效应考虑生育和死亡对人口的影响（生育率和死亡率） 生育率 b(t,t, N)= B(t,t, N)/ N 死亡率 d(t,t, N)= D(t,t, N)/N则有,由于R(t, t,N)|t=0=0，将R(t,t,N) 关于t展开,令 t0 取极限可得,假设5. 群体增长恒定. 则 r(t, N) = r( N),假设6. 个体增长独立. 则 r( N) = r.,在离散时间点k=0, 1, 2, , 上有 N(k+1) = erN(k) = N(k)Marthus 模型“若我的两个假设是成立的,那么,我认为人口繁殖的能量是无限地大于自然
12、界为人类提供资料的能量的。人口如果不受控制，它会以几何比率增长。而生活资料只能以算术比率增长。” Marthus:论人口原理(1798年),给定初值 N(0)=N0，可得,假设1. 人群个体同质。假设2. 群体规模大。假设3. 群体封闭，只考虑生育和死亡对人口的影响。假设4. 从大群体的平均效应考虑生育和死亡对人口的影响。（生育率和死亡率）假设5. 群体增长恒定。假设6. 个体增长独立。,模型的讨论 1.作为人口自然增长， 模型与实际是不一致的。 许多地区人口增长不符合这个模型。 2.只考虑增长和衰减时， 模型是正确的。 3.模型是可以改进的。,草履虫实验室种群的动态(Gause, 1934)
13、,酵母细胞20小时内的动态(Pearl),例4. Logistic 模型 模型: 在有限的资源内生物种群不可能无限增长. 存在有饱和水平, 种群增加接近饱和时增长速度减慢而趋于零. 假设 60 需要修改: r = r(N), r(N) K, N0 K : 环境承载力或饱和水平. K 时，模型退化为 Malthus 模型， r ：内禀增长率 20. 模型的解：分离变量法,递增，有极限 K，依赖于三个参数。,30. 模型的动态特征 平衡解：模型不依赖于时间的解 N0*=0，N1*=K。 平衡解的稳定性：,N0* 不稳定，N1* 稳定,例5 池水含盐 问题 池中有一定体积的盐水，从池的上部向池中注入
14、一定浓度的盐水混合后的盐水将从池的下部流出。建模描述池中盐水浓度的动态。假设： 1. 盐水注入池中后迅速混合 2. 池中盐水浓度均匀。,平衡关系 在时间段t+t内, 池中盐水体积的改变量等于这段时间内流入盐水的体积与流出盐水体积之差； 在时间段t+t内, 池中(纯)盐的改变量等于这段时间内流入的(纯)盐的量与流出的(纯)盐的量之差。变量、参量： 池中盐水体积 V(t), 池中盐水浓度 p(t); 流入盐水速度 rI(t), 流入盐水浓度 pI(t); 流出盐水速度rO(t), 流出盐水浓度 p(t).,模型池中盐水的改变量 V(t+t)-V(t) 流入盐水量流出盐水量池中盐的改变量 p(t+
15、t)V(t+ t)-p(t)V(t) 流入盐量 流出盐量,利用积分中值定理可得,类似地有,模型,连续模型（微分方程）组建的微元法 在自变量的微小的区间内以简单的形式描述有关变量之间的平衡关系, 再利用微积分学的思想进一步处理它,得到以微分方程的形式描述的数学模型。,例5. 录音机的运行 建模分析磁带录音机的运行规律(计数器的读数与运行时间的关系)。 数据：I. 读数与时间t 1 2 3 4 5 10 15 20 25 30 31 - n 9 18 28 37 47 97 151 211 280 362 382 385 数据：II. 读数与转数k 2 4 10 14 18 22 26 31 35
16、 41 60 1 2 5 7 9 11 13 15 17 20 29,背景 1. 磁带盒内有二个磁带轮：送带轮和收带轮放音时送带轮上的磁带减少，缠于收带轮上 2. 计数器只记录某个磁带轮转动的情况。计数器的读数不刚好是磁带轮的转数。 3. 磁带轮在放音时转动不是匀速的，送带轮加速，收带轮减速。,假设 1. 计数器记录了送带轮的转数。 2. 计数器的读数与送带轮的转数成正比。 3. 磁带运行的线速度定常。 4. 磁带厚度均匀，缠绕松紧一致，无空隙 5. 磁带缠绕一圈的周长等于缠绕的圆周长参量、变量 读数：n，带轮转数：k, 运行时间：t(k) 磁带厚度：d，带芯轮半径：r， 磁带速度：v, 磁带
17、最多圈数：N 第 k 圈磁带的半径：Rk， 第k圈磁带长度：Lk,平衡关系 运行k圈磁带的时间等于磁带的长度与运行速度之商。模型：t=0 时 n=0，送带轮缠满磁带并开始转动 由假设1，送带轮计数从外圈数起。 由假设5，Lk=2Rk 由假设4，Rk=r+(N-k+1)d,由假设2, k = c n则有得模型t(n) = a n + b n2, 其中,参数 a, b, c 的估计 1. 最小二乘法估计 a, b 正规方程组,方程 290942 a + 83322472 b=25827 83322472 a + 2.591711010 b=7237447有解 a = 0.11095，b = -7.
18、744510-5模型 t(n) = 0.11095 n 7.4475 10-5 n2检验 n 9 18 28 37 47 97 151 211 280 t 1 2 3 4 5 10 15 20 25 t .99 1.97 3.05 4.00 5.04 10.03 14.99 19.96 24.99,分析 根据假设2：k = c n, 利用数据II可以给出参数 c 的最小二乘估计。可得 c = 2.04。又可测得 r = 1.1 cm，N=3852.04由可以求出d = 0.001628cm，v =2.75m/min，L=85.25m。,问题：大江截流 1997年11月8日电视正在播放十分壮观的
19、长江三峡工程大江截流的实况。截流8：55开始，水面宽40 m，水深60m。 到11：50时， 水面宽34.4m 到13：00时， 水面宽31 m。,这时电视机旁的小明说，现在可以估算下午几点合龙。8：55到11：50，进展的速度为40-34.4=5.6m，平均每小时宽度减少1.9m。从11：50到13：00，宽度减少34.4-31=3.4m，平均每小时减少2.9m小明认为回填速度是越来越快的，近似地每小时速度加快1m。从下午1：00起，大约要5个多小时，即到下午6点多才能合龙。,但到了下午3点28分，电视里传来了振奋人心的消息：大江截流成功！小明后来想明白了，他估算的方法不好现在请你根据上面的数据设计一种合理的估算方法（建立一种合理的数学模型）进行估算，使你的计算结果更切合实际。,
