Computing on Dirty Paper: Interference-Free Integrated Communication and Computing

Abstract

Inspired by Costa’s pioneering work on dirty paper coding (DPC), this paper proposes a novel scheme for integrated communication and computing (ICC), named Computing on Dirty Paper, whereby the transmission of discrete data symbols for communication, and over-the-air computation (AirComp) of nomographic functions can be achieved simultaneously over common multiple-access channels. In particular, the proposed scheme allows for the integration of communication and computation in a manner that is asymptotically interference-free, by pre-canceling the computing symbols at the transmitters (TXs) using DPC principles. A simulation-based assessment of the proposed ICC scheme under a single-input multiple-output (SIMO) setup is also offered, including the evaluation of performance for data detection, and of mean-squared-error (MSE) performance for function computation, over a block of symbols. The results validate the proposed method and demonstrate its ability to significantly outperform state-of-the-art (SotA) ICC schemes in terms of both bit error rate (BER) and MSE.

I Introduction

The rapid proliferation of Internet of Things (IoT) devices and edge computing applications is driving the need for seamless integration of communication and computation in wireless networks, forming the basis for the integrated communication and computing (ICC) paradigm [1, 2]. This approach aims to reduce latency and resource overhead by enabling simultaneous data transmission and the over-the-air computation (AirComp) [3, 4, 5, 6] of nomographic functions over shared channels, particularly in multiple-access scenarios where interference poses significant challenges [7].

Traditional methods treat the communication and computation operations separately, leading to inefficiencies, but recent advancements in ICC have shown that joint designs can significantly enhance performance [8]. However, when discrete data symbols must be communicated alongside computing signals, mutual interference between the two tasks degrades performance, necessitating innovative interference management techniques [9].

Dirty paper coding (DPC), a foundational concept in information theory first proposed in [10], addresses similar issues by allowing a transmitter to communicate reliably in the presence of known interference through precoding, effectively “writing on dirty paper” to achieve the same channel capacity as if the interference were absent.

While many ways to implement DPC have been proposed including approaches leveraging Tomlinson-Harashima precoding [11], nested lattice codes have emerged as a particularly effective method due to their structured nature and ability to approach the theoretical limits of DPC [12, 13, 14, 15, 16]. In addition, DPC has also been successfully applied in various multifunctional communication system aspects, such as in the integrated sensing and communications (ISAC) paradigm [17, 18, 19].

Inspired by these developments, we introduce in this paper Computing on Dirty Paper, a novel ICC approach that leverages DPC to pre-cancel computing symbols at the users in a single-input multiple-output (SIMO) uplink setting, enabling the simultaneous transmission of discrete data and the computation of nomographic functions over the same multiple access channel. In particular, we demonstrate that typical DPC schemes such as the one proposed in [13] can be incorporated at the transmitter (TX), and combined with a classic linear minimum mean square error (LMMSE) estimator at the receiver (RX), to enable ICC with no impact of the AirComp operation over the bit error rate (BER) of the communication stream. A simulation-based assessment of the method in a SIMO setup is performed, which validates our approach, showing substantial improvements over state-of-the-art (SotA) benchmarks in terms of both functionalities.

Notation: Throughout the manuscript, vectors and matrices are represented by lowercase and uppercase boldface letters, respectively; $\mathbf{I}_{M}$ denotes an identity matrix of size $M$ and $\mathbf{1}_{M}$ denotes a column vector composed of $M$ ones; the Euclidean norm and the absolute value of a scalar are respectively given by $\|\cdot\|_{2}$ and $|\cdot|$; the transpose and hermitian operations follow the conventional form $(\cdot)^{\mathsf{T}}$ and $(\cdot)^{\mathsf{H}}$, respectively; $\Re{\{\cdot\}}$, $\Im{\{\cdot\}}$ and $\mathrm{min}(\cdot)$ represents the real part, imaginary part and the minimum operator, respectively. Finally, $\sim\mathcal{CN}(\mu,\sigma^{2})$ denotes the complex Gaussian distribution with mean $\mu$ and variance $\sigma^{2}$, where $\sim$ denotes “is distributed as”.

II System and Signal Models

Consider a typical SIMO uplink setup composed of $K$ single-antenna user equipment (UE)/edge devices (EDs) and one base station (BS)/access point (AP) equipped with $N$ antennas, as illustrated in Fig. 1. Throughout this article we refer to the UEs/EDs as users, unless otherwise specified, and refer to the BS/AP as the RX.

Figure 1: SIMO ICC system consisting of a RX with $N$ antennas and $K$ single antenna users.

II-A Uplink ICC Signal Model

Under the assumption of perfect symbol synchronization amongst users, the received signal $\bm{y}^{(t)}\in\mathbb{C}^{N\times 1}$ at the BS/AP subjected to fading and noise at a $t$-th time slot is given by

$$\vskip-4.30554pt\bm{y}^{(t)}=\sum_{k=1}^{K}\bm{h}_{k}^{(t)}{x}_{k}^{(t)}+\bm{w}^{(t)},\vskip-2.15277pt$$ (1)

where ${x}_{k}^{(t)}\in\mathbb{C}$ is a multifunctional transmit signal from user $k$, incorporating communication and computing signals; $\bm{w}^{(t)}\in\mathbb{C}^{N\times 1}\sim\mathcal{CN}(0,\sigma^{2}_{w}\mathbf{I}_{N})$ is the additive white Gaussian noise (AWGN) vector, and $\bm{h}_{k}^{(t)}\in\mathbb{C}^{N\times 1}$ is the SIMO channel vector of user $k$ to the RX following the uncorrelated block Rayleigh fading model typically assumed in the AirComp literature [4, 20], such that the $(n,k)$-th elements $h_{n,k}\sim\mathcal{CN}(0,1)$ of the channel matrix $\bm{H}^{(t)}$ are assumed to be independent and identically distributed (i.i.d.) circularly symmetric complex Gaussian random variable with zero mean and unit variance, and sufficient coherence time.

II-B SotA Transmit Signal Design

Under the SotA ICC paradigm [21, 8], a trivial way to model the transmit signal is to decompose it as a sum of communication and computing components, $i.e.,$

$$x_{k}^{(t)}\triangleq d_{k}^{(t)}+\psi_{k}(s_{k}^{(t)}),$$ (2)

where $d_{k}^{(t)}\in\mathcal{D}$ and $s_{k}^{(t)}\in\mathcal{S}$ denote the modulated symbol for communication and computing for user $k$, respectively, with $\mathcal{D}$ representing an arbitrarily scaled discrete constellation of cardinality $D$, $e.g.$ quadrature amplitude modulation (QAM); $\mathcal{S}$ representing a discrete set following [22] for the computing symbols; while $\psi_{k}(\cdot)$ denotes a general pre-processing function for AirComp to be elaborated in the following subsection.

The SotA RX design [8] typically consists of a two-stage approach, where the RX first estimates the individual communication symbols $\{d_{k}^{(t)}\}_{k=1}^{K}$ per time slot via a typical LMMSE or a lower complexity Gaussian belief propagation (GaBP) technique under the assumption that the computing signal is noise, followed by the evaluation of a target function $f(\bm{s}^{(t)})$ via a combiner.

For future convenience, the received signal can now be reformulated in terms of matrices as

$$\bm{y}^{(t)}=\bm{H}^{(t)}\bm{x}^{(t)}+\bm{w}^{(t)}=\bm{H}^{(t)}\big(\bm{d}^{(t)}+\bm{s}^{(t)}\big)+\bm{w}^{(t)},$$ (3)

where the complex channel matrix $\bm{H}^{(t)}\triangleq[\bm{h}_{1}^{(t)},\dots,\bm{h}_{K}^{(t)}]\in\mathbb{C}^{N\times K}$, the concatenated transmit signal $\bm{x}^{(t)}\triangleq[x_{1}^{(t)},\dots,x_{K}^{(t)}]^{\mathsf{T}}\in\mathbb{C}^{K\times 1}$, the data signal vector $\bm{d}^{(t)}\triangleq[d_{1}^{(t)},\dots,d_{K}^{(t)}]^{\mathsf{T}}\in\mathcal{D}^{K\times 1}\subset\mathbb{C}^{K\times 1}$ and the computing signal vector $\bm{s}^{(t)}\triangleq[\psi_{1}(s_{1}^{(t)}),\dots,\psi_{K}(s_{K}^{(t)})]^{\mathsf{T}}\in\mathcal{S}^{K\times 1}\subset\mathbb{C}^{K\times 1}$ are explicitly defined.

In addition, concatenating the received signals over $T$ time slots and assuming that the channel $\bm{H}^{(t)}$ (hereafter denoted as $\bm{H}$ for brevity) remains unchanged, we may express the received signal matrix as

$$\bm{Y}=\bm{H}\bm{X}+\bm{W}=\bm{H}(\bm{D}+\bm{S})+\bm{W},$$ (4)

where $\bm{Y}\triangleq[\bm{y}^{(1)},\dots,\bm{y}^{(T)}]\in\mathbb{C}^{N\times T}$, $\bm{X}\triangleq[\bm{x}^{(1)},\dots,\bm{x}^{(T)}]\in\mathbb{C}^{K\times T}$, $\bm{D}\triangleq[\bm{d}^{(1)},\dots,\bm{d}^{(T)}]\in\mathcal{D}^{K\times T}\subset\mathbb{C}^{K\times T}$, $\bm{S}\triangleq[\bm{s}^{(1)},\dots,\bm{s}^{(T)}]\in\mathcal{S}^{K\times T}\subset\mathbb{C}^{K\times T}$ and $\bm{W}\triangleq[\bm{w}^{(1)},\dots,\bm{w}^{(T)}]\in\mathbb{C}^{N\times T}$ are the corresponding matrices.

Finally, for the sake of derivation of the proposed scheme, we assume that the computing symbols $\bm{s}^{(t)}$ across a block of $T$ time slots are identical (repeatedly transmitted), $i.e.,$ $s_{k}^{(t)}=s_{k},\forall k$, and $\bm{s}^{(t)}=\bm{s},\forall t$, and thus $\bm{S}=[\bm{s},\dots,\bm{s}]\in\mathcal{S}^{K\times T}\subset\mathbb{C}^{K\times T}$.

II-C Description of the AirComp Operation

The AirComp operation consists of the evaluation of a target function $f(\bm{s})$ at the RX, which can be described as [4]

$$f(\bm{s})=\phi\bigg(\sum_{k=1}^{K}\psi_{k}(s_{k})\bigg),\vskip-4.30554pt$$ (5)

where $\phi$ represents the AirComp post-processing function for a general nomographic expression.

While many examples of nomographic functions are often considered in the AirComp literature [3, 23], we choose the arithmetic sum operaton for ease of exposition, given by

$$f(\bm{s})=\phi\bigg(\sum_{k=1}^{K}\psi_{k}(s_{k})\bigg)=\sum_{k=1}^{K}s_{k},\vskip-4.30554pt$$ (6)

where the corresponding pre- and post-processing functions are defined as $\psi_{k}(s_{k})\triangleq s_{k}$ and $\phi\big(\sum_{k=1}^{K}\psi_{k}(s_{k})\big)\triangleq\sum_{k=1}^{K}\psi_{k}(s_{k})$.

III Computing on Dirty Paper

III-A Fundamentals on Dirty Paper Coding

Dirty paper coding (DPC), first introduced by Costa in [10], is a powerful coding technique that allows the transmission of information over a noisy channel in the presence of known interference at the TXs, without being affected by the latter. In more detail, the TX– the arbitrary UEs or EDs in our case – needs to transmit a signal $d$ subject to interference $s$. Then, the RX receives signal $y=d+s+w$, where $w$ is random noise. When $s$ is unknown to both the TX and RX, it can only be incorporated into the noise as previously proposed in [8] for the ICC paradigm. However, when $s$ is known to the TX, the RX can completely remove the interference $s$ by leveraging DPC, even though it does not know $s$. In this paper, we introduce a novel scheme for the DPC operation in the context of ICC, termed Computing on Dirty Paper.

III-B Fundamentals on Nested Lattice Codes

For the design of the lattice, consider a standard Gaussian integer lattice $\mathbb{Z}[j]=\{a+bj\mid a,b\in\mathbb{Z}\}$ with generator matrix defined to be $\mathbf{I}_{2}$, as illustration in Fig. 2, which is simple and effective for complex signals (equivalent to a square lattice in $\mathbb{R}^{2}$) in a similar fashion to [13] as follows:

Coarse lattice: Let $\Lambda_{c}=\Delta\mathbb{Z}[j]\subset\mathbb{C}$, where $\Delta>0$ sets the modulo region; $i.e.,$ the “coarse lattice cube” bounded by sides of length $[-\Delta/2,\Delta/2]$. This choice ensures that the transmit power is bounded during the modulo operation.

Fine lattice: Let $\Lambda_{f}=\delta\mathbb{Z}[j]$, with $\delta=\Delta/\sqrt{M}$ and $M=2^{2R}$ where $M$ represents the (integer) number of constellation points inside each coarse lattice cell and rate $R$ denotes the number of bits per complex symbol for $d_{k}$.

Cosets: Each data symbol $d_{k}$ now corresponds to a coset of $\Lambda_{f}$ where

$$\vskip-4.30554ptv_{k}+\Lambda_{c},v_{k}\in\Lambda_{f}\cap\mathcal{V}(\Lambda_{c}),$$ (7)

with $\mathcal{V}(\Lambda_{c})$ denoting the Voronoi region of the coarse lattice $\Lambda_{c}$, which is a square in the complex plane (of area $\Delta^{2}$), with distortion bounded by $\Delta/\sqrt{2}$ in Euclidean norm for nearest-neighbor quantization.

Note that the transmit constellation is the set of all coset representatives inside a single Voronoi cell, and not the entire coset. For a more compresensive overview of nested lattice codes, the reader is referred to [24].

III-C Transmitter Design

Reintroducing the time-varying component for the data, let each transmit user know its own encoded data symbol $v_{k}^{(t)}\in\Lambda_{f}\cap\mathcal{V}(\Lambda_{c})$ and computing symbol picked from the rotated set $\mathcal{S}_{\text{rot}}\triangleq\mathbf{R}\,\mathcal{S}=\bigl\{\mathbf{R}\,s_{k}\,\big|\,s_{k}\in\Lambda_{f}\cap\mathcal{V}(\Lambda_{c})\bigr\}$ to maximize orthogonality, where

$$\vskip-4.30554pt\mathbf{R}\triangleq\frac{1}{\sqrt{2}}\begin{bmatrix}1&-1\\ 1&1\end{bmatrix}.$$ (8)

Nested lattice-based DPC encodes the communications symbols $v_{k}^{(t)}$ leveraging a modulo operation such that the transmit symbol $\tilde{v}_{k}^{(t)}$ is given by

$$\tilde{v}_{k}^{(t)}=\text{modL}(v_{k}^{(t)}-s_{k},\Delta),$$ (9)

where the $\text{modL}(a,b)$ operation is defined as

$$\text{modL}(a,b)\triangleq$$

$$\hskip-43.05542pt\Big(\text{mod}\big(\Re\{{a}\}\!+\!\tfrac{b}{2},b\big)\!-\!\tfrac{b}{2}\Big)\!+\!j\Big(\text{mod}\big(\Im\{{a}\}\!+\!\tfrac{b}{2},b\big)\!-\!\tfrac{b}{2}\Big),$$

with $\text{mod}(a,b)$ denoting the standard modulo operation.

Then, the multifunctional transmit signal $x_{k}^{(t)}$ across $T$ time slots (analogous to the SotA expression in equation (2)) can be expressed as

$$x_{k}^{(t)}=\tilde{v}_{k}^{(t)}+s_{k},$$ (11)

where we emphasize that the computing symbol $s_{k},\forall k$ remains unchanged across $T$ time slots.

Figure 2: Nested lattice structure with coarse lattice $\Lambda_{c}$ and fine lattice $\Lambda_{f}$.

III-D Receiver Design

After transmission of $T$ slots over a typical Rayleigh fading channel, the received signal in matrix form can be expressed as

$$\vskip-4.30554pt\bm{Y}=\bm{H}\bm{X}+\bm{W}=\bm{H}(\tilde{\bm{V}}+\bm{S})\in\mathbb{C}^{N\times T},$$ (12)

where similarly to previous definitions, $\tilde{\bm{V}}\triangleq[\tilde{\bm{v}}^{(1)},\dots,\tilde{\bm{v}}^{(T)}]\in\mathbb{C}^{K\times T}$ with $\tilde{\bm{v}}^{(t)}\triangleq[\tilde{v}_{1}^{(t)},\dots,\tilde{v}_{K}^{(t)}]^{\mathsf{T}}\in\mathbb{C}^{K\times 1}$ and $\bm{S}=[\bm{s},\dots,\bm{s}]\in\mathcal{S}_{\text{rot}}^{K\times T}\subset\mathbb{C}^{K\times T}$.

Note that different time slots contain different data symbols $\tilde{v}_{k}$, retaining the full rate in terms of data transmission identical to typical communications systems, while the computing symbols $s_{k}$ remain unchanged across $T$ time slots to improve the mean-squared-error (MSE) performance of the function evaluation. One can then apply a LMMSE estimator per $t$-th time frame at the receiver to obtain $\hat{\bm{x}}^{(t)}$, which can be expressed as

$$\vskip-4.30554pt\hat{\bm{x}}^{(t)}=(\bm{H}^{\mathsf{H}}\bm{H}+\sigma_{w}^{2}\mathbf{I}_{K})^{-1}\bm{H}^{\mathsf{H}}\bm{y}^{(t)}.$$ (13)

Next, a temporary variable $\hat{\bm{z}}$ can be calculated via the modulo operation

$$\hat{\bm{z}}^{(t)}=\text{modL}(\hat{\bm{x}}^{(t)},\Delta),$$ (14)

which can then be leveraged to estimate $\hat{\bm{v}}^{(t)}$ via

$$\hat{\bm{v}}^{(t)}=\text{modL}\big(\mathcal{Q}_{\Lambda_{f}}(\hat{\bm{z}}^{(t)},\delta),\Delta\big),$$ (15)

where $\mathcal{Q}_{\Lambda_{f}}(\cdot)$ denotes the nearest-neighbor quantizer to the fine lattice $\Lambda_{f}$, and is defined as

$$\mathcal{Q}_{\Lambda_{f}}(a,b)\!\triangleq\!b\Big(\text{round}\big(\tfrac{\Re\{{a}\}}{b}-\epsilon\big)+\epsilon\Big)+jb\Big(\text{round}\big(\tfrac{\Im\{{a}\}}{b}-\epsilon\big)+\epsilon\Big),$$ (16)

with $\text{round}(\cdot)$ denoting the rounding operation to the nearest integer and $\epsilon\triangleq\text{mod}\Big(\frac{\sqrt{M}-1}{2},b\Big)$.

Finally, a hard decision against the constellation $\Lambda_{f}\cap\mathcal{V}(\Lambda_{c})$ can be performed to obtain the final estimate of the encoded data symbols $\hat{\bm{v}}^{(t)}$.

In order to recover the computing symbols, a novel best-case bound  can be derived in the form of a maximum likelihood estimator which can be expressed as

$$\hat{\bm{s}}\!=\!\arg\min_{\bm{s}\in\mathcal{S}_{\text{rot}}^{K\times 1}}\sum_{t=1}^{T}\big\|\bm{y}^{(t)}\!-\!\bm{H}\big(\mathrm{modL}(\hat{\bm{v}}^{(t)}\!-\!\bm{s},\Delta)\!+\!\bm{s}\big)\big\|_{2}^{2},$$ (17)

where each element $\hat{s}_{k}$ of $\hat{\bm{s}}$ is chosen from the discrete set $\mathcal{S}_{\text{rot}}$ that has the lowest error across all time instances $T$, naturally improving the MSE performance as $T$ increases.

The computational complexity of equation (17) is dominated by the number of $K$ users, since an exponential scaling is needed to compare with a constellation of size $M$ such that the dominant complexity is given by $M^{K}$ in order to obtain a maximum likelihood estimate.

Finally, the target function can be computed as

$$\hat{f}(\bm{s})=\sum_{k=1}^{K}\hat{s}_{k}.$$ (18)

IV Performance Analysis

To evaluate the performance of the proposed scheme, let us first visualize the complete nested lattice structure that is used, as shown in Fig. 3. Notice that each coarse lattice cell contains $M=2^{2R}$ constellation points, where $R$ is the number of bits per complex symbol for data and this fine lattice (corresponding to a scaled quadrature phase-shift keying (QPSK)) is used for data symbol transmission. In addition, to maximize orthogonality between data and computing symbols, the computing symbols are chosen from a rotated version of the fine lattice constellation, as previously described.

Figure 3: Nested lattice structure with coarse lattice $\Lambda_{c}$ and fine lattice $\Lambda_{f}$.

(a) $T=5$ (Short Frame Structure)

(b) $T=10$ (Long Frame Structure)

Figure 4: Comparison of BER and MSE performances of SotA [8] and proposed computing on dirty paper ICC schemes.

In order to bound the transmit power for a fair comparison with the SotA, we set the coarse lattice parameter to be $\Delta=2$, which leads to a transmit power of $E[|x_{k}|^{2}]=E_{d}+E_{s}=0.5+0.5=1$ per user, where $E_{d}=E[|d_{k}|^{2}]=0.5$ and $E_{s}=E[|s_{k}|^{2}]=0.5$ are the average powers of data and computing symbols, respectively.

The parameters used in the simulations are $K=2$ users, $N=5$ antennas at the RX, an alternating time slot of $T=5,10$, and $R=2$ bits/symbol (i.e., for QPSK).

Figure 4(a) and 4(b) show the BER and MSE performance of the proposed scheme compared to the SotA ICC method in [8], where the latter uses a two-stage approach with a GaBP detector for data detection followed by a linear combiner for function estimation. As seen, the proposed scheme achieves significant performance gains in both BER and MSE compared to the SotA, with the performance gap increasing with signal-to-noise ratio (SNR). This is because the proposed scheme effectively eliminates the interference between data and computing symbols via DPC, allowing both tasks to be performed more accurately. In addition, increasing the number of time slots $T$ improves the MSE performance for the proposed scheme, as more observations are available for function estimation.

V Conclusion

We proposed Computing on Dirty Paper, a new scheme for ICC inspired by DPC. By pre-canceling computing symbols at the transmitting users, the scheme enables interference-free communication and computation over multiple-access channels. Analysis under a SIMO ICC system and supporting simulations show that the approach achieves superior BER and MSE performance compared to SotA ICC methods.