Complex random variable

From HandWiki

In probability theory and statistics, complex random variables are a generalization of real-valued random variables to complex numbers, i.e. the possible values a complex random variable may take are complex numbers.[1] Complex random variables can always be considered as pairs of real random variables: their real and imaginary parts. Therefore, the distribution of one complex random variable may be interpreted as the joint distribution of two real random variables. Some concepts of real random variables have a straightforward generalization to complex random variables—e.g., the definition of the mean of a complex random variable. Other concepts are unique to complex random variables.

Applications of complex random variables are found in digital signal processing,[2] quadrature amplitude modulation and information theory.

Definition

A complex random variable [math]\displaystyle{ Z }[/math] on the probability space [math]\displaystyle{ (\Omega,\mathcal{F},P) }[/math] is a function [math]\displaystyle{ Z \colon \Omega \rightarrow \mathbb{C} }[/math] such that both its real part [math]\displaystyle{ \Re{(Z)} }[/math] and its imaginary part [math]\displaystyle{ \Im{(Z)} }[/math] are real random variables on [math]\displaystyle{ (\Omega,\mathcal{F},P) }[/math].

Examples

Simple example

Consider a random variable that may take only the three complex values [math]\displaystyle{ 1+i , 1-i , 2 }[/math] with probabilities as specified in the table. This is a simple example of a complex random variable.

Probability [math]\displaystyle{ P(z) }[/math] Value [math]\displaystyle{ z }[/math]
[math]\displaystyle{ \frac{1}{4} }[/math] [math]\displaystyle{ 1+i }[/math]
[math]\displaystyle{ \frac{1}{4} }[/math] [math]\displaystyle{ 1- i }[/math]
[math]\displaystyle{ \frac{1}{2} }[/math] [math]\displaystyle{ 2 }[/math]

The expectation of this random variable may be simply calculated: [math]\displaystyle{ \operatorname{E}[Z] = \frac{1}{4}( 1+i ) + \frac{1}{4}( 1-i ) + \frac{1}{2}2 = \frac{3}{2} . }[/math]

Uniform distribution

Another example of a complex random variable is the uniform distribution over the filled unit circle, i.e. the set [math]\displaystyle{ \{ z \in \mathbb{C} \mid |z| \le 1 \} }[/math]. This random variable is an example of a complex random variable for which the probability density function is defined. The density function is shown as the yellow disk and dark blue base in the following figure.

Probability density function of a complex random variable shich is uniformly distributed inside the unit circle

Complex normal distribution

Main page: Complex normal distribution

Complex Gaussian random variables are often encountered in applications. They are a straightforward generalization of real Gaussian random variables. The following plot shows an example of the distribution of such a variable.

Probability density function of a complex Gaussian random variable

Cumulative distribution function

The generalization of the cumulative distribution function from real to complex random variables is not obvious because expressions of the form [math]\displaystyle{ P(Z \leq 1+3i) }[/math] make no sense. However expressions of the form [math]\displaystyle{ P(\Re{(Z)} \leq 1, \Im{(Z)} \leq 3) }[/math] make sense. Therefore, we define the cumulative distribution [math]\displaystyle{ F_Z : \mathbb{C} \to [0,1] }[/math] of a complex random variables via the joint distribution of their real and imaginary parts:

[math]\displaystyle{ F_Z(z)=F_{\Re{(Z)},\Im{(Z)}}(\Re{(z)},\Im{(z)})=P(\Re{(Z)} \leq \Re{(z)} , \Im{(Z)} \leq \Im{(z)}) }[/math]

 

 

 

 

(Eq.1)

Probability density function

The probability density function of a complex random variable is defined as [math]\displaystyle{ f_Z(z)=f_{\Re{(Z)},\Im{(Z)}}(\Re{(z)},\Im{(z)}) }[/math], i.e. the value of the density function at a point [math]\displaystyle{ z \in \mathbb{C} }[/math] is defined to be equal to the value of the joint density of the real and imaginary parts of the random variable evaluated at the point [math]\displaystyle{ (\Re{(z)},\Im{(z)}) }[/math].

An equivalent definition is given by [math]\displaystyle{ f_Z(z)=\frac{\partial^2 }{\partial x \partial y} P(\Re{(Z)} \leq x , \Im{(Z)} \leq y) }[/math] where [math]\displaystyle{ x=\Re{(z)} }[/math] and [math]\displaystyle{ y=\Im{(z)} }[/math].

As in the real case the density function may not exist.

Expectation

The expectation of a complex random variable is defined based on the definition of the expectation of a real random variable:[3]:p. 112

[math]\displaystyle{ \operatorname{E}[Z] = \operatorname{E}[\Re{(Z)}]+ i \operatorname{E}[\Im{(Z)}] }[/math]

 

 

 

 

(Eq.2)

Note that the expectation of a complex random variable does not exist if [math]\displaystyle{ \operatorname{E}[\Re{(Z)}] }[/math] or [math]\displaystyle{ \operatorname{E}[\Im{(Z)}] }[/math] does not exist.

If the complex random variable [math]\displaystyle{ Z }[/math] has a probability density function [math]\displaystyle{ f_Z(z) }[/math], then the expectation is given by [math]\displaystyle{ \operatorname{E}[Z] = \iint_{\mathbb{C}}z \cdot f_Z(z) \, dx \, dy }[/math].

If the complex random variable [math]\displaystyle{ Z }[/math] has a probability mass function [math]\displaystyle{ p_Z(z) }[/math], then the expectation is given by [math]\displaystyle{ \operatorname{E}[Z] = \sum_{z \in \mathbb{Z}}z \cdot p_Z(z) }[/math].

Properties

Whenever the expectation of a complex random variable exists, taking the expectation and complex conjugation commute:

[math]\displaystyle{ \overline{\operatorname{E}[Z]} = \operatorname{E}[\overline Z] . }[/math]

The expected value operator [math]\displaystyle{ \operatorname{E}[\cdot] }[/math] is linear in the sense that

[math]\displaystyle{ \operatorname{E}[aZ + bW] = a\operatorname{E}[Z] + b\operatorname{E}[W] }[/math]

for any complex coefficients [math]\displaystyle{ a, b }[/math] even if [math]\displaystyle{ Z }[/math] and [math]\displaystyle{ W }[/math] are not independent.

Variance and pseudo-variance

The variance is defined in terms of absolute squares as:[3]:{{{1}}}

[math]\displaystyle{ \operatorname{K}_{ZZ} = \operatorname{Var}[Z] = \operatorname{E}\left[\left|Z-\operatorname{E}[Z] \right|^2\right]=\operatorname{E}[|Z|^2]-\left|\operatorname{E}[Z]\right|^2 }[/math]

 

 

 

 

(Eq.3)

Properties

The variance is always a nonnegative real number. It is equal to the sum of the variances of the real and imaginary part of the complex random variable:

[math]\displaystyle{ \operatorname{Var}[Z]= \operatorname{Var}[\Re{(Z)}]+\operatorname{Var}[\Im{(Z)}] . }[/math]

The variance of a linear combination of complex random variables may be calculated using the following formula:

[math]\displaystyle{ \operatorname{Var}\left[\sum_{k=1}^N a_{k}Z_k \right] = \sum_{i=1}^N \sum_{j=1}^N a_{i}\overline{a_{j}}\operatorname{Cov}[Z_i,Z_j]. }[/math]

Pseudo-variance

The pseudo-variance is a special case of the pseudo-covariance and is defined in terms of ordinary complex squares, given by:

[math]\displaystyle{ \operatorname{J}_{ZZ}=\operatorname{E}[(Z-\operatorname{E}[Z])^2] = \operatorname{E}[Z^2]-(\operatorname{E}[Z])^2 }[/math]

 

 

 

 

(Eq.4)

Unlike the variance of [math]\displaystyle{ Z }[/math], which is always real and positive, the pseudo-variance of [math]\displaystyle{ Z }[/math] is in general complex.

Covariance matrix of real and imaginary parts

For a general complex random variable, the pair [math]\displaystyle{ (\Re{(Z)},\Im{(Z)}) }[/math] has a covariance matrix of the form:

[math]\displaystyle{ \begin{bmatrix} \operatorname{Var}[\Re{(Z)}] & \operatorname{Cov}[\Im{(Z)},\Re{(Z)}] \\ \operatorname{Cov}[\Re{(Z)},\Im{(Z)}] & \operatorname{Var}[\Im{(Z)}] \end{bmatrix} }[/math]

The matrix is symmetric, so [math]\displaystyle{ \operatorname{Cov}[\Re{(Z)},\Im{(Z)}] = \operatorname{Cov}[\Im{(Z)},\Re{(Z)}] }[/math]

Its elements equal:

[math]\displaystyle{ \begin{align} & \operatorname{Var}[\Re{(Z)}] = \tfrac{1}{2}\operatorname{Re}(\operatorname{K}_{ZZ} + \operatorname{J}_{ZZ}) \\ & \operatorname{Var}[\Im{(Z)}] = \tfrac{1}{2}\operatorname{Re}(\operatorname{K}_{ZZ} - \operatorname{J}_{ZZ}) \\ & \operatorname{Cov}[\Re{(Z)},\Im{(Z)}] = \tfrac{1}{2}\operatorname{Im}(\operatorname{J}_{ZZ}) \\ \end{align} }[/math]

Conversely:

[math]\displaystyle{ \begin{align} & \operatorname{K}_{ZZ} = \operatorname{Var}[\Re{(Z)}] + \operatorname{Var}[\Im{(Z)}] \\ & \operatorname{J}_{ZZ} = \operatorname{Var}[\Re{(Z)}] - \operatorname{Var}[\Im{(Z)}] + i2 \operatorname{Cov}[\Re{(Z)},\Im{(Z)}] \end{align} }[/math]

Covariance and pseudo-covariance

The covariance between two complex random variables [math]\displaystyle{ Z,W }[/math] is defined as[3]:{{{1}}}

[math]\displaystyle{ \operatorname{K}_{ZW}=\operatorname{Cov}[Z,W]=\operatorname{E}[(Z-\operatorname{E}[Z])\overline{(W-\operatorname{E}[W])}] = \operatorname{E}[Z\overline{W}]-\operatorname{E}[Z]\operatorname{E}[\overline{W}] }[/math]

 

 

 

 

(Eq.5)

Notice the complex conjugation of the second factor in the definition.

In contrast to real random variables, we also define a pseudo-covariance (also called complementary variance):

[math]\displaystyle{ \operatorname{J}_{ZW}=\operatorname{Cov}[Z,\overline{W}]=\operatorname{E}[(Z-\operatorname{E}[Z])(W-\operatorname{E}[W])] = \operatorname{E}[ZW]-\operatorname{E}[Z]\operatorname{E}[W] }[/math]

 

 

 

 

(Eq.6)

The second order statistics are fully characterized by the covariance and the pseudo-covariance.

Properties

The covariance has the following properties:

  • [math]\displaystyle{ \operatorname{Cov}[Z,W]=\overline{\operatorname{Cov}[W,Z]} }[/math] (Conjugate symmetry)
  • [math]\displaystyle{ \operatorname{Cov}[\alpha Z,W]=\alpha\operatorname{Cov}[Z,W] }[/math] (Sesquilinearity)
  • [math]\displaystyle{ \operatorname{Cov}[Z,\alpha W]=\overline{\alpha}\operatorname{Cov}[Z,W] }[/math]
  • [math]\displaystyle{ \operatorname{Cov}[Z_1+Z_2,W]=\operatorname{Cov}[Z_1,W]+\operatorname{Cov}[Z_2,W] }[/math]
  • [math]\displaystyle{ \operatorname{Cov}[Z,W_1+W_2]=\operatorname{Cov}[Z,W_1]+\operatorname{Cov}[Z,W_2] }[/math]
  • [math]\displaystyle{ \operatorname{Cov}[Z,Z]={\operatorname{Var}[Z]} }[/math]
  • Uncorrelatedness: two complex random variables [math]\displaystyle{ Z }[/math] and [math]\displaystyle{ W }[/math] are called uncorrelated if [math]\displaystyle{ \operatorname{K}_{ZW}=\operatorname{J}_{ZW}=0 }[/math] (see also: uncorrelatedness (probability theory)).
  • Orthogonality: two complex random variables [math]\displaystyle{ Z }[/math] and [math]\displaystyle{ W }[/math] are called orthogonal if [math]\displaystyle{ \operatorname{E}[Z\overline{W}] = 0 }[/math].

Circular symmetry

Circular symmetry of complex random variables is a common assumption used in the field of wireless communication. A typical example of a circular symmetric complex random variable is the complex Gaussian random variable with zero mean and zero pseudo-covariance matrix.

A complex random variable [math]\displaystyle{ Z }[/math] is circularly symmetric if, for any deterministic [math]\displaystyle{ \phi \in [-\pi,\pi] }[/math], the distribution of [math]\displaystyle{ e^{\mathrm i \phi}Z }[/math] equals the distribution of [math]\displaystyle{ Z }[/math].

Properties

By definition, a circularly symmetric complex random variable has [math]\displaystyle{ \operatorname{E}[Z] = \operatorname{E}[e^{\mathrm i \phi} Z] = e^{\mathrm i \phi}\operatorname{E}[Z] }[/math] for any [math]\displaystyle{ \phi }[/math].

Thus the expectation of a circularly symmetric complex random variable can only be either zero or undefined.

Additionally, [math]\displaystyle{ \operatorname{E}[ZZ] = \operatorname{E}[e^{\mathrm i \phi} Z e^{\mathrm i \phi}Z] = e^{\mathrm 2i \phi} \operatorname{E}[ZZ] }[/math] for any [math]\displaystyle{ \phi }[/math].

Thus the pseudo-variance of a circularly symmetric complex random variable can only be zero.

If [math]\displaystyle{ Z }[/math] and [math]\displaystyle{ e^{\mathrm i\phi}Z }[/math] have the same distribution, the phase of [math]\displaystyle{ Z }[/math] must be uniformly distributed over [math]\displaystyle{ [-\pi,\pi] }[/math] and independent of the amplitude of [math]\displaystyle{ Z }[/math].[4]

Proper complex random variables

The concept of proper random variables is unique to complex random variables, and has no correspondent concept with real random variables.

A complex random variable [math]\displaystyle{ Z }[/math] is called proper if the following three conditions are all satisfied:

  • [math]\displaystyle{ \operatorname{E}[Z] = 0 }[/math]
  • [math]\displaystyle{ \operatorname{Var}[Z] \lt \infty }[/math]
  • [math]\displaystyle{ \operatorname{E}[Z^2] = 0 }[/math]

This definition is equivalent to the following conditions. This means that a complex random variable is proper if, and only if:

  • [math]\displaystyle{ \operatorname{E}[Z] = 0 }[/math]
  • [math]\displaystyle{ \operatorname{E}[\Re{(Z)}^2]=\operatorname{E}[\Im{(Z)}^2]\neq\infty }[/math]
  • [math]\displaystyle{ \operatorname{E}[\Re{(Z)}\Im{(Z)}]=0 }[/math]

Theorem — Every circularly symmetric complex random variable with finite variance is proper.

For a proper complex random variable, the covariance matrix of the pair [math]\displaystyle{ (\Re{(Z)},\Im{(Z)}) }[/math] has the following simple form:

[math]\displaystyle{ \begin{bmatrix} \frac{1}{2} \operatorname{Var}[Z] & 0 \\ 0 & \frac{1}{2} \operatorname{Var}[Z] \end{bmatrix} }[/math].

I.e.:

[math]\displaystyle{ \begin{align} & \operatorname{Var}[\Re{(Z)}] = \operatorname{Var}[\Im{(Z)}] = \tfrac{1}{2}\operatorname{Var}[Z] \\ & \operatorname{Cov}[\Re{(Z)},\Im{(Z)}] = 0 \\ \end{align} }[/math]

Cauchy-Schwarz inequality

The Cauchy-Schwarz inequality for complex random variables, which can be derived using the Triangle inequality and Hölder's inequality, is

[math]\displaystyle{ \left| \operatorname{E} \left[ Z\overline{W} \right] \right|^2 \leq \left| \operatorname{E} \left[ \left| Z\overline{W} \right| \right] \right|^2 \leq \operatorname{E} \left[ |Z|^2 \right] \operatorname{E}\left[ |W|^2 \right] }[/math].

Characteristic function

The characteristic function of a complex random variable is a function [math]\displaystyle{ \mathbb{C} \to \mathbb{C} }[/math] defined by

[math]\displaystyle{ \varphi_Z(\omega) = \operatorname{E} \left [ e^{i\Re{(\overline{\omega} Z)}} \right ] = \operatorname{E} \left [ e^{i( \Re{(\omega)}\Re{(Z)} + \Im{(\omega)}\Im{(Z)} )} \right ]. }[/math]

See also

References

  1. Eriksson, Jan; Ollila, Esa; Koivunen, Visa (2009). "Statistics for complex random variables revisited". 2009 IEEE International Conference on Acoustics, Speech and Signal Processing. Taipei, Taiwan: Institute of Electrical and Electronics Engineers. pp. 3565-3568. doi:10.1109/ICASSP.2009.4960396. 
  2. Lapidoth, A. (2009). A Foundation in Digital Communication. Cambridge University Press. ISBN 9780521193955. 
  3. 3.0 3.1 3.2 Park,Kun Il (2018). Fundamentals of Probability and Stochastic Processes with Applications to Communications. Springer. ISBN 978-3-319-68074-3. 
  4. Peter J. Schreier, Louis L. Scharf (2011). Statistical Signal Processing of Complex-Valued Data. Cambridge University Press. ISBN 9780511815911.