matlab 灰色预测 模型（1,1）

把数据改一下应该就行了
clc
clear
x0=[1827 1729.3 1431 1231.7 1109.9 768.2 792.9 805.8];
%道路交通事故财产损失原始数据
for i=2:8
x1(1)=x0(1);
x1(i)=x1(i-1)+x0(i);
end
x1
%对原始数据进行累加
for i=1:7
yn(i)=x0(i+1);
end
yn
for i=1:7
b(i)=(-0.5)*(x1(i)+x1(i+1));
end
b
for i=1:7
B(i,1)=b(i);
B(i,2)=1;
end
B
%求得数据矩阵B
c=inv(B'*B)*B'*yn'
a=c(1,1)
u=c(2,1)
%通过最小二乘法求出变量a,u
for t=1:8
x11(t)=u/a+exp(-a*t)*(-u+1827*a)/exp(-a)/a;
end
x11
%建立时间响应函数
q0=x1-x11
for i=2:8
q1(1)=q0(1);
q1(i)=q1(i-1)+q0(i);
end
q1
%残差第一次累加
for i=2:8
q2(1)=q1(1);
q2(i)=q2(i-1)+q1(i);
end
q2
%残差第二次累加
for i=1:7
b1(i)=(-0.5)*(q2(i)+q2(i+1));
end
b1
for i=1:7
B1(i,1)=b1(i);
B1(i,2)=1;
end
B1
c1=inv(B1'*B1)*B1'*yn'
a1=c1(1,1)
u1=c1(2,1)
syms t
q2t=u1/a1+exp(-a1*t)*(-u1+q2(1)*a1)/exp(-a1)/a1;
D=diff(q2t)
%对二次残差求灰导
X=u/a+exp(-a*t)*(-u+1827*a)/exp(-a)/a+D;
C=1:5;
y=subs(X,t,C)
z2012=y(3)-y(2)
z2013=y(4)-y(3)
z2014=y(5)-y(4)
%求得2012-2014年道路交通事故财产损失预测值