例4. 兔子的繁殖 由一对兔子开始，一年可以繁殖成多少对兔子？ 假设兔子的生殖力是这样的：一对兔子每一个月可以生一对兔子，

Presentation on theme: "例4. 兔子的繁殖 由一对兔子开始，一年可以繁殖成多少对兔子？ 假设兔子的生殖力是这样的：一对兔子每一个月可以生一对兔子，"— Presentation transcript:

1 例4. 兔子的繁殖 由一对兔子开始，一年可以繁殖成多少对兔子？ 假设兔子的生殖力是这样的：一对兔子每一个月可以生一对兔子， 并且兔子在出生两个月以后就具有繁殖后代的能力。 n a(n) a(n+1) = a(n) + a(n-1) 斐波那契数列(黄金数)

2 假设： 1. 每对兔子每一个月定生一对兔子。 2. 兔子出生两个月后都具有繁殖能力。 3. 兔子每经过一个月底就增加一个月令。 平衡关系 本月初(一月令)的幼兔是上月成兔繁殖的后代。 本月的成兔是上月的成兔和上个月(一月令)的幼兔发育结果的总和。

3 月份：n，幼兔：a0(n)，成兔: a1(n)
变量、参量 月份：n，幼兔：a0(n)，成兔: a1(n) 模型 a0(n) = a1(n-1) a1(n) = a0(n-1) + a1(n-1) 令 a(n) = (a0(n), a1(n))’, 则 a(n) = A a(n-1) 其中

4 分析 1. 模拟. a0(1)=1, a1(1)=0 n a0(n) a1(n) a(n) 2. 证明 a(n+1) = a(n) + a(n-1) a0(n+1)=a1(n) a1(n+1)=a(n)

5 a(n)=Ana(0)=RλnL’a(0)+o(λn)
3. 模型的作用机理 a11幼兔的繁殖能力， a12成兔的繁殖能力， a21幼兔的发育为成兔的比例， a22成兔存活的比例。 4. 群体的渐近性质 A 有主特征值 λ=1.618 和右、左特征向量 R= L=( )’,. a(n)=Ana(0)=RλnL’a(0)+o(λn)

6 发展 一对兔子每月可生一对幼兔， 幼兔出生二个月后就具有繁殖能力， 三个月后就离开群体。 问一对幼兔一年后繁殖的群体多大？ 假设： 参量、变量 月份: n, 幼兔: a0(n), 成兔: a1(n), 老兔: a2(n)

7 平衡关系 本月初的幼兔是上月成兔老兔繁殖的后代。 本月初的成兔是上月幼兔发育的结果 本月初的老兔是上月成兔发育的结果

8 模型 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) = A a(n-1) 其中

9 分析 1. 模拟. a0(1)=1, a1(1)=0, a2(1)=0 n a0(n) a1(n) a2(n) a(n) 2. 证明 a(n+1) = a(n-1) + a(n-2) a(n+1) = a(n) + a(n-4) Padovan 数列(塑料数)

10 例5. 录音机的运行 建模分析磁带录音机的运行规律(计数器的读数与运行时间的关系)。 数据：I. 读数与时间 t n 数据：II. 读数与转数 　k ｎ 　　

11

12

13 背景 1. 磁带盒内有二个磁带轮：送带轮和收带轮。放音时送带轮上的磁带减少，缠于收带轮上。 2. 计数器只记录某个磁带轮转动的转动情况。 计数器的读数不刚好是磁带轮的转数。 3. 磁带轮在放音时转动不是匀速的，送带轮加速，收带轮减速。

14 假设 1. 计数器记录了送带轮的转数。 2. 计数器的读数与送带轮的转数成正比。 3. 磁带运行的线速度定常。 4. 磁带厚度均匀，缠绕松紧一致，无空隙。 5. 磁带缠绕一圈的周长等于缠绕的圆周长。

15 参量、变量 计数器读数：n，带轮转数：k, 运行时间：t(k)， 磁带厚度：d，带芯轮半径：r，磁带速度：v, 磁带最多圈数：N 第 k 圈磁带的半径：Rk，第圈磁带长度：Lk

16 平衡关系 运行k圈磁带的时间等于磁带的长度与运行速度之商。 模型：t=0时n=0，送带轮缠满磁带并开始转动。 由假设1，送带轮计数从外圈数起。 由假设5，Lk=2πRk 由假设4，Rk=r+(N-k+1)d

17 平衡原理与机理模型

18 由假设2, k = c n．则有 得模型　t(n) = a n + b n2, 其中

19 参数 a, b, c 的估计 1. 最小二乘法估计 a, b 正规方程组

20 平衡原理与机理模型 方程 a b=25827 a ·1010 b= 有解 a = ，b = ·10-5 模型 t(n) = n – · 10-5 n2 检验 n t t

21 分析 根据假设2：k = c n, 利用数据II可以给出参数 c 的最小二乘估计。 可得 c = 2.04。又可测得 r = 1.1 cm，N=385×2.04 由 可以求出d = cm，v = 2.75m/min，L=85.25m。

22 问题一. 证明PADOVAN数列的递推公式。
问题二. 根据录音机运行的数学模型及观测数据 I 给出模型参数最小二乘估计的正规方程组。 问题三. 在数据 I、II 的基础上，使用MATLAB 给出录音机运行模型的数值分析。