chengxu.txt

来自「我这儿有一组若干年份的数据」· 文本 代码 · 共 58 行

clc;clear 
A=[... 
1996 27000 
1997 30619 
1998 33282 
1999 33689 
2000 41020 
2001 43344 
2002 46030 
2003 53887 
2004 59271 
2005 64507 
2006 70597] 
B=[... 
071 20505 
072 19619 
073 19874 
074 19223 
081 19703] 
t=[A(:,1);2007] 
GNP=[A(:,2);sum(B(1:4,2))] 
pre=inline('a(1)*exp(a(2)*(t-1996))','a','t') 
format short g 
a=nlinfit(t,GNP,pre,[3.305e-81 0.09795]) 
T=1995:2015; 
PreGNP=a(1)*exp(a(2)*(T-1996)); 
temp=[NaN;GNP;NaN*ones(2015-2007,1)]' 
Table=[T;PreGNP;temp]' 
plot(t,GNP,'o',T,PreGNP) 
grid 

结果： 
a = 

26758 0.09795 
Table = 

1995 24261 NaN 
1996 26758 27000 
1997 29512 30619 
1998 32549 33282 
1999 35898 33689 
2000 39593 41020 
2001 43667 43344 
2002 48161 46030 
2003 53117 53887 
2004 58583 59271 
2005 64611 64507 
2006 71260 70597 
2007 78594 79221 
2008 86682 NaN 
2009 95602 NaN 
2010 1.0544e+005 NaN 
2011 1.1629e+005 NaN 
2012 1.2826e+005 NaN 
2013 1.4146e+005 NaN 
2014 1.5601e+005 NaN 
2015 1.7207e+005 NaN