袁 玲，汪 慧，梁 靜

（安徽新華學(xué)院通識(shí)教育部，安徽 合肥 230088）

0 引言

本文研究Stratonovich型標(biāo)量自治的隨機(jī)微分方程[1]如下：

其中，W(t)={Wt,t0≤t≤T}為標(biāo)準(zhǔn)維納過程，f與g是[t0，T]×R上均滿足Lipschitz條件與線性增長條件的可測(cè)函數(shù)。本文運(yùn)用彩色樹理論[2]，構(gòu)造求解該方程的三階隱式型Runge-Kutta算法——IMRK算法。

1 三階隱式型隨機(jī)Runge-Kutta算法IMRK算法的構(gòu)造

為了構(gòu)造3階隱式型隨機(jī)Runge-Kutta算法，先選定一個(gè)三階隱式型的確定性Runge-Kutta方法如下將其作為3階隱式型Runge-Kutta算法的確定性部分，則可以設(shè)在隱式型Runge-Kutta算法的Butcher表[3]中有：

則使上述算法具有1.0階全局收斂性的充要條件[4]是：

不超過1階的樹只有樹1、樹2和樹6；1.5階的樹為樹4、樹5、樹20和樹22，它們的Stratonovich型局部誤差系數(shù)[3]如下：

由此推出的等價(jià)方程為：

樹4、樹5、樹20和樹22對(duì)應(yīng)的階條件為：

再以得到最小主誤差常數(shù)[4]的原則推出相應(yīng)的等價(jià)方程為：

運(yùn)用已知條件aT和A，并且取，解得滿足上述方程組（*）和（**）的一組解為

即得到了具有最小主誤差常數(shù)的強(qiáng)1階收斂的3級(jí)隱式型隨機(jī)Runge-Kutta算法，并將其稱為IMRK算法。

IMRK算法具體如下：

2 IMRK算法的均方穩(wěn)定性

結(jié)論：IMRK2算法的均方穩(wěn)定函數(shù)[5]為R(p,q)，其中

證明：將IMRK2算法（3）應(yīng)用于求解與Ito型線性檢驗(yàn)隨機(jī)微分方程：dy=aydt+bydw(t)

即

則均方穩(wěn)定函數(shù)為：

圖1 Runge-Kutta——IMRK算法與5種算法的均方穩(wěn)定域比較

將隱式型隨機(jī)Runge-Kutta算法——IMRK算法與 Euler算法[1]，Heun算法[2]，Milstein 算法[4]，PL算法[6]，M2算法[4]的均方穩(wěn)定域比較如圖1。由圖像可知，IMRK算法具有更廣的穩(wěn)定區(qū)間，即本文構(gòu)造在穩(wěn)定性方面具有其自身的優(yōu)勢(shì)。

3 IMRK算法精度數(shù)值實(shí)驗(yàn)

選取方程：

取a=-26，b=1，用平均誤差[7]M來表示算法的精度，，其中，和分別表示第次模擬時(shí)在點(diǎn)處的數(shù)值解和準(zhǔn)確解，是模擬次數(shù)取為，比較結(jié)果如表1所示。

表1 不同步長時(shí)誤差精度比較

從表1中可以看出，本文的IMRK算法與現(xiàn)有算法相比，具有更高的精度。