粗大メモ置き場

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

2017-11-27から1日間の記事一覧

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(…