|
为了捕捉这些模型,我们定义了通用的微分算子{·}=-r·+·t+~n(s)·s+σ(s)·s、 然后变分不等式(4.5),(4.6),(4.7),(4.8)和(4.9)的形式与下面的变分不等式问题相同:Lg(t,s)≤ 0,g(t,s)≥ ξ(t,s),(t,s)∈ [0,^T)×R+,(Lg(T,s))(ξ(T,s)- g(t,s))=0(t,s)∈ [0,^T)×R+,g(^T,s)=ξ(^T,s),s∈ R+。(6.1)这里,g(t,s)表示值函数V(t,s),J(t,s),-U(t,s)、K(t,s)或P(t,s)。函数ξ(t,s)表示f(t,s;t)-c、 (V(t,s)-(f(t,s;t)+^c)+,-(f(t,s;t)+^c,(f(t,s;t)-c)-U(t,s))+,或max{A(t,s),B(t,s)}。期货价格f(t,s;t),带^t≤ T分别由OU、CIR和XOU模型下的(2.1)、(2.4)和(2.10)给出。我们现在考虑在时间(δt=^TN)和空间(δs=SmaxM)离散的均匀网格上,偏微分方程Lg(t,s)=0的离散化。我们采用CrankNicolson方法,该方法涉及有限差分方程:-αigi-1,j-1+ (1 - βi)gi,j-1.- γigi+1,j-1=αigi-1,j+(1+βi)gi,j+γigi+1,j,其中gi,j=g(jδt,iδs),ξi,j=ξ(jδt,iδs),ψi=ψ(iδs),σi=σ(iδs)。αi=δt4δsσiδs- ~ni, βi=-δtr+σi(δs), γi=δt4δsσiδs+ψi,对于i=1,2。。。,M- 1和j=1,2。。。,N- 1.需要及时解决的系统isMgj-1=rj,其中右侧isrj=Mgj+αg0,j-1+g0,j。。。+ γM-1....通用汽车,j-1+gM,j,,安德姆=1.- β-γ-α1 - β-γ-α1 - β-γ.........-αM-21- βM-2.-γM-2.-αM-11- βM-1.,M=1 + βγα1 + βγα1 + βγ.........αM-21+βM-2γM-2αM-11+βM-1.,gj=g1,j,g2,j,转基因的-1,j这导致了一系列平稳互补问题。因此,在每个时间点j∈{1, 2, . . .
|