matlab - Matlab: diff(interp1(...'pchip')) は奇妙な結果を生成しますか?

Question

次の任意のデータがあります。

t=[0 1 2 3 4 5 6 7 8 9 10 12 14];
C=[0 1 5 8 10 8 6 4 3.0 2.2 1.5 0.6 0];
area=trapz(t,C);
E=C./area;
F=cumtrapz(t,E);

プロットすると、次のようになります。

今、私はプロットを滑らかにしたい:

t_int=[];
E_int=[];
F_int=[];

xmin=min(t);
xmax=max(t);
points=5000;
stepSize=(xmax-xmin)/points;

for i=xmin:stepSize:xmax
    t_int=[t_int i];
    E_int=[E_int interp1(t, E, i,'pchip')];
    F_int=[F_int interp1(t, F, i,'pchip')];
end

さて、ここで問題のある実際の部分に到達します。実際には、いくつかの F プロットがあり、それに対応する E プロットを取得したいと考えています。F プロットは E の積分であるため、E は微分 F プロットです。

それで：

dy=diff(F_int);
dx=diff(t_int);
dydx = dy./dx;
plot(t_int(:,1:length(dydx)),dydx,t_int,E_int)
legend('E (calculated)','E (actual)')

計算された E プロットが不安定なのはなぜですか?