The Stochastic Crb For Array Processing A Textbook Derivation -

[ \mathbfx(t) \sim \mathcalCN(\mathbf0, \mathbfR) ] [ \mathbfR(\boldsymbol\theta, \mathbfp, \sigma^2) = \mathbfA(\boldsymbol\theta) \mathbfP \mathbfA^H(\boldsymbol\theta) + \sigma^2 \mathbfI ]

[ [\mathbfF(\boldsymbol\eta)]_ij = N \cdot \textTr\left( \mathbfR^-1 \frac\partial \mathbfR\partial \eta_i \mathbfR^-1 \frac\partial \mathbfR\partial \eta_j \right) ]

(from Slepian–Bangs formula): The log-likelihood (ignoring constants) is: [ L = -N \log \det \mathbfR - \sum_t=1^N \mathbfx^H(t) \mathbfR^-1 \mathbfx(t) ] Taking derivatives and expectations yields the above trace formula. 3. Partitioning the Unknown Parameters Let: [ \boldsymbol\eta = [\boldsymbol\theta^T, \ \mathbfp^T, \ \sigma^2]^T ] We want the CRB for ( \boldsymbol\theta ), i.e., the top-left ( d \times d ) block of ( \mathbfF^-1 ). This guide focuses on the derivation — showing

This guide focuses on the derivation — showing the logical steps, assumptions, and mathematical manipulations required to arrive at the closed-form expression for the CRB when signals are modeled as stochastic (Gaussian) processes. We consider an array of ( M ) sensors receiving ( d ) narrowband signals from far-field sources. 1.1 Data Model (Stochastic Assumption) The ( M \times 1 ) snapshot vector at time ( t ) is:

where ( \boldsymbol\eta ) is the real parameter vector. Define the FIM as: [ \mathbfF = \beginbmatrix

Define the FIM as: [ \mathbfF = \beginbmatrix \mathbfF \theta\theta & \mathbfF \theta p & \mathbfF \theta \sigma^2 \ \mathbfF p\theta & \mathbfF pp & \mathbfF p\sigma^2 \ \mathbfF \sigma^2\theta & \mathbfF \sigma^2 p & \mathbfF_\sigma^2\sigma^2 \endbmatrix ]

[ \textCRB(\boldsymbol\theta) = \frac\sigma^22N \left[ \Re \left( \mathbfD^H \mathbf\Pi_A^\perp \mathbfD \odot \mathbfP^T \right) \right]^-1 ] [ \mathbfx(t) \sim \mathcalCN(\mathbf0

The CRB for ( \boldsymbol\theta ) (with nuisance parameters ( \mathbfp, \sigma^2 )) is: [ \textCRB(\boldsymbol\theta) = \left( \mathbfF \theta\theta - [\mathbfF \theta p \ \mathbfF \theta \sigma^2] \beginbmatrix \mathbfF pp & \mathbfF p\sigma^2 \ \mathbfF \sigma^2 p & \mathbfF \sigma^2\sigma^2 \endbmatrix^-1 \beginbmatrix \mathbfF p\theta \ \mathbfF_\sigma^2\theta \endbmatrix \right)^-1 ]

This is the Schur complement of the nuisance parameter block. Let ( \mathbf\Pi_A^\perp = \mathbfI - \mathbfA(\mathbfA^H\mathbfA)^-1\mathbfA^H ) (projector onto noise subspace). 4.1 Derivative w.r.t. ( \theta_k ) [ \frac\partial \mathbfR\partial \theta_k = \mathbfA_k' \mathbfP \mathbfA^H + \mathbfA \mathbfP (\mathbfA_k')^H ] where ( \mathbfA_k' = \frac\partial \mathbfA\partial \theta_k = [\mathbf0, \dots, \mathbfa'(\theta_k), \dots, \mathbf0] ) (derivative of the ( k )-th column). 4.2 Derivative w.r.t. ( p_k ) [ \frac\partial \mathbfR\partial p_k = \mathbfa(\theta_k) \mathbfa^H(\theta_k) ] (because ( \mathbfP ) is diagonal). 4.3 Derivative w.r.t. ( \sigma^2 ) [ \frac\partial \mathbfR\partial \sigma^2 = \mathbfI ] 5. Simplifying the FIM Blocks We use: ( \mathbfR^-1 = \sigma^-2 \mathbf\Pi_A^\perp ) only if ( \mathbfA ) is full rank and ( \mathbfP ) nonsingular? Actually, via Woodbury: [ \mathbfR^-1 = \sigma^-2 \left( \mathbfI - \mathbfA (\mathbfA^H \mathbfA + \sigma^2 \mathbfP^-1)^-1 \mathbfA^H \right) ] But for CRB derivations, a cleaner way: define ( \mathbfR^-1/2 \mathbfA = \mathbfU \Sigma \mathbfV^H ) etc. However, the standard result uses:

Let ( \mathbfB = \mathbfA \mathbfP^1/2 ). Then ( \mathbfR = \mathbfB \mathbfB^H + \sigma^2 \mathbfI ). The projection matrix onto the column space of ( \mathbfB ): [ \mathbfP_B = \mathbfB(\mathbfB^H \mathbfB)^-1 \mathbfB^H ] but ( \mathbfB^H \mathbfB = \mathbfP^1/2 \mathbfA^H \mathbfA \mathbfP^1/2 ).