XN Xnim I K N M: The Mean Value Is
XN Xnim I K N M: The Mean Value Is
It involves three
steps. (1) The N - point sequence is subdivided in to K nonoverlapping segments. Here each
segment has length M. This results in the K data segments.
1
Pxx( i ) ( f )
M
x [n]e
2
j 2 fn
n 0
i = 0,1,...,K 1
We can compute the Bartlett power spectrum estimate by averaging periodogram computed for
each segment as,
1
P (f)
K
B
xx
K 1
(i )
xx
i 0
(f)
1
E[ P ( f )]
K
B
xx
K 1
E[ P
i 0
(i )
xx
( f )]
E[ Pxx(i ) ( f )]
E[ P ( f )]
(i )
xx
M 1
m ( M 1)
m
j 2 fm
xx ( m)e
M
12
sin ( f ) M
)
1 2 xx sin ( f ) d
1 sin fM
WB ( f )
M sin f
WB ( f ) 1
= 0,
m
M
, for m M 1
otherwise
1 K 1
var[P ( f )] 2 var[Pxx(i ) ( f )]
K i 0
B
xx
1
var[Pxx( i ) ( f )]
K
sin 2 fM
1
var[Pxx( B ) ( f )] 2xx ( f ) 1
K
M
sin
2
Welch method.
xi [n] x[n iD ]
n 0,1,..., M 1
i 0,1,..., L 1
1
(i )
P%
(
f
)
xx
MU
1
U
M
(W )
xx
M 1
M 1
x [n]w[n]e
n 0
2
j 2 fn
i = 0,1,...,L 1
w [ n]
2
n 0
1 L1 %
( f ) Pxx( i ) ( f )
L i 0
(W )
xx
E[ P
1 L1 %
( f )] E[ Pxx( i ) ( f )]
L i 0
B & T method
P (f)
BT
xx
M 1
m ( M 1)
E[ P ( f )]
BT
xx
rxx [m]w[m]e j 2 fm
M 1
m ( M 1)
E[rxx [m]]w[m]e j 2 fm
M 1
m ( M 1)
xx [m]wB [m]w[m]e j 2 fm
1
var[ P ( f )] ( f )
N
BT
xx
2
xx
2
w
[
m
]
m ( M 1)
M 1