首页 > 编程语言> > 【数学建模】计划生育政策调整对人口数量、结构及其影响的研究matlab源码

【数学建模】计划生育政策调整对人口数量、结构及其影响的研究matlab源码

2021-04-11 10:01:30 作者：互联网

一、简介

计划生育政策调整对人口数量、结构及其影响的研究

二、源代码

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