我们假设选择网格尺寸δt、δx和δξ,使网格数Nt、NxandNξ为自然数。我们定义了向量φk(ix,iξ)九∈{0,···,Nx},iξ∈{0,····,Nξ}乘以φk(ix,iξ)=φtk(xix,ξiξ),从中我们得到矩阵形式的HJBQVI(12):maxφk+1(ix,iξ)- φk(ix,iξ)δt+supl∈LδnλLhφk(ix-il,iξ)- φk(ix,iξ)+lsio+λΞ(ξiξ)hφk(ix,iξ)-(1)- φk(ix,iξ)+xixΔΞi,(13)supζ∈Zδ(xi-x) (φk(ix-ixζ,iξ+iξζ)- xixΓ(ζ)- φk(ix,iξ))= 0,其中il=dlδxe,ixζ=dζδxe,iξζ=dΓ(ζ)δΞe。我们现在引入常数h∈ R+满足H>δt+2λΞ(g(x))+λL,并将HJBQVI(13)转换为等效的定点问题φk=maxsupl公司∈Lδn'Llφk+'fl,ko,supζ∈Zδ(xix)nMζφk+kζo, (14) 其中,Llij=1.-Hδt+λΞ(ξiξ)+λL, i=j=(ix,iξ)λΞ(ξiξ)h,i=(ix,iξ),j=(ix,iξ)- 1) λLh,i=(ix,iξ),j=(ix- il,iξ)0,否则,\'fl,k(ix,iξ)=hδtφk+1(ix,iξ)+λΞ(ξiξ)xiixδΞ+λL,Mζij=1ni=(ix,iξ),j=(ix-il,iξ+iξζ)o,Kζ(ix,iξ)=xixΓ(ζ)。我们发现,由于因子h,它是以矩阵形式显示的收缩图,而Mζ是非扩张图。我们参考[Ied13,第3.2节]了解方程(14)中问题的求解过程。4数值结果在本节中,我们研究了通过前一节所述方法获得的最优策略的特征。根据恢复强度的类型和限制指令的限制,我们将本节分为四个小节。
|