---
title: 'AMP-A-EC: AMP for Grant-Free OFDM Access'
url: https://www.emergentmind.com/topics/amp-a-ec
type: topic
---

# AMP-A-EC: AMP for Grant-Free OFDM Access

AMP-A-EC is an approximate-message-passing algorithm for OFDM-based massive grant-free access in wideband systems under frequency-selective fading. Introduced as **AMP-based Activity detection and Effective Channel estimation**, it is derived from an exact time-domain Bayesian model and a corresponding factor graph that preserves the Bernoulli-Gaussian coupling between device activity and channel coefficients. Its purpose is joint **MAP-based device activity detection** and **MMSE-based effective-channel estimation**, with an additional best-iterate tracking mechanism to mitigate the non-monotonic convergence behavior that can arise when pilot length is smaller than or comparable to the number of active devices [2508.06540].

## 1. Definition and problem setting

AMP-A-EC is formulated for a single-cell uplink in which a base station has \(M\) antennas and \(N\) single-antenna devices. The channel is block fading and frequency-selective with \(P\) taps. Device activity is binary, \(a_n\in\{0,1\}\), with sparse activity satisfying \(\sum_{n\in\mathcal N} a_n=N_a\ll N\). Large-scale fading coefficients \(\beta_n>0\) are assumed known at the base station, and the small-scale fading taps satisfy
\[
g_{n,p,m}\sim \mathcal{CN}(0,1),
\qquad
h_{n,p,m}=\sqrt{\beta_n}\,g_{n,p,m}.
\]
The effective channel variable is defined by
\[
x_{n,p,m}=a_n h_{n,p,m},
\]
so that \(x_{n,p,m}\) follows a Bernoulli-Gaussian law conditioned on \(a_n\) [2508.06540].

The pilot structure is OFDM-based. The number of subcarriers is \(K\), the pilot length is \(L=KQ\), and each device has a unique length-\(L\) pilot split across \(Q\) OFDM pilot symbols. After DFT-domain manipulation and stacking across OFDM symbols and receive antennas, the received signal is written in compact form as
\[
\mathbf Y=\mathbf A\mathbf X+\mathbf N,
\]
with \(\mathbf Y\in\mathbb C^{L\times M}\), \(\mathbf A\in\mathbb C^{L\times NP}\), and \(\mathbf X\in\mathbb C^{NP\times M}\). In this model, the block rows of \(\mathbf X\) encode the activity-weighted channel taps of each device. For large \(L\), the sensing matrix obeys
\[
\sum_{l\in\mathcal L}|A_{l,n,p}|^2=1,\qquad
A_{l,n,p}=\mathcal O(L^{-1/2}),\qquad
|A_{l,n,p}|^2=\mathcal O(L^{-1}),
\]
which is the asymptotic scaling used by the algorithmic derivation [2508.06540].

The paper places AMP-A-EC alongside a second algorithm, AMP-A-AC, but the two solve different estimation targets. AMP-A-EC is specifically the method for **activity detection and effective channel estimation**, whereas AMP-A-AC addresses activity detection and actual-channel estimation of active devices.

## 2. Bayesian formulation and factor-graph structure

The method is built from three Bayesian estimation problems. The first is **MAP activity detection**:
\[
\max_{a_n\in\{0,1\}} p(a_n\mid \mathbf Y),
\]
with decision rule
\[
\hat a_n^\star(\mathbf Y)=
\begin{cases}
0, & \theta_n<0,\\
1, & \theta_n\ge 0,
\end{cases}
\qquad
\theta_n=\log\frac{p(a_n=1\mid \mathbf Y)}{p(a_n=0\mid \mathbf Y)}.
\]

The second is **MMSE estimation of the effective channel**:
\[
\hat x_{n,p,m}^\star(\mathbf Y)=\int x_{n,p,m}\,p(x_{n,p,m}\mid \mathbf Y)\,dx_{n,p,m}.
\]

The third, introduced for comparison in the same framework, is **MMSE estimation of the actual channel of active devices**:
\[
\hat h_{n,p,m}^\star(\mathbf Y)=\int h_{n,p,m}\,p(h_{n,p,m}\mid \mathbf Y,a_n=1)\,dh_{n,p,m}.
\]
AMP-A-EC targets the first two of these [2508.06540].

The joint distribution factorizes as
\[
p(\mathbf X,\mathbf a,\mathbf Y)
=
\prod_{l\in\mathcal L}\prod_{m\in\mathcal M} p(y_{l,m}\mid \mathbf X_{:,m})
\prod_{n\in\mathcal N} p(a_n)
\prod_{p\in\mathcal P}\prod_{m\in\mathcal M} p(x_{n,p,m}\mid a_n),
\]
with Gaussian likelihood
\[
p(y_{l,m}\mid \mathbf X_{:,m})
=
f_{\mathcal{CN}}\!\left(y_{l,m};\mathbf A_{l,:}\mathbf X_{:,m},\sigma^2\right).
\]

This factorization induces a factor graph that explicitly captures the likelihood factors \(p(y_{l,m}\mid \mathbf X_{:,m})\), the Bernoulli priors \(p(a_n)\), and the Bernoulli-Gaussian coupling \(p(x_{n,p,m}\mid a_n)\). A key structural property is that, for fixed \(n\), all \(x_{n,p,m}\) share the same \(a_n\), so they are dependent unconditionally but independent conditioned on \(a_n\). Rows of \(\mathbf X\) are independent across \(n\). This graph is the basis for message passing toward \(p(a_n\mid\mathbf Y)\) and \(p(x_{n,p,m}\mid\mathbf Y)\) [2508.06540].

## 3. Algorithmic structure of AMP-A-EC

AMP-A-EC applies AMP-style Gaussian approximations to the factor-graph messages and then performs parallel updates of residuals, activity scores, effective-channel estimates, and variance surrogates. The paper introduces the denoising function
\[
\eta_n(x,y,z)
\triangleq
\frac{\beta_n x}{\left(1+\frac{(1-z)f_{\mathcal{CN}}(0;x,y)}{z f_{\mathcal{CN}}(0;x,y+\beta_n)}\right)(y+\beta_n)},
\]
together with its derivative \(\eta_n'(x,y,z)\). This function couples the current pseudo-observation, noise level, and activity probability proxy [2508.06540].

The AMP residual and residual-energy estimator are
\[
z^{(t)}_{f_{l,m}}
\triangleq
y_{l,m}-\sum_{n\in\mathcal N}\sum_{p\in\mathcal P}
A_{l,n,p}\hat x^{(t)}_{x_{n,p,m}\to f_{l,m}},
\qquad
\hat\tau_m^{(t)}
\triangleq
\frac1L\sum_{l\in\mathcal L}|z^{(t)}_{f_{l,m}}|^2.
\]

The activity-related soft factor is
\[
\lambda^{(t+1)}_{n,p,m}
=
\frac{1}{1+\frac{(1-\rho_n)}{\rho_n}
\prod_{k\in\mathcal M,\,k\ne m}\prod_{j\in\mathcal P,\,j\ne p}
\frac{1-\hat\rho^{(t)}_{n,j,k}}{\hat\rho^{(t)}_{n,j,k}}},
\]
where
\[
\hat\rho^{(t)}_{n,p,m}
=
\frac{1}{1+\frac{f_{\mathcal{CN}}(0;\hat r^{(t)}_{x_{n,p,m}},\hat\sigma^{(t)}_{x_{n,p,m}})}
{f_{\mathcal{CN}}(0;\hat r^{(t)}_{x_{n,p,m}},\hat\sigma^{(t)}_{x_{n,p,m}}+\beta_n)}}.
\]

Under the Gaussian-message approximation, the approximate posteriors are written as
\[
p_n^{(t)}(a_n)\propto p(a_n)\prod_{m\in\mathcal M}\prod_{p\in\mathcal P}\mu^{(t)}_{q_{n,p,m}\to a_n}(a_n),
\]
\[
p^{(t)}_{n,p,m}(x_{n,p,m})
\propto
\mu^{(t)}_{q_{n,p,m}\to x_{n,p,m}}(x_{n,p,m})
\mu^{(t)}_{x_{n,p,m}\to q_{n,p,m}}(x_{n,p,m}),
\]
and
\[
\tilde p^{(t)}_{n,p,m}(x_{n,p,m};a_n)
\propto
p(x_{n,p,m}\mid a_n)\,
\mu^{(t)}_{x_{n,p,m}\to q_{n,p,m}}(x_{n,p,m}).
\]
The associated estimators are
\[
\hat a_n^{(t)}=
\begin{cases}
0,& \theta_n^{(t)}<0,\\
1,& \theta_n^{(t)}\ge0,
\end{cases}
\qquad
\hat x^{(t)}_{n,p,m}
=
\int x_{n,p,m}p^{(t)}_{n,p,m}(x_{n,p,m})\,dx_{n,p,m},
\]
with
\[
\theta_n^{(t)}=\log\frac{p_n^{(t)}(a_n=1)}{p_n^{(t)}(a_n=0)}.
\]

For large \(L\) and \(N\), the algorithm further uses
\[
\tau^{(t)}_{l,m}\approx \hat\tau_m^{(t)},\qquad
\hat\gamma^{(t)}_{f_{l,m}\to x_{n,p,m}}\approx \hat\tau_m^{(t)},\qquad
\hat\sigma^{(t)}_{x_{n,p,m}\to f_{l,m}}\approx \hat\tau_m^{(t)}.
\]
This leads to a fully explicit AMP recursion, including
\[
\hat x^{(t+1)}_{n,p,m}
=
\eta_n\!\left(
\hat x^{(t)}_{n,p,m}+\sum_l A^*_{l,n,p}z^{(t)}_{f_{l,m}},
\hat\tau_m^{(t)},
\lambda^{(t+1)}_{n,p,m}
\right),
\]
and the Onsager-corrected residual update
\[
z^{(t+1)}_{f_{l,m}}
=
y_{l,m}
-
\sum_{n,p}A_{l,n,p}\hat x^{(t+1)}_{n,p,m}
+
\frac{1}{L}z^{(t)}_{f_{l,m}}
\sum_{n,p}
\eta_n'\!\left(
\hat x^{(t)}_{n,p,m}+\sum_l A^*_{l,n,p}z^{(t)}_{f_{l,m}},
\hat\tau_m^{(t)},
\lambda^{(t+1)}_{n,p,m}
\right).
\]
This places AMP-A-EC in the standard AMP family in the sense that it uses a residual update with an Onsager-style correction, although its denoising and activity-coupling terms are tailored to the Bernoulli-Gaussian wideband access model [2508.06540].

## 4. Best-iterate tracking and convergence behavior

A distinctive feature of AMP-A-EC is that it does not simply return the final iterate. The paper emphasizes that AMP may fail to converge monotonically when the pilot length \(L\) is smaller than or comparable to the number of active devices. In that regime, the error probability or MSE can decrease first and then increase. AMP-A-EC addresses this by tracking the iterate that minimizes the GROUP-LASSO objective
\[
f(\mathbf X)\triangleq \frac12\|\mathbf Y-\mathbf A\mathbf X\|_F^2+\sum_{i=1}^{NP}\|\mathbf X_{i,:}\|_2,
\]
and defining
\[
f^{(t)}_{\text{best}}
=
\min\{f(\hat{\mathbf X}^{(1)}),\ldots,f(\hat{\mathbf X}^{(t)})\}.
\]

The corresponding best iterate \(\hat{\mathbf X}_{\text{best}}^{(t)}\) and activity scores \(\boldsymbol\theta_{\text{best}}^{(t)}\) are used for the final output:
\[
\hat a_n^{\text{best},(t)}=
\begin{cases}
0,& \theta_{n,\text{best}}^{(t)}<0,\\
1,& \theta_{n,\text{best}}^{(t)}\ge0,
\end{cases}
\qquad
\hat h_{n,p,m}=\hat x_{n,p,m}^{\text{best},(t)} \text{ if active.}
\]

The paper attributes the improvement to three elements: retention of AMP-style low-complexity parallel updates, best-iterate tracking using the GROUP-LASSO objective, and outputting the best past estimate rather than the last one. Reported qualitative consequences are that error probability and MSE first decrease and then remain unchanged, and that the gain is especially visible at \(L=96\) and \(L=128\), where \(L\lesssim N\rho\) or only slightly exceeds it [2508.06540].

A common misconception is to treat AMP-A-EC as a generic name for any expectation-consistent or orthogonalized AMP method. In the current arXiv literature, AMP-A-EC is a specific algorithm for grant-free wideband access, whereas broader AMP papers on first-order cancellation, EP/GAMP equivalence, or rotationally invariant models describe related methodology but do not define this algorithm by name [1907.04235][1907.02247][2412.01574]. This suggests that the suffix “A-EC” should be read operationally within this communications context rather than as a universal AMP taxonomy.

## 5. State evolution and asymptotic characterization

The paper analyzes AMP-A-EC through a scalar state-evolution formalism. The effective observation model is
\[
\mathbf r_t=\mathbf x+\tau^{(t)}\mathbf v,
\]
where \(\mathbf x\sim f_{\mathcal{B-CN}}(\mathbf x;\mathbf 0,\beta\mathbf I,\rho)\), \(\mathbf v\sim\mathcal{CN}(\mathbf 0,\mathbf I)\), and \(\mathbf x\) and \(\mathbf v\) are independent. The recursion is
\[
\tau^{(t+1)}
=
\sigma^2+\frac{NP}{L}
\left(
\frac{\rho\beta \tau^{(t)}}{\beta+\tau^{(t)}}+\phi_t(\tau^{(t)})
\right),
\qquad
\tau^{(0)}=\sigma^2+\frac{NP}{L}\rho\beta.
\]

The term \(\phi_t(\tau^{(t)})\) is
\[
\phi_t(\tau^{(t)})
=
\frac{1}{MP}
\mathbb E_{\mathbf r_t}\!\left[
g(\mathbf r_t,\tau^{(t)}\mathbf I)
\left(1-g(\mathbf r_t,\tau^{(t)}\mathbf I)\right)
\frac{\beta^2}{(\beta+\tau^{(t)})^2}
\mathbf r_t^H\mathbf r_t
\right],
\]
with
\[
g(\mathbf r_t,\tau^{(t)}\mathbf I)
=
\frac{1}{1+\frac{(1-\rho)f_{\mathcal{CN}}(0;\mathbf r_t,\tau^{(t)}\mathbf I)}
{\rho f_{\mathcal{CN}}(0;\mathbf r_t,(\tau^{(t)}+\beta)\mathbf I)}}.
\]

For activity detection, the error probability at iteration \(t\) is
\[
p_n^{\mathrm e,(t)}=\Pr(\hat a_n^{(t)}\ne a_n).
\]
The paper gives a closed form involving incomplete Gamma functions,
\[
p_n^{\mathrm e,(t)}
=
\rho_n \frac{\underline{\Gamma}(PM,\cdot)}{\Gamma(PM)}
+
(1-\rho_n)\frac{\bar{\Gamma}(PM,\cdot)}{\Gamma(PM)},
\]
where the thresholds are explicit functions of \(\tau^{(t)},\beta_n,\rho_n,P,M\). The qualitative consequences stated in the paper are that error probability decreases with \(M\) and \(P\), and that \(p_n^{\mathrm e,(t)}\to 0\) as \(M\to\infty\) [2508.06540].

For active-device channel estimation, the normalized MSE
\[
\frac1{PM}\mathbb E\!\left[\|\mathbf h_n-\hat{\mathbf h}_n^{(t)}\|_2^2\right]
\]
is characterized analytically as a function of \(\tau^{(t)}\), \(\beta_n\), and \(g(\cdot)\). The paper states that this MSE decreases as state evolution improves. It also reports that the analytical state-evolution predictions for AMP-A-EC match simulation closely.

## 6. Computational complexity, numerical results, and positioning

The dominant per-iteration flop counts for AMP-A-EC are reported stepwise. The dominant term is
\[
2LMNP + 2LM^2 + 8NPM^2,
\]
yielding overall complexity
\[
\mathcal O(LNPM).
\]
The paper states that AMP-A-EC is more complex than AMP-A-AC because AMP-A-EC computes \(\eta_n(\cdot)\), \(\eta_n'(\cdot)\), and the more detailed \(\lambda_{n,p,m}\) for every \((n,p,m)\). Comparative orders given in the paper are: ML-MMSE has \(\mathcal O(L^2NP+L^3)\); AMP-FS and OMP-ext. have \(\mathcal O(LNKM)\); and AMP-FL-ext., AMP-A-EC, and AMP-A-AC have \(\mathcal O(LNPM)\) [2508.06540].

The simulation setup includes \(K=32\), subcarrier spacing \(30\) kHz, noise PSD \(-174\) dBm/Hz with \(\sigma^2=10^{-11.418}\), activity probability \(\rho=0.1\), and typical defaults
\[
P_t=10\text{ dBm},\quad L=128,\quad M=64,\quad N=1000,\quad P=3.
\]
Pilots are i.i.d. \(\mathcal{CN}(0,1)\) and normalized so that \(\sum_q\|\tilde{\mathbf s}_{q,n}\|_2^2=K\). The experiments use 500 Monte Carlo realizations.

The paper reports that AMP-A-EC and AMP-A-AC significantly outperform AMP-FL-ext, AMP-FS, and OMP-ext, with up to **94%** reduction in error probability, up to **33%** reduction in MSE, and up to **96%** reduction in computation time versus ML-MMSE. It also reports several monotonic trends: increasing \(L\) improves error probability, false alarm, missed detection, and MSE; increasing \(M\) improves all estimation metrics; increasing \(P\) or \(N\) worsens performance; increasing SNR improves performance; and computation time grows with \(L\), \(M\), and \(N\) but is essentially independent of SNR [2508.06540].

The paper distinguishes the preferable operating regions of the two proposed algorithms. **AMP-A-EC** is preferable at **small \(L\)** and **small \(M\)**. **AMP-A-AC** is preferable at **large \(L\)** and **large \(M\)**, and has lower dominant complexity. AMP-A-EC is described as somewhat more accurate in some small-dimension regimes, while AMP-A-AC can be faster and better in larger regimes. A plausible implication is that the two algorithms are best viewed as complementary solvers within a shared Bayesian and factor-graph framework rather than as strict substitutes.

## 7. Relation to the broader AMP literature

AMP-A-EC belongs to the family of approximate-message-passing algorithms in that it uses Gaussian approximations, scalar denoising structure, and an Onsager-corrected residual update. In the broader AMP literature, standard linear-regression AMP is often motivated by the Onsager term’s cancellation of leading self-interaction, yielding an asymptotically Gaussian effective noise characterized by state evolution [1907.04235]. EP-based derivations likewise connect AMP and GAMP to a unified message-passing rule under large-system approximations, especially for AWGN measurement channels [1907.02247].

However, AMP-A-EC is not a generic statement about AMP theory. It is a communications-specific construction for **OFDM-based grant-free access under frequency-selective fading**, based on an exact time-domain signal model, a factor graph that couples activity and channel taps, and a best-iterate strategy added to address practical non-monotonicity. By contrast, recent rotationally invariant AMP frameworks derive long-memory Onsager terms from free cumulants and OAMP reductions, but do not define an algorithm literally named AMP-A-EC [2412.01574]. The nomenclature therefore identifies a particular algorithmic instance rather than a general AMP subclass.

Within that narrower meaning, AMP-A-EC’s significance lies in combining four elements in a single procedure: exact wideband modeling, MAP activity detection, MMSE effective-channel estimation, and best-iterate tracking under challenging pilot-limited regimes. The paper’s central conclusion is that this combination yields strong accuracy-complexity tradeoffs for massive grant-free access, especially when pilot length is limited and conventional AMP behavior is least reliable [2508.06540].

Source: https://www.emergentmind.com/topics/amp-a-ec