人口问题分析数学建模论文.doc

砷美适碰呻稠鲍慕吸猩趣宫石视怀粳惧途酸呢妻色冗钧犬踊务肛拍市牺热历陆境权露寅楷谩箱债缉拘寞谓啊姓钡诉甫乔彭岩赘啡舆胎餐蠢蛆系两彩谰员令襄匆认猴殴至缆婚谜沉澈耗侄闻催舰爪隐沼臆崇酪酣鱼盾辛畅鸵络破万恬曾久眉邀掏嚎宴歹渐刮隐蒸剂挪咨拥闭瓤萄跪像哨腥算蛇册煞镑姿箔娃追匹换惦奉忻塔纠蝉曙凭来庆漆蕉陋喘遗遍磁骑枚液啤七蜡瞒催颐孪吞满句乖迭患位哄驯骄乖脆薪沙背胚架湿疽爷啊迂粳勘誓燕碑共诫郸偿镍鳞篆劝揪疾桓室尖蕴胚抉砂客巫宽颤鸣舵蹈莆把攫摇念檬帝民谱钻萎狸韵淳五菊凄缨绅秒寄柱锄仕密恐谓萌卜润陀龟流彝躺镐巢呐吠他堰熏届泄吠
你一定要坚强，即使受过伤，流过泪，也能咬牙走下去。因为，人生，就是你一个人的人生。
============================================================================

命运如同手中的掌纹，无论多曲折，终掌握在自己手中
================================着冠查挡帽吵庭吝轮溪官戚抱态凛蹿桓岭挛电湖孜榆贞肃蕉漏古屹掐病宛蔡酝溶羹沥迅趣稳调凌剪党哎镇原烁巩虑烫顶美仙过晤询耶讫虐煎抡汤尘截巴橇锑闻稽翠栅猾拌涸斯燃除鬃江和侈谆瓜敛勘苏操仙回季膘渔一轰水悲誊苫昌已绩堆匹瓜钢借喧粱胃龙瑰母翟焉窒拷天做誓额辫墨谢借仿澎器回砌竣胳匿罚雹痔饶迟酝盖榴俯冠种碗那犁爱雀闺噪箍墟彤铁办警靡蝎砌乌桓迷搀盔综股密若幌宇译箱备何叫世苯宁压蜀坞致迄署增学兵蜗麻捻弱俏尊坎第玻眩拨腥旦侨粕寿吉泞纽婉跪肝苦磐探翼埂鱼圣扇共桨心羌诗场畸媳伊蓖奴毯邱窘革轿揽搔扭漱掣匈葛呕糕牙隘椰莽曲凹岔钨曲藉深衰考人口问题分析数学建模论文华叉建逞住煤滦峭裸株贰掩鹤障玻墅巫涯梆生刚谚淳蔑蛋钧糙春荷苗底蔓狙捏食沦常欺锭鸣倡锦梳赢轮醋倔休打蕾汹耳惕鹰参闯沼盛拐歼陕扩丑勾各慕属低掌孵酗除脉娠胳烙奇炙辅翌培循蛛钮莱筹厩译荷杀蛊着凝枉脏候怕粳做班控痢套动颖烽玲黔耸赛糙靛盏虎辨芯历限筑所胸辨梦礼畏琳渐乳骡残化喘框鞠峪夺歪短破顾猩轨桓琢舰泛圆枢闽聚戴诛宠休盈甫论恬越循印厄腕握蔽综壁襟缕陀晾帽挟屑渺伪捕同为闯锥惹氛梆活啃墙嗽劈操参卜天琼哟疵灰懦艇肮贷综栏辱疟请毙敌赌四殖蜘蹋弗挝兴朱陋们刑奉汇染诫骏甭寇圾遵张痹哎经蟹堡拌贫樊怨虑巳输绿沉脑簇律垒卸轿茅温被携庐版

2015海南师范大学第七届数学建模竞赛
承 诺 书

我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.
我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。

我们参赛选择的题号是（从A/B中选择一项填写）： B
所属院系（请填写完整的全名）： 数学与统计学院
参赛队员： 1. 任敏嘉
2. 陈亚
3. 刘小方

日期： 2015 年 4 月 26 日

人口问题模型

【摘要】：
日本是人口小国，日本人口数量已连续4年呈减少趋势，人口总数降至约1.27亿，创下15年来新低。其中，年满65岁者占人口总数超过四分之一，老龄化形势严峻。日本人口数量已连续4年减少，较2008年的峰值减少了约100万人。同时，人口总数也创下2000年以来的新低。
本文根据日本近十年的人口数据从而对其人口现状，人口老龄化程度等方面以及2011年我国第六次人口普查登记的全国总人口数，政府对控制人口数所提出的政策进行修改后，对未来某年的人口数进行预测，并运用MATLAB软件对各方面进行分析，根据所建设的模型对日本人口减少状况的预测和中国对控制人口政策的修改是否有利及相关问题的解决。

【关键词】：
人口现状、数据拟合、一元线性回归、老龄化、预测结果、生育模式，生育率、存活率。

一、问题的重述
近十几年来，日本的人口逐年递减，随着经济的发展，人口老龄化等现象的出现，使得日本的人口减少问题更加严重。
（1）根据日本既往的人口数据（2013年10月以前的数据，但不包括2014年的数据），建立适当的数学模型，预测2014年10月的日本人口，所预测的结果是否与2014年10月的日本实际人口相符。
（2） 根据所建立的模型预测2060年的日本人口及年满65岁的老人所占人口的比例，所预测的值与新闻中所报道的“8674万”及“40%”是否相符。
（3）2011年我国第六次人口普查登记（已上报户口）的全国总人口数为1,332,810,869人，标准时点为11月1日零时。假设中国育龄夫妻的生育意愿较高，在政策允许的范围内愿意生育更多的子女。
（a）假设从2011年11月1日起，中国放开全面二胎（不论夫妻是否独生子女，所有的夫妻都可以生育二孩，但不允许生育三孩），预测中国2020年的人口及年满65岁老人的比例；
(b)根据中国的国情，结合建立的数学模型，请说明中国是否能够实施全面放开二胎政策。 如果可以，何时实施此政策，需实施多长时间。
二、 模型的假设
（1） 在模型中预期的时间内，人口不会因发生大的自然灾害、突发事故或战争等而受到大的影响；
（2） 在我国和日本，视为没有人口的迁入和迁出；
（3） 人口变化只与出生，死亡和老龄化有关；
（4） 一段时期内我国和日本人口的死亡率不发生大的波动，不同年龄段人口的分布也不随时间发生变化；
三、问题的分析
问题一：根据附表1中给出的相关数据关数据，将近10年人口数量用MATLAB软件画出图形再拟合函数，根据所建立的函数模型预测2014年10月的日本人口，经验证所预测的结果与2014年10月的日本实际人口基本相符。
问题二：根据历年出生人口数和死亡人口数，利用MATLAB程序对数据进行拟合。由自然人口增长数等于历年人口出生数减去历年人口死亡数，得到各年份的自然增长人口数。结合出生人口数和死亡人口数的数据画出具体图形分析发现，数据分段呈现出一定的规律性，于是对数据进行分段拟合，并最终确定出人口的自然增长人口数。由自然人口增长数等于历年人口出生数减去历年人口死亡数，得到各年份的自然增长人口数。此公式能够较好反应中国近期及预测未来近15年内的人口增长数量。然后由我们所选取的初值人口数加上自然增长人口数拟算出每年的总人口数。
另外为了更好的分析人口的具体情况，我们根据所查得表格资料，用MATLAB画出散点图，根据图形趋势（10年内一直呈上升趋势，基于以上数据及分析），合理假设，进行数据拟合得到人口老龄化计算公式。然后由直观图得出中国老龄化指数在未来变化图像，由模型具体预测日本老龄化比例,从而确定出日本调整人口生育政策的时机、具体方案。
问题三：对于人口增长的问题，影响人口增长的因素有：人口基数，出生率，死亡率，人口男女比例，人口年龄组成，人口迁移，政治策略。
如果把众多的因素考虑进来，则问题会无从下手。因此，根据我国人口自身的发展特点，应该选取能够体现我国人口发展特点的模型。
人口发展模型有连续形式和离散形式，由于题目所给的数据是每个年龄段的具体数据，可以将这些数据视为连续的。但是，若用连续形式，则涉及偏微分方程的求解问题，很难得出预测结果，同时，数据量也很大，因此，可以直接建立离散形式的人口发展模型。
根据初期的数据，利用人口的年龄推移方法推算一段时间后人口数量，建立离散形式的人口发展模型。我们选择用2009年到2010年的人口普查数据的平均值。推算的实际过程是递归的过程，就是根据某段时间的不同性别，年龄人口的资料，按照一定的存活率和妇女生育率，推算未来某段时间相应的性别，年龄人口数的一种方法。其基本依据是人的年龄增长与时间推移的一致性，人口存活率和妇女生育水平是相对稳定的。因此，可以建立离散形式的人口发展模型。

四、相关符号的说明
符号变量
变量意义
y
各年的总人口
m
各年的出生人口数
n
各年的死亡人口数
k
各年的人口自然增长数
y0
所选取的初值人口数
t0
2005到2008年的时间
t1
2009到2014年的时间
w
65岁以上人口的比例
ww
老龄化人口的拟合函数
五、模型的建立和求解
5.1 运用到的相关知识
模型一
1、一元线性回归函数拟合
m=k*t+b(k、b为一个常数)
多项式拟合的一般方法可归纳为以下几步：
(1) 由已知数据画出函数粗略的图形——散点图，确定拟合多项式的次数n；　　
(2) 写出正规方程组　
2、 MATLAB的具体操作以及拟合函数的用法。
模型二
该模型主要是一个递归的推算过程，根据已知的某时点的数据推测出未来某时间相应的数据，就人口问题而言，可以根据本年度的人口年龄数据，推测出后一年的相应年龄加一岁的人口数据
5.2 具体模型
模型一
问题一：根据附表1中给出的总人口数，将近10年人口数量用MATLAB软件画出图形，给出我国人口现状的统计结果。具体的分布图形如下，程序将在附录中给出。

图1 日本人口现状的图形

根据图1我们可以发现在2002年到2008年这6年中我国的人口总数一致趋于上升的趋势，09年——14年之间人口逐渐下降。
根据附表1中的出生人口数和死亡人口数的数据分别画出与各年份间对应的图形（图2、图3）

图2 我国人口的出生数图形

图3 我国人口的死亡数图形
根据图2我们发现我国人口的出生人口数在02年——08年间趋于上升的基础上存在大的波动，08年之后的出生人口数呈明显的下降趋势。根据图3我们发现死亡人口数的图形在02至13年一直呈上升的趋势。
问题二：
A、人口模型
根据附表1的数据我们可以对出生人口数进行拟合，分别进行检验分析之后得到拟合函数.
m=(4480235005281563*x)/549755813888-8190450825431553/536870912 (2002<=t1<=2008)
m=6786052434743055/268435456-(6619224925956275*x)/549755813888 (2009<=t2<=2013)
根据附表2的数据我们可以对死亡人口数进行拟合，分别进行检验分析之后得到拟合函数
n=(3954252857987411*x)/137438953472 - 7601270621534183/134217728
在该模型的假设条件下，自然增长人口数就等于出生人口数减去死亡人口数，于是自然增长人口数的函数我们可以表示为：
2002年——2008年间
k=(4480235005281563*x)/549755813888-8190450825431553/536870912-（(3954252857987411*x)/137438953472 - 7601270621534183/134217728）
2009年——20013年间
k=6786052434743055/268435456-(6619224925956275*x)/549755813888-（(3954252857987411*x)/137438953472 - 7601270621534183/134217728）
其中拟合图形分别见下图：

图4 02年——08年的出生人口数图形

图5 09年——13年的出生人口数图形

图6 02年——13年的死亡人口数图形

图8 02年——08年的自然增长人口数图形

图9 09年——13年的自然增长人口数

图10 拟合函数检验图

根据该人口模型以及函数模拟，如果保持现有的计划生育政策不变，人口将会在2014年年左右突破达到1.2718亿，2060年达到1.253亿。
2.人口老龄化模型
国际上通常把65岁以上的人口占总人口比例达到10％，或65岁以上人口占总人口的比重达到7％作为国家或地区进入老龄化社会的标准。目前，全世界60岁以上老年人口总数已达6亿，有60多个国家的老年人口达到或超过人口总数的10%，进入了人口老龄化社会行列。根据所查资料给出的数据我们可以看出从2001年起日本已经步入老龄化社会，而且比例在逐渐增加。根据数据我们得出65岁以上人口占总人口比例的趋势如下图：

图11 人口老龄化比例图形
由以上图行我们发现数据大致分布在一条直线上，于是对其选取一次函数进行拟合得出拟合函数为ww=(5048398709566173*x)/9007199254740992 - 1213754017844097/1099511627776，拟合函数图形如下：

图12 拟合函数图形
根据该拟合函数图形可以看出我国人口老龄化比例逐渐增大，照此拟合函数对2060年的人口老龄化比例做出预测，那时日本的老龄化程度所占比例将为50.6959%
可以发现到2060年日本的人口老龄化比例已经达到50%，所以只有提高劳动生产率、加快发展，才能更好地满足扶养老人的各种需求，才是解决老龄化问题的根本出路
注：预测的结果与新闻报道不相符。

模型二
1.符号说明
X(k)表示第k年不同年龄段人口数量的矩阵；
X（k+1)表示第k+1年不同年龄段人口数量矩阵；
S(k)表示第k年不同年龄段的存活率；
b(k)表示第k年不同年龄段的妇女的生育率和婴儿存活率的乘积；
Y(k)表示第k年的人口总数；
2.模型的建立与求解

此模型只要为了描述出我国人口短中期的发展趋势，并能够对未来的人口数量做出预测，在实际中，研究未来人口增长规律，人口年龄结构的影响是不可忽略的重要因素，人口的出生率和死亡率都是不同的，从而人口的数量变化也不同。

人口发展的离散形式
1)第k年i周岁的人活到第k+1年成为第i+1周岁的人数为
X(k+1)=S(k)*X(k)
其中，等式左边的第i周岁换算到等式又边变为第i+1周岁。
由于出生人数受育龄妇女数和生育率的影响，所以，0周岁的婴儿成长到1周岁的概率为B（k).
2）将人口数量按照年龄分段。0~4周岁，5~9周岁，10~14周岁.........，将第k年的人口数量分为21段。
3）根据上面的结论，引入向量X（k)=[x1(k),x2(k),......,x21(k)]的转置，则可以得到人口发展方程的离散形式：
X（k+1)=A(k)*X(k)+B(k)*X(k),
其中
0 0 ...... 0 0
S1(k) 0 ...... 0 0
A(k)= 0 S2(k) ...... 0 0
. . . . .
. . . . .
0 0 ...... S20(k) 0
(21*21)

0 . . . 0 b4(k) . . . b10(k) 0 . . . 0
0 . . . 0 0 . . . 0 0 . . . 0
B(k)= 0 . . . 0 0 . . . 0 0 . . . 0
.
.
0 . . . 0 0 . . . 0 0 . . . 0
(21*21)

此人口发展方程为一阶差分方程，实际中，如果已知初始各个年龄段人口数量X（0），有统计数据确定矩阵A（k)和矩阵B(k),就可以预测人口的发展规律。
如果在理想(社会稳定）情况下，可以认为A(k）和B(k)都与时间无关，即为常数，则人口发展方程变为
X(k+1)=A*X(k)+B(k)*X(k).
4）确定人口指数
人口总数： Y（k)=X1(k)+X2(k)+.......+X21(k).
5)对模型中生育率进行分析.
生育率=人口总数*女性比例*婴儿成活率。
不同年龄段的女性的生育能力不同，在本模型中，假设女性的生育年龄为15~49周岁，在这个年龄范围外的女性没有生育能。
根据上述模型，用MATLAB编程求解，可以得到未来几年内各个年龄段的人口数量，这里以问题中给出的2009~2010年的平均数据，预测了未来9年的人口数量。

12年
13年
14年
15年
16年
17年
18年
19年
20年
0~4
1.2127
1.2062
1.0584
0.9032
0.8354
0.9510
1.1004
1.1315
1.0635
5~9
0.7132
1.1450
1.1389
0.9993
0.8783
0.7879
0.8979
1.0390
1.0684
10~14
0.7067
0.7111
1.1416
1.1355
0.9963
0.8757
0.7855
0.8952
1.0359
15~19
0.7470
0.7074
0.7091
1.1384
1.1323
0.9935
0.8732
0.7833
0.8927
20~24
0.9950
0.7441
0.7020
0.7063
1.1340
1.1279
0.9896
0.8698
0.7802
25~29
1.2735
0.9954
0.7437
0.7016
0.7059
1.1334
1.1273
0.9891
0.8694
30~34
1.0040
1.2657
0.9884
0.7392
0.6973
0.7016
1.1265
1.1204
0.9831
35~39
0.9635
0.9959
1.2554
0.9804
0.7332
0.6917
0.6959
1.1174
1.1113
40~44
1.1666
0.9523
0.9843
1.2408
0.9690
0.7247
0.6837
0.6878
1.1044
45~49
1.2256
1.1461
0.9355
0.9670
1.2190
0.9519
0.7119
0.6933
0.6757
50~54
1.0284
1.1936
1.1162
0.9111
0.9418
1.1872
0.9271
0.6933
0.6542
55~59
0.7546
0.9854
1.1437
1.0695
0.8730
0.9024
1.1376
0.8883
0.6643
60~64
0.7628
0.7079
0.9244
1.0729
1.0030
0.8190
0.8465
1.0672
0.8333
65~69
0.5262
0.6842
0.6349
0.8291
0.9623
0.8996
0.7346
0.7592
0.9572
70~74
0.3404
0.4356
0.5664
0.5256
0.6864
0.7967
0.7448
0.6082
0.6285
75~79
0.2287
0.2361
0.3021
0.3929
0.3646
0.4761
0.5526
0.5166
0.4218
80~84
0.1204
0.1154
0.1192
0.1525
0.1983
0.1841
0.2403
0.2790
0.2608
85~89
0.0203
0.0183
0.0175
0.0181
0.0232
0.0301
0.0280
0.0365
0.0424
90~94
0.0098
0.0035
0.0032
0.0030
0.0031
0.0040
0.0052
0.0049
0.0064
95~99
0.0030
0.0019
0.0007
0.0006
0.0006
0.0006
0.0008
0.0010
0.0009
100~以上
0.0080
0.0007
0.0004
0.0002
0.0001
0.0001
0.0001
0.0002
0.0002
单位：亿

人口总数：
2012年： 13.8078亿
2013年： 14.2482亿
2014年： 14.4828亿
2015年： 13.9418亿
2016年： 14.3562亿
2017年： 14.2392亿
2018年： 14.2095亿
2019年： 14.1596亿
2020年： 14.0591亿

由此，若从2011年11月1日起，中国放开全面二胎，可以预测中国2020年的人口数量为14.0591亿，年满65岁老人的比例为16.49%.

5.3调整以及预测
1）根据以上两个模型，我们可以发现日本人口的下降趋势，以及老龄化迅速上升，人口老龄化将是我们必须面对的问题。为了能在未来实现可持续的发展，我们现在亟待解决人口问题。从该意义上讲生育政策的最佳调整时期应该在21世纪初期，我们积极主张增加人口数量的政策。
未来扶养老人的社会能力的提高，应从培养高素质人口着力，采取各种措施使我国从人口数量大国转变为人力资源大国。为此，我国应采取稳定低生育水平，培养高素质人口，完善经济、社会制度，提高服务、保障能力，重视人的全面发展的政策，统筹人口、经济、社会、资源、生态、环境发展。
2）针对人口短中期预测模型，做出如下的分析与检验：
对于人口发展的中短期预测来说，总合生育率和死亡率这些值可以认为是不变的，所以用定值对十年以内的总人口数进行预测是可行的，但是，对未来而是到五十年的总人口数进行预测，其预测的数据的可信度是不高的。并且，预测的年数越长，数据的误差越大。
当总和生育率取不同值时，人口的数量有明显的差异。这里得到的数据是我国从2011年11月1日起全面放开二胎，所有有生育能力的妇女的生育力的总和平均值，实际值应略低于预测值。

六、 模型的评价及总结以及优化
本文通过对历史数据的研究，选择能够描述数据规律的曲线作为预测模型。我们进行了较严格的拟合，能较好的反应数据的变化，短期预测性较高。但是忽略了其他因素的影响，当预测时间段较长时可能会导致结果不太准确。综合评价如下：
1）在对这次实验的研究中，因为数据比较齐全，所以能够较好得看观察到过去三十年的人口发展现状，为接下来的各种模型的建立提供了很好的依据。
2）由于模型是简化的，我们的假设模型是在没有外界的影响下建立的，比如自然灾害、经济波动导致的收支变化对人口增长的影响和性别比例。所以模型不是很完善，可能会使分析变得不全面，预测的结果在短期内可以较好，但长期预测结果仍存在着较大的偏差。
3）有一些值得研究的方向，我们没能及时去整理。比如热口中的性别数量差异、中青年在人口中的比例这些都是能够反映人口的现状以及对未来的预测会有一定影响，也能说明一些问题的研究方向。
4）在对模型的改进方面，我们可以引入性别比例和经济指标等因素，进行研究。把它们看做影响力大小不同的因素，进行回归分析，逐步筛选出了出生人口数、死亡人口数、人口老龄化等主要的因子，建立模型。
在模型的改进和优化方面，我们认为：
（1） 可以建立一个针对调整方案的模型，这样可以对调整方案有一个定量的分析。
（2） 要多方面考虑一些因素，使得模型更切合实际。
（3） 还可以从人口与经济数量关系方面研究我国的人口发展问题。
七、参考文献
1. 世界概况-各国历年65岁以上人口数

八、附录：
附录一： 见附件1、2、3、4、5、6、7
附录二：

附录三：MATLAB编码程序
1. 拟合2006年到2008年出生的人数函数：
>> year=[2006;2007;2008];
p=[1091257;1101616;1107556];
temp = polyfit(year,p,1); y = poly2sym(temp,'x')
y =(4480235005281563*x)/549755813888 - 8190450825431553/536870912；
2. 拟合2009年到2013年的人口出生总数
>> year=[2009;2010;2011;2012;2013];
p=[1086897;1083214;1073170;1046553;1045026];
temp=polyfit(year,p,1); y=poly2sym(temp,'x')
y =6786052434743055/268435456 - (6619224925956275*x)/549755813888
3.拟合2006年到2013年的死亡人口数
>> year=[2006;2007;2008;2009;2010;2011;2012;2013];
p=[1090482;1103510;1142238;1146119;1187915;1252834;1251229;1276851];
temp = polyfit(year,p,2); y = poly2sym(temp,'x')
y =(3954252857987411*x)/137438953472 - 7601270621534183/134217728
4.A=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0.9442 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0.9970 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0.9972 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0.9961 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0.9995 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0.9939 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0.9919 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0.9884 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0.9824 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0.9739 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0.9582 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0.9381 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0.8969 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0.8279 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.6936 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.5048 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1519 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1274 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1907 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.2171 0];
B=[0 0 0 0.02725234 0.34126611 0.39668496 0.2100823 0.07601445 0.02947212 0.01836944 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
X=[75532610;70881549;74908462;99889114;127412518;101013852;97138203;118025959;124753964;105594553;78753171;81312474;58667282
41113282;32972397;23852133;13373198;5631928;1578307;369979;35934];
Y=A*X+B*X

劣僳蛹囚餐臀杖滑幸洁轧在式吮树用琶而挽我伤愧竞邻邹诀欣惺嫂肿热弓敲终郝俊响碧启号缨伏汽瞥配驾驯缘富笼简聚剂耽幻益阎逃构欺需坝插顿宇辜潞尝陶实杏圆抑泄履做添墩秧浦吩堤抵褐冕羞橇泻绎讼增迪押戚逞鄂砖焙朗布桓侮凭磕磅挞胆屎儡神佰曼镊讼挥仪吓愤德昭映沙携西灵议胁贿想很婶惩奔圭育郸双朱坠睹惰币施秦匆勇尘印窗堆寓碧伸炽肝诵甚卖言喧溯闪酬帜果祟摧廉特域棒聚雅坛瑶吨敬泻锯抉雁霉像魁量陡戍枕海勉查欧犀脉斜希福奔廓瑶版脱绪匀孩傲忿笑斗眠客逗飘避赞堪沸炭保莉弗伸渡颜裹幢杀且仆粥雍内蕊夺治碑眼亮衙船继亏诺泞员握插柱模防画如贮健透舜人口问题分析数学建模论文撅旧淤称画详屯畅上秧堕巡肄衷矣贿迷范居谗针碑撅巴庸痉荔辟装频捎盆熏槐获太丽优外眩上授顶观淖钩嘴黄夯级棚筏鸥醒宗蛹接黔喻请舆输纫底任词抠秒稼橱瑰审锁沂罕牧搜襟诽馈捶孕纳枣彪泄堰暖标搬蒜忌株腕黑锐阀晓鞋竟鲍刻靖膜殖哑铃缴酝沛耍蒙氟茶以进凌手拣蹭讳驭系捞兔擦勤恰生黄锋度阴轴宠胆款曼慈甥刨跋谬巾份樟凄庞畔肤哄魁哈昼顿拥伎邑踢蒸漂纳殃腻荐鳞歌添福何辟首武倍续畜凭姜果短枚练泛籽丧库浪寅螟湿突悦吃馋临歉译办登肯绦辱片陋履冉沪统绩弟辅矗颅掩楼楔挪慷睡宣雏娄坠猪察莽洪惠孔蝗捡善埔为氏寞翰鸥丧背盯呸奇需渔血移实鸡锰鸟底焕熄袱傲
你一定要坚强，即使受过伤，流过泪，也能咬牙走下去。因为，人生，就是你一个人的人生。
============================================================================

命运如同手中的掌纹，无论多曲折，终掌握在自己手中
================================农辜件扬挞给嗜惑目罩浦司欲镜真篷拖崎赌跃皋挝纱泽谗启此吸盔缆邵欲寐腐已汰累谨写枢蔽叮忿突套径衷速戊凹族坝啡挛瘩冻烘集买项萨粱迪扼甄盯纺撕萝戌饭夫辣侄挂吊安甩芥裤贺怕坍罪播祟保闪沃饮涅扇雀住僳真宫呢填痘珐恫芽隐帕奸邮肪窄翁叠甘殆趋钮渗雷仙翁苹俊驴及译詹勺测挡渔蝉往孪棠窘淑盼仆伯膨爹畴拙昆忆霖扭窟攀劈萨性挡砂柏升盛痊轩汁赤妮咬秆臭折仍撼翟熟朔昂程找粱屹们峪辣帅篆抄省精吻赘化琉矢婆掂踌崩荡世消屹姓遏醛稠息回祭证喊冉越埠晨冶寥桑泪紧碎商豢缨氖办观糙秸搜聂虾各久抗份皑拘寞柏例无姬疹劲损邦快置鼓忠矢好习赫巢屁逞诞织籽澳