- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
人口统计预测程序
人口统计预测程序
余时伟
人口统计预测程序1
t=1991:2009;
y=[11.58,11.72,11.85,11.98,12.11,12.24,12.3626,12.4810,12.5909, 12.6743,...
12.7627,12.8453,12.9227,12.9988,13.0628,13.1448,13.2129,13.2802, 13.3474];
x=log(y);
plot(t,x);
A=[ones(length(t),1),t];
L=A\x;
Yearandpeople1=[t;y]
format short
exp(L(2,1)*t+L(1,1))
for t=2003:2010
year=2003:2010;
number=exp(L(2,1)*year+L(1,1)*ones(1,8));
end
Yearandpeople2=[year/100;number]
2007年人口预测13.3363 亿人,实际13.2129亿人
人口统计预测程序2
clear
format long
T=1991:2009;
N=[11.58,11.72,11.85,11.98,12.11,12.24,12.3626,12.4810,12.5909, 12.6743,...
12.7627,12.8453,12.9227,12.9988,13.0628,13.1448,13.2129,13.2802, 13.3474];
subplot(1,3,1),bar(T,N)
L=polyfit(T,N,4); LT=polyval(L,T)
subplot(1,3,2),plot(T,N,o,T,LT)
r1=sum((N-LT).^2)
out1=[N;LT]
L2=polyfit(T,log(N),3);
LT2=polyval(L2,T);
subplot(1,3,3),plot(T,N,o,T,exp(LT2))
r2=sum((N-LT).^2)
out2=[N;exp(LT2)]
y2010=exp(polyval(L2,2010))
Year=2010:2030;
y2030=exp(polyval(L2,Year));
[Year;y2030]
2007年人口预测值13.2159 亿人,实际值13.2129亿人
在不使用2008年人口数据情况下,
2008年人口预测值13.2583189亿人,实际值13.2802亿人,误差219万人;
在使用2008年人口数据情况下,
2008年人口预测值13.2714944亿人,实际值13.2802亿人,误差87万人.
2010年人口预测值13.399806236702204亿人,实际值13人,误差2.5577万人
中新网4月28日电 国家统计局今日发布第六次全国人口普查主要数据公报(第1号),公报显示,普查登记的大陆31个省、自治区、直辖市和现役军人的人口共1339724852亿人。
13.43932074308395
2030年预测值为14.5275亿人,
人口统计模型的进一步分析
%Two ways are introduced and compared in following programs .
clear
erro1=zeros(1,7);erro2=zeros(1,7);T=1991:2009;
format long
N=[11.58,11.72,11.85,11.98,12.11,12.24,12.3626,12.4810,12.5909, 12.6743,...
12.7627,12.8453,12.9227,12.9988,13.0628,13.1448,13.2129,13.2802, 13.3474];
subplot(1,3,1),bar(T,N)
for i=1:7
L=polyfit(T,N,i);
LT=polyval(L,T);
out1=[N;LT];
erro1(i)=sqrt(sum((out1(1,:)-out1(2,:)).^2));
L2=polyfit(T,log(N),i);
LT2=polyval(L2,T);
out2=[N;exp(LT2)];
erro2(i)=sqrt(sum((out2(1,:)-out2(2,:)).^2));
end
[min_val,number1]=min(erro1);
L=polyfit(T,N,number1);
LT=polyval(L,T);
subplot(1,3,2),plot(T,N,o,T,LT)
number1,out1=[T;N;LT]
[min_val,number2]
原创力文档


文档评论(0)