This preview has intentionally blurred sections. Sign up to view the full version.
View Full DocumentThis preview has intentionally blurred sections. Sign up to view the full version.
View Full DocumentThis preview has intentionally blurred sections. Sign up to view the full version.
View Full DocumentThis preview has intentionally blurred sections. Sign up to view the full version.
View Full Document
Unformatted text preview: 1 Chapter 12 Spectrum Estimation Teacher: SinHorng Chen Office: Engineering Bld. #4, Room 805 Tel: ext. 31822 Email: [email protected] 2 12.1 Ergodicity It is needed to estimate the statistical parameters of a random process from real data. Although most parameters to be estimated can be formulated as the expectation of a function of the random process ( ) t X , the problem lies in the fact that we need to implement the expectation via taking average (ensemble average) over the whole sample space. An example is given in the following. For a specific t , ( ) t X is a random variable. Let { } ( ) ( ) t E t η = X . We can estimate the mean of ( ) t X using n observed samples ( , ) i t ξ X of ( ) t X by 1 ˆ ( ) ( , ) i i t t n η ξ = ∑ X . This is a point estimation. It is a consistent estimate of ( ) t η if its variance → as n → ∞ . In this case, it is a good estimate for large n . 3 But, we usually only have a sample of ( ) t X . Can we use it to estimate { } ( ) E t X ? If { } ( ) E t X depends on t , then this is impossible. If ( ) t X is a regular stationary process, then the time average will approach to { } ( ) E t X as the length → ∞ . This is the meaning of ergodicity : time average → ensemble average Meanergodic process Given a real stationary process ( ) t X , we want to estimate { } ( ) ( ) t E t η = X . We form the time average 1 ( ) 2 T T T t dt T = ∫ η X . 4 Since { } { } 1 ( ) 2 T T T E E t dt T η = = ∫ η X , T η is an unbiased estimator of η . If the variance of T η , 2 T σ → as T →∞ , then T η → η in the MS sense. We say that ( ) t X is meanergodic, i.e., ( ) t X is meanergodic if the time average T η → the ensemble average η , as T →∞ . So, “meanregodic” can be proved via calculating T σ and showing T σ → as T →∞ . 5 Ex. 12.1 Let c be an RV with mean c η and let ( ) t = X c . Then { } {} ( ) c E t E η η = = = X c . A specific sample ( ) t c = X is a line (constant). In this case, ( ) ( ) T ξ ξ = η c (i.e., a specific constant determined by ξ ). Obviously, ( ) T ξ η ≠ η . So ( ) t X is not meanergodic. 6 Ex. 12.2 Given two meanergodic processes 1 ( ) t X and 2 ( ) t X , their means are 1 η and 2 η , respectively. We form 1 2 ( ) ( ) ( ) t t t = + X X cX , where c is a RV independent of 2 ( ) t X and = c and 1 with probabilities 0.5, respectively. Then { } { } { } { } { } 1 2 1 2 1 2 ( ) ( ) ( ) ( ) 0.5 E t E t E t E E t η η η = + = + = + X X cX c X But for a specific sample, if ( ) ξ = c , then 1 ( ) ( ) t t = X X and 1 T η η → as T → ∞ . If ( ) 1 ξ = c , then 1 2 ( ) ( ) ( ) t t t = + X X X . In this case, 1 2 T η η = + η . Hence, ( ) t X is not meanergodic. 7 To calculate 2 T σ of T η for an RP ( ) t X , we can form 1 ( ) ( ) 2 t T t T t d T α α + = ∫ w X which is a moving average of ( ) t X . Obviously, (0) T = η w ....
View
Full Document
 Fall '08
 SinHorngChen
 Variance, spectral density, Estimation theory, Autoregressive model, n Levinson

Click to edit the document details