2018年中青杯数学建模

2018年中青杯研究生组数学建模B题：开放生育限制政策，包括源代码与相关材料

clear
load（‘data.mat‘）;
N = data（:3）;
N0=N/100000000; %第 0 年（2015 年）女性各个年龄（0-100 岁）的人口数（亿）
A=eye（100）;
b = 1-data（:2）‘/1000;
%存活率（1-死亡率）
for i=1:100
A（i:）=A（i:）*b（1i）;
end
A;
c = data（1:end-15）‘/1000; %生育率
c1=0.8*c;p=0.48; %女性占总人口的比例
M=sum（c1‘）; %总合生育率
d=zeros（1011）; B=[c1;A]; L=[Bd]; %构造的 lestie 矩阵
for i=0:1:35
X=L^i*N0; %第 i 年后女性各个年龄段的人口数（亿）
Z（:i+1）=X./p; %第 i 年在各个年龄段的人口总数预测
W（i+11）=sum（Z（:i+1））;
end
W; %2015-2050 的总人数
x = 2015:1:2050;
plot（xW‘-r‘）%成图
xlabel（‘年份‘）;ylabel（‘人口数（亿）‘）;
title（‘2015-2050 年 leslie 人口预测‘）;
set（gca‘ytick‘[13.9:0.15:15.1]）
set（gca‘yticklabel‘{‘13.9‘‘14.0‘‘14.1‘‘14.2‘‘14.3‘‘14.4‘‘14.5‘‘14.6‘‘14.7‘}）
grid on;
for i = 1:1:36
young（i1） = sum（Z（19:61i））;
oldman（i1） = sum（Z（62:101i））;
end
figure

YY = young./W*100;
OO = oldman./W*100;

[AXH1H2] = plotyy（x[youngoldman]x[YYOO]）;
set（H1（1）‘linestyle‘‘-‘‘Marker‘‘x‘‘color‘‘k‘）
set（H1（2）‘linestyle‘‘-‘‘Marker‘‘x‘‘color‘‘r‘）
set（H2（1）‘linestyle‘‘-‘‘Marker‘‘s‘‘color‘‘b‘）
set（H2（2）‘linestyle‘‘-‘‘Marker‘‘s‘‘color‘‘m‘）
set（AX（1）‘ylim‘[210]‘ytick‘[2:2:10]）%设置左y轴数据范围  
set（AX（2）‘ylim‘[080]‘ytick‘[0:10:80]）%设置右y轴数据范围  
set（get（AX（1）‘Ylabel‘）‘String‘‘人口数（亿）‘）
set（get（AX（2）‘Ylabel‘）‘String‘‘人口比例（%）‘）
set（get（AX（1）‘Xlabel‘）‘String‘‘年份‘）
legend（‘劳动人口数量‘‘老年人口数量‘‘劳动人口占比‘‘老年人口占比‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件      351486  2018-06-10 22:55  A3318.doc
     目录           0  2018-06-10 23:00  A3318材料\
     目录           0  2018-06-10 23:00  A3318材料\中国人口和就业统计年鉴2016\
     文件       36864  2017-02-21 09:34  A3318材料\中国人口和就业统计年鉴2016\0101.xls
     文件       33792  2017-02-21 09:34  A3318材料\中国人口和就业统计年鉴2016\0102.xls
     文件       37376  2017-02-21 09:32  A3318材料\中国人口和就业统计年鉴2016\0103.xls
     文件       38912  2017-02-21 09:50  A3318材料\中国人口和就业统计年鉴2016\0104.xls
     文件       43008  2017-02-21 09:51  A3318材料\中国人口和就业统计年鉴2016\0105.xls
     文件       34816  2017-02-21 09:32  A3318材料\中国人口和就业统计年鉴2016\0106.xls
     文件       54272  2017-02-21 09:31  A3318材料\中国人口和就业统计年鉴2016\0107.xls
     文件       26624  2017-02-21 09:31  A3318材料\中国人口和就业统计年鉴2016\0108.xls
     文件       31744  2017-02-21 09:30  A3318材料\中国人口和就业统计年鉴2016\0109.xls
     文件       41984  2017-02-21 09:30  A3318材料\中国人口和就业统计年鉴2016\0110.xls
     文件       35840  2017-02-21 09:29  A3318材料\中国人口和就业统计年鉴2016\0111.xls
     文件       36864  2017-02-21 09:29  A3318材料\中国人口和就业统计年鉴2016\0112.xls
     文件       38912  2017-02-21 09:28  A3318材料\中国人口和就业统计年鉴2016\0113.xls
     文件       30208  2017-02-21 09:28  A3318材料\中国人口和就业统计年鉴2016\0114.xls
     文件       30720  2017-02-21 09:28  A3318材料\中国人口和就业统计年鉴2016\0115.xls
     文件       36352  2017-02-21 09:27  A3318材料\中国人口和就业统计年鉴2016\0116.xls
     文件       43008  2017-02-21 09:27  A3318材料\中国人口和就业统计年鉴2016\0117.xls
     文件       40960  2017-02-21 09:26  A3318材料\中国人口和就业统计年鉴2016\0118.xls
     文件       30720  2017-02-21 09:26  A3318材料\中国人口和就业统计年鉴2016\0119.xls
     文件       33792  2017-02-21 09:26  A3318材料\中国人口和就业统计年鉴2016\0120.xls
     文件       45568  2017-02-21 09:25  A3318材料\中国人口和就业统计年鉴2016\0121.xls
     文件       45056  2017-02-21 09:25  A3318材料\中国人口和就业统计年鉴2016\0122.xls
     文件       40960  2017-02-21 09:24  A3318材料\中国人口和就业统计年鉴2016\0123.xls
     文件       40960  2017-02-21 09:23  A3318材料\中国人口和就业统计年鉴2016\0124.xls
     文件       28160  2017-02-21 09:23  A3318材料\中国人口和就业统计年鉴2016\0125a.xls
     文件       38400  2017-02-21 09:23  A3318材料\中国人口和就业统计年鉴2016\0125b.xls
     文件       33792  2017-02-21 09:22  A3318材料\中国人口和就业统计年鉴2016\0126.xls
     文件       30720  2017-02-21 09:22  A3318材料\中国人口和就业统计年鉴2016\0127.xls
............此处省略199个文件信息