%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
axis=
[1999.00,2000.00,2001.00,2002.00,2003.00,2004.00,2005.00,2006.00,2007.00,2008.00,200
9.00,2010.00];%年份
a=
[3422.00,3565.00,3866.00,4134.00,5118.00,5855.00,6842.00,7196.00,10320.00,13411,1580
0,21600];%房价
a1=leija(a,12);%得到X1(k)
f=[0];
bf=cat(2,f,a1);
bf2=cat(2,a1,f);
bf3=ones(11,1);
B=-0.5*(bf+bf2);
B=B(2:12)';
B=cat(2,B,bf3);%得到矩阵B
Y=a(2:12)';%得到矩阵Y
A=(inv(B'*B))*B'*Y;%得到矩阵A
anum=A(1)
bnum=A(2)%两个系数
a1yuce=zeros(1,11);
for i=0:11
a1yuce(i+1)=(a(1)-bnum/anum)*exp(-anum*i)+bnum/anum;
end %生成X1(k)预测值
ayuce=leijian(a1yuce,12);%生成X(k)预测值