預估-校正方法

數值分析中,預估-校正方法是一類求解常微分方程的算法 - 找到一個未知的函數以滿足一定微分方程。 所有這類算法以如下兩個步驟進行:

  1. 首先,"預估"步,基於之前若干步的一組函數值及導數值擬合出的函數出發,進而外插此函數在後續點的值。
  2. 其次,"校正"步,通過使用函數的 預估 值和 另一種方法 改進初始近似,以內插這一未知的函數在相同後續點的值。

預估-校正方法求解常微分方程

編輯

對於常微分方程(ODE)的數值解,預估–校正方法通常使用一個顯式方法作為預估步和一個隱式方法作為校正步。

例如:歐拉方法與梯形規則

編輯

一個簡單的預估–校正方法(即Heun方法)可以由歐拉法 (一個顯式方法)和梯形規則 (一個隱式方法)構成。

考慮如下微分方程

 

並將步長大小記作  .

首先,預估步:從當前的值  開始,由歐拉方法計算初步估算值  

 

然後,校正步:使用梯形規則改善初步估算值,

 

這一值在下一步使用。

PEC模式和PECE模式

編輯

取決於應用校正步的頻繁程度,預估-校正方法可以有多種不同的變體。 預估–評價–校正–評價(PECE)模式指的是,以上述例子為例:

 

也可以運用預估–評價–校正(PEC)模式,該模式每步只評價函數 f 一次:

 

此外,校正步可以多次重複,以希望實現一個更好的真實解的逼近值。 如果校正步運行兩次即為PECECE模式:

 

PECEC模式比PECECE模式少了一次函數評價過程。

更一般地說,如果校正步運行 k 次,此方法即為P(EC)k 或P(EC)kE模式。 如果校正步是迭代直到收斂,這可以被稱為PE(CE).

註解

編輯

參考文獻

編輯

外部連結

編輯