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
Web文献
ヤホー
1階微分の4次精度中心差分の導出を教えて下さい df/dx=の形でお願いします ... - Yahoo!知恵袋
丁寧なサイトだが正負が違う。訂正をお願いしたいが連絡が取れないのがな。
4次精度・6次精度の微分公式の導出とFORTRAN90を用いた数値微分