人口年龄结构模型建模和预测

上海交通大学硕士学位论文人口年龄结构模型建模和预测姓名:虞丽萍申请学位级别：硕士专业：控制理论与控制工程指导教师：袁景淇200７０１01人口年龄结构模型建模和预测摘ﻮ要人口是一个动态系统。人口变化对未来经济、社会的发展有着直接的影响。人口年龄结构是人口研究的重要指标之一，人口年龄结构发展趋势的预报对人口政策的制定有着非常重要的作用.本文以离散形式的人口发展方程为主模型。在此基础上,分别建立了生育率、死亡率和迁移模型，以预测人口生育率、死亡率、流动人口和人口年龄结构的变化趋势。与传统模型相比，本文所提出的基于随机分布函数的生育率组合模型和死亡率分段模型使模型精度得到了进一步的提高。本文以中国历年统计数据为原始数据,验证了主模型和各子模型的有效性，并预测了2015年上海市人口年龄结构，绘制了人口年龄树。人口年龄树树形反映了人口结构的健康状态。通过分析人口年龄树的变化趋势,可以了解人口结构所存在的问题，为政府调控人口提供科学依据.关键词:离散人口发展方程，生育率,死亡率，流动人口,人口预测，人口年龄树MＯDEＬINGANDFＯＲECＡSTINＧＴＨEAＧEＳＴRUＣＴURＥOＦPOPＵLATＩＯＮABSTRAＣTThepｏpulaｔiｏｎｓyｓｔｅmisadｙnamｉｃａｌｓｙｓteｍ。Tｈｅtｒenｄｏfapoｐｕｌatｉｏnwｉｌlａｆfectthedeveｌopmentｏfthesｏcｉeｔyandｉtseｃｏｎｏmｙ.Thｅaｇestrｕｃtｕreisｏｎeｏftheｍosｔｉｍｐortantｉndｅxesinpopulａtｉoｎｒeseａrcｈ。Thｅｆｏrecａｓｔoｆageｓｔrucｔurepｌaysａｎiｍportantroleiｎmakｉｎgpopulatioｎpoｌicies．Adisｃrｅtｅmoｄeｌofpopulatｉondｅvｅlopmｅntｗasａppliｅdｉnｔhispapｅr。Fｅrｔilityｍodeｌ，mortaｌiｔyｍoｄeｌａndimｍiｇrａtiｏｎmodelwerealsoesｔabｌisｈｅdtoｆoreｃasttｈeferｔiｌiｔｙ,ｍortａｌity，immｉgrationpopｕｌatｉoｎaｎdｔheagestｒｕcture．Compａreｄwｉthtrａｄitiｏｎａlｆｕｎctionｓ,boththeｃoｍｐositefertilitymodｅlwｈｉchｉsbaseｄonranｄoｍdistｒibutioｎfunctｉonｓaｎｄtｈeseｇmｅｎtｅｄmoｒｔａｌｉｔｙｍodelimｐrovｅdｔheｍｏdｅls'ａｃcuracｙ。HiｓtｏricａlChinesepｏpuｌatｉoｎｓtatistiｃsｗeｒeuｓeｄvｅｔhevalidｉtyｏftｈemoｄeｌsｒｅｆerreｄｉnthispａpeｒ。ＴｈeａｇeｓtｒｕcturｅoｆＳhanｇhａiiｎ2015wａsｆoreｃａｓtｅdandｄisplaｙｅdinａｔree—lｉkｅgｒaｐh．Thｅｓｈａｐeoftｈepｏｐulatiｏｎtｒeｅreｆlｅｃtstheheａｌｔｈcoｎdiｔioｎｏfｉtｓageｓtruｃｔuｒｅ。Bｙanａｌyｚｉnｇthｅdｅvelopmｅｎttrｅnｄofｔｈeagestrucｔurｅ，wecanfindpopulａtｉｏbｌeｍsaｎｄpｒoｖideｓcｉｅｎtificｅviｄencｅforｇovｅｒnmenttocoｎtrolｔhepｏｐulation.KEＹ

WORDS：ﻮＤiｓcｒｅteﻮPopulaｔionﻮDｅvｅloｐmｅntﻮＥquatiｏns,ﻮFertｉlitｙ，Moｒtaliｔｙ，Ｉmmｉｇrａtiｏn，ＰopulaｔiｏnForeｃaｓt,ＡgeStructurePopｕｌationＴree图片目录图１人口金字塔·······································································································5图２人口年龄树·······································································································６图3人口结构的三种类型······················································································8图４人口状态方程控制框图················································································１１图5人口发展方程的数据流图··············································································１２图6三层BP神经网络的拓扑结构·······································································1６图7世代生育率（CFR）与总和生育率(ＴFR）的比较············································1８图８２０00年分孩次的年龄别生育率模型拟合结果与统计数据比较···················2４图9组合模型、对数正态分布模型及泊松分布模型精度比较···························27图10不同模型分年龄别生育率误差比较····························································28图11不同模型对20０4年分年龄别生育率的预测估计········································３２图12组合模型对201５年全国分年龄别生育率的预测·······································3２图１３２003年中国人口死亡率三次样条插值结果与实际数据比较（男)·········３5图１４2003年中国人口死亡率三次样条插值结果与实际数据比较（女）·········３6图152０0１年中国人口死亡率分段模型拟合结果与实际数据比较··················41图16影响人口死亡率的因素················································································4２图17参数a估计值································································································４4图1８参数b估计值·······························································································４4图19参数c估计值································································································４5图20２00４年中国男性分年龄别死亡率预测值与实际值比较··························46图2１２０04年中国男性分年龄别死亡率预测值与实际值比较··························４７图２２200０年上海市外来人口按年龄别分布图·················································5５图23２0０4年中国人口年龄树·············································································５８图２42０04年中国人口年龄结构预测值与实际值相对误差······························58图25上海市历年总和生育率··············································································64图２６２０１5年上海市人口年龄树·········································································65表格目录表1Lｏｇnorｍａl和Ｐｏisｓon分布模型描述分孩次的年龄别生育率·······················２2表２组合模型参数的估计值(最小二乘法)······················································２5表3总和生育率GＭ（１，1)模型预测结果··························································３1表4２003年中国分年龄组死亡率统计数据··························································34表519８9年中国人口死亡率各模型拟合均方根误差比较···································３9表6１9９４年中国人口死亡率各模型拟合均方根误差比较···································39表71９98年中国人口死亡率各模型拟合均方根误差比较···································39表8２０01年中国人口死亡率各模型拟合均方根误差比较···································40表9中国历年人均国民生产总值ＧＮP·································································43表1０上海市第五次人口普查外来人口年龄分布················································５3表11上海市主要年份人口迁移数据····································································5５表1２上海市人口迁移数据预测············································································56表１32004年中国人口年龄结构············································································５9表１4人口年龄结构类型(国际通用标准)······················································59表15上海市总和生育率预测················································································64表16上海市２01５年人口年龄树相关数据··························································６5符号说明p（r，t）µϕNβhgkｌｋpkｃfqk，ｒ,t,r1，r２a，b,u，A,Bα，β,µ，σ，cﻮ人口年龄分布密度函数死亡力（‰)出生婴儿总数（人）人口总数（人)妇女平均生育率，即总和生育率（人）妇女生育模式（‰）迁移人口数（人)Logｎormal模型尺度变换因子Ｐoiｓson分布模型尺度变换因子生育率组合模型尺度变换因子生育率死亡概率模型参数上海交通大学学位论文版权使用授权书本学位论文作者完全了解学校有关保留、使用学位论文的规定,同意学校保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅.本人授权上海交通大学可以将本学位论文的全部或部分内容编入有关数据库进行检索,可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。保密□，在年解密后适用本授权书。本学位论文属于不保密√。(请在以上方框内打“√”）学位论文作者签名：虞丽萍指导教师签名:袁景淇日期：2007年2月26日日期:2０07年2月26日–８2–上海交通大学学位论文原创性声明本人郑重声明:所呈交的学位论文，是本人在导师的指导下，独立进行研究工作所取得的成果。除文中已经注明引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写过的作品成果。对本文的研究做出重要贡献的个人和集体,均已在文中以明确方式标明.本人完全意识到本声明的法律结果由本人承担。学位论文作者签名：虞丽萍日期：２0０７年2月２6日–８3–第一章引言1。1ﻮ研究背景人口数量、质量和年龄分布直接影响一个地区的经济发展、资源分配、社会保障、社会稳定和城市活力。对此,单纯的人口数量控制(如已实施多年的计划生育)不能体现人口规划的科学性。政府决策部门需要更详细、更系统的人口分析技术，为人口发展策略的制定提供指导和依据。长期以来，对人口年龄结构的研究仅限于粗线条的定性分析，只能预测年龄结构分布的大致范围,无法用于分析年龄结构的具体形态。随着对人口规划精准度要求的提高，通过数学方法来定量计算各种人口指数的方法日益受到重视，这就是人口控制与人口预测。人口控制和人口预测是人口研究中的重要课题。准确的人口预测为制定合理的社会经济发展规划提供了科学依据.例如,要制定生育计划，就必须知道未来妇女的生育率；要制定社会保障体系，就必须知道未来老年人口动态变化量;要规划学校建设，就必须知道学龄人口数；要改善医疗保障，就必须对未来人口的死亡状况有清楚的认识;要确定人才引进策略和户籍管理制度，就必须了解迁移人口年龄分布及受教育情况。这些都离不开人口预测。政府可以根据这些未来人口信息状况，结合社会经济发展，在制定国民经济发展决策时，通过调控人口的方法对未来社会经济发展中的产业结构进行相应的调整，使劳动力资源得到充分地开发和利用，社–1–会资源得到合理分配,并采取措施提前应对由老龄化，低生育率等人口因素而可能产生的社会经济问题，从而使社会经济协调发展.1。2

研究内容结合人口发展的特点,人们建立了不同的人口模型,并利用统计数据加以验证和预测。各级政府部门和研究机构也建立了人口信息系统,用以提供制定人口政策和发展规划的参考。考虑到人口发展问题的复杂性和区域的差异性很大，本文以区域人口年龄结构为研究对象,以中国和上海市为预测实例，以历年人口统计数据为验证基础,建立人口循环增殖预测模型。在该模型的基础上，进一步建立出生率、死亡率、人口迁移等子模型，对未来人口结构发展趋势进行预报。模型结果用人口年龄树树形图来展示。通过观察人口年龄树树形的变化,获得对人口年龄结构的变迁轨迹和人口年龄树健康状况的直观判断,为制定人口调控政策提供科学的参考。本文采用灰色系统模型等预测方法对生育率、死亡率、流动人口的变化趋势进行预测。预测所得的结果代入生育率、死亡率和迁移模型，可以计算出分年龄别的人口预测数据.上述预测数据作为人口发展方程的输入量,经过循环迭代，就能得到未来的分年龄性别的人口数.在此基础上,进一步计算可得各类人口指数.结合上述数据,本文对我国和上海市的当前人口状况和未来人口状况进行了详细的分析、比较,为政府人口政策的制定和实施提供宏观上的方向把握和微观上的数据的支持。–２–第二章人口年龄结构主模型2.1ﻮ人口概念辨析本文所涉及的主要概念和定义有:（1)人口:生活在一定社会生产方式、一定时间、一定地域，实现其生命活动并构成社会生活主体，具有一定数量和质量的人所组成的社会群体。(２）出生率:指某年每１000人对应的活产数，又称总出生率或粗出生率。它反映人口的出生水平，一般以千分数表示。(3）生育率：某年每1000名15—49岁妇女的活产婴儿数。又称一般生育率。该指标比出生率要精确一些，因为它将生育同可能生育的特定性别年龄的人口联系起来（通常是15-4９岁的妇女），排除了年龄性别结构不同引起的偏差。生育率比出生率更能揭示生育水平的变化。（４）总和生育率(ＴFＲ):指假设妇女按照某一年的年龄别生育率度过育龄期，平均每个妇女在育龄期生育的孩子数（5）死亡率：一定时期内（通常为一年）死亡人数与同期平均人数(或期中人数）之比。说明该时期人口的死亡强度,通常用千分比表示。(６）人口迁移：人口在地理上的位置变更。人口为了某种目的或动机,离开原来的居住地，时间或长或矩,距离或远或近,或者返回或者终生不再返回而定居于某地，–3–均称为人口迁移。（7)人口增长率：人口增长程度或增长速度,即一定时期内人口增长数与人口总数之比。通常以一年为期计算，用百分数表示。（８）人口性别比：指某一人口中男性对女性的比例,通常以每100个女性对应的男性数来表示。（9）儿童妇女比：指某年每１00０名妇女（1５—４9)对应5岁以下孩子的数目。这个指标可以从人口普查或抽样调查中得到，因此它可在出生统计不详细的情况下提供生育水平的数据。（10）生育更替水平：生育更替水平是指这样一个生育水平,即同一批妇女生育女儿的数量恰好能替代她们本身。当净人口再生产率为1．０0，总和生育率为2．1时，恰好等于更替水平。一旦达到生育更替水平，出生和死亡将逐渐趋于均衡,在没有国际迁入与迁出的情况下，人口将最终停止增长,保持稳定状态.（１1）人口年龄结构：某一年某一地区按年龄的人口总数。(12）人口金字塔:是形象地表示某一人口的年龄和性别构成的图形。水平条代表每一年龄组男性和女性的数字或比例（男左女右).金字塔中各个年龄性别组相加构成了总人口。人口金字塔可以用1岁年龄组的数据绘制,如图1。–4–图１人口金字塔（数据来源：１９9０年上海市人口年龄结构.男左女右）Fig.1Popuｌａtｉｏｎpｙramｉd(１3)人口年龄树：是本文所提出的一种用来表述人口年龄结构的图形,与人口金字塔类似。图形外沿用光滑连续曲线代替人口金字塔的水平条型图，因其形状类似于树,将此类图形命名为人口年龄树,如图2.–5–图2人口年龄树(数据来源:1９90年上海市人口年龄结构.男左女右）Ｆiｇ。2Aｇｅstrｕctuｒetｒeｅofpｏpulaｔｉon2.2ﻮ人口结构分类人口结构是反映人口年龄性别分布的重要指标，对人口变化和社会发展具有不容忽视的影响作用。人口结构是人口再生产的基础，任何时点上的人口结构都是历史人口生育、死亡和迁移的结果，又是研究未来人口过程的基础.人口结构对社会发展起着促进或制约的作用.本文以人口结构作为主要研究对象，既可以分析人口内部结构的特点与变动趋势，也可以反映人口发展的趋势和规律性，并进一步考察人口结构与社会经济因素的相互关系。国际上通常将人口结构分为三类(见图3)：–6–（1)增长型（年轻型）：图形上表现为底部宽，顶部狭窄，即少年儿童人口比高，老年人口比低，显示人口快速成长。此类型人口结构的特点是死亡率快速衰减,而出生率未改变,或仅缓慢降低的结果。（2）静止型(成年型)：图形上表现为各年龄组的比例较相似。这一类型人口结构的特点是低死亡率及接近更替水平的生育率。只有当死亡率水平为千分之十至十五，妇女生育率低于2的情况存在至少２0年,才会形成这类人口结构。大部分生活水准高，预期寿命长,及成长率低的发达国家属于此类型。(3）缩减型（老年型):图形表现为顶部宽，底部相对较窄,显示一种负的人口成长结构。通常发生在长期死亡率超过出生率时。这种类型的人口通常面临低生育率和老龄化的问题。（a）增长型(ｂ）静止型–7–(c）缩减型图3人口结构的三种类型Fｉｇ。3Tｒｅｅtypｅsofpｏpulaｔｉonstrｕｃｔure与人口金字塔相比,人口年龄树有明确的物理意义。它通过年龄树的生长变化来模拟人口结构的发展过程；用年龄树枝叶的生长来表征人口的年龄增长;用年龄树枝叶的凋零来表征人口的死亡；用人口年龄树树形的健康状况来反映人口结构的健康状况。对照图1中的上海市人口年龄树可以发现，上海市目前的人口结构属于缩减型，顶端枝叶茂盛，而底端枝叶稀少，这样的人口年龄树已经与发达国家的形状相同。因此，上海市的人口年龄结构也面临与其它发达国家一样的问题，也就是人口负增长以及人口老龄化。2.3

人口发展方程按人口控制论的观点,人口状态指某一区域内人口按年龄分布的状况，人口状态随时间变化的过程叫人口发展过程.要对人口结构的变化趋势进行研究，就必须–8–[1］［2[1］［2][3][4]［5］的人口数进行预测，并进一步得到各项人口指标.通过将这些人口指标与决策者关注的问题相联系，可以对人口问题进行定量与定性分析,为政策措施提供相应的参考方案.根据不同时期不同地区人口发展的特点，专家学者建立了各种人口预测方法来模拟人口发展过程，如一元线性回归法、自回归法、指数函数法、幂函数法、多元回归模型法、灰色系统GM（１，１)法、系统动力学法等，上世纪7０年代末80年代初,宋健、于景元等人建立了人口发展的偏微分方程，将我国的人口研究从定性分析引入定量分析，对人口数量和出生率、死亡率等人口指数进行了预测，因此应用最为广泛和成功。由于人口政策的变动和生育习惯的改变，现今人口的生育模式、死亡率等参数已经发生了巨大的变化，当时建立模型的结论和条件已经不适用于现在的情况,本文在宋健等人工作的基础上,重新考虑现在的人口状况对人口动力系统的影响，计算并给出各种人口指数。2．3．1

人口发展方程的连续形式引起人口年龄结构变化的三大要素是出生、死亡、迁移,人口发展方程描述了人口年龄结构与出生率、死亡率和迁移人口数的关系，从动力学的角度反映了人口发展的变迁过程，人口发展方程的连续模型为:–9–∂p(ｒ,ｔ）∂r

ﻮ∂p（r，t）∂ｔﻮ−µ（r，ｔ)ｐ（r，t)g（ｒ,ｔ）p（ｒ,t)pﻮ0ﻮ（r）ｒϕ（t）β（t)ｈ（r，ｔ)k(r，ｒϕ（t）β（t)ｈ（r，ｔ)k(r，ｔ）p（r，t）dr∫１式中,p（ｒ，t):人口年龄分布密度函数，ﻮ······························································（2．１）ｐ（r,ｔ)∂N（r，t）/∂ｒ。N（r，ｔ）表示ｔ时µ（ｒ，t）N（ｔ）µ（ｒ，t）N（ｔ）ｈ(r,ｔ)k(r,t)g(r，t）：某地区t时刻年龄为r的人的死亡率ϕ（t)：t时刻某地区单位时间内出生婴儿总数：t时刻某地区人口总数β（ｔ)：妇女平均生育率，即总和生育率：妇女生育模式函数：性别比例函数：某地区ｔ时刻迁移人口数r1,r２

：妇女育龄区间，通常ﻮr１ﻮ=15，ﻮr2ﻮ=49从控制论的角度来看,人口状态β（t）ϕ(t）ﻮp（r,t）ﻮ为输出量，通过可控变量妇女生育率可调控出生人口数ﻮ,进而影响人口状态,形成反馈控制(图4）。–10–２.3.2ﻮ图4人口状态方程控制框图Ｆiｇ。４Cｏnｔroｌｃｈaｒtofpopulａtioｎｄｅvelopmenteｑuations人口发展方程的矩阵形式［6］根据年龄移算理论,可以从某一时点的某年龄组人数推算一年（或n年）后年［6］龄相应增长一岁（或增长ｎ岁）的人口数.在这个人口数的基础上减去相应年龄的死亡人数，并考虑人口迁移，就可以得到未来某年龄组的实际人口数.对于0岁的新生人口,则需要通过生育率作重新计算.当社会经济条件变化不大时，各年龄组死亡率比较稳定，相应活到下一年龄组的比例即存活率也基本上稳定不变。因而可以根据现有的分性别年龄组存活率推算未来各相应年龄组的人数.即，若某t年年初有r岁人口数xr（ｔ)人，次年(t＋1)年年初这些人长了一岁为(r+1）岁.若µｒ(t）为这批人在一年内的死亡率，ｇr(t)为迁移进来的人口数，则（t+１）年年初(ｒ＋1）岁的人口数为xr（ｔ)×（１−µr(t)）gr（t）。0岁人口数x0（t)需通过妇女生育情况另行计算，人口发展方程数据流图见图5。–11–图5人口发展方程的数据流图根据上述推导，以一岁一组的人口统计数据为观测值，将人口发展方程离散化,用矩阵形式表示，可以得到人口发展方程的另一种形式：ｘ(ｔ1）H（t)ｘ(t）β（t)Ｂ(t)ｘ（t）G(t）x0(t1)（１−µ00（t）)Ψ（t)Ψ（ｔ)β（t）ﻮr∑rｉr1

ｋ（ｔ）hｋ（ｔ）h（t）x(t)ﻮＴ

mx（x（t）（ｘG（t)（gﻮ(t),x（(t),x（t),……x０(ｔ），g1(t),……ｇｍ

（t）)（ｔ））ﻮTﻮ为最大存活年龄–１２–ﻮ０ﻮ0ﻮΛ

0ﻮ００１−µ(ｔ)１−µ(ｔ)1H(t)0ﻮ１１−µ２ﻮ（ｔ）ﻮΛΛΟﻮ００1−µｍﻮ（t)

000000

Λﻮ0

b１bﻮ(t)ﻮΛﻮｂ2ｂ

（t）ﻮ0ﻮΛﻮ００B(ｔ）ﻮ０Μ００ΜﻮΛΛﻮ０Μ0ﻮ0Μ0ﻮΛΛﻮ0Μ0

０Μ０ﻮΛΛ

０0０ﻮ·········································(2.2)bｉ(t)(1−µb（t)ﻮ00ﻮ（t)）（1−µﻮ0ﻮ（t)）kﻮｉ

（t）ｈi（t),iｒ1…ｒ2ﻮ，iﻮ：某地区t年i岁妇女所对应的生育率本文以人口发展矩阵方程为主模型，并在此基础上进一步建立生育率、死亡率、迁移率的子模型.２。４ﻮ人口预测方法概述２。4。1ﻮ自回归法假定人口发展过程近似于直线状,由于当前人口数量受到历史变化状况的影响，可以用自回归模型来预测未来人口数.X(t１)ａbＸ(t）···························································（2．3)其中,Ｘ（t）为人口数,ａ、ｂ为参数。２.４.２ﻮ指数函数法假定人口发展过程近似于指数状态,前一段时间内发展缓慢，越往后人口增长–１３–越快,可以用指数模型预测人口数。X(t)ae其中,ａ、ｂ为参数.ﻮbtﻮ····································································（２．４)２．4。3ﻮ幂函数法b假设人口随时间变化曲线前部分斜率大，后部分斜率逐渐减小,可以选用幂函b数来预测。X（ｔ)ａｔ·····································································(2.5)其中,a、b为参数。2。４.4

多元回归模型法人口系统除了人口本身，还受经济、政策等各种要素的影响.人口发展是人口与各要素之间相互关联的结果。如果将这些因素都考虑在内，预测未来人口，称为多元回归。Ｙbｂ1x１b2x２Λbnxn·····································（２。6)其中x1、x2…xn为影响人口的各种因素，ｂ、b１…bn为参数.2.4。5灰色系统GＭ（１，1）法在没有规律可寻或资料不全的情况下，可以用灰色系统型为−aｋﻮＧM（1,1)进行预测,模ｘ（ｋ1)(x（1)⋅u/ａ）⋅ｅ式中，ｘ（k)对应一组人口数据序列ﻮｕ/a·······························(２。7）2.4。6ﻮ［２９，［２９，3０］–1４–时间序列预测方法的基本思想是:预测一个现象的未来变化时,用该现象的过去行为来预测未来。即通过时间序列的历史数据揭示现象随时间变化的规律，将这种规律延伸到未来，从而对该现象的未来做出预测。时间序列预测法是一种重要的预测方法，对资料的要求比较单一，只需变量本身的历史数据.但对数据的完整性及样本数量要求较高。用于人口数据预测的时间序列模型主要有自回归移动平均模型ARＭＡ（ｐ,q）、自回归综合移动平均模型AＲＩMA（ｐ，d,q）等几类。2.４．７ﻮ［２8]神经网络法［２8]人口增长具有非线性动力学特性，BＰ神经网络有很强的自学习、自适应能力，可以克服某些人口预测方法中的人为随机因素。它通过对连续或断续的输入作为状态响应而进行信息处理，采用逼近的方式来解决问题。神经网络是一种由多个神经元以某种规则连接而形成的层次网络结构,通过向环境学习获取知识并改进自身性能,即通过对输入样本的学习训练建立输入与输出之间的非线性映射关系。利用神经网络预测就是把已知的一个序列作为输入值,把要预测的值作为输出值,对该序列进行学习训练，构造网络建立非线性映射来逼近该序列的实际变化规律,对输出结果进行预测.ａ1ﻮk

Λａﻮnﻮkﻮ为观测到的人口时序数，作为输入层，输出层c1ﻮｋﻮ为要预测的数据,建立输入节点，输出节点数分别为n、１的３层BＰ神经网络模型（图6)。选取一组输入输出数据作为训练样本输入网络,经过正向传播输出和反向传播修正权值完成网络学习的过程。然后进一步输入下一组人口时序数据得到预测值。–15–Ｗﻮ11

Ｗﻮc１cﻮWｐ１Wｂﻮ1ﻮVﻮn1ﻮ…

Ｖﻮ１iﻮi1bＶﻮｉ

Vﻮniﻮ…ﻮVﻮ1ｐﻮｂpｂﻮWV

11ａ

Ｖﻮh1…ﻮｈiａﻮ…ﻮＶﻮhpﻮaﻮVﻮｎｐﻮV1aﻮｋﻮhaﻮkﻮnａﻮk１ﻮh图６三层BP神经网络的拓扑结构ﻮｎFig.6TｈetopｏlｏgｙｓtｒucturｅofＴｈreelａｙｅrBPNeuraｌNｅtｗｏｒｋ自回归、指数函数、幂函数都要求人口数据具有明显的规律性，只有对于特定地区的人口状况才适用。多元回归模型涉及较多影响因素，而如何对这些因素进行定量是研究的难点。灰色系统、时间序列和神经网络法适用于不同类型的人口数据,但后两种方法需要大量的历史数据,相对操作也比较复杂。因此本文的预测主要采用灰色系统GM（1,1）模型。–1６–第三章生育率模型3.1ﻮ生育率指标说明［７］生育率反映了育龄妇女的生育能力和生育水平，按参照标准的不同，对生育率［７］的定义也不同。人口学上用来衡量生育率的指标主要有年龄别/组生育率，世代生育率和总和生育率。年龄别/组生育率指某个年龄(或年龄组)的育龄妇女在某一年的平均生育率，可以通过人口统计资料直接获得，能够反映生育率的年龄分布趋势.世代生育率(ＣomｐletｅdCｏhortＦertilitｙRaｔｅ）以每一世代各年龄层的方式,记录一个妇女一生中生育的子女总数,反映各世代实际的生育率变化,是生育率的队列指标.总和生育率(TｏｔａlＦeｒtｉlｉtｙRaｔe）假设妇女按照某一年的年龄别生育率度过育龄期，平均每个妇女在育龄期生育的孩子数被称为总和生育率。实际上,它就是假设一个妇女在整个育龄期都按照某一年的年龄别生育率生育的情况下，一生所生育孩子的总数,是生育率的时期指标.–１7–图７世代生育率（ＣＦR)与总和生育率（TFR)的比较Fｉｇ。7CoｍpaｒｉsｏnofＣFRａndＴFRﻮ[8］要统计世代生育率必须对一代人的整个生育年龄区间的生育行为进行跟踪，因此总和生育率使用更为普遍.总和生育率和世代生育率的主要区别在于参照的时间体系不同，总和生育率可以视为对世代生育率的一种近似。但实质上，两者的内涵有显著差异（图７）.生育率的分析方法也可相应地分为两大类，即队列分析法（以世代生育率CFR为指标）和时期分析法（以总和生育率TFR为指标）。从理论上来说，队列分析法更能反映妇女一生的生育水平，但总和生育率更容易获得。因此，时期分析法使用更广泛.在本课题研究中,以总和生育率来表征生育水平的大小。３．2ﻮ生育率模型生育率模型可以分为两大类

[9］.第一类为演绎模型,即从某些理论假设出发来[9］,推导出的模型。例如，１９7４年A.Ｊ．Coale和T．J.Trｕssell提出用模型生育率表来研,究妇女生育规律。19７８年,法国人口学家ＪohnBonｇａaｒｔｓ根据生育率和生育率–１８–[12-１６]直接决定因素(如结婚、避孕、人工流产和产后不孕概率等)之间的关系建立了综[12-１６]合生育率模型.第二类模型为归纳模型,其特点是通过对大量统计数据的分析,找出共同点来。此类模型往往有明确的数学形式，便于计算和仿真，也是本文研究的重点，以下是几个使用较广的归纳模型.3。2.1ﻮ[17－２0］Cｏmｐeｒtz模型[17－２0］BxxCｏmｐｅrtz函数早先用于人口死亡率的研究，随后被推广到了生育率的研究上.BxxWunsch(1966）和Mａrｔin(19６7)评估了将Compeｒｔｚ模型运用在生育率研究上的可行性。Compｅrｔz模型的数学形式如下:Ｆ(x)FA·····························································(3。1)F（x）∑f（i)·····························································(3．２）ir1其中，x：年龄;F(x）：累计生育率；F：总和生育率;f（i)：年龄为ｉ的妇女的生育率;Ａ、B：描述生育率趋势的参数.3.２.2

[20—22］Gammａ[20—22］αβｆ（ｘ)xeＫｘｅΓ（ααβｆ（ｘ)xeＫｘｅΓ（α）时间序列模型估计Gammａ函数的参数以预测年龄别生育率。宋健亦将Ｇammａ函数应用于中国的生育率数据，得到较好的拟合效果。Ｇammａ函数数学形式如下:···································(3。３)其中，α—1表示函数转折高度系数；β为函数转折收敛系数，代表函数尾端的收敛速度。３。２．３ﻮ［23］Loｇnormａl模型［23］–19–[2４］[2５,2６]2［ln（x−ｘ）−µ[2４］[2５,2６]2［ln（x−ｘ）−µ］０ｆ(ｘ）k××ｅxp｛−｝22σ２σｅxp（µ）x基础上构造生育率模型.Lｏgｎｏrmal（对数正态分布模型）最早由黄荣清提出,其数学表述为：·················(3．4）式中，x０为起始生育年龄，本文令x01４岁；kl为对数正态分布函数所对应的尺度变换因子,与生育水平相关。由模型计算得到的期望生育年龄为：2参数µ和σ决定了分年龄别生育率曲线的形状:在µ相同的条件下，σ越小，到达期望生育年龄的时间就越短;在σ相同的条件下，µ越小,期望生育年龄周围生育率越高.3．2。４ﻮPｏｉsｓon分布模型（x−x（x−x)−λλｅｆ（x）ｋ×·················································（3.5)描述分年龄别生育率。模型的数学表述为：(ｘ−x0)！该模型的期望生育年龄为:λｘ0式中，起始生育年龄x0也取14岁；kp是泊松分布函数所对应的尺度变换因子;λ表征了从起始生育年龄x0至平均生育年龄的时间长度。若λi为第i孩生育率模型的参数,则λi１−λｉ可以表征生育第i孩至生育第i1孩的平均间隔时间的大小.３。２。５生育率组合模型–20–2（x−x）2（x−x）−λ［ｌｎ（x−ｘ）−µ]0f(ｘ）ｋ×{αeｘｐ{−｝（１−α）}····························（3.6）2（x−x）!［27］[23]法满足不同生育率数据对模型的精度要求。为了改进模型精度,扩大模型的适用范围，本文构造了分年龄别生育率组合模型(ｃoｍｐoｓiteferｔｉlityratｅmodel）：σ２π０除了上述五个模型以外,国内外一些学者还研究了双指数曲线、多项式等数学函数对生育率的拟合效果.研究发现，在一般情况下，某一地区的生育率年龄分布模式是不会随时间变迁而发生明显改变的，因此可以使用一些简单的函数（包含较少的参数，一般为４—５个）来模拟生育率的年龄分布。另一方面，模型的精度主要取决于实际生育率统计数据与所用模型函数的曲线形状的接近程度。用不同模型拟合不同地区的生育率数据，往往会得到不同的效果.研究表明，由于中国当前的生育率模式与以往的生育模式有较大不同，Ｃｏｍｐeｒtz模型和Gamma模型对年龄别生育率的拟合精度不如Logｎormal模型。因此，本文将重点比较Lognormal模型、Pｏｉsson分布模型和生育率组合模型对中国及上海市分年龄别生育率的拟合情况以及这三个模型在生育率预测应用上的复杂度.3。３ﻮ模型验证和比较3。３。1ﻮLognｏrｍal模型与Ｐoisson模型验证[３１，3２］[３１，3２］４９最小二乘法,优化目标函数为:min∑(yx−ｆｘ１５

２对上述模型进行验证，参数辨识使用２x）即要求实际生育率ｙﻮｘ

与模型输出fﻮx

的误差平方和最小。验证结果如表1所示，表中–2１–k年份孩次k年份孩次期望生育年龄/岁ＲMSE表１Lｏｇnｏｒmal和Ｐoisson分布模型描述分孩次的年龄别生育率Tabｌe1ComｐarisonｏfPoissonｍｏdｅｌａｎdＬｏｇnormalmoｄelforage-sｐeｃificfｅｒtｉlitｙrateslLogｎorｍａｌﻮPoiｓsonﻮLｏgnormalﻮＰoiｓｓｏn１9９０ﻮ1ﻮ１.02４ﻮ1.１１４ﻮ２3.４4ﻮ２３.20ﻮ0.００4９ﻮ0。0１０62ﻮ０．71６ﻮ０．6９9ﻮ26．７8

26．16ﻮ0．0030ﻮ0.０0443ﻮ0。284ﻮ0．282ﻮ29。２2ﻮ28．6７

0.0０１8ﻮ0。０0２４1９9４ﻮ１ﻮ0.92７ﻮ1.038

２4。16ﻮ2４．００ﻮ0。0０3８ﻮ０.０１１32ﻮ0。４43ﻮ0.430

27.5７ﻮ2６．9４ﻮ0。０02９

０。００293ﻮ0.114ﻮ0.107

29。２７ﻮ28。５2ﻮ０．０0１３ﻮ0。０0132０００ﻮ１ﻮ0。８６0ﻮ0．９１５ﻮ24。０7

23。76

０.0０22ﻮ0．0０5９2ﻮ0．28９ﻮ0。２７8ﻮ2８．96

２8．３０ﻮ０。00１9

0。00１１３ﻮ0．06５ﻮ0.0５7ﻮ３0．55ﻮ２９。33

0．００0１

０。000４RMＳE

指标反映了优化目标函数值的大小，本文将其作为衡量模型精度的指标。图８所示为不同模型对20０0年一孩、二孩、三孩生育率统计数据的拟合结果。就中国历年生育数据而言，在描述分孩次的年龄别生育率时，对数正态分布模型的精度一般高于泊松分布模型（也有例外，见图８(ｂ））。–22–（a)一孩–23–（b)二孩(c)三孩图82000年分孩次的年龄别生育率模型拟合结果与统计数据比较Fｉg．8Comparisonｏfｔhｅmoｄeｌwithsｔaｔistｉｃdaｔafoｒage—sｐeｃificfertilｉty由图8可知，对数正态分布模型和泊松分布模型反映的生育率变化趋势是一致的。在１９90年至2０0０年间,生育水平呈下降趋势，二孩和三孩生育率下降幅度更为明显。就期望生育年龄来看，妇女生育一孩的平均年龄介于23、24岁间,随年份的变化不大。生育二孩和三孩的平均年龄随着年份的增加有明显的上升趋势.这种变化趋势完全符合我国的实际情况，说明运用对数正态分布模型和泊松分布模型来描述分年龄别生育率分布是可行的.３。３．２

组合模型参数辨识–２４–2组合模型使用非线性最小二乘法辨识参数。非线性最小二乘法是否收敛，很大2程度上依赖于待定参数起始值设置的质量。参数搜索的起始值越接近真值，计算结果精度越高，收敛越快。对于复杂的非线性模型,使用最小二乘法时往往会有多个局部极小点，设定的初值不同,可能会收敛到不同的局部极小点.因此，对于组合模型而言,待定参数初始值的选取是非常重要的.由于生育率组合模型是描述妇女分年龄别生育率分布的模型,各个参数都有明确的物理意义，在设定初值时必须予以充分考虑。参数µ、σ、λ的初值可由妇女期望生育年龄推算。分析８0年代以来的妇女生育统计数据,期望生育年龄一般为23至26岁（由妇女生育峰值年龄近似得到)。而对于Lｏgｎormａｌ模型和Poｉsson分布模型，期望生育年龄分别为exｐ（µσ/2）１4以及λ14。因此，λ初值可取9～1２.由于σ的范围一般为０~1，则µ初值可取1。7～２。5。尺度变换因子kc是一个与生育水平高低密切相关的参数，而总和生育率是衡量生育水平最常用的指标,因此可以用总和生育率作为kｃ的初值。通过调节α值,能使组合生育率模型更贴近分年龄别生育率分布的实际情况。由于对数正态分布模型在描述一般生育率时，精度要比泊松分布模型高(见图8），故将α的取值范围设定为0。5～２.α>1时，加权和变为加权差的形式。综上所述，可将组合模型的初值设定为ｋc=总和生育率，µ=2.4，σ=0．５，λ=11，通过调节权值α，即可得到模型参数的估计值。经反复验证后发现,选择不同的k值和α值对算法的收敛影响最大，k的初值为总和生育率时，有很好的收敛性。对α初值的选择则直接影响到计算精度。表2组合模型参数的估计值（最小二乘法)Taｂle２Estimateparameｔeｒvaｌｕeｓofcｏmpｏsｉteｍodelbyuｓｉnｇleasｔsquarｅｍethodｓ–25–年份

αﻮkﻮcﻮµﻮσ

λ1986年

1.21５０

２．3９８8ﻮ2.4１09ﻮ0.3７49ﻮ１2.５75019８9年ﻮ1。24９3ﻮ２.24２9ﻮ2.４02２ﻮ0。３8４5ﻮ12.5３１0１990年ﻮ1.4040ﻮ2．1997ﻮ2。4２０２

０．3534ﻮ1１.７3９519９4年ﻮ1。6６６7ﻮ1．54１８

2．449７ﻮ0．3１4５ﻮ12．6３5８1９９5年ﻮ1.371２ﻮ１。４１８６ﻮ2。35０6ﻮ0。3４3１ﻮ12．０79819９7年

1。60９3

1。4532ﻮ２.4４06ﻮ０。３０７1ﻮ12。275０１9９８年

1．5474

1.4４9４ﻮ2.4372ﻮ０.307８ﻮ1２．０4362000年

1。27８2ﻮ1。２1５４ﻮ2．３８66ﻮ０.３3８2ﻮ11。8３192003年ﻮ１．3２96ﻮ1。４1３５

2.４0６0ﻮ0．3４33ﻮ１1.４4673．3.３

模型精度比较–2６–图９组合模型、对数正态分布模型及泊松分布模型精度比较(均方根误差)Fig.９Coｍpariｓonｏfaccｕｒacyofｔhｅｔhｒeｅｍodels相比对数正态分布模型和泊松分布模型，组合模型的均方根误差更小,精度得到明显提高（图９）。另外，利用分年龄别生育率计算出生人数及出生率时，要分别乘以各年龄别平均育龄妇女人数或各年龄别育龄妇女人数占育龄妇女总人数的比重，因此生育高峰年龄段内（22—３2岁）生育率的准确性显得尤为重要.而组合模型在生育高峰年龄段内的分年龄别生育率误差在三个模型中一般是最小的（图10)，进一步显示了组合模型的优点。–27–图10不同模型分年龄别生育率误差比较(200３年抽样数据ﻮ［3１］

）Fiｇ.10Cｏmparisonｏｆaｇｅ—ｓpecificfｅrtiｌｉtｙｒatｅserrorｏfｔheｔhreｅmodｅlsfｏr20０3Logｎｏrmａｌ模型和Pｏissoｎ分布模型数学形式简单,物理意义明确。Pｏissoｎ分布模型参数最少，易于辨识参数，大大减少了参数预测的工作量,只需预测两组参数.组合模型精度最高，考虑到实际生育率不一定满足特定分布,该模型适用范围最广,但参数个数较多。Ｌogｎorｍaｌ模型参数个数适中，拟合精度较高。这三个模型可以根据其优缺点及对实际数据的拟合情况选择最合适的模型进行应用。３。4ﻮ生育率预测每一年的分年龄别生育率均有35组数据，即使将其压缩成以每五岁年龄为间隔的年龄组生育率,也有7组数据。如果以不同年龄别或年龄组的生育率为预测单位，–28–［3３－36］预测的工作量会很大。198７年Ｂoｚik和Bｅｌl最早提出将主成分分析法应用于生育率［3３－36］，其基本思想是将所有的分年龄别生育率视为多个变量,透过主成分分析缩减维度简化计算.由于各主成分间相互独立，可使用回归分析或时间序列等方法对每一主成分进行预测，将多变量预测问题转化为若干个单变量的预测。基于此原理,可以利用合适的生育率模型拟合历史统计数据,以辨识得到的各个参数作为观测序列,即要预测的主成分。分别对这些参数进行预测。然后，只需将各参数的预测值代入生育率模型，即可得到分年龄别生育率的预测值。无论是Loｇnｏｒmaｌ模型、Pｏｉsson分布模型或组合模型都是基于随机分布函数的生育率模型，具有类似的数学形式，可以概括为:f(x)k×h（x）·····························································(3.7）其中，x指妇女生育年龄;f（ｘ)为生育率;h(ｘ）是特定的生育模式近似函数；k为尺度变换因子，与生育水平高低有关。假设妇女的生育意愿和期望生育年龄在20年内不发生显著变化,则可认为生育模式在这段时间内保持恒定，只需要对ｋ值进行预测即可。3.4.1ﻮ总和生育率预测f(x）TFR×ｈ（x)如果将分年龄别生育率ff(x）TFR×ｈ（x)························································(3.8）即分年龄别生育率为总和生育率TFR和归一化后的生育模式的乘积.基于之前的假设，只需将原始生育率先做归一化处理，则对年龄别生育率的预测就可以简化为对总和生育率的预测。由于总和生育率的历史数据比分年龄别生育率的历史数据完备得多，在缺乏统计数据的情况下，用这种方法可以方便地得到分年龄别生育率–２9–ｋｄＹ（ｔ)aＹ（ｔ）u·······················································（３．1０）−a(k−ｋｄＹ（ｔ)aＹ（ｔ）u·······················································（３．1０）−a(k−１）a令p１−（Y（1)Ｙ(2））1２1−(Y（2)Y(3））１B2ΟΜ1−（Y（ｎ−1)Y（n)）１′ａ−1p(B’B）B＇X·············································（3.1２)以总和生育率TFＲ为观测序列，应用灰色模型GＭ（1,1)进行预测。(1)建立GM（1，１)模型原始数列X(ｋ）｛x(1）,x（2)，Λ,x（ｎ）｝一次累加生成向量Y（k），以弱化数列的随机性,强化其规律性.Y（k）∑ｘ（i)k1，2,Λ,ｎ·································(３．9）i1GＭ（１，1）模型的一般形式为:ｄｔ其解为：Y（k）［Ｘ（1）−u／ａ］eu／a·······························（3．1１）其中a、u为需要辨识的参数。（２)参数辨识u2ＸNx（２）x(3)Λx（ｎ）u–３0–将辨识得到的参数a、u代入式3。１1，得到Ｙ(k１）的估计值Y(k１)ˆ则待预测值X(k1）的估计值Ｘ(ｋ1）为ˆˆ(k１）Ｙ(k1)−Y（k)···········································（3。1３）ˆˆ(３)预测实例以１994年至200４年的中国总和生育率作为观测数列进行预测。表３总和生育率GＭ（１,1)模型预测结果Ｔａｂle3ＦorecaｓtedvａｌuesｏfTFRbyGM(1，１)年份199４１９９５19９6199７１9981９9９20０020012002２０0３2０04原始值１．５６１。４31．551.46１.46１.4５１.40１.39１．391。4１1.４5估计值１。56１.481。471。4６１．４5１.4４1.43１.431．42１。411。40年份预测值ﻮ2０0５１．３9ﻮ20061.３8ﻮ2０071．38ﻮ2００81.37ﻮ20091．36ﻮ２０１01。３5ﻮ2011１.3４ﻮ201２1．34ﻮ2013１．33ﻮ2014１.３2

2０１５1.313.4．2

分年龄别生育率预测根据前文的假设,以２003年的生育模式为基本生育模式，分别预测２004年及2０1５年分年龄别生育率。结合表２和表3的各项参数估计值,计算所得的２００4年分年龄别生育率及2０15年分年龄别生育率如图11、１2。–31–图１1不同模型对2004年分年龄别生育率的预测估计Fig.１１Forecａｓtedvaｌuｅsｏf2０0４ａｇe－speｃｉfiｃferｔilitｉｓｂyｄｉｆfｅrｅnｔmｏdels图１２组合模型对20１５年全国分年龄别生育率的预测Fig。1２Fｏｒecasteｄvaｌuesｏｆ2０15ａge-specificｆertilｉtieｓbyｃomposiｔｅmｏdel–3２–第四章死亡率模型４。1ﻮ人口死亡风险主要度量指标−nＭ−nＭ(ｘ,xn）µ(x)ｌimq（x）·························································（4。3）都表征了个体死亡的可能性,但含义各不相同。死亡力µ(x)(Fｏrｃeofmortality)：又称瞬间死亡率,指某一确切年龄的死亡概率水平。分年龄死亡率Ｍ（ｘ，xn)（ａｇespｅｃifｉcｍoｒtａｌitｙ）:某年度从x岁至ｘｎ岁之间的死亡人数，可以直接由统计数据获得.分年龄死亡概率（agｅsｐｅciｆｉｃprｏbabiｌityofｍｏrtality）：已活到ｘ岁的人有多大比例将在到达ｘn岁之前死亡.死亡力µ（x）数据难以得到，实际应用中一般使用平均死亡力µ(x）或累计死亡力.µ（x），M(x,xn)和ｑ(x,xｎ）虽然都是描述死亡率的指标,但三者的物理意义有明显区别。若假定活到x岁的存活人数随年龄x的增加呈指数递减，则有以下关系：q(x，xn)1−e·············································(4。1）µ(x）M(ｘ,xn)························································（４.2)∆x−0–３3–4。２

数据修匀[3７］所谓数据修匀，即根据一个连续变量的不规则观察序列,通过修匀的方法得到[3７］一个光滑的有规则的修正序列,使其与观察值序列相和谐。当要修匀的数据范围比较大,用单个函数形式不能得到满意的修匀时,可以将数据分段。在不同的子范围上，用不同的函数拟合它，这种修匀方式通常称为样条修匀.１980年Bｅnjａmin和Pｏlｌaｒd将样条函数引入修匀理论中，首次用自然三次样条解决了人口死亡率的修匀问题。在很多情况下,死亡率的统计数据以分年龄组的数据形式存在,数据修匀可以把这些被压缩的数据恢复成分年龄别数据.并且，数据修匀可以调整和修改异常值，以保持死亡率曲线的光滑性和必要的连续性。表４２003年中国分年龄组死亡率统计数据Ｔaｂｌe４Chｉnesemortａｌitysｔａtｉｓticｓｆor2003年龄男女0—43．74。1２５-90．60．49１0－１4０.3６０。5215－19１。02０．4320-２41．170.９32５－2９１。1３０．81３0－341．４０．8735－３92。１4１。２4０—44２.881。６84５-４９3.722.６2５0-546。31３。3２55—5９9．875。756０—6４1５。8６8.3４65－６92３。61１8．42–3４–70－7475—798０-84８5—8９９０＋ﻮ４4．187３。５2115.11７４.９２９0.7８ﻮ30.514９．2６9２．931１７。４4289．5７以每一年龄组所对应的死亡率数据作为三次样条插值的连接点，用Mａtlab的插值函数ｉｎterｐ１可以方便地将分年龄组数据还原为分年龄别数据.以2003年中国人口死亡率数据为例，插值结果见图１3、１4。图１３20０３年中国人口死亡率三次样条插值结果与实际数据比较（男）Fｉg．13Ｃｏmpariｓｏnｏfｔheｉnterpoｌａtevａluesaｎｄrealmｏrtaliｔｙfoｒ200３(ｍale）–35–图142０03年中国人口死亡率三次样条插值结果与实际数据比较（女)Ｆｉｇ。14Ｃompaｒｉsonｏｆtheinｔerpoｌａｔevａluesａndreａｌmortalityfor200３(female)由于连接点处的数据使用分年龄组数据，即该年龄组平均死亡率,而非连接点本身所对应的实际死亡率数据，因此在低龄段尤其是高龄段有一定的误差，而在5～80岁期间,通过插值得到的数据与实际数据非常接近.4。3ﻮ死亡率模型４.3。1

死亡率模型综述人口死亡模型主要分为两大类，即间接度量模型和直接度量模型。６0年代，英国人口学者W.Brass提出了logiｔ模型生命表。该方法认为不同生命表上的存活概率经过lｏgiｔ变换后存在一种近似的线性关系，可以利用一个已知人口的生命过程来度–36–［3８，３9]［４０]βxβxµ（ｘ)αeγ［3８，３9]［４０]βxβxµ（ｘ)αeγβµ(x)αｘβxαeµ（ｘ)···················································(４.7）βｘ在loｇit体系基础上又发展为Bａsiａ和Eｗbank的4个参数模型。Ｂｒaｓs的lｏgit体系作为人口死亡分析的方法，被中国人口学者广泛使用在人口分析上。我国学者黄荣清经研究发现，在一定条件下，同样是2个参数的l（x)的双对数模型的稳定性要比lｏgiｔ体系更好一些.死亡的直接度量模型就是以年龄为自变量、死亡风险（通常为死亡力或死亡概率）为因变量的函数模型.其中以Gomｐeｒｔz(１825）的研究影响最为深远。１８６0年麦可海姆(Mａkeｈaｍ)对该法则进行了补充。在此基础上，人口学家提出了各种不同的死亡模型。在以下各模型中,x为年龄,µ（ｘ)为死亡力函数（1）Ｃoｍpeｒtz模型:µ(x)αe·························································（4。4)（２）Maｋehａｍ模型·························································(4.5）（３)Ｗｅibuｌl模型:···························································（４。６)此模型不考虑０岁婴儿死亡率.(4）Kannisｔo模型:１α(e−１）观察死亡率曲线可以发现，死亡率的年龄分布通常为U型(死亡水平较高时）–37–[４1]［4２]［43]［44］或J型(死亡水平较低,如图[４1]［4２]［43]［44］老年期死亡率随年龄的上升而上升，而青壮年期的死亡率变化则比较平稳。经研究发现，以中国的死亡率数据为例,在儿童少年期,Weｉbull模型的精度要高于Cｏmｐertz模型,在高龄期，Kannｉsto模型的拟合效果好于其余模型.Cｏmpertz模型低估低龄段死亡率，高估高龄段死亡率。Mａkeｈａｍ模型改进了对低龄段死亡率估计的准确性.因此，用上述模型对全年龄区间的死亡率进行描述并不合适,针对这一问题,一些学者提出了全年龄区间的死亡率模型。根据人从出生开始到生命结束的死亡风险，按不同年龄阶段，一般概括为三种不同形式的作用力.在死亡率模型中,从Thｉeｌｅ（187２)到Rogers＆Ｐlank（１９83）的模型都把死亡力设计成三种不同函数之和：即一个递减函数（儿童少年期）、一个递增函数（老年期)和一个先增后减的函数(青壮年期).（５)Ｔ.Ｎ.Tｈile模型：µ（ｘ）A1eﻮ−B1xﻮA2ｅﻮ−0．5B２（x−C）

2ﻮA3eﻮB3xﻮ·························(4．8）Ａ(xＢ）Ａ(xＢ）q（x)eﻮcﻮDeﻮ−Ｅ(lnｘ−lｎＦ)

2ﻮﻮGHﻮｘﻮｘﻮ·····················（4.9）−Cx其中,q（x）−Cx(７)Rｏgｅr＆Ｐlank模型ｑ(x）AＢｅ

Deﻮ−Ｆ（x−E）ﻮ−eﻮ−Ｇ（ｘ−Ｅ）ﻮＬｘHe1GＨＬｘHe··························（４。1０）从本质上说，死亡间接度量和直接度量两类模型之间并没有根本区别。所谓直接和间接不过是在度量这个过程时的比较参照系不同而已。我国学者就这两种模型–38–[4５]的内在联系进行了研究。[4５]４.3．2ﻮ分段死亡率模型以1９８９年、19９4年、1998年、２001年的中国死亡率人口统计数据为原始数据，通过最小二乘法或非线性最小二乘法分别对Comperｔz模型、Mａkｅhａm模型、Wｅｉbulｌ模型以及Kａｎnistｏ模型进行验证，并计算各年龄区间的均方根误差，结果列于表5－表8。儿童少年期青壮年期老年期ﻮ表5１98９年中国人口死亡率各模型拟合均方根误差比较Tabｌe５CoｍｐaｒｉsonoftｈemoｒｔａliｔｙRMSEbｙｄifferenｔｍoｄｅｌsCompertz模型Maｋeham模型Wｅibull模型０。０058０．005８0。00130.0０２２０．00110．００43０.0４1２0．00560。0８2３ﻮＫannｉｓto模型0。０0595．４1Ｅ－040。０05８儿童少年期青壮年期老年期ﻮ表61994年中国人口死亡率各模型拟合均方根误差比较Taｂle6ＣomparｉsｏｎofthemｏrtaｌiｔyRMＳＥbydifferenｔmｏｄelsＣｏｍpeｒｔz模型Mａkehａm模型Wｅibulｌ模型0。0１０。０10.0０１5０.00230.00１60．00４2０．0４480。0０880.0854ﻮKannisto模型0．０1０１9.4０Ｅ—０40.0083儿童少年期青壮年期老年期ﻮ表7１9９8年中国人口死亡率各模型拟合均方根误差比较Tａble７ＣomｐariｓonoftheｍｏrtalｉｔｙRMSEbｙdｉffｅrenｔmodeｌｓCoｍpｅrｔz模型Mａｋｅｈam模型Wｅｉbull模型０。008５０。0084０。００120.00240.０0140。０0４０。０4430．01０40。0785ﻮKａnnｉsto模型０.０0８６8。74Ｅ-０40．0１05–3９–儿童少年期青壮年期老年期

表8２00１年中国人口死亡率各模型拟合均方根误差比较Tabｌe8ＣｏmparisonoｆthemoｒtalityＲMSEｂydｉfｆｅrentmoｄｅｌｓCｏmpertz模型Makehaｍ模型Weｉｂuｌｌ模型0。00530．０0518．92E-０40.0020.0０1１０。0０340.03