Virtual Signal Processing-Based Integrated Multi-User Detection

TL;DR

This paper introduces a new multi-user detection method that reduces errors and improves performance in communication systems.

Contribution

The novel VSP-IMUD method integrates virtual signal processing to reduce SIC error propagation.

Findings

VSP-IMUD reduces the frequency of SIC applications.

The method improves system bit-error rate performance.

Theoretical analysis and simulations confirm the effectiveness of VSP-IMUD.

Abstract

The demand for high data rates and large system capacity has posed significant challenges for medium access control (MAC) methods. Successive interference cancellation (SIC) is a classical multi-user detection (MUD) method; however, it suffers from an error propagation problem. To address this deficiency, we propose a method called Virtual Signal Processing-Based Integrated Multi-User Detection (VSP-IMUD). In VSP-IMUD, the received mixed multi-user signals are treated as an equivalent signal. The channel ambiguity corresponding to each user’s signal is then examined. For channels with non-zero ambiguity values, the signal components are detected using zero-forcing (ZF) reception. Next, the detected ambiguous signal components are reconstructed and subtracted from the received mixed signal using SIC. Once all the ambiguous signals are detected, the remaining signal components with zero ambiguity values are equated to a virtual integrated signal, to which a matched filter (MF) is applied. Finally, by selecting the signal with the highest channel gain and adopting its data as the reference symbol, the remaining signals’ dataset can be determined. Our theoretical analysis and simulation results demonstrate that VSP-IMUD effectively reduces the frequency of SIC applications and mitigates its error propagation effects, thereby improving the system’s bit-error rate (BER) performance.

Figures

Variation of the probability of direct applying IMUD with K under NT=1, NR=2, and various εs.

Full text

1. Introduction

The ever-increasing demand for wireless broadband services has driven the development of 5G technology. Compared to previous mobile communication systems, 5G is expected to offer larger system capacity, higher data rates, lower latency, and improved reliability [1]. As a key application in the 5G era, the Internet of Things (IoT) has been rapidly developing in recent years, leading to explosive growth in the number of IoT devices. It is estimated that by 2025, more than 41.6 billion IoT devices will be connected to networks worldwide [2,3]. Therefore, efficiently supporting the data transmissions of more users under limited communication resources has become a critical issue warranting thorough investigation.

By allocating different types of communication resources—such as frequency, time slots, code words, and spatial characteristics—to multiple users in an orthogonal manner, interference among concurrent data transmissions can be effectively avoided, thereby enabling resource sharing among users [4]. However, since the frequency resources available for wireless communication are limited, the emergence of various communication systems and the increase in mobile users pose significant challenges to traditional orthogonal multiple access (OMA) methods. As a result, dynamic spectrum sharing methods have been proposed to facilitate multi-user communications. For example, cognitive radio (CR) [5,6] enables unlicensed or cognitive users to access the available authorized spectrum for data transmission. However, timely sensing of the licensed user’s signal and withdrawal is critical for CR applications, as inappropriate spectrum access may interfere with authorized communications. Due to advancements in microelectronic technology, the processing capability of hardware has improved while costs have decreased. A receiver (Rx) can now decode data carried in multiple concurrent signals through complex signal processing. In recent years, non-orthogonal multiple access (NOMA) [7] has emerged as a promising medium access control (MAC) method for the 5G system. Uplink NOMA allows multiple transmitters (Txs) to transmit to a common Rx over the same frequency channel in a non-orthogonal manner, enabling the Rx to employ successive interference cancellation (SIC) [8] to mitigate co-channel interference (CCI).

As mentioned above, in addition to resource sharing, receiver-side multi-user detection (MUD) is essential for achieving high system capacity and data rates. Existing MUD methods can be classified into linear detection and non-linear detection based on whether their structures involve feedback mechanisms. Linear MUD methods include linear de-correlation detection [7] and minimum mean squared error (MMSE) detection [8]. However, linear MUD [9] requires computing the inverse of the cross-correlation matrix of all users’ channels, which incurs high computational complexity. Non-linear MUD, on the other hand, includes decision feedback detection methods such as SIC. SIC first recovers data information carried in certain signal components from the received mixed signal, then uses the decoded data and estimated channel state information (CSI) to reconstruct the previously detected signals. These reconstructed signals are subtracted from the received mixed signal, resulting in a reduced mixed signal. By performing this iterative processing on the reduced mixed signal, multiple users’ signals can be detected. By combining MMSE and SIC, MMSE-SIC detection has been proposed and applied within the V-BLAST framework [10,11]. However, SIC has an error propagation problem [12], leading to a high bit-error rate (BER) for later-decoded users’ data, which hampers its practical application.

Although increasing the number of receiving antennas can enhance the Rx’s signal processing capability and support more users’ data transmissions, practical limitations such as device size and hardware costs make it unfeasible to increase the number of receiving antennas indefinitely. Therefore, utilizing a limited number of receiving antennas for multi-user detection (MUD) holds significant research importance. By exploiting interactions among multiple wireless signals, the authors of [12] proposed a MUD scheme named ICom/SIC-ZF, which incorporates interference combination (ICom) with zero-forcing (ZF) reception and SIC. This method treats all signals other than the desired signals (i.e., those to be detected) as a structured equivalent interference, and uses ZF to cancel this interference while recovering the desired signals. Subsequently, SIC is employed to subtract the detected signals from the received mixed signal. By iterating this process on the remaining mixed signal, all users’ signal components can be detected sequentially. In [13], the authors utilized maximum likelihood (ML) to decode data from two overlapping wireless signals simultaneously under a single receiving antenna setting, based on the characteristics of various combinations of symbols carried in the mixed signal. However, with high-order modulation or a large number of users’ signals, multiple symbol combinations can produce identical constellation points, leading to decoding errors. In such cases, the BER performance of this method diminishes significantly. The authors of [14] utilized intelligent reflecting surfaces (IRSs) to reflect a transmitted signal towards the intended Rx. By intelligently adjusting the reflecting parameters of the IRS, the reflected signal can interact with the directly arriving component to produce a desired signal at the Rx. Since neither the directly transmitted nor the reflected signal component carries legitimate data information, it becomes difficult or even impossible for an eavesdropper to intercept the legitimate user’s information. In [15], an IRS was employed to reflect interference towards the affected Rx. The reflected interference serves as a recycling signal that interacts with the direct interference to generate a desired signal, thereby achieving interference recycling and enhancing data transmission.

Based on the above discussions, this paper proposes a method called Virtual Signal Processing-Based Integrated Multi-User Detection (VSP-IMUD). This approach achieves integrated MUD (IMUD) by leveraging the spatial features of the equivalent desired signals that consist of multiple users’ signals. First, we design a channel ambiguity detection (CAD) module to avoid similarity or ambiguity among different equivalent (virtual) desired signal structures. In scenarios where similarity or ambiguity exists, directly applying IMUD may result in decoding errors and poor BER [13]. The CAD module identifies channels that yield ambiguity and then employs ICom/SIC-ZF to detect signals transmitted through these ambiguous channels. Once all signals transmitted through the ambiguous channels are decoded, the remaining signal components are treated as an equivalent desired signal through virtual signal combination. Based on the spatial characteristics of this equivalent desired signal, we can derive possible symbol combinations of the various desired signal components. Finally, we use ZF to recover the data sent from the Tx with the largest channel gain relative to the common Rx, and adopt it as the reference data to determine the exact symbol set of the equivalent signal. In this way, VSP-IMUD is realized.

The contributions of this paper are two-fold:

• Proposal of VSP-IMUD. This method allows for the integral decoding of data information carried in multiple signal components at once, based on the overall spatial characteristics of the received mixed signal. This enables the common Rx to perform multi-user detection (MUD) with fewer receiving antennas.
• Proposal of CAD. This module first identifies channels that cause the structure of one equivalent desired signal to be ambiguous relative to another. By detecting signals transmitted through these ambiguous channels, the uniqueness of the remaining signals’ equivalent structure is ensured, allowing for IMUD to effectively decode multiple users’ data simultaneously.

The rest of this paper is organized as follows. Section 2 describes the system model, while Section 3 details the design of VSP-IMUD. Then, Section 4 evaluates the performance of the proposed scheme. Finally, Section 5 concludes the paper.

Throughout this paper, we will use the following notations. Let $[eqn]$ denote the set of complex numbers. Vectors and matrices are represented by bold letters. $[eqn]$ , $[eqn]$ and $[eqn]$ denote the transpose, Hermitian, and inverse of matrix $[eqn]$ , respectively. $[eqn]$ represents the Euclidean norm of vectors and matrices, and $[eqn]$ denotes the inner product of two vectors.

2. System Model

We consider a multi-user uplink communication system consisting of K Txs and one common Rx, as shown in Figure 1. This figure illustrates a general multi-user uplink communication model. In practice, the Rx can function as either a base station (BS) or an access point (AP), and the Txs are mobile stations. In the system model, each Tx and Rx are equipped with $[eqn]$ and $[eqn]$ antennas, respectively. Assuming $[eqn]$ , a Tx can employ either beamforming (BF) or spatial multiplexing (SM) to send one or multiple directional beams to the Rx. For simplicity, we assume that all Txs transmit to the Rx simultaneously, and their signals can arrive at the Rx synchronously. The transmit power of Tx is denoted as $[eqn]$ . We use $[eqn]$ (where $[eqn]$ ) to represent the data sent from $[eqn]$ to Rx. $[eqn]$ holds where $[eqn]$ denotes statistical expectation. We denote the composite channel matrix of the K Txs as $[eqn]$ , where $[eqn]$ ( $[eqn]$ ) represents the channel matrix from $[eqn]$ to Rx. We employ a spatially uncorrelated Rayleigh flat fading channel model (in this commonly used channel model, “spatially uncorrelated” is based on the assumption that the inter-antenna spacing is sufficiently large to prevent correlations among the channel fading experienced by various transmit–receive antenna pairs. Furthermore, “Rayleigh flat fading” indicates that the channel varies relatively slowly compared to the processing time of the proposed scheme. Typically, we assume that the relative speed between the Tx and the Rx is not high), where the channel elements are modeled as independent and identically distributed complex Gaussian random variables with zero mean and unit variance. We assume that each Tx, say Tx_k_, can accurately estimate $[eqn]$ and report it to the Rx via a low-rate, error-free control link. The delivery delay of CSI and signaling is negligible compared to the time scale on which the channel state varies [16,17].

As Figure 1 shows, at each Tx, the bit stream is first divided into multiple sets, and each set is converted into a data symbol through constellation mapping. This symbol is then fed to the transmit antenna for transmission. We denote the modulation symbol set used by $[eqn]$ as $[eqn]$ where $[eqn]$ and $[eqn]$ denotes the number of elements in the set. Here, L represents the modulation order, and all Txs adopt the same modulation scheme. Since we focus on the design of Rx-side processing, we set $[eqn]$ for simplicity. However, it is important to note that the proposed scheme can be extended to the case where $[eqn]$ . In such instances, the Tx would need to preprocess the data before transmission.

3. Design of VSP-IMUD

In this section, we present the design of VSP-IMUD and describe the receiver structure employing VSP-IMUD. By exploiting interactions among wireless signals, the detection of multiple users’ signals can be achieved with fewer receiving antennas.

3.1. Integrated Multi-User Detection

As Figure 1 plots, when multiple Txs simultaneously send signals to the common Rx, the mixed signals received by the Rx can be expressed as

$$[eqn]$$

where $[eqn]$ denotes the Additive White Gaussian Noise (AWGN) vector whose elements have zero mean and variance $[eqn]$ . For clarity, we denote the mixed signal $[eqn]$ as $[eqn]$ .

Under IMUD, the received mixed signal $[eqn]$ is treated as a structured equivalent signal. Based on the characteristics of such an equivalent signal, the data carried in multiple signal components can be decoded simultaneously, thereby realizing MUD. Given K and L, there will be $[eqn]$ forms of $[eqn]$ . Assuming that all Txs adopt L-order modulation, we can represent the modulation symbol set as $[eqn]$ , with each Tx_k_’s desired data $[eqn]$ . Based on the composite channel matrix $[eqn]$ and all possible combinations of data from the K Txs, the Rx can derive $[eqn]$ possible forms of $[eqn]$ , as presented in Equation (2).

$$[eqn]$$

As Equation (2) shows, there are $[eqn]$ possible combinations of the symbol set $[eqn]$ . Superscripts $[eqn]$ are used to index these symbol combinations. We denote the data combination indexed by $[eqn]$ as $[eqn]$ , where $[eqn]$ and $[eqn]$ ( $[eqn]$ ) represents the data symbol sent from $[eqn]$ . Consequently, the m-th form of the received mixed signal $[eqn]$ can be expressed as $[eqn]$ . If $[eqn]$ is uniquely related to $[eqn]$ , we can obtain $[eqn]$ by detecting the structured $[eqn]$ in an integrated manner.

We define the spatial feature of $[eqn]$ as

$$[eqn]$$

According to Equations (2) and (3), given $[eqn]$ and $[eqn]$ , the structure and spatial features of $[eqn]$ are determined. Based on Equation (3), we can employ a matched filter (MF) $[eqn]$ to process $[eqn]$ . Since $[eqn]$ is the target that the Rx aims to detect, we need $[eqn]$ filter chains to attempt to decode the data vector $[eqn]$ from $[eqn]$ . For each filter chain, we calculate the modulus of its output and select the one with the largest modulus as the correct chain. Then, we can determine the receiving filter as $[eqn]$ where M is computed in terms of $[eqn]$ . In this way, $[eqn]$ is decoded integrally. In practical applications, as K and/or L increase, the structure of SVP-IMUD becomes more complicated. To mitigate this complexity, we can divide the K Txs into multiple smaller groups, effectively reducing K before applying VSP-IMUD. Furthermore, we can increase $[eqn]$ to enhance the hardware processing capabilities at the Rx, achieving a trade-off between the number of matched filter chains and the number of receiving antennas. Additionally, it is important to note that the matched filter chains are implemented in the digital domain, eliminating the need for radio frequency chains, which makes them more cost-effective compared to receiving antennas.

We denote the decoded data vector using $[eqn]$ as $[eqn]$ . It is important to note that two possible forms of $[eqn]$ , namely $[eqn]$ and $[eqn]$ , may arise due to the symmetric nature of the constellation points with respect to the origin. However, only one of $[eqn]$ and $[eqn]$ corresponds to the correct data vector sent from the K Txs. To determine the exact data vector, we use ZF to recover the data from the Tx with the largest channel gain relative to the RX, and adopt this as referential data for determining the exact data vector. In Figure 2, we illustrate the realization of VSP-IMUD using binary phase shift keying (BPSK) as an example.

As the figure shows, with K Txs and $[eqn]$ , two data vectors $[eqn]$ and $[eqn]$ correspond to the same MF vector $[eqn]$ . We denote these two symbol sets as $[eqn]$ and $[eqn]$ , respectively. Since the elements in $[eqn]$ and $[eqn]$ are opposite, we use ZF to decode the data sent from the Tx with the largest channel gain to obtain the correct symbol set. We denote this referential symbol as $[eqn]$ (where $[eqn]$ ). We then compare the referential symbol $[eqn]$ with the symbol under the same index in sets $[eqn]$ and $[eqn]$ . If there is an $[eqn]$ in $[eqn]$ that matches $[eqn]$ , we can determine that $[eqn]$ is the correct symbol set. Otherwise, $[eqn]$ becomes the final decoding output.

In practice, due to the randomness of the Tx channels, the Euclidean distance between two $[eqn]$ s may be very small, i.e., $[eqn]$ (for $[eqn]$ and $[eqn]$ ), where $[eqn]$ is a small positive real number. This situation can lead to incorrect outputs from the Rx structure shown in Figure 2, resulting in a high bit-error rate (BER). We refer to this condition as channel ambiguity, and we will design an ambiguity detection and elimination method in the following subsection.

3.2. Channel Ambiguity Detection and Elimination

In the previous subsection, we assumed the structure of the received mixed signal to be unique and employed an MF to recover the data/symbol set sent from multiple TXs. In practice, however, the structure of the mixed signal may arise from various data and channel combinations, resulting in similarities and ambiguities among the characteristics of the equivalent signals. In such cases, the MF cannot achieve accurate data recovery [13].

Therefore, in this subsection, we will design a channel ambiguity detection (CAD) module to identify the channels that contribute to the similarity of the mixed signal’s structure, and then decode the signals transmitted via these ambiguous channels successively using SIC. Once the channel ambiguity is completely mitigated, the structure of the remaining mixed signal becomes unique, allowing for exclusive recovery of the dataset from multiple Txs using the IMUD method presented in Section III-A. Below, we will detail the detection and elimination of channel ambiguity.

We denote the set of symbols to be detected as $[eqn]$ , and initialize it to $[eqn]$ , where $[eqn]$ (for $[eqn]$ ) is the data symbol sent from $[eqn]$ . According to Equation (2), the mixed signal perceived by the Rx is $[eqn]$ . The Rx can then obtain all possible combinations of the aggregated channel matrix $[eqn]$ and the symbol set $[eqn]$ . Note that although Rx does not know the desired data sent by the Txs at this time, it can still derive $[eqn]$ based on the known modulation scheme, where each element indicates a possible structure or characteristic of the received mixed signal. Next, we use the Euclidean distance between $[eqn]$ and $[eqn]$ as the ambiguity coefficient $[eqn]$ of $[eqn]$ and $[eqn]$ , defined in Equation (4):

$$[eqn]$$

where $[eqn]$ and $[eqn]$ .

We employ a real value $[eqn]$ as the ambiguity threshold. If $[eqn]$ , the structures of $[eqn]$ and $[eqn]$ are similar (ambiguous), indicating detection of ambiguity between $[eqn]$ and $[eqn]$ . It should be noted that the value of $[eqn]$ affects the processing complexity and BER performance of the Rx. When $[eqn]$ is large, more pairs $[eqn]$ will be classified as similar, leading to a greater number of Tx signals that need to be decoded to eliminate ambiguity, which increases the processing burden. Conversely, if $[eqn]$ is small, the similarities in features of the mixed signals arising from various datasets may go undetected. In this situation, while IMUD may still be applicable, it could yield an incorrect symbol set, resulting in poor BER performance. It is important to note that $[eqn]$ serves to balance decoding complexity and accuracy. Moreover, the degree of ambiguity depends on the K channel matrices associated with the K Txs, the K data symbols, and their modulation order L. Among these factors, the K channel matrices and data symbols are random and independent, making a quantitative analysis of the relationship among $[eqn]$ , K and L difficult. For this reason, we utilize simulations to investigate the impact of different values of $[eqn]$ on the probability of direct decoding with IMUD, as illustrated in Figure 4. Consequently, we can empirically determine the value of $[eqn]$ to achieve an optimal balance between signal processing complexity and system performance.

In order to detect the ambiguity among various pairs $[eqn]$ , we define an $[eqn]$ ambiguity matrix $[eqn]$ as shown in Equation (5). As depicted, the element in the i-th row and j-th column of $[eqn]$ is either 1 or 0, indicating whether the spatial features of $[eqn]$ and $[eqn]$ are ambiguous. Specifically, if the Euclidean distance between $[eqn]$ and $[eqn]$ is less than $[eqn]$ , then $[eqn]$ ’s element $[eqn]$ is set to 1; otherwise, it is set to 0. It should be noted that under $[eqn]$ , we set $[eqn]$ .

$$[eqn]$$

Next, we calculate the sum of the elements in each row of $[eqn]$ , and record the result of the i-th row as $[eqn]$ . The value $[eqn]$ reflects the number of mixed signals that are similar to $[eqn]$ . The larger $[eqn]$ is, the more ambiguous the structure of $[eqn]$ is in relation to other mixed signals. In this case, if one performs IMUD on the received signal, s/he will have to choose among multiple possible symbol sets, resulting in a high BER due to this ambiguity.

The ambiguity arises from the fact that the same aggregated channel matrix $[eqn]$ and various symbol sets $[eqn]$ where $[eqn]$ can produce identical structures for the mixed signals. To avoid ambiguity, we identify the channels in $[eqn]$ that contribute to the ambiguity of the signal structure and subsequently remove them from $[eqn]$ . In what follows, we will establish the criterion for determining the contribution of channel $[eqn]$ ( $[eqn]$ ) to the similarity between different mixed signals. We use $[eqn]$ to denote the ambiguity coefficient of $[eqn]$ and initialize $[eqn]$ to 0, indicating that $[eqn]$ does not contribute to ambiguity.

We assume that the $[eqn]$ -th row of $[eqn]$ that outputs the largest $[eqn]$ has been determined. If $[eqn]$ , the structures of $[eqn]$ and $[eqn]$ are similar. Then, according to Equation (4), we have $[eqn]$ . Under $[eqn]$ , if $[eqn]$ ( $[eqn]$ ), we can have $[eqn]$ . In this case, we cannot determine the contribution of $[eqn]$ to the ambiguity of the mixed signal structures; thus, the value of $[eqn]$ remains unchanged. Conversely, if $[eqn]$ ( $[eqn]$ ), $[eqn]$ holds, indicating that $[eqn]$ contributes to the structure ambiguity of the received signals. Therefore, we increment $[eqn]$ as $[eqn]$ . If $[eqn]$ , $[eqn]$ does not contribute to the ambiguity of the received signals’ structures, and the value of $[eqn]$ remains unchanged.

After checking each $[eqn]$ in the $[eqn]$ - $[eqn]$ row, we can obtain the ambiguity coefficient $[eqn]$ for each $[eqn]$ . Similarly, we can examine all rows of $[eqn]$ to determine the ambiguity coefficients for the channels of all K Txs. We then sort $[eqn]$ ( $[eqn]$ ) in descending order. This process establishes the decoding order for the signals that generate ambiguity. Specifically, the signal transmitted via $[eqn]$ with the largest $[eqn]$ , denoted as $[eqn]$ , should be detected first.

Since Rx has $[eqn]$ receiving antennas, it can utilize ZF reception to concurrently decode up to $[eqn]$ signal components associated with the channels with the largest ambiguity coefficients in the spatial domain. For the remaining $[eqn]$ signals, they are treated as equivalent interference [11]. The Rx then employs SIC to eliminate the detected ambiguous signal components from the original mixed signal and updates the symbol set $[eqn]$ by removing the decoded symbols.

Next, the Rx recalculates the ambiguity matrix and the channel ambiguity coefficients $[eqn]$ for the remaining signals and decodes any signal components that still exhibit structural ambiguity. This processing is repeated iteratively until $[eqn]$ , at which point the IMUD method introduced in the previous subsection can be applied to the remaining mixed signals, as there will be no similarities among these components. In this way, the desired data carried in the remaining mixed signal can be uniquely determined using the spatial features of the mixed signal and the decoded referential data symbol.

It is important to note that if $[eqn]$ exceeds the number of signals to be detected, we can directly apply ZF to detect the signals. In other words, the proposed VSP-IMUD is primarily focused on scenarios where the number of signals to be detected is greater than $[eqn]$ , where ZF is infeasible due to the RX’s limited number of antennas.

3.3. VSP-IMUD-Based Rx Structure

Based on the designs of IMUD and CAD in the previous two subsections, we now present the VSP-IMUD-based Rx structure. As Figure 3 shows, there are three processing steps in the Rx structure.

Step 1: At Switch A, we check the relationship between the number of symbols to be decoded, denoted as $[eqn]$ , and the total number of Txs, denoted as K. If $[eqn]$ , we switch to A1, and the received mixed signal is fed into the CAD module for ambiguity detection (i.e., proceed to Step 2). If $[eqn]$ , we switch A to A2 and continue with ambiguity detection for the remaining mixed signal (i.e., proceed to Step 2).

Step 2: For the mixed signal that has passed through the CAD module, we check $[eqn]$ at switch B. If $[eqn]$ , we switch B to B1, and the mixed signal is fed into the IMUD module. As Figure 3 shows, there are $[eqn]$ filter chains in this structure. Each chain uses a filter vector that matches the spatial feature of a possible equivalent/structured mixed signal.

We then calculate the modulus of the output from each filter chain and determine the index of the filter chain yielding the largest modulus. Based on this index information, we can obtain two potential desired symbol sets, namely $[eqn]$ and $[eqn]$ . To finalize the symbol set, the Rx employs a ZF filter vector $[eqn]$ to recover the signal from the Tx with the largest channel gain relative to the Rx, and adopts the decoded data $[eqn]$ as the referential symbol. If there is an $[eqn]$ in $[eqn]$ such that $[eqn]$ , we switch C to C1 and output $[eqn]$ as the final decoded symbol set. Conversely, if $[eqn]$ , we switch C to C2 and output $[eqn]$ . In this way, signals from all TXs are detected.

At switch B, if $[eqn]$ , it indicates that the ambiguity of the mixed signal has not been fully resolved. Therefore, we switch B to B2 to allow the Rx to sort the channels by their ambiguity values in descending order and to detect the signals propagating through these channels. If the number of signals with $[eqn]$ exceeds $[eqn]$ , the Rx employs ICom/SIC-ZF to detect the signals from the Txs corresponding to the first $[eqn]$ channels with the highest $[eqn]$ values. This results in obtaining $[eqn]$ data symbols, which we denote as $[eqn]$ . We can then reconstruct the detected signals and utilize SIC to subtract them from the mixed signal. After this, we can update the symbol set $[eqn]$ as $[eqn]$ . Finally, we compare $[eqn]$ with $[eqn]$ and proceed to Step 3.

Step 3: At switch D, if $[eqn]$ , we switch D to D1 and concurrently switch E to E1. The Rx can then apply the ZF matrix $[eqn]$ to directly detect the remaining mixed signals and obtain the detected symbol set $[eqn]$ . This completes the MUD process. If $[eqn]$ , it indicates that the Rx does not have enough antennas to decode $[eqn]$ data. Moreover, it remains uncertain whether the ambiguity in the symbol set $[eqn]$ has been mitigated. Consequently, we switch D to D2, switch F to F1, and return to Step 1, sending the mixed signal back to switch A. At switch A, since $[eqn]$ , we switch A to A2 and perform channel ambiguity detection. Based on the output from the CAD module, we switch B to either B1 or B2 (as discussed in Step 2). These steps are repeated until all K signals are successfully detected.

4. Simulation Results

In this section, we use MATLAB simulations to demonstrate the performance of VSP-IMUD. The system consists of K Txs and one Rx. Since we focus on a multi-user communication scenario, we set $[eqn]$ in our simulations. Although reducing K to 1 can eliminate CCI among multiple Txs, it also limits the system’s ability to accommodate only one user’s data transmission at a time, thereby degrading the system’s spectral efficiency (SE). Therefore, we omit the study of the single-user scenario in this work. Each Tx is equipped with $[eqn]$ antennas and employs BPSK modulation to transmit a single data stream to the Rx with power $[eqn]$ . The common Rx has $[eqn]$ antennas. If $[eqn]$ , the Tx transmits omnidirectionally; otherwise, when $[eqn]$ , the Tx employs singular value decomposition (SVD)-based precoding to obtain the transmitted signal. In what follows, we will evaluate both the probability of direct decoding using IMUD and BER performance under VSP-IMUD. In our simulations, we assume perfect CSI for simplicity. When imperfect CSI is considered, $[eqn]$ can be re-modeled as $[eqn]$ [18], where $[eqn]$ and $[eqn]$ denote accurate and inaccurate channel matrices, respectively. The coefficient $[eqn]$ indicates the degree of CSI imperfection, with $[eqn]$ representing perfect CSI. The matrix $[eqn]$ is an $[eqn]$ complex diagonal Gaussian matrix with zero mean and unit variance. It is important to note that all of the investigated methods are susceptible to inaccuracies in CSI. However, since the modeling of CSI imperfections and their influence on data transmission, as well as efforts to improve the algorithms’ robustness to these imperfections, have been extensively studied [19], we have omitted this discussion in the present work.

Figure 4 plots the variation in the probability of directly applying IMUD with K under $[eqn]$ , $[eqn]$ , and various $[eqn]$ s. As the figure shows, for a fixed K, the probability of direct decoding with IMUD increases as $[eqn]$ grows. When $[eqn]$ is held constant, the probability increases as K reduces. This trend is due to the fact that a reduction in either $[eqn]$ or K decreases the detection probability of channel ambiguity, which in turn increases the likelihood of directly applying IMUD. Conversely, if $[eqn]$ or K increases, the detection probability of channel ambiguity rises, thereby reducing the probability of direct decoding with IMUD. The probability curve plotted in Figure 4 indicates that IMUD can be directly applied when no channel ambiguity is detected. Under very small $[eqn]$ (e.g., $[eqn]$ ), the probability of directly using IMUD remains unchanged as K varies, as no channel ambiguity is detected.

Variation of the probability of direct applying IMUD with K under NT=1, NR=2, and various εs.

We define the transmit power normalized by noise power $[eqn]$ as $[eqn]$ , and set $[eqn]$ in the simulation. Figure 5 simulates the impact of SIC on the BER of signal components under $[eqn]$ , $[eqn]$ , and $[eqn]$ , where the Rx does not apply IMUD to decode the Txs’ data. Given $[eqn]$ , the Rx must perform SIC five times to decode all $[eqn]$ signals until the last two signal components can be directly recovered. As illustrated in the figure, due to the error propagation issue associated with SIC, the BER for later decoded data is worse than that for previously decoded data.

Figure 6 displays the BER of VSP-IMUD under $[eqn]$ using various referential symbol selection criteria. We use $[eqn]$ to indicate that the data transmitted via the channel with highest channel gain is selected as the referential symbol to determine the final symbol set. For comparison, we denote random referential symbol selection as random. As the figure shows, for fixed K, the BER using the $[eqn]$ criterion outperforms the random referential symbol selection. This is because random selection is more likely to produce referential symbol errors, leading to an incorrect final symbol set.

Figure 7 simulates the system BER of VSP-IMUD under $[eqn]$ , $[eqn]$ , $[eqn]$ , and various $[eqn]$ s. The figure shows that the BER performance improves as $[eqn]$ increases. This improvement is due to the enhanced signal processing gain achieved by the Tx-side array when equipped with more antennas, thereby improving the overall system BER.

Figure 8 compares the system BER of VSP-IMUD, ICom/SIC-ZF, and minimum mean square error (MMSE)-SIC under $[eqn]$ , $[eqn]$ , $[eqn]$ , and $[eqn]$ . We assume $[eqn]$ in the simulations. In this scenario, conventional ZF and MMSE become inapplicable. To make MMSE applicable when $[eqn]$ , we treat $[eqn]$ signal components as an equivalent signal, similar to ICom/SIC-ZF. As the figure shows, under VSP-IMUD, the system BER decreases as K increases. This occurs because, as K rises, the probability of channel ambiguity occurring also slightly increases. In this situation, although the Rx needs to employ ZF to decode the signals transmitted via the ambiguous channels and subsequently use SIC to remove the detected signals from the mixed signal, the number of SIC applications remains small. That is, the increase in the probability of needing SIC is limited as K rises, which minimally contributes to the degradation of system BER. Furthermore, as K increases, the referential symbol can be selected from a larger set of candidate signals, improving the correctness of the referential symbol due to selective diversity gain. Thus, the BER performance of VSP-IMUD improves as K increases. In contrast, ICom/SIC-ZF and MMSE-SIC require more SIC applications as K grows, resulting in degraded BER performance as K increases. Furthermore, VSP-IMUD achieves the lowest system BER, followed by MMSE-SIC, while ICom/ZF-SIC exhibits the worst system BER. This is because the reduced reliance on SIC in VSP-IMUD enhances the decoding of multiple user signals, resulting in superior system BER compared to the other SIC-based schemes. Moreover, MMSE-SIC outperforms ICom/ZF-SIC in the medium- $[eqn]$ regime and converges to the performance of ICom/ZF-SIC at low- and high- $[eqn]$ regimes. This is because, when $[eqn]$ is low, noise is the dominant factor affecting the system BER; thereby, even MMSE outperforms ZF in balancing interference suppression and desired signal preservation, and both ICom/ZF-SIC and MMSE-SIC yield similar BER performance. As $[eqn]$ increases, CCI dominates the system BER; therefore, MMSE-SIC outperforms ICom/ZF-SIC in system BER. When $[eqn]$ becomes too large, ZF and MMSE exhibit similar interference suppression performance; thus, ICom/ZF-SIC and MMSE-SIC yield approximately the same system BER.

Figure 9 plots the comparison of system SE of VSP-IMUD, ICom/ZF-SIC and MMSE-SIC under $[eqn]$ , $[eqn]$ , $[eqn]$ , and $[eqn]$ . As the figures show, VSP-IMUD achieves the highest system SE, followed by MMSE-SIC, while ICom/ZF-SIC exhibits the lowest system SE. The analysis is shown in Figure 7. Moreover, MMSE-SIC outperforms ICom/ZF-SIC in the low- $[eqn]$ regime and converges to the performance of ICom/ZF-SIC as $[eqn]$ increases. This observation is consistent with the characteristics of MMSE and ZF.

5. Conclusions

In this paper, we propose a VSP-IMUD-based reception scheme by exploiting the interactions among multiple wireless signals. By virtually combining multiple signal components, we can obtain an equivalent structured signal, and then employ IMUD to decode this effective signal based on its spatial characteristics. Before applying IMUD, we must detect and mitigate channel ambiguity for each signal component to ensure the uniqueness of the equivalent signal. Ultimately, we select the signal with the highest channel gain to decode a referential symbol. With this referential symbol, the Rx can output the final symbol set. Compared to other existing MUD methods, such as ICom/SIC-ZF, VSP-IMUD can effectively avoid the error propagation problem by reducing the number of SIC applications, thereby achieving improved BER performance.