南京财经大学本科毕业论文（设计）
for j=1:n-1 %计算数据矩阵B的第二列数据 e(j,:)=1; end
for i=1:n-1 %构造数据矩阵B B(i,1)=c(i,:); B(i,2)=e(i,:); end
alpha=inv(B'*B)*B'*y; %计算参数?,u矩阵
for i=1:n+1 %计算数据估计值的累加数列，如改为n+1为n+m可预测后m-1个值
ago(i,:)=(x(1,:)-alpha(2,:)/alpha(1,:))*exp(-alpha(1,:)*(i-1))+alpha(2,:)/alpha(1,:); end
var(1,:)=ago(1,:)
for i=1:n %如改n为n+m-1，可预测后m-1个值
var(i+1,:)=ago(i+1,:)-ago(i,:); %估计值的累加数列的还原，并计算出下一预测值 end for i=1:n
error(i,:)=var(i,:)-x(i,:); %计算残差 end
c=std(error)/std(x); %调用统计工具箱的标准差函数计算后验差的比值c ago %显示输出预测值的累加数列 alpha %显示输出参数?,u数列 var %显示输出预测值 error %显示输出误差
20
南京财经大学本科毕业论文（设计）
c %显示后验差的比值c
2．中国人口的基本数据（1950-2006） 总人口出生率死亡率自然增长年份 (万) (‰) (‰) 率(‰) 1950 55196 37 18 19 1951 56300 37.8 17.8 20 1952 57482 37 17 20 1953 58796 37 14 23 1954 60266 37.97 13.18 24.79 1955 61456 32.6 12.28 20.32 1956 62828 31.9 11.4 20.5 1957 64653 34.03 10.8 23.23 1958 65994 29.22 11.98 17.24 1959 67207 24.78 14.59 10.19 1960 66207 20.86 25.43 -4.57 1961 65859 18.02 14.24 3.78 1962 67295 37.01 10.02 26.99 1963 69172 43.37 10.04 33.33 1964 70499 39.14 11.5 27.64 1965 72538 37.88 9.5 28.38 1966 74542 35.05 8.83 26.22 1967 76368 33.96 8.43 25.53 1968 78534 35.59 8.21 27.38 1969 80671 34.11 8.03 26.08 1970 82992 33.43 7.6 25.83 1971 85229 30.65 7.32 23.33 1972 87177 29.77 7.61 22.16 1973 89211 27.93 7.04 20.89 1974 90859 24.82 7.34 17.48 1975 92420 23.01 7.32 15.69 1976 93717 19.91 7.25 12.66 1977 94974 18.93 6.87 12.06 1978 96259 18.25 6.25 12 1979 97542 17.82 6.21 11.61 1980 98705 18.21 6.34 11.87 1981 100072 20.91 6.36 14.55 1982 101654 22.28 6.6 15.68 1983 103008 20.19 6.9 13.29 1984 104357 19.9 6.82 13.08
21
净增人口(万) 1029 1104 1182 1314 1470 1190 1372 1825 1341 1213 -1000 -348 1436 1877 1327 2039 2004 1826 2166 2137 2321 2237 1948 2034 1648 1561 1297 1257 1285 1283 1163 1367 1582 1354 1349 总人口增长率(‰) 19 20 20.99 22.86 25 19.75 22.32 29.05 20.74 18.38 -14.88 -5.26 21.8 27.89 19.18 28.92 27.63 24.5 28.36 27.21 28.77 26.95 22.86 23.33 18.47 17.18 14.03 13.41 13.53 13.33 11.92 13.85 15.81 13.32 13.1 南京财经大学本科毕业论文（设计）
1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 105851 107507 109300 111026 112704 114333 115823 117171 118517 119850 121121 122389 123626 124761 125786 126743 127627 128453 129227 129988 130756 131448 21.04 22.43 23.33 22.37 21.58 21.06 19.68 18.24 18.09 17.7 17.12 16.98 16.57 15.64 14.64 14.03 13.38 12.86 12.41 12.29 12.4 12.09 6.78 6.86 6.72 6.64 6.54 6.67 6.7 6.64 6.64 6.49 6.57 6.56 6.51 6.5 6.46 6.45 6.43 6.41 6.4 6.42 6.51 6.81 14.26 15.57 16.61 15.73 15.04 14.39 12.98 11.6 11.45 11.21 10.55 10.42 10.06 9.14 8.18 7.58 6.95 6.45 6.01 5.87 5.89 5.28 1494 1656 1793 1726 1678 1629 1490 1348 1346 1333 1271 1268 1237 1135 1025 957 884 826 774 761 768 692 14.32 15.64 16.68 15.79 15.11 14.45 13.03 11.64 11.49 11.25 10.6 10.47 10.11 9.18 8.22 7.61 6.97 6.47 6.03 5.89 5.91 5.29 3．MATLAB程序求解长期人口模型（即通过该程序求的中国每年净增加的人口数目来预测其未来的人口数） for i=1:50
x(i)=-17159.2214*exp(-0.0524*i)+18116.2214; x(i+1)=-17159.2214*exp(-0.0524*(i+1))+18116.2214; n= x(i+1)-x(i) end
22