Let
X
i
X
i
(
i=1
i
1
to NN) be the
NN independent random
processes, each will zero mean and variance
σ2
σ
2
, which are combined to give
X=1N∑i=1N
Xi
X
1
N
i
1
N
Xi
(9)
Then, if the characteristic function of each input process
before scaling is
Φu
Φ
u
and we use
Equation 5 to
include the scaling by
1N
1
N
, the characteristic function of
XX is
ΦX
u=∏i=1N
Φ
X
i
uN=ΦNuN
ΦX
u
i
1
N
Φ
X
i
u
N
Φ
u
N
N
(10)
Taking logs:
log
ΦX
u=NlogΦuN
ΦX
u
N
Φ
u
N
(11)
Using Taylor's theorem to expand
ΦuN
Φ
u
N
in terms of its derivatives at
u=0
u
0
(and hence its moments) gives
ΦuN=Φ0+uNdΦ0d+12uN2dΦ0d+16uN3dΦ0d+124uN4dΦ0d+…
Φ
u
N
Φ
0
u
N
Φ
0
1
2
u
N
2
2
Φ
0
1
6
u
N
3
3
Φ
0
1
24
u
N
4
4
Φ
0
…
(12)
From the
Moments property of characteristic
functions with zero mean:
-
valid pdf
Φ0=E
Xi
0=1
Φ
0
Xi
0
1
-
zero mean
Φ′0=iE
Xi
=0
Φ
0
Xi
0
-
variance
Φ0′′=i2E
Xi
2=−σ2
2
Φ
0
2
Xi
2
σ
2
-
scaled skewness
Φ0′′′=i3E
Xi
3=−(iγσ3)
3
Φ
0
3
Xi
3
γ
σ
3
-
scaled kurtosis
Φ04=i4E
Xi
4=(κ+3)σ4
4
Φ
0
4
Xi
4
κ
3
σ
4
These are all constants, independent of
NN, and dependent only on the
shape of the pdfs
f
Xi
f
Xi
.
Substituting these moments into Equation 11 and Equation 12 and
using the series expansion,
log(1+x)=x
1
x
x
+ (terms of order
x2
x
2
or smaller), gives
log
ΦX
u=NlogΦuN=Nlog(1−u22Nσ2+**)=N(−u2σ22N+**)=−u2σ22+##
ΦX
u
N
Φ
u
N
N
1
u
2
2
N
σ
2
**
N
u
2
σ
2
2
N
**
u
2
σ
2
2
##
(13)
where ** represents the terms of order
N−32
N
3
2
or smaller and ## represents the terms of order
N−12
N
1
2
or smaller. As
N→∞
N
,
log
ΦX
u→−u2σ22
ΦX
u
u
2
σ
2
2
Therefore, as
N→∞
N
ΦX
u→e−u2σ22
ΦX
u
u
2
σ
2
2
(14)
Note that, if the input pdfs are symmetric, the skewness will
be zero and the error terms will decay as
N-1
N
-1
rather than
N−12
N
1
2
; and so convergence to a Gaussian characteristic
function will be more rapid.
Hence we may now infer from Equation 6, Equation 7 and Equation 14 that the pdf of
XX as
N→∞
N
will be given by
fX
x=12πσ2e−x22σ2
fX
x
1
2
σ
2
x
2
2
σ
2
(15)
Thus we have proved the required central limit
result.