matlab 灰色GM（1，1）預測模型 預測房價

 哈哈。我又來了，這一次寫的依然是校數學建模裏面的比賽的代碼。灰色預測模型。具體代碼就是下面的了。。。沒有啥玩意了。

預測廣州十年內房價matlab程序：
clc,clear;
syms a b;
c=[a b]’;
A=[6795,7455.6,8677.6,11534.9,14605.5,15472.1,18163.9,18757.1,19711.5,21163.5];
B=cumsum(A); %原始數據累加
n=length(A);
for i=1:(n-1)
C(i)=(B(i)+B(i+1))/2; %生成累加矩陣
end
%計算待定參數的值
D=A;D(1)=[];
D=D’;
E=[-C;ones(1,n-1)];
c=inv(E*E’)*E*D;
c=c’;
a=c(1);b=c(2);
%預測後續數據
F=[];F(1)=A(1);
for i=2:(n+10) %只推測後10個數據，可以從此修改
F(i)=(A(1)-b/a)/exp(a*(i-1))+b/a;
end
G=[];G(1)=A(1);
for i=2:(n+10) %只推測後10個數據，可以從此修改
G(i)=F(i)-F(i-1); %得到預測出來的數據
end
t1=2007:2016;
t2=2007:2019; %多10組數據
G
h=plot(t1,A,’o’,t2,G,’-‘); %原始數據與預測數據的比較
set(h,’LineWidth’,1.5);