Start with \( \hat{P}_k(0) = P(0) \).

Real and non-negative coefficients \( c_{\ell,k} \) and \( a_{\ell,k} \) are determined.

For each time instant \( n \geq 1 \), perform:

For each cluster \( k = 1, 2, ..., K \), perform:

1. Power estimation

\[ \hat{\psi}_k(n+1) = \hat{P}_k(n) + \mu_k \sum_{\ell \in \mathcal{N}_k} c_{\ell,k} \left| x_\ell(n) \right|^2 - \hat{P}_k(n) \]

\[ \hat{P}_k(n+1) = \sum_{\ell \in \mathcal{N}_k} a_{\ell,k} \hat{\psi}_\ell(n+1) \]

2. Decision making for detection

\[ H_0 : \hat{P}_k(n+1) < \gamma \quad \text{or} \quad H_1 : \hat{P}_k(n+1) > \gamma \]

Refer to equation (32) in [42] for selecting the threshold value.

End