灰色预测程序哪里出错

一路追逐

a=[390.6,412,320,559.2, 380.8,542.4,553,310,561,...
300,632,540,406.2,313.8,576,587.6,318.5]';
t0=find(a<=320);n=length(t0);
- J& J4 [  C. B* ~( o/ d" w/ k t1=cumsum(t0); %累加运算
9 ?5 q+ M6 g, D# O  K, ^9 L9 `* y1 xB=[-0.5*(t1(1:end-1)+t1(2:end)),ones(n-1,1)];Y=t0(2:end);
r=B\Y
4 t% W8 V- t& S y=dsolve('Dy+a*y=b','y(0)=y0');
, n( J( c( m8 o) `0 S& \ y=subs(y,{'a','b','y0'},{r(1),r(2),t1(1)});
& M4 s4 A  l% C6 l6 Y# N2 [) m yuce1=subs(y,'t',[0:n+1])
%为提高预测精度，先计算预测值，再显示微分方程的解
( L  t9 e' |8 P; n6 d/ M& z/ _y=vpa(y,6) %其中的6 表示显示6 位数字
. O4 E* V0 l4 P4 V$ fyuce=diff(yuce1); %作差分运算，进行数据还原
yuce=[t0(1),yuce]
yuce_new=yuce(n+1:end) %求得的两个预测值
  N- ~1 U0 |+ s2 E  F  T
! V! ~" `" F5 p 为什么在做差分运算时yuce=diff(yuce1）得到的结果是yuce =[ 0, 0, 0, 0, 0, 0, 0]
1 C. n7 v6 T9 k& r2 F求赐教