一、简介
计划生育政策调整对人口数量、结构及其影响的研究
二、源代码
function pre_popilation()
clear all
close all
clc
%根据人口、年龄、性别、生育策略等估计未来人口数模型。
r1=15; %生育起始年龄
r2=45; %生育终止年龄
n=2; %晚育年龄
st=2; %生育模型参数1
rf=n/2; %生育模型参数2
r=r1:r2; %生育年龄
h=(r-r1).^(rf-1).*exp(-(r-r1)/st)/st^rf/gamma(rf); %生育模式
bt=2.5 %设妇女生育数为常数，即妇女在生育年龄生育概率为常数，且生育总数为常数
k=0.5; %设女性比例为常数0.5
p_sz=1249200 %2000年深圳人口数
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fil='人口年龄分布.xls';
p0=xlsread(fil,'b2:b102');
p0=p0/sum(p0)*p_sz;
fil='各年龄死亡人口数.xls';
u=xlsread(fil,'e3:e103');
pn=p0; %为第n年的人口分布设置初始分布
px=pn; %为第n+1年估计分配初值
year=input('请输入需要预测的年份(>2000)');
year=year-2000;
for num=1:year
t=pn(r1:r2);
f=bt*sum(k*h.*t');%一年期新生儿数量。
px(1)=f;
for n=1:length(pn)-1
px(n+1)=pn(n)*(1-u(n)*0.001);
end
plot(num+2000,sum(px),'.');
hold on
pn=px;
hold on
end
P_sum=sum(pn)
P0_sum=sum(p0)
figure
plot(pn);
hold on
plot(p0,'r--');
hold off
grid off
title('各年龄段的人口分布图');
xlabel('年龄/岁'),ylabel('人口数/人');
%预测未来人口数量及分布的函数，可直接运行，后面几个名字类似的是相关的修改程序
clear all
close all
clc
%根据人口、年龄、性别、生育策略等估计未来人口数模型。
r1=15; %生育起始年龄
r2=45; %生育终止年龄
n=2; %晚育年龄
st=2; %生育模型参数1
rf=n/2; %生育模型参数2
r=r1:r2; %生育年龄
h=(r-r1).^(rf-1).*exp(-(r-r1)/st)/st^rf/gamma(rf); %生育模式
bt=1.22 %设妇女生育数为常数，即妇女在生育年龄生育概率为常数，且生育总数为常数
k=0.5; %设女性比例为常数0.5
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fil='人口年龄分布.xls';
p0=xlsread(fil,'b2:b102');
fil='各年龄死亡人口数.xls';
u=xlsread(fil,'e3:e103');
pn=p0; %为第n年的人口分布设置初始分布
px=pn; %为第n+1年估计分配初值
year=input('请输入需要预测的年份(>2000)');
year=year-2000;
for num=1:year
t=pn(r1:r2);
f=bt*sum(k*h.*t');%一年期新生儿数量。
px(1)=f;
for n=1:length(pn)-1
px(n+1)=pn(n)*(1-u(n)*0.001);
end
plot(num+2000,sum(px),'.');
hold on
pn=px;
hold on
end
三、运行结果
四、备注
版本：2014a
