人口结构模型.docx

谈研究年龄结构之数学模型 Leslies Model 许世壁 对外搜寻关键词.人口动力学. 对外搜寻关键词 .人口动力学 .固有值 .固有向量 .Primary decomposition Theorem .特征多项式 .Frobenius .basis 关年龄结构(Age Strucure)之问题。举个例来说明：目前台湾人口中0?5岁有xi 人，6?10岁有X2人，……,76?80岁有 刘6人。试问20年或100年后，这些人 口的变化如何？还有每一个年龄分类 (Age class)在总人口上的比例会不会很稳定 地趋向某一个固定值？如果是的话，多快？就社会学、经济学而言，这是一个很实 际而且值得研究的重要问题。下面我们就要导出有关这个问题之数学模型 Leslies Model 。它同时也可以应用到其他生物，如鱼类及昆虫等。 首先，假设从现有的统计数据，我们能选择出一个适当的单位时间 T,而后将人口 Vat 分成A类。令向量 代表在第N期时(时间为NT)人口里的女性年龄结构(在 此我们假设男性，女性人口数目相等)，简而言之令 在此，分量Vkn，k= 1,…，A代表年龄介于(k-1)T及kT中间之女性总人数。譬如 说应用到实际人口时，我们通常取 T=5年，而且将人口分成16类，即A=16， V1,n =在第N期时，介于0?5岁之女性总人数。 V2,n =在第N期时，介于6?10岁之女性总人数。 V16,n =在第N期时，介于76?80岁之女性总人数。 如果年龄超过80岁时，则我们不予讨论。 下一步我们要做的工作是如何找出第 N+1期的年龄结构向量(Age structure vector) ，与在第N期之年龄结构向量 之关系。假设下列的bk及 mk, k=1,…，A为已知， 加=在軍位琦間個年常分IB之女雀平均每人所生出的文夜敷目叫二第k個年爺分類能活過軍位時間T而變成第k + 1年齡分類之百分比 利用bn, mk及Vk,N之定义，我们可以导出下列关系式： 口｝%“十 1 = fc = 2g 3j 1 9A 门)=虹旳小+…+ bmVA,N 其中(2)式说明了从第N期到N+1期所出生之女孩总数。所以，如果我们假设 虹昇「)亦 TJ11，…，TTLA 及 这些非负之实数均可由人口之统计资料得到， 则我们有下 列式子 VS3w+i = + *1 - + ^aVa^n 幻小十h = (3); Va,n^i = 或者用矩阵之表示法，(3)式可改写为(4)(匚西Be旬 或者用矩阵之表示法，(3)式可改写为 (4) (匚西Be旬Mo dal) 771] 、、 亦_1 … 0 (5) M = ntz [Rffpioduction Mitrix) 0 m4-l 如果，我们假设现在的年龄结构向量。为已知，则由⑷式，我们得到: 如果，我们假设现在的年龄结构向量 。为已知，则由⑷式，我们得到: 1^n=MnV0 所以，我们将年龄结构的问题变成一个线性代数的问题: 为了解决这个问题，必须利用线性代数中有关固有值 (eige nvalue) 及固有向量 (eigenvector) 之观念及其重要定理 Primary decomposition Theorem (参考[1])。 首先我们考虑 A x A矩阵M之固有值 入。从固有值之定义，入为M之特征多 f(A)=detf Af- XZ) = 0 项式(characteristic poly no mial) 之根。通常特征多项 式很难算，但在这里矩阵 M有其特殊形式(5)，所以利用降级展开行列式det M _ XI ( )，得到 f(A)=屮— f(A)=屮—虹沪一1 — — X0 fw o /(A)=o 因为f(0) 0而且当 够大时 ，所以 必有一正实根。事实 上,因为矩阵M里之系数皆大于或等于零，根据有名的Frobenius定理(参考[2]) An 0 |A| Ao TOC \o 1-5 \h \z 我们得知存在一固有值— ，而且其他A-1个固有值 入，满足 。现 在，我们假设矩阵 M满足下列性质： Ao 0 | A| An (H)存在一固有值 而且其他A-1个固有值 人满足 。 、 一 Ao 在假设(H)下，我们可以用数值分析之方法 Power Method实际地算出 (参考 Ao Ag 妙口 [3])。有了 ，因为 满足⑺ 式，我们可以检查一下下列 为一对应于 入 之固有向量； 矶 Eo c R? 令Eo为由向量 所产生之一维子空间， 。由线性代数的Primary decomposition Theorem，我们可将 RA写成RA ，其中Ei是以对于固 ，… 人*入0、一，，— ，、 一，，「 有值 之固有向量(eigenvectors) 或广义固有向量(generali