|
时间倒向求解的系统isMj-1Vj-1=rj,(4.4),其中右侧isrj=MjVj+α1,j-1V0,j-1+α1,jV0,j。。。+...γM-1,j-1VM,j-1+γM-1,jVM,j, (4.5)和MJ=1.- β -γ1,j-α2,j1- β -γ2,j-α3,j1- β -γ3,j。。。。。。。。。-αM-2,j1- β -γM-2,j-αM-1,j1- β, (4.6)Mj=1+βγ1,jα2,j1+βγ2,jα3,j1+βγ3,j。。。。。。。。。αM-2,j1+βγM-2,jαM-1,j1+β, (4.7)Vj=V1,j,V2,j,虚拟机-1,jT、 (4.8)这导致了一系列平稳互补问题。由于交易者可以在到期前的任何时候平仓,价值函数V(t,x)必须满足约束条件V(t,x)≥ f(t,x),Xmin≤ x个≤ Xmax,0≤ t型≤\'\'T。(4.9)相应地,离散方案可以写成VI,j≥ fi,j,0≤ 我≤ M、 0个≤ j≤ N、 (4.10)因此,在每个时间步骤j∈ {1,2,…,N- 1} ,我们需要解决线性互补问题Mj公司-1Vj-1.≥ rj,Vj-1.≥ fj公司-1,(Mj-1gj-1.- rj)T(fj-1.- gj公司-1) = 0.(4.11)特别是,我们的算法强制执行支配支付函数的约束,如下所示:VNEW,j-1=最大值沃迪,j-1,fi,j-1.. (4.12)投影SOR方法用于求解线性系统。请注意,约束与V(k+1)i,j同时被加强-1在每次迭代中计算。因此,在每个时间步j,SOR算法是迭代(在k上)方程sv(k+1)1,j-1=最大值f1,j-1,V(k)1,j-1+ω1 - β[r1,j- (1 - β) V(k)1,j-1+γ1,j-1V(k)2,j-1],V(k+1)2,j-1=最大值f2,j-1,V(k)2,j-1+ω1 - β[r2,j+α2,j-1V(k+1)1,j-1.- (1 - β) V(k)2,j-1+γ2,j-1V(k)3,j-1],...V(k+1)M-1,j-1=最大值fM公司-1,j-1,V(k)M-1,j-1+ω1 - β[rM-1,j+αM-1,j-1V(k+1)M-2,j-1.- (1 - β) V(k)M-1,j-1],其中k是迭代计数器,ω是超松弛参数。迭代格式从初始点V(0)j开始-1并继续,直到满足收敛标准,如| | V(k+1)j-1.-V(k)j-1 | |<,w这里的是一个公差参数。
|