2013-09-01から1ヶ月間の記事一覧

数値微分 その3

前回は1回微分近似を扱いましたが,今回は2回微分近似公式を確認してみます やってることはほとんど同じですが・・・ num_differential3.wxm f(x) : 評価関数 h : 摂動(≠0) 今f(x)をx0点周りで4次の項までテイラー展開します(%o1) 上式をx = x0 + hとして…

数値微分 その2

前回の差分公式より精度が高い3点近似公式を確認してみます num_differential2.wxm f(x) : 評価関数 h : 摂動(≠0) 今f(x)をx0点周りで3次の項までテイラー展開します(%o1) 上式をx = x0 + hとして%o2式にまとめます 同様に,x = x0 + 2*hとして%o3式にま…

数値微分 その1

前回の数値積分からの流れで,今回は数値微分(Numerical differentiation)のお話です 数値微分は解析的に微分することが難しい問題を近似的に解くための大変便利な手法ですヽ( ´ー`)ノ 微分方程式の数値解法も広義の意味で数値微分に含まれます 今回は数値…