灰色预测程序哪里出错

一路追逐

a=[390.6,412,320,559.2, 380.8,542.4,553,310,561,...
* v" w" n$ q' |! E 300,632,540,406.2,313.8,576,587.6,318.5]';
t0=find(a<=320);n=length(t0);
/ J/ m7 k  H# W. G4 r" S( d t1=cumsum(t0); %累加运算
( ~- @" n9 E% g. A9 Q& NB=[-0.5*(t1(1:end-1)+t1(2:end)),ones(n-1,1)];Y=t0(2:end);
# @$ L2 ?; ]5 }8 b2 a" O) o0 J r=B\Y
y=dsolve('Dy+a*y=b','y(0)=y0');
" W& w& ?! Q1 X' T3 g& r- H y=subs(y,{'a','b','y0'},{r(1),r(2),t1(1)});
7 [; V6 E' h  J yuce1=subs(y,'t',[0:n+1])
2 u6 n5 z  K0 K5 u' r %为提高预测精度，先计算预测值，再显示微分方程的解
y=vpa(y,6) %其中的6 表示显示6 位数字
yuce=diff(yuce1); %作差分运算，进行数据还原
yuce=[t0(1),yuce]
2 o7 Y* R* ~4 a( [# E' J2 o/ a' q( X yuce_new=yuce(n+1:end) %求得的两个预测值
1 @+ F( |+ g& j  j0 ]' U" W
为什么在做差分运算时yuce=diff(yuce1）得到的结果是yuce =[ 0, 0, 0, 0, 0, 0, 0]
求赐教
6 G, d) k6 h% b+ _0 _. g