|
(1995)对预计SOR方法进行了详细讨论。然后,变分不等式(2.12)–(2.16)与下面的变分不等式问题具有相同的形式:Ligi(t,s)≤ 0,gi(t,s)≥ 嗨(t,s),(t,s)∈ [0,^T)×R+,(Ligi(T,s))(hi(T,s)- gi(t,s))=0(t,s)∈ [0,^T)×R+,gi(^T,s)=hi(^T,s),s∈ R+。(3.2)这里,gi(t,s)表示函数Vi(t,s),Ji(t,s)的值,-Ui(t,s)、Ki(t,s)或Pi(t,s)。函数hi(t,s)表示fi(t,s)-c、 (六)(t,s)-(fi(t,s)+^c))+,-(fi(t,s)+^c),(fi(t,s)-c)-Ui(t,s))+,或max{Ai(t,s),Bi(t,s)}。期货价格fi(t,s),含^t≤ T,由(2.4)给出。我们现在考虑在时间(δt=^TN)和空间(δs=SmaxM)离散的非均匀网格上,偏微分方程Ligi(t,s)=0的离散化。将s-导数的CrankNicolson方法和t-导数的后向差分应用于求积方程,得出了有限差分方程:-r(gm,ni+gm,n-1i)+gm,ni- gm,n-1iδt+Дmi(gm+1,ni- 转基因的-1,ni2δs+gm+1,n-1i-转基因的-1,n-1i2δs)+(σmi)(gm+1,ni- 2gm,镍+gm-1,ni2(δs)+gm+1,n-1i- 2gm,n-1i+gm-1,n-1i2(δs))+qii(gm,ni+gm,n-1i)+Xj6=iqij(gm,nj+gm,n-1j)=0,(3.3),其中qii=-Pj6=使用iqijis。为方便起见,我们可以使用这些常数的下标符号,例如gm、ni≡ gi(nδt,mδs),hm,ni≡ hi(nδt,mδs),Дmi=Дi(mδs),σmi=σi(mδs)。通过明确对待政权政变条款来实现。更换gm,n-1J与gm,njfor j 6=获得- αmigm-1,n-1i+(1- βmi)gm,n-1i- γmigm+1,n-1i=αmigm-1,ni+(1+βmi)gm,ni+γmigm+1,n+δtXj6=iqijgm,nj,(3.4)其中αmi=δt4δs(σmi)δs- ^1英里,βmi=-δt(r)-qii)+(σmi)(δs),γmi=δt4δs(σmi)δs+Дmi,(3.5)对于i,j∈ E、 m=1,2。。。,M- 1和n=1,2。。。,N- 1.
|