灰色预测程序哪里出错

一路追逐

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]';
: \: G8 e3 w  E t0=find(a<=320);n=length(t0);
t1=cumsum(t0); %累加运算
& k/ j: f$ n2 @$ {B=[-0.5*(t1(1:end-1)+t1(2:end)),ones(n-1,1)];Y=t0(2:end);
( S. {2 N; B) I1 @7 T9 u& z; D  S r=B\Y
; e- D; T) |4 U4 j. h  a/ Z4 h y=dsolve('Dy+a*y=b','y(0)=y0');
y=subs(y,{'a','b','y0'},{r(1),r(2),t1(1)});
( I0 q- ~+ ?/ O7 a yuce1=subs(y,'t',[0:n+1])
%为提高预测精度，先计算预测值，再显示微分方程的解
y=vpa(y,6) %其中的6 表示显示6 位数字
! B$ f7 i( t# r4 j+ M9 D' Kyuce=diff(yuce1); %作差分运算，进行数据还原
1 C* K6 g9 t! k1 W" I% Ayuce=[t0(1),yuce]
yuce_new=yuce(n+1:end) %求得的两个预测值
5 U7 J( H' F* N
为什么在做差分运算时yuce=diff(yuce1）得到的结果是yuce =[ 0, 0, 0, 0, 0, 0, 0]
求赐教
) I/ o1 J# i/ ~; }! c: ?9 C* v