灰色预测编程谁会？数据有30个，一共4组。C++ matlab 都行

syms a b

c=[a,b]'

A=[128 129 130 132]%这里换成你的数据

B=cumsum(A)

n=length(A)

for i=1:(n-1)

C(i)=(B(i)+B(i+1))/2

end

%计算待定参数的值

D=AD(1)=[]

D=D'

E=[-Cones(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)

F(i)=(A(1)-b/a)/exp(a*(i-1))+b/a

end

G=[]G(1)=A(1)

for i=2:(n+10)%n后面加几就是预测几个数据，这里是预测是10个数据

G(i)=F(i)-F(i-1)

end

t1=1:4

t2=1:14

G%G是预测结果，想显示结果可以去掉分号

plot(t1,A,'o',t2,G)%这是改坦画图检验，可以换成其告皮他的方法检验

之前我毕含想帮你答这个题的，后蠢李来想了想，不太可能，四个状态随机出现，在出现十次的前提下，预测第十一次将会出现什么状态，并且精准度要求 99.9%； 如果真有这个算法，估计这个世带数迟界已经不是现在这个模样了

2014.4.11~4.20 预测价格为 5.9535

>>%灰色预测模型键搏程序

clear

syms a b

c=[a b]'

A=[5.76 5.92 5.91 5.92 5.93 5.93 5.93 5.94 5.94 5.95 5.95 5.94 ]% 原始序列

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=[-Cones(1,n-1)]

c=inv(E*E')*E*D

c=c'

a=c(1)

b=c(2)

%预测 往后预测1个数据（搏亮胡即4.11-4.20价格）

F=[]F(1)=A(1)

for i=2:(n+1)

F(i)=(A(1)-b/a)/exp(a*(i-1))+b/a

end

G=[]G(1)=A(1)

for i=2:(n+1)

G(i)=F(i)-F(i-1)

end

t1=2002:2013

t2=2002:2014

G

plot(t1,A,'基拦o',t2,G)

G =

5.76005.91555.91895.92245.92585.92935.93275.93625.93965.94315.94655.95005.9535

---