Skip to content
homeaboutworkprojectsthesiswritingresume
Loading
~/blog/gaussian-vectors-and-processes0%dark
  1. home/
  2. writing/
  3. Gaussian Vectors and Processes

06 June 2026 · 6 min read · updated 09 June 2026

Gaussian Vectors and Processes

A Gaussian law is the unique distribution whose linear images are again Gaussian, and it is pinned down entirely by a mean and a covariance. We compute the Gaussian characteristic function, define a Gaussian vector through the Gaussianity of all its linear combinations, prove its law is determined by mean and covariance, prove uncorrelated Gaussian coordinates are independent, and construct a Gaussian process from any mean function and positive semidefinite covariance function. This is the class of processes the Karhunen-Loeve expansion diagonalises and Brownian motion belongs to.

  • 13 results
  • 15 connections
  • probability
  • gaussian
  • stochastic-processes
On this page▾
  • The Gaussian characteristic function
  • Gaussian vectors
  • Gaussian processes

6 min left

  • The Gaussian characteristic function1m
  • Gaussian vectors3m
  • Gaussian processes2m

The Gaussian distribution is the fixed point of linear operations. A linear image of a Gaussian is Gaussian, a sum of independent Gaussians is Gaussian, and the limit of normalised sums is Gaussian, which is why it is the universal law of aggregated randomness. Measure-theoretically a Gaussian is determined by two numbers, its mean and its variance, and a Gaussian vector by a mean vector and a covariance matrix, the minimal second-order data. This post builds the Gaussian vector and the Gaussian process from the characteristic function of the previous post [1], [2].

#The Gaussian characteristic function

Definition1

A random variable is standard normal, written Z∼N(0,1)Z\sim N(0,1)Z∼N(0,1), when it has density 12πe−x2/2\frac{1}{\sqrt{2\pi}} e^{-x^2/2}2π​1​e−x2/2. A variable X=μ+σZX=\mu+\sigma ZX=μ+σZ is Gaussian N(μ,σ2)N(\mu,\sigma^2)N(μ,σ2), including the degenerate constant μ\muμ when σ=0\sigma=0σ=0.

Proposition2

The standard normal has characteristic function φZ(t)=e−t2/2\varphi_Z(t)=e^{-t^2/2}φZ​(t)=e−t2/2, and N(μ,σ2)N(\mu,\sigma^2)N(μ,σ2) has φX(t)=eiμt−σ2t2/2\varphi_X(t)=e^{i\mu t-\sigma^2 t^2/2}φX​(t)=eiμt−σ2t2/2.

Proof

Differentiating φZ(t)=∫eitx12πe−x2/2 dx\varphi_Z(t)=\int e^{itx}\frac{1}{\sqrt{2\pi}}e^{-x^2/2}\,dxφZ​(t)=∫eitx2π​1​e−x2/2dx under the integral, licensed because E∣Z∣<∞\E\abs Z<\inftyE∣Z∣<∞, gives φZ′(t)=∫ix eitx12πe−x2/2 dx\varphi_Z'(t)=\int ix\,e^{itx}\frac{1}{\sqrt{2\pi}}e^{-x^2/2}\,dxφZ′​(t)=∫ixeitx2π​1​e−x2/2dx. Integrating by parts with ddxe−x2/2=−xe−x2/2\frac{d}{dx}e^{-x^2/2}=-xe^{-x^2/2}dxd​e−x2/2=−xe−x2/2 moves the xxx onto the exponential. The boundary term [−i2πeitxe−x2/2]x=−∞x=+∞[-\tfrac{i}{\sqrt{2\pi}}e^{itx}e^{-x^2/2}]_{x=-\infty}^{x=+\infty}[−2π​i​eitxe−x2/2]x=−∞x=+∞​ vanishes because ∣eitx∣=1\abs{e^{itx}}=1​eitx​=1 and e−x2/2→0e^{-x^2/2}\to0e−x2/2→0 at ±∞\pm\infty±∞, leaving φZ′(t)=i2π∫(it)eitxe−x2/2 dx=−t φZ(t)\varphi_Z'(t)=\frac{i}{\sqrt{2\pi}}\int(it) e^{itx}e^{-x^2/2}\,dx=-t\,\varphi_Z(t)φZ′​(t)=2π​i​∫(it)eitxe−x2/2dx=−tφZ​(t). With φZ(0)=1\varphi_Z(0)=1φZ​(0)=1 this linear equation has the unique solution φZ(t)=e−t2/2\varphi_Z(t)=e^{-t^2/2}φZ​(t)=e−t2/2. For X=μ+σZX=\mu+\sigma ZX=μ+σZ, φX(t)=E[eit(μ+σZ)]=eiμtφZ(σt)=eiμt−σ2t2/2\varphi_X(t)=\E[e^{it(\mu+\sigma Z)}]=e^{i\mu t} \varphi_Z(\sigma t)=e^{i\mu t-\sigma^2t^2/2}φX​(t)=E[eit(μ+σZ)]=eiμtφZ​(σt)=eiμt−σ2t2/2.

The quadratic exponent is the signature of the Gaussian. The transform of N(μ,σ2)N(\mu,\sigma^2)N(μ,σ2) is an exponential of a quadratic in ttt, and reading off the coefficients recovers the mean and variance.

#Gaussian vectors

The right definition of a Gaussian vector asks that the distribution survive every projection to a line.

Definition3

A random vector X∈RnX\in\R^nX∈Rn is Gaussian when every linear combination a⊤X=∑iaiXia^\top X=\sum_i a_iX_ia⊤X=∑i​ai​Xi​ is a univariate Gaussian. Its mean is μ=E[X]\mu=\E[X]μ=E[X] and its covariance is the matrix Σij=Cov⁡(Xi,Xj)=E[(Xi−μi)(Xj−μj)]\Sigma_{ij}=\Cov(X_i,X_j) =\E[(X_i-\mu_i)(X_j-\mu_j)]Σij​=Cov(Xi​,Xj​)=E[(Xi​−μi​)(Xj​−μj​)].

The covariance is symmetric and positive semidefinite, since a⊤Σa=Var⁡(a⊤X)≥0a^\top\Sigma a=\Var(a^\top X)\ge 0a⊤Σa=Var(a⊤X)≥0. Mean and covariance are all the data there is.

Theorem4

A Gaussian vector has characteristic function φX(t)=exp⁡(i t⊤μ−12t⊤Σt)\varphi_X(t)=\exp(i\,t^\top\mu-\tfrac12 t^\top\Sigma t)φX​(t)=exp(it⊤μ−21​t⊤Σt), so its law is determined by μ\muμ and Σ\SigmaΣ alone.

Proof

For fixed t∈Rnt\in\R^nt∈Rn the scalar t⊤Xt^\top Xt⊤X is Gaussian by definition, with mean t⊤μt^\top\mut⊤μ and variance E[(t⊤X−t⊤μ)2]=t⊤Σt\E[(t^\top X-t^\top\mu)^2]=t^\top\Sigma tE[(t⊤X−t⊤μ)2]=t⊤Σt. Evaluating its scalar characteristic function at argument 111, φX(t)=E[ei t⊤X]=E[ei(t⊤X)]=exp⁡(i t⊤μ−12t⊤Σt)\varphi_X(t)=\E[e^{i\,t^\top X}]=\E[e^{i(t^\top X)}]=\exp(i\,t^\top\mu-\tfrac12 t^\top\Sigma t)φX​(t)=E[eit⊤X]=E[ei(t⊤X)]=exp(it⊤μ−21​t⊤Σt) by Proposition 2. The characteristic function depends only on μ\muμ and Σ\SigmaΣ. If two laws on Rn\R^nRn share a characteristic function, then for every ttt the projections x↦t⊤xx\mapsto t^\top xx↦t⊤x have equal scalar characteristic functions along the ray, so by the uniqueness theorem all one-dimensional projections coincide, and by the Cramer-Wold device the projections determine the joint law. Hence the law depends only on μ\muμ and Σ\SigmaΣ.

In the Gaussian world the second moment is the whole story, and that collapses the distinction between independence and zero correlation.

Corollary5

The coordinates of a Gaussian vector are independent if and only if they are uncorrelated. More generally, subvectors X(1)X_{(1)}X(1)​ and X(2)X_{(2)}X(2)​ are independent if and only if their cross-covariance block vanishes.

Proof

Independence always implies zero covariance. Conversely, if the cross-covariance block is zero, then Σ\SigmaΣ is block diagonal, so the quadratic form splits, t⊤Σt=t(1)⊤Σ(1)t(1)+t(2)⊤Σ(2)t(2)t^\top\Sigma t=t_{(1)}^\top\Sigma_{(1)}t_{(1)}+ t_{(2)}^\top\Sigma_{(2)}t_{(2)}t⊤Σt=t(1)⊤​Σ(1)​t(1)​+t(2)⊤​Σ(2)​t(2)​, and the characteristic function factors, φX(t)=φX(1)(t(1)) φX(2)(t(2))\varphi_X(t)=\varphi_{X_{(1)}}(t_{(1)})\,\varphi_{X_{(2)}}(t_{(2)})φX​(t)=φX(1)​​(t(1)​)φX(2)​​(t(2)​). Each factor is a genuine marginal characteristic function, because a subvector X(1)X_{(1)}X(1)​ is itself Gaussian (any a⊤X(1)=(a,0)⊤Xa^\top X_{(1)}=(a,0)^\top Xa⊤X(1)​=(a,0)⊤X is a linear combination of all coordinates of XXX, hence univariate Gaussian by Definition 3) with covariance the leading block Σ(1)\Sigma_{(1)}Σ(1)​, so by Theorem 4 its characteristic function is exp⁡(i t(1)⊤μ(1)−12t(1)⊤Σ(1)t(1))\exp(i\,t_{(1)}^\top\mu_{(1)}-\tfrac12 t_{(1)}^\top\Sigma_{(1)}t_{(1)})exp(it(1)⊤​μ(1)​−21​t(1)⊤​Σ(1)​t(1)​), the first factor, and likewise for X(2)X_{(2)}X(2)​. A joint characteristic function that factors into the marginals is the product law, so the subvectors are independent.

In general zero correlation is weaker than independence; the Gaussian collapses the two, which is what makes it the tractable model of dependence.

Conversely, any mean and any positive semidefinite covariance are realised by some Gaussian vector.

Proposition6

For every μ∈Rn\mu\in\R^nμ∈Rn and symmetric positive semidefinite Σ\SigmaΣ, there is a Gaussian vector with mean μ\muμ and covariance Σ\SigmaΣ.

Proof

Being symmetric positive semidefinite, Σ\SigmaΣ factors as Σ=AA⊤\Sigma=AA^\topΣ=AA⊤, for instance through its spectral decomposition Σ=QΛQ⊤\Sigma=Q\Lambda Q^\topΣ=QΛQ⊤ with A=QΛ1/2A=Q\Lambda ^{1/2}A=QΛ1/2. Let Z=(Z1,…,Zn)Z=(Z_1,\dots,Z_n)Z=(Z1​,…,Zn​) have independent standard normal coordinates and set X=μ+AZX=\mu+AZX=μ+AZ. Then E[X]=μ\E[X]=\muE[X]=μ and Cov⁡(X)=A Cov⁡(Z) A⊤=AA⊤=Σ\Cov(X)=A\,\Cov(Z)\,A^\top=AA^\top=\SigmaCov(X)=ACov(Z)A⊤=AA⊤=Σ, and every a⊤X=a⊤μ+(A⊤a)⊤Za^\top X=a^\top\mu+(A^\top a)^\top Za⊤X=a⊤μ+(A⊤a)⊤Z is a linear combination of independent Gaussians, hence Gaussian, so XXX is a Gaussian vector with the required mean and covariance.

#Gaussian processes

A process is a family of random variables indexed by a parameter, usually time. The Gaussian property is imposed on every finite subfamily at once.

Definition7

A Gaussian process (Xt)t∈T(X_t)_{t\in T}(Xt​)t∈T​ is a family of random variables for which every finite vector (Xt1,…,Xtn)(X_{t_1},\dots,X_{t_n})(Xt1​​,…,Xtn​​) is a Gaussian vector. It is described by its mean function m(t)=E[Xt]m(t)=\E[X_t]m(t)=E[Xt​] and its covariance function K(s,t)=Cov⁡(Xs,Xt)K(s,t)=\Cov(X_s,X_t)K(s,t)=Cov(Xs​,Xt​).

The covariance function is symmetric and positive semidefinite, meaning every matrix (K(ti,tj))i,j(K(t_i,t_j))_{i,j}(K(ti​,tj​))i,j​ is positive semidefinite, being the covariance of (Xt1,…,Xtn)(X_{t_1},\dots,X_{t_n})(Xt1​​,…,Xtn​​). These two functions determine the process, and conversely any such pair is realised.

Theorem8

For every function mmm on TTT and every symmetric positive semidefinite kernel KKK on T×TT\times TT×T, there is a Gaussian process with mean function mmm and covariance function KKK, unique in law.

Proof

For each finite set t1,…,tnt_1,\dots,t_nt1​,…,tn​ define μ(n)=(m(ti))i\mu^{(n)}=(m(t_i))_iμ(n)=(m(ti​))i​ and Σ(n)=(K(ti,tj))ij\Sigma^{(n)}=(K(t_i,t_j))_{ij}Σ(n)=(K(ti​,tj​))ij​, which is symmetric positive semidefinite by hypothesis, and let Pt1,…,tnP_{t_1,\dots,t_n}Pt1​,…,tn​​ be the Gaussian law of mean μ(n)\mu^{(n)}μ(n) and covariance Σ(n)\Sigma^{(n)}Σ(n) from Proposition 6. This family is consistent, since a marginal of a Gaussian over a subset of coordinates is the Gaussian with the corresponding subvector mean and submatrix covariance, which is exactly PPP on the smaller index set, and permuting indices permutes the law correspondingly. The Kolmogorov extension theorem, which builds a measure on the product space RT\R^TRT from a consistent family of finite-dimensional laws by Caratheodory extension from the algebra of cylinder sets, then produces a process with these finite-dimensional distributions, Gaussian by construction. Uniqueness in law holds because the finite-dimensional distributions determine the law on the cylinder sigma-algebra, again an intersection-closed generating system.

The covariance function carries the entire second-order structure, and when the index set is an interval and the kernel is continuous, it is exactly the Mercer kernel of the previous track. The Karhunen-Loeve expansion applies the spectral decomposition of that kernel to the process itself, writing Xt=m(t)+∑nλn ξnφn(t)X_t=m(t)+\sum_n\sqrt{\lambda_n}\,\xi_n \varphi_n(t)Xt​=m(t)+∑n​λn​​ξn​φn​(t) with independent standard Gaussian coefficients ξn\xi_nξn​, where the eigenvalues are the variances of the modes and the eigenfunctions their shapes. Brownian motion is the Gaussian process with m=0m=0m=0 and K(s,t)=min⁡(s,t)K(s,t)=\min(s,t)K(s,t)=min(s,t), and its construction, the foundation of stochastic calculus, is the next step.

[1]
R. Durrett, Probability: Theory and Examples, 5th ed. in Cambridge Series in Statistical and Probabilistic Mathematics. Cambridge University Press, 2019.
[2]
O. Kallenberg, Foundations of Modern Probability, 3rd ed. Springer, 2021.

Part 4 of 9 in Probability

← previousCharacteristic Functionsnext →Convergence and Limit Theorems

Explore connections

see in the atlas →

related

  • The Karhunen-Loeve Expansion
  • The Ornstein-Uhlenbeck Process
  • Independence

referenced by (6)

  • Eigenvalues and the Spectral Theorem
  • Positive Definite Matrices
  • Quadratic Variation
  • Second-Order Processes and Mean-Square Calculus
  • The Black-Scholes Equation
  • The Construction of Brownian Motion
cite
@misc{gaussian-vectors-and-processes,
  author = {Zac Kienzle},
  title  = {Gaussian Vectors and Processes},
  year   = {2026},
  month  = {06},
  url    = {https://zackienzle.com/blog/gaussian-vectors-and-processes}
}