[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
Multi-Task Visual Perception for Object Detection and Semantic Segmentation in Intelligent Driving
Previous Article in Journal
3D Point Cloud Shape Generation with Collaborative Learning of Generative Adversarial Network and Auto-Encoder
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sparse Reconstruction-Based Joint Signal Processing for MIMO-OFDM-IM Integrated Radar and Communication Systems

1
National Key Laboratory of Radar Signal Processing, Xidian University, Xi’an 710071, China
2
Department of Electrical and Computer Engineering, National University of Singapore, Singapore 119077, Singapore
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(10), 1773; https://doi.org/10.3390/rs16101773
Submission received: 29 March 2024 / Revised: 4 May 2024 / Accepted: 13 May 2024 / Published: 16 May 2024
Graphical abstract
">
Figure 1
<p>Schematic diagram of the proposed MIMO-OFDM integrated system model.</p> ">
Figure 2
<p>Schematic diagram of MIMO-OFDM-IM integrated signal design in transmitter.</p> ">
Figure 3
<p>Flow chart of the joint signal processing schemes in the proposed system.</p> ">
Figure 4
<p>Flowchart of the conventional processing.</p> ">
Figure 5
<p>The PSLR performance. (<b>a</b>) The variation in PSLR with SNR. (<b>b</b>) The variation in PSLR with the number of activated subcarriers (SNR = 10 dB).</p> ">
Figure 6
<p>The variation in PSLR with SNR under different number of activated subcarriers.</p> ">
Figure 7
<p>The variation in MSEs of range and velocity estimation with SNR under different methods (<math display="inline"><semantics> <msub> <mi>N</mi> <mi>r</mi> </msub> </semantics></math> = <math display="inline"><semantics> <msub> <mi>N</mi> <mi>t</mi> </msub> </semantics></math> = 4, <math display="inline"><semantics> <msub> <mi>N</mi> <mrow> <mi>c</mi> <mo>_</mo> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> </semantics></math> = 8, <span class="html-italic">K</span> = 4, <span class="html-italic">H</span> = 4). (<b>a</b>) Range estimation. (<b>b</b>) Velocity estimation.</p> ">
Figure 8
<p>Complexity comparison. (<b>a</b>) The variation in number of complex multiplications with the number of activated subcarriers. (<b>b</b>) The variation in number of complex multiplications with the number of transmitting antennas.</p> ">
Figure 9
<p>The BER performance of different information detection methods under the proposed integrated MIMO-OFDM-IM system (<math display="inline"><semantics> <msub> <mi>N</mi> <mi>r</mi> </msub> </semantics></math> = <math display="inline"><semantics> <msub> <mi>N</mi> <mi>t</mi> </msub> </semantics></math> = 4). (<b>a</b>) <math display="inline"><semantics> <msub> <mi>N</mi> <mrow> <mi>c</mi> <mo>_</mo> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> </semantics></math> = 8, <span class="html-italic">K</span> = 3; (<b>b</b>) <math display="inline"><semantics> <msub> <mi>N</mi> <mrow> <mi>c</mi> <mo>_</mo> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> </semantics></math> = 8, <span class="html-italic">K</span> = 6.</p> ">
Figure 10
<p>The BER performance of different information detection methods under the proposed integrated MIMO-OFDM-IM system (<math display="inline"><semantics> <msub> <mi>N</mi> <mi>r</mi> </msub> </semantics></math> = <math display="inline"><semantics> <msub> <mi>N</mi> <mi>t</mi> </msub> </semantics></math> = 8). (<b>a</b>) <math display="inline"><semantics> <msub> <mi>N</mi> <mrow> <mi>c</mi> <mo>_</mo> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> </semantics></math> = 8, <span class="html-italic">K</span> = 3; (<b>b</b>) <math display="inline"><semantics> <msub> <mi>N</mi> <mrow> <mi>c</mi> <mo>_</mo> <mi>s</mi> <mi>u</mi> <mi>b</mi> </mrow> </msub> </semantics></math> = 8, <span class="html-italic">K</span> = 6.</p> ">
Figure 11
<p>The variation in radar–communication trade-off curve with the number of activated subcarriers under different SNRs (<math display="inline"><semantics> <msub> <mi>N</mi> <mi>r</mi> </msub> </semantics></math> = <math display="inline"><semantics> <msub> <mi>N</mi> <mi>t</mi> </msub> </semantics></math> = 4). (<b>a</b>) SNR = 6 dB; (<b>b</b>) SNR = 12 dB.</p> ">
Versions Notes

Abstract

:
Multiple-input multiple-output orthogonal frequency division multiplexing (MIMO-OFDM) technology is widely used in integrated radar and communication systems (IRCSs). Moreover, index modulation (IM) is a reliable OFDM transmission scheme in the field of communication, which transmits information by arranging several distinguishable constellations. In this paper, we propose a sparse reconstruction-based joint signal processing scheme for integrated MIMO-OFDM-IM systems. Combining the advantages of MIMO and OFDM-IM technologies, the integrated MIMO-OFDM-IM signal design is realized through the reasonable allocation of bits and subcarriers, resulting in better intercarrier interference (ICI) resistance and a higher transmission efficiency. Taking advantage of the sparseness of OFDM-IM, an improved target parameter estimation method based on sparse signal reconstruction is explored to eliminate the influence of empty subcarriers on the matched filtering at the receiver side. In addition, an improved sequential Monte Carlo signal detection method is introduced to realize the efficient detection of communication signals. The simulation results show that the proposed integrated system is 5 dB lower in the peak sidelobe ratio (PSLR) and 1.5 × 10 5 lower in the number of complex multiplications than the latest MIMO-OFDM system and can achieve almost the same parameter estimation performance. With the same spectral efficiency, it has a lower bit error rate (BER) than existing methods.

Graphical Abstract">

Graphical Abstract

1. Introduction

Integrated radar and communication technology has recently seen rapid development and applications in both military and civilian arenas. This technology uses a single transmitted signal to simultaneously perform radar detection and communication transmission functions, which reduces the hardware, energy consumption, and electromagnetic interference of the mobile platform and also offers higher energy efficiency and spectral efficiency compared to a single system [1,2,3]. Currently, there are four main application scenarios for integrated radar and communication systems (IRCSs), including unmanned aerial vehicle (UAV) patrol and inspection systems to ensure the safety of the natural environment [4], civil aviation air traffic control systems with automatic real-time tracking [5], military combat systems with real-time target tracking and information sharing [6], and collaborative intelligent transportation systems with environmental awareness and real-time communication [7]. It is worth noting that IRCSs have significant information security issues in different application areas. However, communication information needs to be processed by source coding, channel coding, and information encryption before transmission, which can ensure the accuracy and security of information transmission. It is also very important to optimize transmitted waveforms and further improve received signal processing for better interoperability and efficiency [8,9,10,11].
Orthogonal frequency division multiplexing (OFDM) is a special multicarrier modulation technique that has gained wide attention in the field of digital communications and radar sensing [12,13,14]. This technology can realize high-speed data transmission by virtue of effectively overcoming intersymbol interference, strong resistance to multipath fading, and high spectrum utilization. Several in-depth studies on integrated waveform design and received signal processing schemes based on OFDM signals can be found in the literature [15,16,17,18,19,20]. Index modulation, on the other hand, is a reliable modulation technology with high spectrum and energy efficiencies. It transmits information using embedded indexes of building blocks of antennas, subcarriers, spreading codes, and time slots [21,22]. OFDM with index modulation (OFDM-IM) is introduced to take advantage of both technologies. It can flexibly adjust the number and location of activated subcarriers to achieve a favorable compromise between system performance, spectrum utilization, and energy efficiency [23,24]. Because of its flexible design, OFDM-IM has attracted great attention [25,26]. In addition, the multiple-input multiple-output (MIMO) technology is one of the key technologies for future mobile and wireless communication systems, which is characterized by obtaining both reliability and effectiveness gains through the use of spatial resources on the basis of the full utilization of existing spectrum resources. Therefore, MIMO-OFDM technology is also considered as one of the key technologies in 5G or even 6G digital communication [27,28,29].
By combining the advantages of MIMO technology and OFDM-IM technology, the system proposed in [30] adopts a hybrid analog–digital (HAD) beamforming structure based on MIMO-OFDM-IM, which can be effectively applied to the millimeter wave field. In [31], a feasible solution is proposed by using a Pareto waveform optimization design to approximate the global Pareto optimal compromise between spectrum efficiency and energy efficiency, and it can effectively solve the variational constraint problem in combinatorial optimization models. Based on the signal sparsity characteristics of the MIMO-OFDM-IM system, a low-complexity signal detector based on deep learning (DL) is proposed in [32], which firstly performs the feature information extraction by the acquired signals of the two subnetworks and the DL method, and then the optimized least squares method is used to demodulate the required communication information. It improves the bit error performance and reduces the computational complexity. Paper [33] demonstrated the effect of different signal detection methods on the transmission error bit rate of MIMO-OFDM-IM systems through the derivation of theoretical equations and puts forward its adaptability to visible light communication systems in another reference. In the proposed scheme, the spectral efficiency of the system is greatly improved due to the fact that the OFDM-IM waveforms emitted by each transmitting antenna do not interfere with each other. The method in [34] designs a novel MIMO-OFDM-IM waveform by deactivating a subset of subcarriers. Compared with the conventional MIMO-OFDM waveforms, the subcarrier design of this MIMO-OFDM-IM waveform is more flexible and can achieve a better BER performance, which leads to higher energy efficiency for practical systems. However, the detection of the MIMO-OFDM-IM communication signal becomes a very challenging problem due to the mutual interference between multiple channels and the correlation between different subcarriers of the same subblock [35,36]. In addition, when MIMO-OFDM-IM is applied to radar, it is found that the traditional signal processing method has the problem of dynamic range decline [37,38]. Therefore, under the integrated MIMO-OFDM-IM system, it is urgent to develop more waveform design schemes and corresponding signal processing schemes to improve the system performance.
In this paper, MIMO-OFDM-IM technology is applied to the designed IRCS. After realizing the integrated signal design based on MIMO-OFDM-IM, an improved signal processing scheme is proposed. Firstly, according to the technical characteristics of MIMO-OFDM and IM, an integrated pulse signal design scheme based on MIMO-OFDM-IM is proposed. At the radar receiver, the communication information and phase difference information of the multipulse echo signal are compensated, and the sparse signal model is obtained by screening the effective data. In addition, a segmented coherent accumulation method based on sparse signal reconstruction is adopted to achieve target parameter estimation within the detection range. At the communication receiver, an improved communication signal detection method based on sequential Monte Carlo is explored. To this end, we summarize our contributions as follows:
(1) The integrated signal design adopts the mode of multiple OFDM symbols in one pulse, which improves the communication transmission rate under the same bandwidth, and IM technology increases the subcarrier frequency spacing, thus having a better intercarrier interference (ICI) resistance capability.
(2) In the integrated MIMO-OFDM-IM systems, an improved method based on information compensation and sparse signal reconstruction is proposed to estimate the range and velocity of radar targets. Compared with the existing methods, it has a better parameter estimation performance.
(3) According to the subcarrier index, we introduce an improved communication signal detection method based on SMC in the proposed system. Compared with the existing methods, the improved method can achieve a lower communication error rate with the same spectral efficiency while ensuring a low computational complexity.
The rest of this paper is organized as follows. In Section 2, we introduce the system model and signal design. In Section 3, the optimized radar and communication signal processing schemes are proposed. In Section 4, we use simulation results to verify the performance of the system. Finally, Section 5 summarizes the main work of this paper.

2. System Model and Signal Design

2.1. System Model

The application scenario we consider is an integrated MIMO-OFDM system with N t transmitting antennas and N r receiving antennas, which simultaneously transmit integrated signals to multiple targets and multiple communication users. Each antenna transmits integrated signals in the form of pulses, and each pulse contains multiple OFDM symbols. This transmission mode improves the communication data rate under the same bandwidth. The model of the proposed integrated platform is shown in Figure 1.
The transmitted signals of the different antennas in Figure 1 are separated in the frequency domain and only occupy the pre-set continuous frequency spectrum of the whole system’s working frequency band. Therefore, the waveforms transmitted by each antenna are orthogonal, which meets the requirements of the MIMO radar for the orthogonality of transmitted waveforms. It is worth mentioning that in order to obtain the subcarrier complex valued weight sequence with the nonoverlapping frequency domain on each transmitting antenna, we design the subcarrier complex weight sequence by interleaving the structure [39,40].
Next, we apply IM technology to the above integrated MIMO-OFDM system and explore the integrated MIMO-OFDM-IM signal design method. At the transmitter side, the integrated signal design scheme is divided into two stages: the first is bit and subcarrier allocation and the second is the integrated signal design, which is shown in Figure 2.

2.2. Bit and Subcarrier Allocation

The main principle of OFDM-IM is to selectively activate some OFDM subcarriers according to the bit stream information and only carry out traditional constellation modulation on the activated subcarriers while the unactivated subcarriers are directly set to zero. After the receiver demodulates the OFDM signal, it can obtain the constellation information on the activated subcarrier, as well as obtain additional index information according to the index position of the activated subcarrier [41]. It is worth noting that the IM equivalently increases the subcarrier frequency spacing, resulting in better ICI resistance.
The allocations of bits and subcarriers are as follows. Assume that the number of bits to be transmitted by the transmitter is N and the bits carried by each antenna are N 1 = N / N t . In addition, to decrease mutual interference between subcarriers, we divide all subcarriers into I subblocks. It is assumed that each antenna has N c available subcarriers, and then each subblock contains N c _ s u b = N c / I subcarriers and N 2 = N 1 / I bits.
If K activations are selected from each subblock, there are a total of C N c _ s u b K combinations of activated subcarriers, and C is the number of combinations. Therefore, each subblock can transmit p 1 = log 2 C N c _ s u b K bits of index information, where · means rounding down. Meanwhile, the indexes of K activated subcarriers are stored in matrix A , and the rest are saved in matrix B . In addition, the constellation modulation of M points is performed on the activated subcarriers so that the number of bits of modulation information in each subblock is p 2 = K log 2 M . Suppose that the number of OFDM symbols contained in each pulse signal is N s , and the modulation information on the k th subcarrier in the m th OFDM symbol is denoted as d M ( m , k ) , m = 1 , 2 , , N s ; then, the data symbol on the k th subcarrier, m th OFDM symbol, ith subblock, and gth antenna can be formulated as
d g i ( m , k ) = d M ( m , k ) k A 0 k B
The symbol matrix generated by the i th subblock of g th antennas can be written as
d g i = d g i ( 1 , 1 ) d g i ( 1 , 2 ) d g i 1 , N c _ s u b d g i ( 2 , 1 ) d g i ( 2 , 2 ) d g i 2 , N c _ s u b d g i ( N s , 1 ) d g i ( N s , 2 ) d g i ( N s , N c _ s u b )
where i = 1 , 2 , , I and g = 1 , 2 , , N t . Each subblock is processed in the same way as described above, and then all the subblocks are concatenated to obtain the OFDM-IM signal block for each antenna, and the obtained bit data matrix on the gth antenna can be represented as
d g = d g 1 , d g 2 , , d g I , g = 1 , 2 , , N t
Assume that the duration of OFDM symbols is T b . The complex envelope of a pulse signal transmitted on the g th antenna [15] can be written as
s g ( t ) = m = 0 N s 1 k = 0 N c 1 d g ( m , k ) e j 2 π k Δ f t m T b rect t m T b T b
where Δ f = 1 / T b is the frequency interval of neighboring subcarriers, d g ( m , k ) is the element of the mth row and kth column of d g , and rect ( · ) represents the rectangular window function defined on [0,1].

2.3. Integrated Signal Design

Next, we perform a I × N c _ s u b block interleaver and inverse fast Fourier transform (IFFT) on all transmit antennas. In addition, in order to eliminate intersymbol interference (ISI), the guard intervals (GI) sequence should be inserted before all OFDM signal blocks [42]. The traditional way is to directly add a cyclic prefix (CP), but the introduction of a cyclic prefix will produce a high CP sidelobe in radar signal processing, which will affect the target detection performance [43]. Therefore, according to the performance standard of the proposed system, this paper adopts the nonlinear optimization method in [44] to obtain the optimal GI sequences of all OFDM symbols.
Then, after conventional OFDM operations such as serial-to-parallel conversion, the signal is simultaneously transmitted from N t transmitting antennas via the frequency-selective Rayleigh fading MIMO channel.
According to the above method, all OFDM signals of P pulses are generated, and the total number of activated subcarriers is K = K · I · P . Then, the indexes of the activated subcarriers are arranged in the vector ξ in ascending order, which is expressed as
ξ = ξ 1 ξ 2 ξ K
where ξ 1 , ξ 2 , ⋯, ξ K denote the first to K activated subcarrier indexes, respectively. According to the signal generation flow in Figure 2, the spectral efficiency (in bps/Hz) of the proposed system is given by the following formula:
E = N t p 1 + p 2 I N c + N c p = N t log 2 C N c _ s u b K + K log 2 M I N c + N c p
It is worth mentioning that the balance among spectrum utilization, transmission reliability, and complexity can be achieved by adjusting relevant parameters in practical applications.

3. Optimized Joint Signal Processing Schemes

In this section, we propose the optimized radar signal processing and communication signal processing schemes. At the radar receiving end, an improved target range and velocity estimation method based on information compensation and sparse signal reconstruction is introduced. At the communication receiving end, an improved low-complexity detection method based on sequential Monte Carlo is explored. Figure 3 shows the process of the proposed schemes. This figure displays the progressive relationship between different notations, which will be introduced in the next subsection.

3.1. Sparsity Reconstruction-Based Improved Parameter Estimation

For target echo signals, compared with the MIMO-OFDM system, the subcarriers in the proposed integrated MIMO-OFDM-IM system are sparse, so the variance of the ambiguity function of the integrated signal will increase, which leads to the increase in the sidelobe level obtained by the traditional radar signal processing method. In the theoretical system of compressive sensing, in order to express a signal more accurately, the signal can usually be transformed into a new basis or framework. As long as we know in advance that the signal has sparse expression under this framework, we can reconstruct the transmitted signal. Therefore, we propose a target parameter estimation method based on information compensation and sparse signal reconstruction.

3.1.1. Conventional Processing

The conventional processing of the radar echo signal in the proposed integrated system is shown in Figure 4.
Supposing that there are Q detectable targets in the application scenario, the distance between the radar transmitter and the q th target is R q and the relative velocity of the q th target moving toward the radar is v q . It can be obtained that the time delay of the q th target is τ q = 2 R q / c ; the Doppler frequency shift is f d q = 2 f c v q / c , where f c is the carrier frequency; and c is the propagation velocity. For simplicity, we suppose that the direction of arrival at the receiving end of the integrated platform is accurately known, and the echo is processed in units of pulses. Then, we can perform digital beamforming; that is, use the aperture of the array antenna to form the receiving beam in the desired direction through digital signal processing so as to form the equivalent single-transmit waveform for the colocated MIMO radar at the receiver [4]. The echo signal of the pth pulse can be formulated as
y p ( t ) = a r H ( θ r ) [ y p , 1 ( t ) , y p , 2 ( t ) , , y p , N r ( t ) ] T = N r q = 1 Q β q e j 2 π f c τ q e j 2 π f d q p T p a t T ( θ t ) s p ( t τ q ) + ω p ( t )
where a r ( θ r ) denotes the receive steering vector, β q is the scattering factor of the q th target, T p is the pulse repetition period of the transmitted waveform, s p ( t ) is the transmitted signal vector of the p th pulse, ω p ( t ) is Gaussian noise, and a t θ t indicates the transmit steering vector, and it can be written as
a t θ t = 1 , e j 2 π d t sin θ t / λ , , e j 2 π N t 1 d t sin θ t / λ T
where d t is the transmitting antenna spacing and λ = c / f c is the signal wavelength.
Then, the radar-received signal in (7) is sampled at the sampling frequency f s = N c Δ f , and the discrete format of the signal in (7) can be denoted as
y p = N r q = 1 Q β q e j 2 π f d q p T p F N c 1 × ψ 1 s ˜ p τ q , ψ 2 s ˜ p τ q , , ψ N t s ˜ p τ q × a t θ t + ω p
where ψ g indicates the diagonal matrix of the g th channel set for realizing the subcarrier interleaving structure, ⊙ indicates the Hadamard product, s ˜ p indicates the frequency domain weight vector in the p th pulse, and F N c 1 indicates an inverse discrete Fourier transform (IDFT) matrix, which can be written as
F N c 1 = 1 N c e j 2 π ( i 1 ) ( j 1 ) N c N c × N c , i , j = 1 , 2 , , N c
In addition, τ q is the delay vector of the q th target, which is given by
τ q = e j 2 π f c τ q , e j 2 π f c + Δ f τ q , , e j 2 π f c + N c 1 Δ f τ q T
To minimize the influence of the phase difference information of different channels on the signal processing results, we need to compensate for the received signals of different channels. However, the phase difference information of different channels is different, so it is necessary to separate the channels before compensation. Therefore, according to the structural characteristics of the MIMO-OFDM-IM system, we use the diagonal matrix group in (9) to extract data from different channels through staggered allocation. The received signal of the pth pulse on the g th channel after separation can be represented as
y p , g = N r q = 1 Q β q e j 2 π f d q p T p e j 2 π ( g 1 ) d t sin θ t / λ F N c 1 ψ g s ˜ p τ q + ω ˜ p , g
Next, we perform N c -points DFT on the signal in (12) to obtain the following equation
y ˜ p , g = N r q = 1 Q β q e j 2 π f d q p T p e j 2 π ( g 1 ) d t sin θ t / λ ψ g s ˜ p τ q + ω ^ p , g
where ω ^ p , g = ω ˜ p , g F N c .

3.1.2. Information Compensation

According to the ambiguity function theory of integrated signal in [42], it is known that communication information will affect the echo matching filtering effect, thus reducing the radar detection performance. Based on this, it is necessary to compensate for the signal in (13) in the frequency domain, and the signal compensated by s ˜ p 1 can be rewritten as
y ^ p , g = N r q = 1 Q β q e j 2 π f d q p T p e j 2 π ( g 1 ) d t sin θ t / λ ψ g τ q + W p , g
By stacking all OFDM signals into a matrix, the following formula can be given as
Y g = N r e j 2 π ( g 1 ) d t s i n θ t / λ ψ g × q = 1 Q β q e j 2 π f d q T p τ q , q = 1 Q β q e j 2 π f d q 2 T p τ q , , q = 1 Q β q e j 2 π f d q P T p τ q + [ W 1 , g , W 2 , g , , W P , g ] = N r e j 2 π ( g 1 ) d t s i n θ t / λ ψ g q = 1 Q β q τ q D q T + W g
where W g = W 1 , g , W 2 , g , , W P , g , ⊗ denotes the Kronecker product of a matrix, and D q is the Doppler frequency shift vector of the q th target, which can be written as
D q = [ e j 2 π f d q T p , e j 2 π f d q 2 T p , , e j 2 π f d q P T p ] T
Next, each channel is compensated for phase difference information, and (15) can be rewritten as
Y = g = 1 N t Y g e j 2 π ( g 1 ) d t s i n θ t / λ = N r q = 1 Q β q τ q D q T + W ˜

3.1.3. Improved Parameter Estimation

Generally speaking, the traditional method based on the modulation symbol domain can be used to realize target parameter estimation; that is, we can perform DFT on the signal in (17) along the pulse symbol dimension and perform an IDFT along the subcarrier dimension. Based on this, (17) can be updated to the following formula:
Y ^ = F N c 1 Y F P = N r q = 1 Q β q F N c 1 τ q D q T F P + W ^
where F P is a discrete Fourier transform (DFT) matrix and its size is P × P . However, due to the use of IM, there may be zero elements in the transmission symbols, which will cause some elements in the matrix Y to converge to infinity. Therefore, it is very difficult to obtain a good parameter estimation performance by directly adopting the method in (18). We need to extract valid data through data filtering and use new methods to solve it.
According to the relationship between the vectorization of a matrix product and Kronecker product [45], (17) can be transformed into the following formula:
vec ( Y ) = vec N r q = 1 Q β q τ q D q T + vec ( W ˜ )
where vec ( · ) indicates that the matrix is stacked in a vector form by columns. Suppose that R q = F N c 1 τ q and V q = D q T F P . According to the properties of the Kronecker product, (19) can be rewritten as
vec ( Y ) = vec N r q = 1 Q β q F N c F N c 1 τ q D q T F P F P 1 + vec ( W ˜ ) = F P H F N c vec N r q = 1 Q β q R q V q + vec ( W ˜ )
where R q and V q contain the range and velocity information of the q th target, respectively. According to the index vector ξ of the activated subcarrier in (5), the position of valid data in vec ( Y ) can be determined. Based on this, we construct a matrix C with the size of K × N c P to filter the data in (20), and its elements can be obtained by the following formula
C i , j = 1 j = ξ ( i ) 0 otherwise
According to the characteristics of the matrix and vector, we can filter out invalid data by multiplying (20) by C and extract valid data to realize the sparse reconstruction of the echo signal. After the sparse signal reconstruction of (20), the echo signal model can be reformulated as
Y ˜ = C ˜ vec N r q = 1 Q β q R q V q + W ˜
where
Y ˜ = C · vec ( Y )
C ˜ = C · ( F P H F N c )
W ˜ = C · vec ( W ˜ )
It can be seen that (22) is a sparse model of received signals corresponding to all targets. In (22), a large matrix dimension will lead to high algorithm complexity. To reduce the complexity of the algorithm, we adopt a segmented coherent accumulation method based on the sparse signal model, and the specific implementation process is as follows:
Assuming that the total number of pulses is P, and it is divided into H equal parts on average, then the number of pulses in each pulse group is P ^ = P / H . For the h th pulse group data, (22) can be transformed into the following formula:
Y ˜ h = C ˜ h χ h , R , v + W ˜ h
where C ˜ h = C h ( F P H F N c ) , C h is similar to C , its size is K × ( N c P / H ) , and χ h , R , v is a two-dimensional matrix containing range and velocity information, which can be written as
χ h , R , v = N r q = 1 Q β q R q ( V q , h )
where
V q , h = [ e j 2 π f d q T p , e j 2 π f d q 2 T p , , e j 2 π f d q P ^ T p ] F P ^
According to the sparsity of χ h , R , v , (26) can be further equivalently transformed into an optimization problem, which can be described as
f o p t = arg min χ h 1 2 Y ˜ h C ˜ h χ h , R , v 2 2 + μ χ h , R , v 1
A large amount of literature points out that there are many methods for solving the unconstrained optimization problem in (29). In this paper, an improved method with low complexity proposed in [46] is used to solve this problem.
After solving the above problem, the matrix χ h , R , v corresponding to the h th pulse group can be obtained, and then the range dimension information of the pulse group can be obtained by performing an IDFT on this matrix along the pulse group dimension, which can be represented as
ψ h , R = χ h , R , v F P ^ 1
The H pulse groups are processed according to the above method, and all the results are then recombined. Then, DFT is performed on the reconstructed data, and finally the two-dimensional image information of multiple targets is obtained, which can be denoted as
χ ^ R , v = ψ 1 , R , ψ 2 , R , , ψ H , R F P = χ 1 , R , v F P ^ 1 , χ 2 , R , v F P ^ 1 , χ H , R , v F P ^ 1 F P
To measure the estimation accuracy of the proposed method based on MIMO-OFDM-IM, we use the Cramér–Rao bound (CRB) as a standard to evaluate all unbiased estimators. Combined with the detailed derivation process in [15], the CRBs of range and velocity under the proposed system are obtained, which can be formulated as
CRB R = I v v c 2 4 I v v I τ τ I v τ I τ v
CRB v = I τ τ I v v I τ τ I v τ I τ v
where I v v , I τ τ , I v τ , and I τ v are the corresponding elements in the Fisher information matrix.
In addition, for the proposed signal design scheme, the range resolution can be expressed as
R r e s = c 2 B = c 2 N c Δ f
where c denotes the signal propagation velocity and B denotes the integrated signal bandwidth. In addition, we can effectively determine the maximum detectable range of the proposed system in the environment through the spacing of subcarriers, which can be written as
R max = c 2 Δ f
Similarly, the velocity resolution and the maximum detectable velocity of the system are expressed as
v r e s = c 2 f c P T p
v max = c 2 f c T p
According to (35), we can realize a larger detection range of the system by reducing the subcarrier spacing. However, (37) shows that the subcarrier spacing should be increased for a higher detectable velocity. Therefore, it is necessary to choose the appropriate subcarrier spacing according to different requirements in practical applications.
Next, we discuss the computational complexity of the proposed parameter estimation method. First, define one complex multiplication operation as a unit of computational complexity; it is certain that the greater the number of complex multiplications, the higher the computational complexity. According to the processing flow of the above algorithm, the computational complexity of the direct sparse reconstruction of P pulse data is O ( Q K I N c P 2 ) , while the computational complexity of the proposed segmented coherent accumulation method is H · O ( Q K I N c P 2 / H 2 ) + P 2 · O ( N c lg N c ) .

3.2. Improved Communication Signal Detection

In communication applications, the received signal is the superposition of the transmitted signal polluted by Gaussian noise. The conventional processing of posterior distribution is to search all possible transmitted signal blocks exhaustively, but its complexity is high in the case of large signal blocks. Sequential Monte Carlo (SMC) is an inference method for estimating nonlinear and non-Gaussian state space models, which has been widely used in wireless communication [47,48]. The SMC method approximates the posterior distribution of a state variable by generating high-dimensional Monte Carlo samples dimension by dimension. Specifically, symbol samples are extracted from a defined set in a certain order, and new sequential particles are created by the extracted samples. Then, the importance weight associated with each particle is updated by the probability distribution characteristics. The research shows that the detector based on SMC can achieve a near-maximum likelihood performance while the cost of the receiver is lower. However, the basic SMC can not fully meet the application requirements, so it is necessary to study improved methods suitable for different scenarios. In this section, to reduce the complexity and obtain a better communication performance, we propose an improved SMC signal detection method.
For the communication receiver, we first need to achieve high-precision communication synchronization according to the cross-correlation between all GI sequences and the received communication signal. Then, the interference cancellation method in [49] is adopted to effectively reduce the interference between signals. Meanwhile, the position information of the OFDM signal block carrying modulation symbols is extracted by using the relative position relationship between sequences. All GI sequences will be removed after synchronization is completed. In addition, N c -point fast Fourier transform (FFT) is adopted at each receiving antenna, and the signal block in the frequency domain is obtained by the I × N c _ s u b block deinterleaver. Specifically, for the subcarriers of each OFDM-IM subblock, the received signal vectors of the ith subblock can be stacked as follows:
f i m , k = α H i m , k T d i m , k + ω i m , k
where f i m , k = f 1 i m , k f 2 i m , k f N r i m , k T is the received signal vector, α is the balance factor, ω i m , k is the Gaussian noise vector with zeros mean and unit variance, d i m , k = d 1 i m , k d 2 i m , k d N t i m , k T is the modulated transmission signal vector, and H i m , k is the channel frequency responses matrix, which can be written as
H i m , k = h 1 , 1 i m , k h 1 , 2 i m , k h 1 , N r i m , k h 2 , 1 i m , k h 2 , 2 i m , k h 2 , N r i m , k h N t , 1 i m , k h N t , N r i m , k h N t , N r i m , k
To reduce the complexity, we propose an improved communication signal detection method based on SMC. Without changing the noise-whitening characteristics, we need to perform matched filtering and noise whitening on the received signal vector to maximize the signal-to-noise ratio (SNR). The traditional deterministic SMC algorithm can be used to obtain the noise-whitening matrix [47], which can be written as
U i m , k = H i m , k H i m , k T 1 2
According to the theoretical basis of the continuous interference cancellation method, we perform QR decomposition on the noise whitening matrix obtained in (40), which can be specifically represented as
U i m , k = Q i m , k R i m , k
where Q i m , k denotes a unitary matrix and R i m , k denotes an upper triangular matrix. After matched filtering, noise whitening, and triangulation of the received signal vector, (38) can be transformed into the following formula:
f ^ i m , k = Q i m , k H U i m , k 1 H i m , k f i m , k = Q i m , k H H i m , k H i m , k T 1 2 1 H i m , k × α H i m , k T d i m , k + ω i m , k = α Q i m , k H H i m , k H i m , k T 1 2 d i m , k + ω ^ i m , k = α R i m , k d i m , k + ω ^ i m , k
Due to the high complexity of directly calculating the posterior probability of each subblock, we adopt a sequential iterative method to numerically approach the true posterior probability infinitely. According to the structure in (42), the symbolic information detection is performed in the order of the antenna index. Since the states of the subcarriers corresponding to each transmitting antenna are independent, the SMC method can iteratively enumerate the symbolic information in the order of the antenna index.
The traversal process of symbol information starts from the first subcarrier, so the number of subcarrier state combinations l s u b and the number of samples ε will increase continuously. We set the lower limit value of the number of state combinations as L s u b and the upper limit of the number of samples as ε . The sample extraction process is realized by the important weight of each piece of symbol information.
Combining the results of (42) with the characteristics of maximum likelihood detection, we can obtain the sequence of probability distributions d i m , k conditioned on f ^ i m , k as
p f ^ i m , k | d i m , k = 1 π e f ^ i m , k α R i m , k d i m , k 2
Suppose that subcarriers in the same subblock are uniformly activated; the prior probability of symbols in any subblock can be expressed as
p d g i m , k = K M N c s u b d g i m , k d M m , k N c s u b K N c s u b otherwise
Combining the signal model in (42) with the probability distribution function of (43) and (44), the importance weight of each symbol sample for SMC [48] can be represented as
ρ g i m , k j = ρ g 1 i m , k j p f ^ g i m , k | κ g 1 i m , k j = ρ g 1 i m , k j p f ^ g i m , k | d g i m , k , κ g 1 i m , k j × p d g i m , k , j = 1 , 2 , , ε
where κ g 1 i m , k j is the sample extracted based on subcarriers, and p f ^ g i m , k | κ g 1 i m , k j is the prediction distribution of f g i m , k under the condition of the particles κ g 1 i m , k j .
According to the comparison results of the importance weight, the L s u b state combinations and symbol information corresponding to all subcarriers of each OFDM symbol can be obtained. Then, the posteriori probability corresponding to l th symbol vector can be expressed as
P m , k l = p f ^ i m , k | d ˜ i m , k l p d ˜ i m , k l
where d ˜ i m , k l is the l th symbol vector. We arrange them in descending order by the posteriori probability P m , k l . In addition, after N c _ s u b consecutive detections, the validity of the combination with the largest posterior probability on each subcarrier is also checked sequentially.
To reduce the complexity, we use a combination method to check the legality of all combinations without using a lookup table, which can be written as
γ = k = 1 K C A k 1 k
where A denotes the vector storing the indexes of all activated subcarriers. According to (47), this paper regards the state combination of γ < C N c _ s u b K as legal. Therefore, we can effectively check the legality of all combinations and discard all illegal status combinations. Finally, the legal state combination with the highest posterior probability is selected as the final MIMO-OFDM-IM subblock estimation.

4. Simulation and Numerical Results

In this section, some simulation results are presented to verify the performance of the proposed joint signal processing schemes based on sparse reconstruction in an integrated MIMO-OFDM-IM system. Some parameters required in the simulation are shown in Table 1. In the experiment, we consider an application scenario in which the integrated platform simultaneously sends integrated signals to two targets at a range of 50 m and 100 m, as well as to two communication users, and the velocity of the targets is 20 m/s and 10 m/s, respectively. The integrated platform performs matched filtering and parameter estimation after receiving the echo signal. Meanwhile, the communication users receive and demodulate the communication signals.

4.1. Radar Performance

In this subsection, we verify the radar performance of different radar signal processing methods. To ensure the validity of the comparison results, we set the spectral efficiency of both the MIMO-OFDM system and the MIMO-OFDM-IM system to 4 bps/Hz by adjusting the relevant parameters.

4.1.1. PSLR Performance

First, we pay attention to the peak sidelobe ratio (PSLR) of range and velocity two-dimensional images. The calculation method of the PSLR can be expressed as follows:
PSLR = 10 log 10 P s l P m l
where P s l denotes the peak level of the highest sidelobe and P m l denotes the peak level of the main lobe.
The PSLR performance of the proposed scheme is shown in Figure 5. In the simulation, N t = N r = 4 , the standard Monte Carlo technique with 10 4 independent trials is used for different SNR levels, and K = 4 in Figure 5a. In Figure 5a, the PSLR of the improved method under a different number of pulse groups and the PSLR of the traditional processing method [27] are plotted as a function of the SNR. It can be observed that the PSLR of the traditional FFT method with information compensation gradually decreases with the SNR, but tends to a constant value at a high SNR. This is because this method does not filter the valid data, which leads to the deterioration of the results with the uncertainty of index modulation. However, the PSLR of the proposed method decreases almost linearly with the SNR. Moreover, because a smaller number of pulse groups will result in more pulses in the same pulse group, the PSLR will decrease with the decrease in H. In addition, the PSLR performance of the multichannel is slightly improved compared with that of a single channel. Figure 5b shows the variation in the PSLR of different methods with K under SNR = 10 dB. It can be observed that the PSLR of the traditional FFT method based on information compensation decreases with the increase in K. When K < 4 , the PSLR drops rapidly with increasing K. On the other hand, when K > 4 , the PSLR drops very slowly. However, the PSLR of the proposed method is much lower and hardly changes with K under the same SNR, which shows that the uncertainty of index modulation has little influence on the proposed method and the overall performance is better.
Figure 6 shows the variation in the PSLR with the SNR for different methods when K is 3, 5, and 7. It can be found that the PSLR of FFT methods decreases with the increase in the SNR, but when the SNR reaches a certain value, it remains almost unchanged. When all subcarriers in a subblock are activated, the MIMO-OFDM-IM signal reverts to the traditional MIMO-OFDM signal, and the PSLR decreases linearly with an increasing SNR. The PSLR performance of the improved parameter estimation method (IPEM) based on the MIMO-OFDM system in [15] is better compared with the MUSIC algorithm. In contrast, the proposed method has a lower PSLR under the same conditions and decreases with the increase in the SNR. Also, Figure 6 validates that the PSLR performance of the proposed scheme hardly varies with K at different SNRs, which is consistent with the conclusion drawn in Figure 5. It is also noted that, while the IPEM [15] exhibits the same advantage, its performance is consistently about 6dB worse than the proposed method.

4.1.2. Range and Velocity Estimation Performance

In order to quantify the target range and velocity estimation performance of the proposed method, Monte Carlo experiments with 10 4 repetitions were carried out under different SNR settings. In the simulation, K = 4 . The other parameters are as given in Table 1. The four range and velocity estimation methods used for simulation comparisons are as follows: (1) the MUSIC algorithm under the MIMO-OFDM system, (2) the improved parameter estimation method (IPEM) under the MIMO-OFDM system in [15], (3) the FFT algorithm under the MIMO-OFDM-IM system, and (4) the proposed method under the MIMO-OFDM-IM system. Figure 7 shows the variation in the mean-square errors (MSEs) of range and velocity estimation with the SNR under different methods. It can be seen that the MSE of range and velocity estimations both decrease with the increasing SNR. The estimation accuracies of the MUSIC algorithm under the MIMO-OFDM system are basically similar to the published IPEM method in [15], and both are higher than that of the FFT algorithm. In contrast, the proposed method can obtain better estimation accuracy. This is because the proposed method screens the effective data under the designed integrated waveform condition based on MIMO-OFDM-IM. In addition, although the velocity estimation accuracy of the proposed method is better than the methods in the literature and the MUSIC algorithm under the MIMO-OFDM system, its improvement is not as significant as the range estimation improvement, mainly because the proposed method performs FFT on the velocity dimension after dimensionality reduction.

4.1.3. The Complexity of Parameter Estimation

Figure 8 shows the complexity of the proposed parameter estimation method. In the simulation, 10 4 Monte Carlo experiments were carried out. According to Figure 8a, when H = 1, the algorithm does not reduce the dimension, and the number of complex multiplications is large. Having said that, the complexity of the proposed method can be effectively reduced by increasing H, as seen in Figure 8a. In addition, Figure 8b shows that when H = 4, the complexity of the proposed system will be slightly higher than that of the MIMO-OFDM system. However, when H is increased to eight, the complexity of the proposed system is greatly reduced, which is lower than the complexity of the improved algorithm in the traditional MIMO-OFDM system. Ultimately, the H to be chosen should be a reasonable compromise between the algorithm complexity and parameter estimation performance.

4.2. Communication Performance

4.2.1. BER Performance

In this subsection, the bit error rate (BER) performance of the proposed signal detection method for the integrated MIMO-OFDM-IM system is evaluated via Monte Carlo simulations with 10 4 trials. In this simulation, the communication symbol is modulated by BPSK modulation. Similar to the previous section, the spectral efficiency is uniformly set to 4 bps/Hz. Five detection methods are used for comparative analysis, including the maximum likelihood (ML) detection method, minimum mean square error (MMSE) detection method, improved MMSE detection method [50], subblock-wise detection method [51], and the proposed method. In Figure 9, the BER of different methods for different SNR levels is depicted with N r = N t = 4. According to the result of Figure 9a, it can be observed that the performance loss of the MMSE detection method is substantial, whilst the performance of the improved MMSE detection method has improved, albeit not significantly. In contrast, the performance of the subblock-wise detection method has greatly improved. Finally, the proposed detection method achieves a better BER performance and has only a small performance loss compared with the ML detection method. A further comparison shows that, with an increasing K (Figure 9b), the BER performance of all detection methods decreases slightly, mainly because the modulation symbols to be detected will increase with the increase in K. In this case, the performance of the proposed method is still slightly better than the subblock-wise detection method, and both methods are still very close to the ML detection method.
Figure 10 shows the BER performances of different detection methods for different SNR levels with N r = N t = 8. As seen from Figure 10, the BER performance of the ML detection method, subblock-wise detection method, and the proposed detection method are obviously superior to other detection methods. Comparing Figure 10 with Figure 9, it can be observed that the BER performances of all information detection methods have obviously improved with the increase in the number of transmitting and receiving antennas. In particular, the performance of the subblock-wise detection method is closer to that of the proposed information detection method. However, the algorithm complexity of the subblock-wise detection method is significantly higher than that of the proposed method, which will be shown in the statistical results of the algorithm complexity later. It is worth mentioning that the ML detection method has a better BER performance than other methods, but its algorithm complexity is very high. Comparing Figure 10a and Figure 10b, it can be observed that all information detection methods are sensitive to K.

4.2.2. The Complexity of Communication Detection Methods

In the following analysis, the number of complex multiplications is used as an evaluation indicator of algorithm complexity. The number of complex multiplications of different detection methods under the same parameter conditions is shown in Table 2. Obviously, the complexity of the ML detection method is very high, and with the increasing N t , its complexity increases exponentially. An MMSE-based detection method can effectively reduce the computational complexity, but its BER performance is not great, as has been observed earlier. In contrast, the improved detection method based on SMC can effectively improve the information detection performance while ensuring the algorithm complexity. Especially, with the increase in N t , the complexity of the proposed method is lower than that of the subblock-wise detection method. Generally speaking, the signal detection method under the integrated MIMO-OFDM-IM system can ensure stable synchronization and a relatively good communication information detection performance.

4.3. Trade-Off Performance

In the integrated MIMO-OFDM-IM system, to evaluate the influence of different activated subcarrier combinations on system performance, we constructed a trade-off curve between radar and communication through simulations employing different activated subcarrier combinations, and the trade-off parameters include radar target detection probability and communication BER. In the simulations, we carried out 10 4 Monte Carlo experiments, and in each experiment, the target was detected in the radar monitoring range. Figure 11 shows the variation in the radar–communication trade-off curve with the number of activated subcarriers under different SNRs. It can be found that the radar target detection probability increases monotonically with K and becomes stable when the number equals or is more than five. On the other hand, when K = 5 , the communication BER performance is the best. Therefore, the appropriate activated subcarrier combinations can be selected according to different requirements of the actual environment, while a general rule-of-thumb approach would be using five activated subcarriers since it gives the best communication performance and relatively high target detection probability.

5. Conclusions

In this paper, a sparse reconstruction-based joint signal processing scheme for integrated MIMO-OFDM-IM systems is proposed, which can achieve a higher target range and velocity estimation accuracy and lower communication BER under the same spectral efficiency. According to the combination of MIMO-OFDM and IM technology, the waveform design based on the MIMO-OFDM-IM integrated system is realized. For radar, after conventional processing and information compensation, a segmented coherent accumulation method based on sparse signal reconstruction is used to realize parameter estimation. The numerical results show that the proposed method can achieve a better PSLR performance and higher accuracy of range and velocity estimation. For communication, we use the subcarriers corresponding to the OFDM symbols of each subblock as a benchmark; the information samples are detected iteratively through an antenna index; and the information symbols are obtained through a series of processes such as weight updating, screening, and sorting. The simulation experiment proved that the proposed scheme has a high information transmission efficiency and good communication information detection performance. In the future, sn MIMO-OFDM waveform design based on adaptive, robust, and PAPR reduction will be our research focus. The optimized waveform design can further solve the problems of radar–communication performance trade-off and algorithm practicability.

Author Contributions

Y.W. and Y.C. (Yunhe Cao) developed the theory and signal model. Y.W. and Y.C. (Yunhe Cao) performed and analyzed the numerical simulations. Y.W., T.-S.Y., Y.C. (Yuanhao Cheng) and Y.Z. wrote and edited the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the National Nature Science Foundation of China under grant 61771367.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author. The data are not publicly available due to restrictions, e.g., privacy or ethics.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
OFDMorthogonal frequency division multiplexing
IMindex modulation
IRCSintegrated radar and communication systems
UAVunmanned aerial vehicle
MIMOmultiple-input multiple-output
HADhybrid analog–digital
DLdeep learning
ICIintercarrier interference
IFFTinverse fast Fourier transform
CPcyclic prefix
GIguard interval
ISIintersymbol interference
DFTdiscrete Fourier transform
CRBCramér–Rao bound
SMCsequential Monte Carlo
SNRsignal-to-noise ratio
PSLRpeak sidelobe ratio
MSEmean-square error
BERbit error rate
MLmaximum likelihood
MMSEminimum mean square error

References

  1. Kumari, P.; Choi, J.; González-Prelcic, N.; Heath, R.W. IEEE 802.11ad-Based Radar: An Approach to Joint Vehicular Communication-Radar System. IEEE Trans. Veh. Technol. 2018, 67, 3012–3027. [Google Scholar] [CrossRef]
  2. Grossi, E.; Lops, M.; Venturino, L. Joint Design of Surveillance Radar and MIMO Communication in Cluttered Environments. IEEE Trans. Signal Process. 2020, 68, 1544–1557. [Google Scholar] [CrossRef]
  3. Gu, M.-X.; Lee, M.-C.; Liu, Y.-S.; Lee, T.-S. Design and Analysis of Frequency Hopping-Aided FMCW-Based Integrated Radar and Communication Systems. IEEE Trans. Commun. 2022, 70, 8416–8432. [Google Scholar] [CrossRef]
  4. Zhang, J.; Liu, F.; Masouros, C.; Heath, R.; Feng, Z.; Zheng, L.; Petropulu, A. An overview of signal processing techniques for joint communication and radar sensing. IEEE J. Sel. Topics Signal Process. 2021, 15, 1295–1315. [Google Scholar] [CrossRef]
  5. Wang, J.; Liang, X.-D.; Chen, L.-Y.; Li, Y.-L. Waveform Designs for Joint Wireless Communication and Radar Sensing: Pitfalls and Opportunities. IEEE Internet Things J. 2023, 10, 15252–15265. [Google Scholar] [CrossRef]
  6. Xiao, Z.; Zeng, Y. Waveform design and performance analysis for full-duplex integrated sensing and communication. IEEE J. Sel. Areas Commun. 2022, 40, 1823–1837. [Google Scholar] [CrossRef]
  7. Liu, F.; Masouros, C.; Petropulu, A.P.; Griffiths, H.; Hanzo, L. Joint Radar and Communication Design: Applications, State-of-the-Art, and the Road Ahead. IEEE Trans. Commun. 2020, 68, 3834–3862. [Google Scholar] [CrossRef]
  8. Wang, Y.; Cao, Y.; Yeo, T.-S.; Cheng, Y.; Fu, J. Robust Multiplexing Waveform Design for Integrated OFDM Radar and Communication via Complex Weight Optimization. Remote Sens. 2023, 15, 4438. [Google Scholar] [CrossRef]
  9. Liu, F.; Zhou, L.; Masouros, C.; Li, A.; Luo, W.; Petropulu, A. Toward Dual-functional Radar-Communication Systems: Optimal Waveform Design. IEEE Trans. Signal Process. 2018, 66, 4264–4279. [Google Scholar] [CrossRef]
  10. Yu, X.; Yao, X.; Yang, J.; Zhang, L.; Kong, L.; Cui, G. Integrated Waveform Design for MIMO Radar and Communication via Spatio-Spectral Modulation. IEEE Trans. Signal Process. 2022, 70, 2293–2305. [Google Scholar] [CrossRef]
  11. Zheng, L.; Lops, M.; Eldar, Y.C.; Wang, X. Radar and Communication Coexistence: An Overview: A Review of Recent Methods. IEEE Signal Process. Mag. 2019, 36, 85–99. [Google Scholar] [CrossRef]
  12. Hwang, S.; Seo, J.; Park, J.; Kim, H.; Jeong, B.J. Compressive Sensing-Based Radar Imaging and Subcarrier Allocation for Joint MIMO OFDM Radar and Communication System. Sensors 2021, 21, 2382. [Google Scholar] [CrossRef] [PubMed]
  13. Rong, J.; Liu, F.; Miao, Y. High-Efficiency Optimization Algorithm of PMEPR for OFDM Integrated Radar and Communication Waveform Based on Conjugate Gradient. Remote Sens. 2022, 14, 1715. [Google Scholar] [CrossRef]
  14. Yamamoto, H.; Nakamura, A.; Itami, M. A study on LDM-BST-OFDM transmission for the next-generation terrestrial broadcasting. IEEE Trans. Broadcast. 2020, 66, 205–215. [Google Scholar] [CrossRef]
  15. Liu, Y.; Liao, G.; Chen, Y.; Xu, J.; Yin, Y. Super-resolution range and velocity estimations with OFDM integrated radar and communications waveform. IEEE Trans. Veh. Technol. 2020, 69, 11659–11672. [Google Scholar] [CrossRef]
  16. Liu, Y.; Liao, G.; Xu, J.; Yang, Z.; Huang, L.; Zhang, Y. Transmit power adaptation for orthogonal frequency division multiplexing integrated radar and communication systems. J. Appl. Remote Sens. 2017, 11, 035017. [Google Scholar] [CrossRef]
  17. Shi, C.; Wang, F.; Salous, S.; Zhou, J. Joint subcarrier assignment and power allocation strategy for integrated radar and communications system based on power minimization. IEEE Sens. J. 2019, 19, 11167–11179. [Google Scholar] [CrossRef]
  18. Temiz, M.; Alsusa, E.; Baidas, M.W. A dual-functional massive MIMO OFDM communication and radar transmitter architecture. IEEE Trans. Veh. Technol. 2020, 69, 14974–14988. [Google Scholar] [CrossRef]
  19. Keskin, M.F.; Koivunen, V.; Wymeersch, H. Limited feedforward waveform design for OFDM dual-functional radar-communications. IEEE Trans. Signal Process. 2021, 69, 2955–2970. [Google Scholar] [CrossRef]
  20. Ahmed, A.; Zhang, Y.D.; Hassanien, A. Joint Radar-Communications Exploiting Optimized OFDM Waveforms. Remote Sens. 2021, 13, 4376. [Google Scholar] [CrossRef]
  21. Saad, M.; Palicot, J.; Bader, F.; Ghouwayel, A.C.A.; Hijazi, H. A Novel Index Modulation Dimension Based on Filter Domain: Filter Shapes Index Modulation. IEEE Trans. Commun. 2021, 69, 1445–1461. [Google Scholar] [CrossRef]
  22. Yu, P.; Huang, F. Spatial Index Modulation with the Combination of Both Antenna Indexes and Two Constellations. IEEE Access 2022, 10, 133778–133788. [Google Scholar] [CrossRef]
  23. Choi, J. Coded OFDM-IM With Transmit Diversity. IEEE Trans. Commun. 2017, 65, 3164–3171. [Google Scholar] [CrossRef]
  24. Zhu, Y.; Wang, B.; Li, J.; Zhang, Y.; Xie, F. Y-Shaped Net-Based Signal Detection for OFDM-IM Systems. IEEE Commun. Lett. 2022, 26, 2661–2664. [Google Scholar] [CrossRef]
  25. Li, J.; Dang, S.; Huang, Y.; Chen, P.; Qi, X.; Wen, M.; Arslan, H. Composite Multiple-Mode Orthogonal Frequency Division Multiplexing with Index Modulation. IEEE Trans. Wireless Commun. 2023, 22, 3748–3761. [Google Scholar] [CrossRef]
  26. Sui, Z.; Yan, S.; Zhang, H.; Yang, L.-L.; Hanzo, L. Approximate Message Passing Algorithms for Low Complexity OFDM-IM Detection. IEEE Trans. Veh. Technol. 2021, 70, 9607–9612. [Google Scholar] [CrossRef]
  27. Xu, Z.; Petropulu, A. A Bandwidth Efficient Dual-Function Radar Communication System Based on a MIMO Radar Using OFDM Waveforms. IEEE Trans. Signal Process. 2023, 71, 401–416. [Google Scholar] [CrossRef]
  28. Xu, J.; Zhou, Z.; Li, L.; Zheng, L.; Liu, L. RC-Struct: A Structure-Based Neural Network Approach for MIMO-OFDM Detection. IEEE Trans. Wirel. Commun. 2022, 21, 7181–7193. [Google Scholar] [CrossRef]
  29. Qian, J.; Lops, M.; Zheng, L.; Wang, X.; He, Z. Joint System Design for Coexistence of MIMO Radar and MIMO Communication. IEEE Trans. Signal Process. 2018, 66, 3504–3519. [Google Scholar] [CrossRef]
  30. Basar, E. Multiple-input multiple-output OFDM with index modulation. IEEE Signal Process. Lett. 2015, 22, 2259–2263. [Google Scholar] [CrossRef]
  31. Yang, Y.; Dang, S.; Wen, M.; Guizani, M. Millimeter Wave MIMO-OFDM with Index Modulation: A Pareto Paradigm on Spectral- Energy Efficiency Trade-Off. IEEE Trans. Wirel. Commun. 2021, 20, 6371–6386. [Google Scholar] [CrossRef]
  32. Zhao, T.; Li, F. Variational-autoencoder signal detection for MIMO-OFDM-IM. Digit. Signal Process. 2021, 118, 103230. [Google Scholar] [CrossRef]
  33. Basar, E. On multiple-input multiple-output OFDM with index modulation for next generation wireless networks. IEEE Trans. Signal Process. 2016, 64, 3868–3878. [Google Scholar] [CrossRef]
  34. Wang, L. Generalized Quadrature Space-Frequency Index Modulation for MIMO-OFDM Systems. IEEE Trans. Commun. 2021, 69, 6375–6389. [Google Scholar] [CrossRef]
  35. Hodge, J.A.; Mishra, K.V.; Zaghloul, A.I. Intelligent Time-Varying Metasurface Transceiver for Index Modulation in 6G Wireless Networks. IEEE Antennas Wirel. Propag. Lett. 2020, 19, 1891–1895. [Google Scholar] [CrossRef]
  36. Chen, R.; Zheng, J. Index-Modulated MIMO-OFDM: Joint Space-Frequency Signal Design and Linear Precoding in Rapidly Time-Varying Channels. IEEE Trans. Wirel. Commun. 2018, 17, 7067–7079. [Google Scholar] [CrossRef]
  37. Noh, H.; Lee, H.; Yang, H.J. IICI-Robust Transceiver Design for Integration of MIMO-OFDM Radar and MU-MIMO Communication. IEEE Trans. Veh. Technol. 2023, 72, 821–838. [Google Scholar] [CrossRef]
  38. Schindler, D.; Schweizer, B.; Knill, C.; Hasch, J.; Waldschmidt, C. Synthetization of Virtual Transmit Antennas for MIMO OFDM Radar by Space-Time Coding. IEEE Trans. Aerosp. Electron. Syst. 2021, 57, 1964–1971. [Google Scholar] [CrossRef]
  39. Hakobyan, G.; Ulrich, M.; Yang, B. OFDM-MIMO Radar With Optimized Nonequidistant Subcarrier Interleaving. IEEE Trans. Aerosp. Electron. Syst. 2020, 56, 572–584. [Google Scholar] [CrossRef]
  40. Multerer, T.; Prechtel, U.; Vossiek, M.; Ziegler, V. Systematic Phase Correction for Direction-of-Arrival Estimation in Spectrally Interleaved OFDM MIMO Radar. IEEE Trans. Microw. Theory Technol. 2019, 67, 4570–4577. [Google Scholar] [CrossRef]
  41. Zhang, S.-Y.; Shahrrava, B. A SLM Scheme for PAPR Reduction in Polar Coded OFDM-IM Systems Without Using Side Information. IEEE Trans. Broadcast. 2021, 67, 463–472. [Google Scholar] [CrossRef]
  42. Deng, S. Reduced-Guard-Interval OFDM Using Digital Sub-Band-Demultiplexing. IEEE Photon. Technol. Lett. 2013, 25, 2174–2177. [Google Scholar] [CrossRef]
  43. Liu, Y.; Yi, J.; Wan, X.; Zhang, X.; Ke, H. Evaluation of Clutter Suppression in CP-OFDM-Based Passive Radar. IEEE Sens. J. 2019, 19, 5572–5586. [Google Scholar] [CrossRef]
  44. Wang, Y.; Cao, Y.; Yeo, T.-S.; Han, J.; Peng, Z. Joint Sequence Optimization-Based OFDM Waveform Design for Integrated Radar and Communication Systems. IEEE Trans. Veh. Technol. 2022, 71, 12734–12748. [Google Scholar] [CrossRef]
  45. Roth, I.; Kliesch, M.; Flinth, A.; Wunder, G.; Eisert, J. Reliable Recovery of Hierarchically Sparse Signals for Gaussian and Kronecker Product Measurements. IEEE Trans. Signal Process. 2020, 68, 4002–4016. [Google Scholar] [CrossRef]
  46. Chen, Y.; Liu, Y.; Chen, J.; Ma, J. A Novel Identification Method for a Class of Closed-Loop Systems Based on Basis Pursuit De-Noising. IEEE Access 2020, 8, 99648–99654. [Google Scholar] [CrossRef]
  47. Vesa, J.; Rasilo, P. Permeability Estimations of SMC Material Particles. IEEE Trans. Magn. 2020, 56, 2000807. [Google Scholar] [CrossRef]
  48. Wang, Y.; Xia, Y.; Shen, H.; Zhou, P. SMC Design for Robust Stabilization of Nonlinear Markovian Jump Singular Systems. IEEE Trans. Autom. Control 2018, 63, 219–224. [Google Scholar] [CrossRef]
  49. Yih, C. Iterative interference cancellation for OFDM signals with blanking nonlinearity in impulsive noise channels. IEEE Signal Process. Lett. 2012, 19, 147–150. [Google Scholar] [CrossRef]
  50. Fang, L.; Xu, L.; Huang, D.D. Low Complexity Iterative MMSE-PIC Detection for Medium-Size Massive MIMO. IEEE Wirel. Commun. Lett. 2016, 5, 108–111. [Google Scholar] [CrossRef]
  51. Zheng, B.; Wen, M.; Basar, E.; Chen, F. Multiple-Input Multiple-Output OFDM With Index Modulation: Low-Complexity Detector Design. IEEE Trans. Signal Process. 2017, 65, 2758–2772. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the proposed MIMO-OFDM integrated system model.
Figure 1. Schematic diagram of the proposed MIMO-OFDM integrated system model.
Remotesensing 16 01773 g001
Figure 2. Schematic diagram of MIMO-OFDM-IM integrated signal design in transmitter.
Figure 2. Schematic diagram of MIMO-OFDM-IM integrated signal design in transmitter.
Remotesensing 16 01773 g002
Figure 3. Flow chart of the joint signal processing schemes in the proposed system.
Figure 3. Flow chart of the joint signal processing schemes in the proposed system.
Remotesensing 16 01773 g003
Figure 4. Flowchart of the conventional processing.
Figure 4. Flowchart of the conventional processing.
Remotesensing 16 01773 g004
Figure 5. The PSLR performance. (a) The variation in PSLR with SNR. (b) The variation in PSLR with the number of activated subcarriers (SNR = 10 dB).
Figure 5. The PSLR performance. (a) The variation in PSLR with SNR. (b) The variation in PSLR with the number of activated subcarriers (SNR = 10 dB).
Remotesensing 16 01773 g005
Figure 6. The variation in PSLR with SNR under different number of activated subcarriers.
Figure 6. The variation in PSLR with SNR under different number of activated subcarriers.
Remotesensing 16 01773 g006
Figure 7. The variation in MSEs of range and velocity estimation with SNR under different methods ( N r = N t = 4, N c _ s u b = 8, K = 4, H = 4). (a) Range estimation. (b) Velocity estimation.
Figure 7. The variation in MSEs of range and velocity estimation with SNR under different methods ( N r = N t = 4, N c _ s u b = 8, K = 4, H = 4). (a) Range estimation. (b) Velocity estimation.
Remotesensing 16 01773 g007
Figure 8. Complexity comparison. (a) The variation in number of complex multiplications with the number of activated subcarriers. (b) The variation in number of complex multiplications with the number of transmitting antennas.
Figure 8. Complexity comparison. (a) The variation in number of complex multiplications with the number of activated subcarriers. (b) The variation in number of complex multiplications with the number of transmitting antennas.
Remotesensing 16 01773 g008
Figure 9. The BER performance of different information detection methods under the proposed integrated MIMO-OFDM-IM system ( N r = N t = 4). (a) N c _ s u b = 8, K = 3; (b) N c _ s u b = 8, K = 6.
Figure 9. The BER performance of different information detection methods under the proposed integrated MIMO-OFDM-IM system ( N r = N t = 4). (a) N c _ s u b = 8, K = 3; (b) N c _ s u b = 8, K = 6.
Remotesensing 16 01773 g009
Figure 10. The BER performance of different information detection methods under the proposed integrated MIMO-OFDM-IM system ( N r = N t = 8). (a) N c _ s u b = 8, K = 3; (b) N c _ s u b = 8, K = 6.
Figure 10. The BER performance of different information detection methods under the proposed integrated MIMO-OFDM-IM system ( N r = N t = 8). (a) N c _ s u b = 8, K = 3; (b) N c _ s u b = 8, K = 6.
Remotesensing 16 01773 g010
Figure 11. The variation in radar–communication trade-off curve with the number of activated subcarriers under different SNRs ( N r = N t = 4). (a) SNR = 6 dB; (b) SNR = 12 dB.
Figure 11. The variation in radar–communication trade-off curve with the number of activated subcarriers under different SNRs ( N r = N t = 4). (a) SNR = 6 dB; (b) SNR = 12 dB.
Remotesensing 16 01773 g011
Table 1. The parameters of simulation.
Table 1. The parameters of simulation.
ParametersValue
Bandwidth20 MHz
Carrier frequency5 GHz
Duration of pulses32 μ s
Number of OFDM symbols8
Number of subcarriers256
Number of subblocks32
Total number of pulses32
Lower limit of the number of state combinations6
Upper limit of the number of samples10
Modulation mode of communication informationBPSK
Table 2. Algorithm complexity comparison of different detection methods.
Table 2. Algorithm complexity comparison of different detection methods.
ParametersNumber of Complex Multiplications
N t = 4 N t = 8 N t = 12
ML detection3.26 × 10 6
MMSE detection6.62 × 10 2 5.28 × 10 3 9.35 × 10 3
Improved MMSE detection2.07 × 10 3 9.23 × 10 3 2.84 × 10 4
Subcarrier-wise detection2.23 × 10 3 8.51 × 10 3 1.5 × 10 4
Proposed improved detection2.27 × 10 3 6.04 × 10 3 9.61 × 10 3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, Y.; Cao, Y.; Yeo, T.-S.; Cheng, Y.; Zhang, Y. Sparse Reconstruction-Based Joint Signal Processing for MIMO-OFDM-IM Integrated Radar and Communication Systems. Remote Sens. 2024, 16, 1773. https://doi.org/10.3390/rs16101773

AMA Style

Wang Y, Cao Y, Yeo T-S, Cheng Y, Zhang Y. Sparse Reconstruction-Based Joint Signal Processing for MIMO-OFDM-IM Integrated Radar and Communication Systems. Remote Sensing. 2024; 16(10):1773. https://doi.org/10.3390/rs16101773

Chicago/Turabian Style

Wang, Yang, Yunhe Cao, Tat-Soon Yeo, Yuanhao Cheng, and Yulin Zhang. 2024. "Sparse Reconstruction-Based Joint Signal Processing for MIMO-OFDM-IM Integrated Radar and Communication Systems" Remote Sensing 16, no. 10: 1773. https://doi.org/10.3390/rs16101773

APA Style

Wang, Y., Cao, Y., Yeo, T. -S., Cheng, Y., & Zhang, Y. (2024). Sparse Reconstruction-Based Joint Signal Processing for MIMO-OFDM-IM Integrated Radar and Communication Systems. Remote Sensing, 16(10), 1773. https://doi.org/10.3390/rs16101773

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop