Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

#### Normal Samples

If we want 10 samples from a Gaussian or normal random process with variance 4 can use rnorm(10,sd=2). Remember the standard deviation (sd) is the square root of the variance.

`set.seed(2)x <- rnorm(10,sd=2)x`
`##   -1.7938291  0.3696984  3.1756907 -2.2607513 -0.1605035  0.2648406##    1.4159095 -0.4793960  3.9689479 -0.2775740`
`var(x)`
`##  3.880803`

The var() function produces an estimate of the variance, if we want a better estimate we need more samples.

`var(rnorm(1000,sd=2))`
`##  4.105966`

#### Complex Normal Samples

If we are using base R and want complex normal (CN) samples, we need to write our own function. When the signal processing literature refers to CN they are usually referring to circularly-symmetric CN. Circularly-symmetric means the samples are independent and their mean is 0.

The function produces N CN samples with variance v. The real and imaginary parts are independent, because they are produced by different calls to rnorm(). Let x,y be independent. The var(ax) = a^2var(x) and var(x+y)=var(x)+var(y). So, if we want a variance of 1 would have to start a variance of sqrt(1/2).

`makeCN <- function(N,v=1) {(sqrt(v/2))*rnorm(N) + (sqrt(v/2))*1i*rnorm(N)}makeCN(10)`
`##    0.0023376-0.2079938i  0.7613032+0.6053620i  0.3946671-0.4049715i##    0.4892950-0.1207824i  0.4651165-0.2871364i -0.2312504+0.9408834i##   -0.2153405-0.9648887i -1.0994866+1.0119199i  1.0396552+0.7824796i##   0.1147878+0.9059002i`

If we want to check the variance, we can’t use var() directly.

`var(makeCN(10))`
`## Warning in var(makeCN(10)): imaginary parts discarded in coercion`
`##  0.6039204`

But the real and imaginary parts are independent, so we can calculate the variance separately.

`z <- makeCN(10)var(Re(z)) + var(Im(z))`
`##  0.5623849`

To make this easier, we can create a function to find the variance.

`varComplex <- function(z) var(Re(z)) + var(Im(z))`

To get a good estimate we need a-lot of samples.

`varComplex(makeCN(1000))`
`##  1.016615`

Let’s set the variance to 2 and then estimate the variance of the samples.

`varComplex(makeCN(1000,v=2))`
`##  1.925119`

Success!