为了保持自然资料的合理开发与利用，人类必须保持并控制生态平衡，甚至必须控制人类自身的增长。
本节将建立几个简单的单种群增长模型，以简略分析一下这方面的问题。一般生态系统的分析可以通过一些简单模型的复合来研究，大家若有兴趣可以根据生态系统的特征自行建立相应的模型。 美丽的大自然种群的数量本应取离散值，但由于种群数量一般较大，为建立微分方程模型，可将种群数量看作连续变量，甚至允许它为可微变量，
由此引起的误差将是十分微小的。
离散化为连续，方便研究
§ 3.2 Malthus模型与 Logistic模型模型 1 马尔萨斯（ Malthus）模型马尔萨斯在分析人口出生与死亡情况的资料后发现，
人口净增长率 r基本上是一常数，（ r=b-d,b为出生率，
d为死亡率）,既：
1 dN r
N dt?
dN rN
dt?
或 （ 3.5）
0()0() r t tN t N e （ 3.6）
（ 3.1） 的解为：
其中 N0=N(t0)为初始时刻 t0时的种群数。
马尔萨斯模型的一个显著特点,种群数量翻一番所需的时间是固定的 。
令种群数量翻一番所需的时间为 T，则有：
002 rTN N e?
ln 2T
r?
故模型检验比较历年的人口统计资料，可发现人口增长的实际情况与马尔萨斯模型的预报结果基本相符，例如,1961年世界人口数为 30.6 （即 3.06× 109），人口增长率约为 2%，人口数大约每 35年增加一倍。检查 1700年至 1961的 260年人口实际数量，发现两者几乎完全一致，且按马氏模型计算，人口数量每 34.6年增加一倍，两者也几乎相同。
1950 2000 2050 2100 2150 2200
0
0,5
1
1,5
2
2,5
3
3,5
x 1 0
11
t/ 年
N/
人马尔萨斯模型人口预测模型预测假如人口数真能保持每 34.6年增加一倍，那么人口数将以几何级数的方式增长。例如，到 2510年，人口达 2× 1014个，
即使海洋全部变成陆地，每人也只有 9.3平方英尺的活动范围，
而到 2670年，人口达 36× 1015个，只好一个人站在另一人的肩上排成二层了。 故 马尔萨斯模型是不完善的。
几何级数的增长
Malthus模型 实际上只有在群体总数不太大时才合理，到总数增大时，
生物群体的各成员之间由于有限的生存空间，有限的自然资源及食物等原因，就可能发生生存竞争等现象。
所以 Malthus模型假设的人口 净增长率不可能始终保持常数，
它应当与人口数量有关。
模型 2 Logistic模型人口净增长率应当与人口数量有关，即,r=r(N)
从而有：
()dN r N Ndt?
（ 3.7）
r(N)是未知函数，但根据实际背景，它无法用拟合方法来求 。
为了得出一个有实际意义的模型，我们不妨采用一下工程师原则。工程师们在建立实际问题的数学模型时，总是采用尽可能简单的方法。
r(N)最简单的形式是常数，此时得到的就是马尔萨斯模型。
对马尔萨斯模型的最简单的改进就是引进一次项（竞争项）
对马尔萨斯模型引入一次项（竞争项），令 r(N)=r-aN
此时得到微分方程：
()dN r a N Ndt (1 )d N NrNd t K或 （ 3.8）
（ 3.8） 被称为 Logistic模型或生物总数增长的统计筹算律，是由荷兰数学生物学家弗赫斯特（ Verhulst）首先提出的。一次项系数是负的，因为当种群数量很大时，会对自身增大产生抑制性，故一次项又被称为竞争项。
（ 3.8） 可改写成：
()dN k K N Ndt （ 3.9）
(3.9)式还有另一解释，由于空间和资源都是有限的，不可能供养无限增长的种群个体，当种群数量过多时，由于人均资源占有率的下降及环境恶化、疾病增多等原因，出生率将降低而死亡率却会提高。设环境能供养的种群数量的上界为 K（近似地将 K看成常数）,N表示当前的种群数量，
K-N恰为环境还能供养的种群数量，（ 3.9）指出，种群增长率与两者的乘积成正比，正好符合统计规律，得到了实验结果的支持，这就是（ 3.9）
也被称为统计筹算律的原因。
图 3-5
对 （ 3.9） 分离变量：
11 d N k K d t
N K N
两边积分并整理得：
1 kK t
KN
Ce
令 N(0)=N0，求得：
0
0
KNC
N
故 （ 3.9） 的满足初始条件 N(0)=N0的解为：
0
00
() () k KtNKNt N K N e （ 3.10）
易见：
N(0)=N0,lim ( )
t N t K
N(t)的图形请看图 3.5
模型检验用 Logistic模型来描述种群增长的规律效果如何呢？ 1945
年克朗皮克（ Crombic）做了一个人工饲养小谷虫的实验，数学生物学家高斯（ E·F·Gauss）也做了一个原生物草履虫实验，
实验结果都和 Logistic曲线十分吻合。
大量实验资料表明用 Logistic模型来描述种群的增长，效果还是相当不错的。例如，高斯 把 5只草履虫放进一个盛有
0.5cm3营养液的小试管，他发现，开始时草履虫以每天 230.9%
的速率增长，此后增长速度不断减慢，到第五天达到最大量
375个，实验数据与 r=2.309,a=0.006157,N(0)=5的 Logistic
曲线：
几乎完全吻合，见图 3.6。
2,3 0 9
375()
1 7 4 tNt e
图 3-6
Malthus模型和 Logistic模型的总结
Malthus模型和 Logistic模型 均为对微分方程（ 3.7）
所作的模拟近似方程。前一模型假设了种群增长率 r为一常数，（ r被称为该种群的内禀增长率）。后一模型则假设环境只能供养一定数量的种群，从而引入了一个竞争项。
用模拟近似法建立微分方程来研究实际问题时必须对求得的解进行检验，看其是否与实际情况相符或基本相符。
相符性越好则模拟得越好，否则就得找出不相符的主要原因，对模型进行修改。
Malthus模型与 Logistic模型虽然都是为了研究种群数量的增长情况而建立的，但它们也可用来研究其他实际问题，只要这些实际问题的数学模型有相同的微分方程即可，下面我们来看两个较为有趣的实例。
历史背景,
例 5 赝品的鉴定在第二次世界大战比利时解放以后，荷兰野战军保安机关开始搜捕纳粹同谋犯。他们从一家曾向纳粹德国出卖过艺术品的公司中发现线索，于 1945
年 5月 29日以通敌罪逮捕了三流画家范 ·梅格伦（ H·A·Vanmeegren），此人曾将 17世纪荷兰名画家扬 ·弗米尔（ Jan Veermeer）的油画,捉奸,等卖给纳粹德国戈林的中间人。可是，范 ·梅格伦在同年 7月 12日在牢里宣称：他从未把,捉奸,卖给戈林，而且他还说，这一幅画和众所周知的油画,在埃牟斯的门徒,以及其他四幅冒充弗米尔的油画和两幅德胡斯（ 17世纪荷兰画家）
的油画，都是他自己的作品，这件事在当时震惊了全世界，为了证明自己是一个伪造者，他在监狱里开始伪造弗米尔的油画,耶稣在门徒们中间,,当这项工作接近完成时，范 ·梅格伦获悉自己的通敌罪已被改为伪造罪，因此他拒绝将这幅画变陈，以免留下罪证。
为了审理这一案件，法庭组织了一个由著名化学家、物理学家和艺术史学家组成的国际专门小组查究这一事件。他们用 X射线检验画布上是否曾经有过别的画。此外，他们分析了油彩中的拌料（色粉），检验油画中有没有历经岁月的迹象。科学家们终于在其中的几幅画中发现了现代颜料钴兰的痕迹，还在几幅画中检验出了 20世纪初才发明的酚醛类人工树脂。根据这些证据，范 ·梅格伦于 1947年 10月 12日被宣告犯有伪造罪，被判刑一年。可是他在监狱中只待了两个多月就因心脏病发作，于 1947年 12月 30日死去。
然而，事情到此并未结束，许多人还是不肯相信著名的,在埃牟斯的门徒,是范 ·梅格伦伪造的。事实上，在此之前这幅画已经被文物鉴定家认定为真迹，并以 17万美元的高价被伦布兰特学会买下。专家小组对于怀疑者的回答是：由于范 ·梅格伦曾因他在艺术界中没有地位而十分懊恼，他下决心绘制
,在埃牟斯的门徒,,来证明他高于三流画家。当创造出这样的杰作后，他的志气消退了。而且，当他看到这幅,在埃牟斯的门徒,多么容易卖掉以后，
他在炮制后来的伪制品时就不太用心了 。这种解释不能使怀疑者感到满意，
他们要求完全科学地、确定地证明,在埃牟斯的门徒,的确是一个伪造品。
这一问题一直拖了 20年，直到 1967年，才被卡内基 ·梅伦（ Carnegie-Mellon）
大学的科学家们 基本上解决。
原理与模型测定油画和其他岩石类材料的年龄的关键是本世纪初发现的放射性现象。
放射性现象,著名物理学家卢瑟夫在本世纪初发现，某些,放射性,元素的原子是不稳定的，并且在已知的一段时间内，有一定比例的原子自然蜕变而形成新元素的原子，且物质的放射性与所存在的物质的原子数成正比。
用 N(t)表示时间 t时存在的原子数,则：
dN N
dt
常数 λ是正的，称为该物质的衰变常数用 λ来计算半衰期 T:
00()
dN N
dt
N t N
与负增长的 Malthus模型完全一样其解为,0()0() ttN t N e
0
1
2
N
N?
令
0
ln 2T t t
则有,
许多物质的半衰期已被测定，如碳 14，其 T=5568；
轴 238，其 T=45亿年。
与本问题相关的其他知识,
(1)艺术家们应用白铅作为颜料之一，已达两千年以上。白铅中含有微量的放射铅 210，白铅是从铅矿中提炼出来的，而铅又属于铀系，其演变简图如下（删去了许多中间环节）
(2)地壳里几乎所有的岩石中均含有微量的铀。一方面，铀系中的各种放射性物质均在不断衰减，而另一方面，铀又不断地衰减，补充着其后继元素。从而，各种放射性物质（除铀以外）在岩石中处于放射性平衡中。根据世界各地抽样测量的资料，地壳中的铀在铀系中所占平均重量比约为百万分之 2.7（一般含量极微）。各地采集的岩石中铀的含量差异很大，但从未发现含量高于 2— 3%的。
(3)从铅矿中提炼铅时，铅 210与铅 206一起被作为铅留下，
而其余物质则有 90— 95%被留在矿渣里，因而打破了原有的放射性平衡。（注：这些有关物理、地质方面的知识在建模时可向相应的专家请教。）
简化假定：
本问题建模是为了鉴定几幅不超过 300年的古画，为了使模型尽可能简单，可作如下假设：
(1)由于镭的半衰期为 1600年，经过 300年左右，应用微分方程方法不难计算出白铅中的镭至少还有原量的 90%，故可以假定，每克白铅中的镭在每分钟里的分解数是一个常数。
(2)铅 210的衰变为：
铅 210
T=22年钋 210 铅 206
T=138天若画为真品，颜料应有 300年左右或 300年以上的历史，容易证明：每克白铅中钋 210的分解数等于铅 210的分解数（相差极微，
已无法区别）。可用前者代替后者，因钋的半衰期较短，易于测量 。
建模：
(1)记提炼白铅的时刻为 t=0，当时每克白铅中铅 210的分子数为 y0，由于提炼前岩石中的铀系是处于放射性平衡的，故铀与铅的单位时间分解数相同。由此容易推算出每克白铅中铅
210每分钟分解数不能大于 30000个，否则铀的含量将超过 4%，
而这是不可能的。因为：
00 30000uUy若
200 3 0 0 0 0 6 0 2 4 3 6 5 1,0 2 1 0
u
U则 （个）
这些铀约重 20
23
1,02 10 23 8 0,04
6,02 10
（克）
即每克白铅约含 0.04克铀，含量为 4%
以上确定了每克白铅中铅分解数的上界，若画上的铅分解数大于该值，说明画是赝品；但若是小于不能断定画一定是真品。
(2)设 t时刻 1克白铅中铅 210含量为 y(t)，而镭的单位时间分解数为 r（常数），则 y(t)满足微分方程：
dy yr
dt
由此解得：
00( ) ( )0( ) [1 ]t t t try t e y e
00( ) ( )0 ( ) [ 1 ]t t t ty y t e r e故：
若此画是真品,t-t0≈300（年）。画中每克白铅所含铅 210
目前的分解数 λy(t)及目前镭的分解数 r均可用仪器测出，从而可求出 λy0的近似值，并利用（ 1）判断这样的分解数是否合理。
若判断结果为不合理，则可以确定此画必是赝品，但反之不一定说明画是真品（因为估计仍是十分保守的且只能证明画的
“年龄”）。
Carnegie-Mellon大学的科学家们利用上述模型对部分有疑问的油画作了鉴定，测得数据如下（见表 3-1）。
油画名称 210分解数 （ 个 /分 ） 镭 226分解数 （ 个 /分 ）
1,在埃牟斯的门徒 8.5 0.8
2,濯足 12.6 0.26
3,看乐谱的女人 10.3 0.3
4,演奏曼陀琳的女人
8.2 0.17
5,花边织工 1.5 1.4
6,笑女 5.2 6.0
计算 λy0 （ 个 /分 ）
98050
157130
127340
102250
1274.8
-10181
表 3-1
对,在埃牟斯的门徒,,λ y0≈98050 （个 /每克每分钟），它必定是一幅伪造品。类似可以判定（ 2），（ 3），（ 4）也是赝品。而（ 5）和（ 6）
都不会是几十年内伪制品，因为放射性物质已处于接近平衡的状态，这样的平衡不可能发生在十九世纪和二十世纪的任何作品中。
判定 结果：
利用放射原理，还可以对其他文物的年代进行测定。
例如对有机物（动、植物）遗体，考古学上目前流行的测定方法是放射性碳 14测定法，这种方法具有较高的精确度，
其基本原理是：由于大气层受到宇宙线的连续照射，空气中含有微量的中微子，它们和空气中的氮结合，形成放射性碳 14（ C14）。有机物存活时，它们通过新陈代谢与外界进行物质交换，使体内的 C14处于放射性平衡中。一旦有机物死亡，新陈代谢终止，放射性平衡即被破坏。因而，通过对比测定，可以估计出它们生存的年代。例如,1950年在巴比伦发现一根刻有 Hammurabi王朝字样的木炭，经测定，其 C14衰减数为 4.09个 /每克每分钟，而新砍伐烧成的木炭中 C14衰减数为 6.68个 /每克每分钟,C14的半衰期为
5568年，由此可以推算出该王朝约存在于 3900-4000年前。
例 6 新产品的推广经济学家和社会学家一直很关心新产品的推销速度问题。怎样建立一个数学模型来描述它，并由此析出一些有用的结果以指导生产呢？以下是第二次世界大战后日本家电业界建立的电饭包销售模型。
设需求量有一个上界，并记此上界为 K，记 t时刻已销售出的电饭包数量为 x(t)，则尚未使用的人数大致为 K－ x(t)，于是由统计筹算律：
()dx x K xdt
记比例系数为 k,则 x(t)满足：
()dx kx K xdt
此方程即 Logistic模型，解为：
() 1 K ktKxt Ce
还有两个奇解,
x=0和 x=K
对 x(t)求一阶、两阶导数：
2
2'( ) (1 )
K k t
K k t
c K k ext
Ce
32
3
( 1 )()
(1 )
K k t K k t
K k t
C K k e C ext
Ce
容易看出,x’(t)>0，即 x(t)单调增加。
由 x’’(t0)=0，可以得出 =1，此时,。
0RKtCe? 2)( 0 Ktx?
当 t<t0时,x’’(t)>0,x’(t)单调增加，而当 t>t0时,x’’(t)<0，
x’(t)单调减小。 实际调查表明，销售曲线与 Logistic曲线十分接近，尤其是在销售后期，两者几乎完全吻合。 在销出量小于最大需求量的一半时，销售速度是不断增大的，
销出量达到最大需求量的一半时，该产品最为畅销，接着销售速度将开始下降。
所以初期应采取小批量生产并加以广告宣传；从有 20%用户到有
80%用户这段时期，应该大批量生产；后期则应适时转产，这样做可以取得较高的经济效果。
本节将建立几个简单的单种群增长模型，以简略分析一下这方面的问题。一般生态系统的分析可以通过一些简单模型的复合来研究，大家若有兴趣可以根据生态系统的特征自行建立相应的模型。 美丽的大自然种群的数量本应取离散值，但由于种群数量一般较大，为建立微分方程模型，可将种群数量看作连续变量，甚至允许它为可微变量，
由此引起的误差将是十分微小的。
离散化为连续，方便研究
§ 3.2 Malthus模型与 Logistic模型模型 1 马尔萨斯（ Malthus）模型马尔萨斯在分析人口出生与死亡情况的资料后发现，
人口净增长率 r基本上是一常数，（ r=b-d,b为出生率，
d为死亡率）,既：
1 dN r
N dt?
dN rN
dt?
或 （ 3.5）
0()0() r t tN t N e （ 3.6）
（ 3.1） 的解为：
其中 N0=N(t0)为初始时刻 t0时的种群数。
马尔萨斯模型的一个显著特点,种群数量翻一番所需的时间是固定的 。
令种群数量翻一番所需的时间为 T，则有：
002 rTN N e?
ln 2T
r?
故模型检验比较历年的人口统计资料，可发现人口增长的实际情况与马尔萨斯模型的预报结果基本相符，例如,1961年世界人口数为 30.6 （即 3.06× 109），人口增长率约为 2%，人口数大约每 35年增加一倍。检查 1700年至 1961的 260年人口实际数量，发现两者几乎完全一致，且按马氏模型计算，人口数量每 34.6年增加一倍，两者也几乎相同。
1950 2000 2050 2100 2150 2200
0
0,5
1
1,5
2
2,5
3
3,5
x 1 0
11
t/ 年
N/
人马尔萨斯模型人口预测模型预测假如人口数真能保持每 34.6年增加一倍，那么人口数将以几何级数的方式增长。例如，到 2510年，人口达 2× 1014个，
即使海洋全部变成陆地，每人也只有 9.3平方英尺的活动范围，
而到 2670年，人口达 36× 1015个，只好一个人站在另一人的肩上排成二层了。 故 马尔萨斯模型是不完善的。
几何级数的增长
Malthus模型 实际上只有在群体总数不太大时才合理，到总数增大时，
生物群体的各成员之间由于有限的生存空间，有限的自然资源及食物等原因，就可能发生生存竞争等现象。
所以 Malthus模型假设的人口 净增长率不可能始终保持常数，
它应当与人口数量有关。
模型 2 Logistic模型人口净增长率应当与人口数量有关，即,r=r(N)
从而有：
()dN r N Ndt?
（ 3.7）
r(N)是未知函数，但根据实际背景，它无法用拟合方法来求 。
为了得出一个有实际意义的模型，我们不妨采用一下工程师原则。工程师们在建立实际问题的数学模型时，总是采用尽可能简单的方法。
r(N)最简单的形式是常数，此时得到的就是马尔萨斯模型。
对马尔萨斯模型的最简单的改进就是引进一次项（竞争项）
对马尔萨斯模型引入一次项（竞争项），令 r(N)=r-aN
此时得到微分方程：
()dN r a N Ndt (1 )d N NrNd t K或 （ 3.8）
（ 3.8） 被称为 Logistic模型或生物总数增长的统计筹算律，是由荷兰数学生物学家弗赫斯特（ Verhulst）首先提出的。一次项系数是负的，因为当种群数量很大时，会对自身增大产生抑制性，故一次项又被称为竞争项。
（ 3.8） 可改写成：
()dN k K N Ndt （ 3.9）
(3.9)式还有另一解释，由于空间和资源都是有限的，不可能供养无限增长的种群个体，当种群数量过多时，由于人均资源占有率的下降及环境恶化、疾病增多等原因，出生率将降低而死亡率却会提高。设环境能供养的种群数量的上界为 K（近似地将 K看成常数）,N表示当前的种群数量，
K-N恰为环境还能供养的种群数量，（ 3.9）指出，种群增长率与两者的乘积成正比，正好符合统计规律，得到了实验结果的支持，这就是（ 3.9）
也被称为统计筹算律的原因。
图 3-5
对 （ 3.9） 分离变量：
11 d N k K d t
N K N
两边积分并整理得：
1 kK t
KN
Ce
令 N(0)=N0，求得：
0
0
KNC
N
故 （ 3.9） 的满足初始条件 N(0)=N0的解为：
0
00
() () k KtNKNt N K N e （ 3.10）
易见：
N(0)=N0,lim ( )
t N t K
N(t)的图形请看图 3.5
模型检验用 Logistic模型来描述种群增长的规律效果如何呢？ 1945
年克朗皮克（ Crombic）做了一个人工饲养小谷虫的实验，数学生物学家高斯（ E·F·Gauss）也做了一个原生物草履虫实验，
实验结果都和 Logistic曲线十分吻合。
大量实验资料表明用 Logistic模型来描述种群的增长，效果还是相当不错的。例如，高斯 把 5只草履虫放进一个盛有
0.5cm3营养液的小试管，他发现，开始时草履虫以每天 230.9%
的速率增长，此后增长速度不断减慢，到第五天达到最大量
375个，实验数据与 r=2.309,a=0.006157,N(0)=5的 Logistic
曲线：
几乎完全吻合，见图 3.6。
2,3 0 9
375()
1 7 4 tNt e
图 3-6
Malthus模型和 Logistic模型的总结
Malthus模型和 Logistic模型 均为对微分方程（ 3.7）
所作的模拟近似方程。前一模型假设了种群增长率 r为一常数，（ r被称为该种群的内禀增长率）。后一模型则假设环境只能供养一定数量的种群，从而引入了一个竞争项。
用模拟近似法建立微分方程来研究实际问题时必须对求得的解进行检验，看其是否与实际情况相符或基本相符。
相符性越好则模拟得越好，否则就得找出不相符的主要原因，对模型进行修改。
Malthus模型与 Logistic模型虽然都是为了研究种群数量的增长情况而建立的，但它们也可用来研究其他实际问题，只要这些实际问题的数学模型有相同的微分方程即可，下面我们来看两个较为有趣的实例。
历史背景,
例 5 赝品的鉴定在第二次世界大战比利时解放以后，荷兰野战军保安机关开始搜捕纳粹同谋犯。他们从一家曾向纳粹德国出卖过艺术品的公司中发现线索，于 1945
年 5月 29日以通敌罪逮捕了三流画家范 ·梅格伦（ H·A·Vanmeegren），此人曾将 17世纪荷兰名画家扬 ·弗米尔（ Jan Veermeer）的油画,捉奸,等卖给纳粹德国戈林的中间人。可是，范 ·梅格伦在同年 7月 12日在牢里宣称：他从未把,捉奸,卖给戈林，而且他还说，这一幅画和众所周知的油画,在埃牟斯的门徒,以及其他四幅冒充弗米尔的油画和两幅德胡斯（ 17世纪荷兰画家）
的油画，都是他自己的作品，这件事在当时震惊了全世界，为了证明自己是一个伪造者，他在监狱里开始伪造弗米尔的油画,耶稣在门徒们中间,,当这项工作接近完成时，范 ·梅格伦获悉自己的通敌罪已被改为伪造罪，因此他拒绝将这幅画变陈，以免留下罪证。
为了审理这一案件，法庭组织了一个由著名化学家、物理学家和艺术史学家组成的国际专门小组查究这一事件。他们用 X射线检验画布上是否曾经有过别的画。此外，他们分析了油彩中的拌料（色粉），检验油画中有没有历经岁月的迹象。科学家们终于在其中的几幅画中发现了现代颜料钴兰的痕迹，还在几幅画中检验出了 20世纪初才发明的酚醛类人工树脂。根据这些证据，范 ·梅格伦于 1947年 10月 12日被宣告犯有伪造罪，被判刑一年。可是他在监狱中只待了两个多月就因心脏病发作，于 1947年 12月 30日死去。
然而，事情到此并未结束，许多人还是不肯相信著名的,在埃牟斯的门徒,是范 ·梅格伦伪造的。事实上，在此之前这幅画已经被文物鉴定家认定为真迹，并以 17万美元的高价被伦布兰特学会买下。专家小组对于怀疑者的回答是：由于范 ·梅格伦曾因他在艺术界中没有地位而十分懊恼，他下决心绘制
,在埃牟斯的门徒,,来证明他高于三流画家。当创造出这样的杰作后，他的志气消退了。而且，当他看到这幅,在埃牟斯的门徒,多么容易卖掉以后，
他在炮制后来的伪制品时就不太用心了 。这种解释不能使怀疑者感到满意，
他们要求完全科学地、确定地证明,在埃牟斯的门徒,的确是一个伪造品。
这一问题一直拖了 20年，直到 1967年，才被卡内基 ·梅伦（ Carnegie-Mellon）
大学的科学家们 基本上解决。
原理与模型测定油画和其他岩石类材料的年龄的关键是本世纪初发现的放射性现象。
放射性现象,著名物理学家卢瑟夫在本世纪初发现，某些,放射性,元素的原子是不稳定的，并且在已知的一段时间内，有一定比例的原子自然蜕变而形成新元素的原子，且物质的放射性与所存在的物质的原子数成正比。
用 N(t)表示时间 t时存在的原子数,则：
dN N
dt
常数 λ是正的，称为该物质的衰变常数用 λ来计算半衰期 T:
00()
dN N
dt
N t N
与负增长的 Malthus模型完全一样其解为,0()0() ttN t N e
0
1
2
N
N?
令
0
ln 2T t t
则有,
许多物质的半衰期已被测定，如碳 14，其 T=5568；
轴 238，其 T=45亿年。
与本问题相关的其他知识,
(1)艺术家们应用白铅作为颜料之一，已达两千年以上。白铅中含有微量的放射铅 210，白铅是从铅矿中提炼出来的，而铅又属于铀系，其演变简图如下（删去了许多中间环节）
(2)地壳里几乎所有的岩石中均含有微量的铀。一方面，铀系中的各种放射性物质均在不断衰减，而另一方面，铀又不断地衰减，补充着其后继元素。从而，各种放射性物质（除铀以外）在岩石中处于放射性平衡中。根据世界各地抽样测量的资料，地壳中的铀在铀系中所占平均重量比约为百万分之 2.7（一般含量极微）。各地采集的岩石中铀的含量差异很大，但从未发现含量高于 2— 3%的。
(3)从铅矿中提炼铅时，铅 210与铅 206一起被作为铅留下，
而其余物质则有 90— 95%被留在矿渣里，因而打破了原有的放射性平衡。（注：这些有关物理、地质方面的知识在建模时可向相应的专家请教。）
简化假定：
本问题建模是为了鉴定几幅不超过 300年的古画，为了使模型尽可能简单，可作如下假设：
(1)由于镭的半衰期为 1600年，经过 300年左右，应用微分方程方法不难计算出白铅中的镭至少还有原量的 90%，故可以假定，每克白铅中的镭在每分钟里的分解数是一个常数。
(2)铅 210的衰变为：
铅 210
T=22年钋 210 铅 206
T=138天若画为真品，颜料应有 300年左右或 300年以上的历史，容易证明：每克白铅中钋 210的分解数等于铅 210的分解数（相差极微，
已无法区别）。可用前者代替后者，因钋的半衰期较短，易于测量 。
建模：
(1)记提炼白铅的时刻为 t=0，当时每克白铅中铅 210的分子数为 y0，由于提炼前岩石中的铀系是处于放射性平衡的，故铀与铅的单位时间分解数相同。由此容易推算出每克白铅中铅
210每分钟分解数不能大于 30000个，否则铀的含量将超过 4%，
而这是不可能的。因为：
00 30000uUy若
200 3 0 0 0 0 6 0 2 4 3 6 5 1,0 2 1 0
u
U则 （个）
这些铀约重 20
23
1,02 10 23 8 0,04
6,02 10
（克）
即每克白铅约含 0.04克铀，含量为 4%
以上确定了每克白铅中铅分解数的上界，若画上的铅分解数大于该值，说明画是赝品；但若是小于不能断定画一定是真品。
(2)设 t时刻 1克白铅中铅 210含量为 y(t)，而镭的单位时间分解数为 r（常数），则 y(t)满足微分方程：
dy yr
dt
由此解得：
00( ) ( )0( ) [1 ]t t t try t e y e
00( ) ( )0 ( ) [ 1 ]t t t ty y t e r e故：
若此画是真品,t-t0≈300（年）。画中每克白铅所含铅 210
目前的分解数 λy(t)及目前镭的分解数 r均可用仪器测出，从而可求出 λy0的近似值，并利用（ 1）判断这样的分解数是否合理。
若判断结果为不合理，则可以确定此画必是赝品，但反之不一定说明画是真品（因为估计仍是十分保守的且只能证明画的
“年龄”）。
Carnegie-Mellon大学的科学家们利用上述模型对部分有疑问的油画作了鉴定，测得数据如下（见表 3-1）。
油画名称 210分解数 （ 个 /分 ） 镭 226分解数 （ 个 /分 ）
1,在埃牟斯的门徒 8.5 0.8
2,濯足 12.6 0.26
3,看乐谱的女人 10.3 0.3
4,演奏曼陀琳的女人
8.2 0.17
5,花边织工 1.5 1.4
6,笑女 5.2 6.0
计算 λy0 （ 个 /分 ）
98050
157130
127340
102250
1274.8
-10181
表 3-1
对,在埃牟斯的门徒,,λ y0≈98050 （个 /每克每分钟），它必定是一幅伪造品。类似可以判定（ 2），（ 3），（ 4）也是赝品。而（ 5）和（ 6）
都不会是几十年内伪制品，因为放射性物质已处于接近平衡的状态，这样的平衡不可能发生在十九世纪和二十世纪的任何作品中。
判定 结果：
利用放射原理，还可以对其他文物的年代进行测定。
例如对有机物（动、植物）遗体，考古学上目前流行的测定方法是放射性碳 14测定法，这种方法具有较高的精确度，
其基本原理是：由于大气层受到宇宙线的连续照射，空气中含有微量的中微子，它们和空气中的氮结合，形成放射性碳 14（ C14）。有机物存活时，它们通过新陈代谢与外界进行物质交换，使体内的 C14处于放射性平衡中。一旦有机物死亡，新陈代谢终止，放射性平衡即被破坏。因而，通过对比测定，可以估计出它们生存的年代。例如,1950年在巴比伦发现一根刻有 Hammurabi王朝字样的木炭，经测定，其 C14衰减数为 4.09个 /每克每分钟，而新砍伐烧成的木炭中 C14衰减数为 6.68个 /每克每分钟,C14的半衰期为
5568年，由此可以推算出该王朝约存在于 3900-4000年前。
例 6 新产品的推广经济学家和社会学家一直很关心新产品的推销速度问题。怎样建立一个数学模型来描述它，并由此析出一些有用的结果以指导生产呢？以下是第二次世界大战后日本家电业界建立的电饭包销售模型。
设需求量有一个上界，并记此上界为 K，记 t时刻已销售出的电饭包数量为 x(t)，则尚未使用的人数大致为 K－ x(t)，于是由统计筹算律：
()dx x K xdt
记比例系数为 k,则 x(t)满足：
()dx kx K xdt
此方程即 Logistic模型，解为：
() 1 K ktKxt Ce
还有两个奇解,
x=0和 x=K
对 x(t)求一阶、两阶导数：
2
2'( ) (1 )
K k t
K k t
c K k ext
Ce
32
3
( 1 )()
(1 )
K k t K k t
K k t
C K k e C ext
Ce
容易看出,x’(t)>0，即 x(t)单调增加。
由 x’’(t0)=0，可以得出 =1，此时,。
0RKtCe? 2)( 0 Ktx?
当 t<t0时,x’’(t)>0,x’(t)单调增加，而当 t>t0时,x’’(t)<0，
x’(t)单调减小。 实际调查表明，销售曲线与 Logistic曲线十分接近，尤其是在销售后期，两者几乎完全吻合。 在销出量小于最大需求量的一半时，销售速度是不断增大的，
销出量达到最大需求量的一半时，该产品最为畅销，接着销售速度将开始下降。
所以初期应采取小批量生产并加以广告宣传；从有 20%用户到有
80%用户这段时期，应该大批量生产；后期则应适时转产，这样做可以取得较高的经济效果。
