|
我们假设网格大小δt、δx和Δξ的选择使得网格数Nt、NxandNξ是自然数。我们定义了向量φk(ix,iξ)九∈{0,··,Nx},iξ∈{0,···,Nξ}φ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ξ))= 其中il=dlδxe,ixζ=dζδxe和iζζ=dΓ(ζ)ΔΞe。我们现在引入常数h∈ R+满足h>δt+2λΞ(g(x))+λL,并将HJBKVI(13)转化为等价的定点问题φk=max苏普勒∈Lδn\'Llφk+\'fl,ko,supζ∈Zδ(xix)nMζφk+kζo, (14) 在哪里=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数值结果在本节中,我们研究了通过前一节所述方法获得的最优策略的特征。根据恢复强度的类型和限制指令的限制,我们将本节分为四个小节。
|