中国人口预测数学建模论文(4)

附录二：MATLAB程序

1、78年——89年的出生率程序 t1=1978:1:1989

y1=[0.01825 0.01782 0.01821 0.02091 0.02228 0.02019 0.0199 0.02104 0.02243 0.02333 0.02237 0.02158] [p,s]=polyfit(t1,y1,2) t2=1989:1:2008

[y delta]=polyval(p,t1,s)

yy=-0.000039777722*t1.^2+0.158207350160*t1-157.286723589417 plot(t1,yy)

2、89年——08年的出生率程序 t2=1989:1:2008

[p,s]=polyfit(t2,y2,2) [y delta]=polyval(p,t2,s)

yy=0.00001948279790*t2.^2-0.07839101390114*t2+78.86500481137304 plot(t2,yy)

3、78年——89年的死亡率程序

z1=[0.00625 0.00621 0.00634 0.00636 0.0066 0.0069 0.00682 0.00678 0.00686 0.00672 0.00664 0.00654]

[p,s]=polyfit(t1,z1,2) [y delta]=polyval(p,t1,s)

zz=-0.00001378621379*t1.^2+0.05473424575582*t1-54.31990883773020 plot(t1,zz)

4、89年——08年的死亡率程序

z2=[0.00654 0.00667 0.0067 0.00664 0.00664 0.00649 0.00657 0.00656 0.00651 0.0065 0.00646 0.00645 0.00643 0.00641 0.0064 0.00642 0.00651 0.00681 0.00693 0.00706] t2=1989:1:2008

[p,s]=polyfit(t2,z2,2) [y delta]=polyval(p,t2,s)

zz=0.00000428229665*t2.^2-0.01710954272160*t2+17.09636128346106 plot(t2,zz)

5、78年——89年的自然增长率程序

r1=-0.000039777722*t1.^2+0.158207350160*t1-157.286723589417-(-0.00001378621379*t1.^2+0.05473424575582*t1-54.31990883773020) plot(t1,r1)

6、89年——08年的自然增长率程序

r2=0.00001948279790*t2.^2-0.07839101390114*t2+78.86500481137304-(0.00000428229665*t2.^2-0.01710954272160*t2+17.09636128346106)

15

plot(t2,r2)

7、对我国人口数量自2006年后的近15年预测程序 t3=2005:1:2020 t4=2006:1:2021

r2=0.00001948279790*t3.^2-0.07839101390114*t3+78.86500481137304-(0.00000428229665*t3.^2-0.01710954272160*t3+17.09636128346106) yy1=130756*(r2+1).^(t4-2005) plot(t4,yy1)

8、人口老龄化比例的拟合程序 [p,s]=polyfit(t3,w,1) [y delta]=polyval(p,t3,s)

ww=0.154835164835*t3-302.673296703281 plot(t3,w,'bo',t5,ww,'r*')

16