Abstract:We develop a efficient, easy-to-implement, and strictly monotone numerical integration method for Mean-Variance (MV) portfolio optimization in realistic contexts, which involve jump-diffusion dynamics of the underlying controlled processes, discrete rebalancing, and the application of investment constraints, namely no-bankruptcy and leverage. A crucial element of the MV portfolio optimization formulation over each rebalancing interval is a convolution integral, which involves a conditional density of the logarithm of the amount invested in the risky asset. Using a known closed-form expression for the Fourier transform of this density, we derive an infinite series representation for the conditional density where each term is strictly positive and explicitly computable. As a result, the convolution integral can be readily approximated through a monotone integration scheme, such as a composite quadrature rule typically available in most programming languages. The computational complexity of our method is an order of magnitude lower than that of existing monotone finite difference methods. To further enhance efficiency, we propose an implementation of the scheme via Fast Fourier Transforms, exploiting the Toeplitz matrix structure. The proposed monotone scheme is proven to be both $\ell_{\infty}$-stable and pointwise consistent, and we rigorously establish its pointwise convergence to the unique solution of the MV portfolio optimization problem. We also intuitively demonstrate that, as the rebalancing time interval approaches zero, the proposed scheme converges to a continuously observed impulse control formulation for MV optimization expressed as a Hamilton-Jacobi-Bellman Quasi-Variational Inequality. Numerical results show remarkable agreement with benchmark solutions obtained through finite differences and Monte Carlo simulation, underscoring the effectiveness of our approach.

Convergence and non-negativity preserving of the solution of balanced method for the delay CIR model with jump

Positivity and convergence of the balanced implicit method for the nonlinear jump-extended CIR model

On the Convergence of the Accelerated Riccati Iteration Method.

Double Smoothed Volatility Estimation of Potentially Non‐stationary Jump‐diffusion Model of Shibor

First order strong approximation of Ait-Sahalia-type interest rate model with Poisson jumps

On the backward Euler method for a generalized Ait-Sahalia-type rate model with Poisson jumps

Unconditionally positivity-preserving approximations of the Ait-Sahalia type model: Explicit Milstein-type schemes

Consistency Problems for Jump-Diffusion Models

A monotone numerical integration method for mean-variance portfolio optimization under jump-diffusion models

Balanced implicit Patankar–Euler methods for positive solutions of stochastic differential equations of biological regulatory systems

Unconditionally positivity-preserving explicit Euler-type schemes for a generalized Ait-Sahalia model

On De Finetti's control under Poisson observations: optimality of a double barrier strategy in a Markov additive model

Jump-Diffusion Risk-Sensitive Asset Management II: Jump-Diffusion Factor Model

Precautionary measures for credit risk management in jump models

On Positive Solutions of a Delay Equation Arising When Trading in Financial Markets

Impulse Control in Finance: Numerical Methods and Viscosity Solutions

Adaptive Control Barrier Functions with Vanishing Conservativeness Under Persistency of Excitation

The role of adaptivity in a numerical method for the Cox-Ingersoll-Ross model

Convergence of Gradient-Based Iterative Solution of Coupled Markovian Jump Lyapunov Equations.

Balanced truncation approach to model reduction of Markovian jump time-varying delay systems

Identifying the covariation between the diffusion parts and the co-jumps given discrete observations