research38 min studyLesson 842 in Mathematics

Stochastic integrals in Hilbert spaces

Use predictable Hilbert–Schmidt integrands and the covariance-weighted Itô isometry to construct infinite-dimensional stochastic integrals.

On this page

What you’ll understand

  • Check the Hilbert–Schmidt condition for a cylindrical-noise integrand.
  • Apply the Hilbert-valued Itô isometry.
  • Weight an integrand norm correctly for a trace-class Q-Wiener input.

Stochastic integration can turn cylindrical noise into a genuine state-space random vector. The integrand must make the combined energy of all input directions summable, and it must be predictable so that it does not use future increments. The relevant size is a Hilbert–Schmidt norm, not an operator determinant or merely a bounded-operator norm. When the driving Wiener process already has covariance Q, the same calculation incorporates its square root. Keeping this covariance weight visible prevents double-counting or omitting the input's spatial correlation.

Sum scalar integrals with a square-summable operator

Let U and H be real separable Hilbert spaces, and let W be cylindrical Brownian motion on U relative to a filtration satisfying the usual completeness and right-continuity conditions. For a predictable integrand Φs\Phi_s taking values in the Hilbert–Schmidt operators from U to H, assume E∫0T∥Φs∥HS2ds<∞\mathbb E\int_0^T\|\Phi_s\|_{\mathrm{HS}}^2ds<\infty. The stochastic integral is the H-valued L² limit of sums of scalar Brownian integrals along an orthonormal input basis. Its value does not depend on that basis.

For simple predictable processes, coefficients are known at the beginning of each increment interval. Orthogonality of independent Brownian increments gives the isometry, which then completes the construction by density. Adaptedness alone needs suitable measurability before it supplies an admissible integrand; predictability is a clean sufficient formulation. A bounded operator can fail the Hilbert–Schmidt requirement in infinite dimension, as the identity does. This is the precise point where finite-dimensional intuition about arbitrary matrices stops applying.

∥Φ∥HS(U,H)2=∑k∥Φek∥H2,E∥∫0TΦs dWs∥H2=E∫0T∥Φs∥HS(U,H)2ds\|\Phi\|_{\mathrm{HS}(U,H)}^2=\sum_k\|\Phi e_k\|_H^2,\qquad \mathbb E\left\|\int_0^T\Phi_s\,dW_s\right\|_H^2=\mathbb E\int_0^T\|\Phi_s\|_{\mathrm{HS}(U,H)}^2ds
The predictable integrand has finite expected time-integrated Hilbert–Schmidt norm. This gives a centered square-integrable H-valued martingale.↔ Scroll wide equations horizontally.

A covariance square root belongs inside the norm

If the input is a U-valued Q-Wiener process, then its eigen-directions already have amplitudes qk\sqrt{q_k}. For a bounded predictable operator G from U to H, the isometry therefore uses ∥GQ1/2∥HS(U,H)2=∑kqk∥Gek∥H2\|GQ^{1/2}\|_{\mathrm{HS}(U,H)}^2=\sum_kq_k\|Ge_k\|_H^2. Equivalently, one can formulate integrands on the covariance space U0=Q1/2UU_0=Q^{1/2}U with its canonical inverse-square-root norm, treating the kernel of Q appropriately. These are two descriptions of the same noise weighting.

An integrand defined only on the covariance space need not extend to a bounded operator on all of U, so the bounded-G formula is a convenient sufficient setup rather than the most general theorem. Conversely, applying Q again after using the covariance-space norm would weight the noise twice. In calculations, identify whether the symbol denotes covariance itself or a noise-amplitude operator. Some SPDE references call the amplitude Q and hence write QQ adjoint in the covariance; that notation differs from the covariance-Q convention here.

E∥∫0TGs dWsQ∥H2=E∫0T∥GsQ1/2∥HS2ds=E∫0TTr⁡(GsQGs∗)ds\mathbb E\left\|\int_0^TG_s\,dW_s^Q\right\|_H^2=\mathbb E\int_0^T\|G_sQ^{1/2}\|_{\mathrm{HS}}^2ds=\mathbb E\int_0^T\operatorname{Tr}(G_sQG_s^*)ds
Here Q is the input covariance, not its amplitude. The displayed bounded-G formulation assumes the integrand is predictable and the right side is finite.↔ Scroll wide equations horizontally.

Energy control gives convergence, not classical time derivatives

If integrands converge in their expected integrated Hilbert–Schmidt norm, the isometry gives convergence of the terminal integrals in mean square. Martingale inequalities also control the expected supremum norm of the differences, yielding continuous H-valued versions in this square-integrable setting. With only almost-sure local square integrability, localization constructs a local martingale instead; a global finite second moment cannot be asserted without the stronger expectation condition. These are distinct levels of integrability.

A stochastic integral is not an ordinary time integral of an H-valued derivative of W. Its increments retain quadratic variation, and their energy is governed by the isometry. This makes integration useful even when the noise has no pointwise state-space realization. Later, a semigroup can appear inside the integrand and provide the missing Hilbert–Schmidt smoothing. That time-dependent regularization must be integrable near the current time, where the semigroup approaches the identity and the high-frequency noise is least damped.

An infinite diagonal integral

Let U=H=ell², W be identity-covariance cylindrical noise, and Phi ek=(1/k)ek be constant in time. At T=3 calculate the full integral's expected squared norm, the first two-coordinate contribution, and the mean squared tail error.

  1. The squared Hilbert–Schmidt norm of Phi is the sum of 1/k², namely pi²/6. Phi is therefore an admissible deterministic integrand despite the input not being ell²-valued.

  2. The isometry gives expected output energy 3pi²/6=pi²/2, approximately 4.934802. The first two directions contribute 3(1+1/4)=15/4=3.75.

  3. The remaining mean squared error is pi²/2−15/4, approximately 1.184802, when the same input Brownian coordinates define both the full and truncated integrals. Summing output energies rather than input energies is what makes this construction converge.

The result

A Hilbert–Schmidt operator can turn cylindrical input into an H-valued integral without changing the input's cylindrical nature.

Explore the model

Interactive lab

Interactive lab · matrix

An integrand norm is not its determinant

Pause and Reset, then set both controls to scale a=1 and shear b=1 and take eight Steps to time 4. The displayed A=[[1,1],[0,1]] has squared Frobenius norm 3 and determinant 1. Reset, set a=2 and b=1, and again take eight Steps to time 4; predict squared Frobenius norm 6 and determinant 2. As two-dimensional unit-covariance integrands, their one-time-unit output energies are 3 and 6. Readouts describe the selected final matrix; playback only interpolates geometry and does not generate stochastic integrals.

An integrand norm is not its determinantThe selected matrix is A = [[a,b],[0,1]]. Playback interpolates geometry between the identity and A; readouts always describe A. Its determinant is a and its singular values come from AᵀA. Positive a makes every supported matrix invertible. This two-dimensional family is not a general eigensolver or a learned network. Current readings: Area scale · det A: 1.4; Singular values · σmax, σmin: 1.597, 0.876; Condition number · κ₂(A): 1.823. All parameters are adjustable with the labeled controls below.Matrix A with first row 1.4, 0.6 and second row zero, one; display interpolates between identity and A while final determinant remains 1.4.−4−2024−2−1012x · equal geometric scaleyA = [ 1.4 0.6 ; 0 1 ]Display blend: 0% · final values below
Interpolated geometryOriginal unit squareSelected final transform
Area scale · det A
1.4
Singular values · σmax, σmin
1.597, 0.876
Condition number · κ₂(A)
1.823
t = 0 s

The selected matrix is A = [[a,b],[0,1]]. Playback interpolates geometry between the identity and A; readouts always describe A. Its determinant is a and its singular values come from AᵀA. Positive a makes every supported matrix invertible. This two-dimensional family is not a general eigensolver or a learned network. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forItô energy sums squared column lengths; a volume determinant answers a different question.

Check your understanding

For a bounded integrand G driven by covariance-Q Wiener noise, which quantity belongs in the Itô isometry?

Try it yourself

Take covariance Q=diag(1,1/4) and constant G=[[1,2],[0,1]] on R². Compute the squared covariance-weighted integrand norm and the expected squared integral norm over time two.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Multiply G by Q's positive square root before summing squared matrix entries.

Compare with the solution

The square root is diag(1,1/2), giving GQ^(1/2)=[[1,1],[0,1/2]]. Its squared Frobenius norm is 1+1+1/4=9/4. Over time two the expected squared integral norm is 9/2. Using G's unweighted squared Frobenius norm six would incorrectly give twelve, because it would treat the second noise direction as having unit variance instead of one-quarter.

Take it with you

  • Predictability and integrated Hilbert–Schmidt energy construct the integral.
  • The covariance square root determines the correct noise-weighted norm.
  • Terminal and pathwise convergence follow from stochastic norm estimates, not ordinary differentiation.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.