中心差分法詳見:
由泰勒公式可得到:
由(1)可得
或者同理,由(2)可得
(6)稱為求 的一階向前差分公式。(7)稱為求 的一階向後差分公式。
由(1)(3)可得求 的一階向前差分公式:
一階向前差分法的係數見下表。
一階向後差分法的係數見下表。
由(1)(3)消去可得即或者
(10)稱為求 的二階向前差分公式。二階向前差分法的係數見下表。
二階向後差分法的係數見下表。
★★★★往期相關★★★★
數值微分|多項式的導數計算c++版
數值微分|多項式的導數計算
通過案例學python之定義函式類
數值積分|第二類反常積分
數值積分|第一類反常積分
數值積分|中點法則(midpoint rule)
數值積分|龍貝格公式
數值積分|自適應辛普森積分公式
數值積分|自適應梯形積分
數值積分|牛頓-柯特斯公式
數值積分|高斯積分
數值積分|泰勒(taylor)公式求積分
數值積分| 辛普森公式
python實現辛普森公式
Python求解常微分數值解
sympy是乙個python的科學計算庫,用一套強大的符號計算體系完成諸如多項式求值 求極限 解方程 求積分 微分方程 級數展開 矩陣運算等等計算問題。雖然matlab的類似科學計算能力也很強大,但是python以其語法簡單 易上手 異常豐富的三方庫生態,個人認為可以更優雅地解決日常遇到的各種計算問...
數值微分與數值積分
注意 diff函式計算的是向量元素間的差分,故差分向量元素的個數比原向量少了 乙個。同樣,對於矩陣來說,差分後的矩陣比原矩陣少了一行或一列。另外,計算差分之後,可以用f x 在某點處的差商作為其導數的近似值。matlab提供了求向前差分的函式diff,其呼叫格式有三種 dx diff x 計算向量x...
數值微分C 實現
取h 0.1時用euler法,改進euler法,rung kutta方法求其數值解並與精確解進行比較。輸入 求解區間,初值,數值解個數 輸出 數值解 改進尤拉公式 四階龍格 庫塔方法 改進的尤拉方程 四階龍格 庫塔方法 定義係數c double c 儲存生成的節點橫座標 double datax 1...