Continuous-time stochastic processes underlie many natural and engineered systems. In healthcare, autonomous driving, and industrial control, direct interaction with the environment is often unsafe or impractical, motivating offline reinforcement learning from historical data. However, there is limited statistical understanding of the approximation errors inherent in learning policies from offline datasets. We address this by linking reinforcement learning to the Hamilton-Jacobi-Bellman equation and proposing an operator-theoretic algorithm based on a simple dynamic programming recursion. Specifically, we represent our world model in terms of the infinitesimal generator of controlled diffusion processes learned in a reproducing kernel Hilbert space. By integrating statistical learning methods and operator theory, we establish global convergence of the value function under sufficiently small generator-approximation error and derive explicit finite-sample guarantees with bounds tied directly to system properties such as smoothness and stability. Our theoretical and proof-of-concept results indicate that operator-based approaches may hold promise for solving offline reinforcement learning using continuous-time optimal control.
Figures & tables
Fig. 1: Overview of the Operator-Based Continuous-Time Offline Reinforcement Learning ( O-CTRL ) algorithm: a generator world model based on a Reproducing Kernel Hilbert Space (RKHS) representation of state–action data enables dynamic programming for optimal value functions, illustrated on the swing-up pendulum task from Gymnasium [ 2 ] .
Algorithm 1 O perator-Based C ontinuous- T ime offline R einforcement L earning ( O-CTRL )
Fig. 2: End-to-end error decomposition in ( 26 ) underlying Theorem VII.5 for bounding the term ∥V⋆−Vk∥Lμ2 .
Task
Complexity
World Model Construction
O(N3+N2(ns+na))
Dynamic Programming
O(N3+kmax(N2na+NC(na)))
TABLE I: Computational complexity of Algorithm 1 . C(na) denotes the complexity of evaluating one sampled Fenchel-conjugate term Da(λ(s(i)),s(i)) , e.g. C(na)=na for separable action penalties or na2 for quadratic (coupled) action penalties.
Fig. 3: Value function learning for one-dimensional benchmark systems. The error convergence confirms our worst-case analysis.
Fig. 4: Learned value function and policy for the linear SDE with additive action using 200 training samples.
Fig. 5: Learned value function and policy for the nonlinear SDE with affine action using 200 training samples.
Fig. 6: Nonlinear policy obtained by O-CTRL on the Pendulum-v1 replay dataset using 8000 offline transition samples converted to infinitesimal tuples ( 11 ).
Fig. 7: AUV depth and control trajectories for the two objectives. Solid curves show the pooled mean over five learned models and 50 trials per model. Shaded regions show one population standard deviation. Dashed lines show the target or input limits.
Mean Episode Reward ( ± std)
Algorithm
Random
Replay
Ours
−141.71±84.23
−104.36±83.00
CQL
−293.09±206.68
−451.00±530.91
IQL
−313.00±222.91
−284.89±196.59
TD3+BC
−1181.08±347.40
−748.11±420.18
Expert
−130.98±79.40
TABLE II: Pendulum results on d3rlpy offline datasets (Random, Replay). Higher (less negative) is better. The expert policy is a Stable-Baselines3 Soft Actor-Critic (SAC) agent (Pendulum-v1) trained online (dataset not applicable). 4 4 4 https://huggingface.co/sb3/sac-Pendulum-v1