Key focus: Discuss scalar parameter estimation using CRLB. Estimate DC component from observed data in the presence of AWGN noise.
Consider a set of observed data samples
If the of PDF
As seen in the previous section, the curvature of the likelihood function (Fisher Information) is related to the concentration of PDF. More the curvature, more is the concentration of PDF, more will be accuracy of estimates. The Fisher Information is calculated from log likelihood function as,
Under the regularity condition that the score of the log likelihood function is zero,
The inverse of the Fisher Information gives the Cramér-Rao Lower Bound (CRLB).
Theoretical method to find CRLB:
1) Given a model for observed data samples –
2) Keep
3) If the result depends on
4) If the result depends on
5) Take the reciprocal of the result and negate it.
Let’s see an example for scalar parameter estimation using CRLB.
Derivation of CRLB for an embedded DC component in AWGN Noise:
Here
Given the fact that the samples are influenced by the AWGN noise with zero mean and variance=
The log likelihood function is formed as,
Taking the first partial derivative of log likelihood function with respect to A,
Computing the second partial derivative of log likelihood function by differentiating one more time,
The Fisher Information is given by taking the expectation and negating it.
The Cramér-Rao Lower Bound is the reciprocal of Fisher Information I(A)
The variance of any estimator that estimates the DC component
Tweaking the CRLB:
Now that we have found an expression for CRLB for the estimation of the DC component, we can look for schemes that may affect the CRLB. From the expression of CRLB, following points can be inferred.
1) The CRLB does not depend on the parameter to be estimated (
2) The CRLB increases linearly with
3) The CRLB decreases inversely with
Hi, Please help me how simulate it with matlab