粗大メモ置き場

個人用,たまーに来訪者を意識する雑記メモ

MATLAB 位相遅れなし微分 (4次精度)

メモ。

導出はTaylor展開してtにh,2hを代入して解きます。


4次精度の微分はこう。

(f(x - 2h) - 8 f(x - h) + 8 f(x + h) - f(x + 2h))/(12h) + O(h⁴)

function V = deriv4(Etime,Epos)
em2 = Epos(1:end-4);
em1= Epos(2:end-3);
%Epos(3:end-2)
ep1 = Epos(4:end-1);
ep2 = Epos(5:end);

V = ( 2/3*(ep1-em1)+1/12*(em2-ep2))/(Etime(2)-Etime(1));
end