Next Article in Journal
Displacement Study of a Large-Scale Freeform Timber Plate Structure Using a Total Station and a Terrestrial Laser Scanner
Next Article in Special Issue
Compressive Sensing Spectroscopy Using a Residual Convolutional Neural Network
Previous Article in Journal
An Ultrasonic Object Detection Applying the ID Based on Spread Spectrum Technique for a Vehicle
Previous Article in Special Issue
A Neuron-Based Kalman Filter with Nonlinear Autoregressive Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cognitive Frequency-Hopping Waveform Design for Dual-Function MIMO Radar-Communications System

1
School of Information Engineering, East China Jiaotong University, Nanchang 330031, China
2
School of Information Science and Engineering, Southeast University, Nanjing 210096, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(2), 415; https://doi.org/10.3390/s20020415
Submission received: 30 November 2019 / Revised: 3 January 2020 / Accepted: 8 January 2020 / Published: 11 January 2020
(This article belongs to the Special Issue Sensor Signal and Information Processing III)

Abstract

:
A frequency-hopping (FH)-based dual-function multiple-input multiple-output (MIMO) radar communications system enables implementation of a primary radar operation and a secondary communication function simultaneously. The set of transmit waveforms employed to perform the MIMO radar task is generated using FH codes. For each transmit antenna, the communication operation can be realized by embedding one phase symbol during each FH interval. However, as the radar channel is time-variant, it is necessary for a successive waveform optimization scheme to continually obtain target feature information. This research work aims at enhancing the target detection and feature estimation performance by maximizing the mutual information (MI) between the target response and the target returns, and then minimizing the MI between successive target-scattering signals. The two-step cognitive waveform design strategy is based upon continuous learning from the radar scene. The dynamic information about the target feature is utilized to design FH codes. Simulation results show an improvement in target response extraction, target detection probability and delay-Doppler resolution as the number of iterations increases, while still maintaining high data rate with low bit error rates between the proposed system nodes.

1. Introduction

The multiple-input multiple-output (MIMO) radar [1,2,3] is a significant technique owing to the enhancement in performance it provides over the traditional radar systems with a single transmit antenna. The transmit waveforms are detected by a matched filter bank in the MIMO radar receiver. Making use of the knowledge of the propagation channels, a superior spatial resolution can be achieved. Furthermore, a MIMO radar has several advantages such as prominent interference rejection capability, enhanced parameter estimation performance, and better flexibility for transmit waveform optimization [4,5,6]. Several works in the literature discussed the subject of MIMO waveform optimization [7,8,9,10,11]. In [7], the authors considered the joint design of both transmit waveforms and receive filters for a collocated MIMO radar with the existence of signal-dependent interference and white noise. The design problem is formulated into a maximization of the signal-to-interference-plus-noise ratio (SINR), including various constraints on the transmit waveforms. Stoica and Li designed the covariance matrix of the transmit waveforms to control the spatial power [8]. However, the cross correlation between the transmit waveforms at a specific target range is minimized. Several research works [12,13,14,15] designed the transmit waveform directly instead of just the covariance matrix. The work of [13] considered prior knowledge of the target impulse response (TIR) and utilized the information to select the transmit waveforms which optimize the mutual information (MI) between the target echoes and the TIR. The work of [15] investigated the problem of the spectrally compatible waveform design for MIMO radar in the presence of multiple targets and signal-dependent interference. A new method was proposed to deal with a more general problem, i.e., designing a spectrally compatible waveform for multiple targets, by minimizing the waveform energy of the overlayed space-frequency bands. The waveform optimization which employs prior knowledge of TIR is also implemented in the single-input multiple-output (SIMO) radar system [16].
The work of [17] derived MIMO radar ambiguity functions. San Antonio and Fuhrmann discussed some properties of the MIMO radar ambiguity function, which provide several ideas for MIMO waveform optimization [17]. The frequency-hopping (FH) sequences presented in [18] were used in MIMO radar configuration [19,20]. The FH orthogonal transmit waveforms discussed in [18] are initially considered for multiuser radar system. Furthermore, [18] employed the FH codes to minimize the peaks in the cross correlation functions of the transmit waveforms as much as possible. However, in the multiuser radar scenario, each operator activates its individual system. This is different from the MIMO radar system where the receiver antennas can cooperate to resolve the target responses. The FH sequences, which were designed by Chen and Yang, optimized the MIMO ambiguity function.
Based on their characteristics of being easily produced and modulus constant, FH codes are considered a good choice for the MIMO radar waveforms. A new scheme for optimizing the MIMO radar waveforms was provided in [19]. The scheme makes the energy of the MIMO ambiguity function spread in the range and angular dimensions evenly, as well as decreases the sidelobes in the MIMO radar ambiguity function. The work of [19] also designed optimal FH waveforms, which had separate FH codes and amplitudes, for a collocated MIMO radar system. The joint design problem can be solved by using game theory, provided by [20]. The authors considered the two objective functions, corresponding to FH codes and amplitudes separately, as two interacting players. By this concept, the joint optimization scheme obtained better integrated code and amplitude matrices that can improve performance much better than separate designs.
The works of [21,22,23,24,25,26] indicated a possibility of employing the radar-communications integration concept to solve the lack of radio frequency (RF) spectrum. Efficient utilization of shared bandwidth between wireless communications and radar can be achieved by using dynamic frequency allocation. For example, using dynamic frequency allocation is a way to make shared bandwidth between wireless communications and radar possess more efficient utilization. The work of [23] proposed a novel dual-function radar-communications (DFRC) strategy to embed quadrature amplitude modulation (QAM) based communication information in the radar waveforms by exploiting sidelobe control and waveform diversity. In [27,28], the authors proposed a novel distributed DFRC MIMO system capable of simultaneously performing radar and communication tasks. The distributed DFRC MIMO system performs both objectives by optimizing the power allocation of the different transmitters in the DFRC system. The proposed strategy can serve multiple communication receivers located in the vicinity of the distributed DFRC MIMO system. Numerous recent studies [29,30,31] considered that the developing concept of DFRC is secondary to the main radar task. Communication source embedding into the illumination of MIMO radar system is realized using waveform diversity, sidelobe control, or the time-modulated array technique, which was studied in [30]. Hassanien and Himed presented a signaling strategy for communication source embedding into the illumination of FH-based MIMO radar system [32]. The main principle behind the signaling strategy is to embed phase modulation (PM)-based symbols by using phase rotating the FH pulses. The phase shift is implemented to each transmit FH waveform of the MIMO radar system. The PM-based symbol embedding does not influence the function of the MIMO radar system, which uses the FH waveforms.
The investigation on cognitive radar waveform optimization has received a lot of interest [33,34,35]. To further enhance the performance of the TIR estimation in a time-varying target scene, the transmitted signal parameters should be constantly adjusted. Then, updated knowledge about the time-varying target scene is employed to allocate fundamental resources like transmitted signal parameters in a cognitive mode [34]. A new strategy for optimizing the waveforms of a cognitive radar was presented in [34]. The aim is to enhance the performance of target estimation by minimizing the mean-square error (MSE) of the estimates of target scattering coefficients (TSC) based on Kalman filtering and then minimizing MI between the radar target echoes at successive time instants. However, there is also an increase in the computational load due to the Kalman filtering step in the waveform optimization. Such a cognitive radar system cannot be used in real applications to address the environmental sensing issues. To improve the performance of the target parameters estimation and classification, the pioneering study by Bell in [36] developed an information-theoretic method for the radar waveforms optimization. The authors in [37,38,39] extended the information-theoretic method by maximizing the MI between the target response and the target-scattering signals as a waveform design criterion in the MIMO radar system. The work of [40] proposed an innovative method to designing the transmit signal of cognitive MIMO radar system, which combines the signal optimization and selection processes. The works of [41,42,43,44] present a signaling scheme for information embedding into the illumination of the radar using FH pulses. An FH-based joint radar-communication system enables implementing a primary radar operation and a secondary communication function simultaneously. Then, the authors consider the problem of radar codes optimization under a peak-to-average-power ratio (PAR) and an energy constraint. However, a time-variant radar scenario is not considered.
Based on the points discussed above, it is interesting to discuss the performance of an adaptive dual-function MIMO radar-communications system that combines the adaptive FH waveform optimization scheme stated in [38] and PM-based information embedding strategy stated in [31]. To further adapt to the dynamic radar environment, we consider the problem of adaptive waveforms design and propose a two-step waveform optimization scheme, which provides better target detection performance and high data-rate communication capability between the proposed dual-function systems. The proposed scheme is summarized as follows:
Step 1. Waveform Design: this part includes the design of the FH waveform for the individual MIMO transmit antennas. The primary goal is to maximize the MI between the target-scattering signal and the estimated target response. The method ensures that the target-scattering signal at each time instant is dependent on the target response.
Step 2. Waveform Selection: after the best waveform ensemble is gained, part two is to choice the suitable PM-based FH sequences for emission. The principle of this module is to minimize the MI between consecutive target scattering signals. This part ensures that we constantly obtain target returns that become independent of each other in time, with the purpose of achieving more information about the target characteristics at each time instant of reception.
The premise of the FH waveform optimization scheme is channel estimation. The target feature estimation can be performed by the MIMO radar receiver through observations implemented in the previous time instant. A feedback loop enables the delivery of the estimates to the dual-function MMO radar-communications transmitter. At a result, the optimization strategy allows the MIMO radar transmitter to constantly adjust FH codes to suit the time-varying channel scene.
We choose the FH waveform for the following reasons:
(1)
The channel environment is complicated to wireless communications due to densely populated scatterers. However, FH waveforms are immune to multipath channel fading under the circumstances.
(2)
The FH waveforms are robust to antagonistic environments by offering low interception probability. Furthermore, FH waveforms are immune to clutter interference.
(3)
The constant-modulus waveforms have the property of high transmission power efficiency. FH waveforms enjoy the constant-modulus feature and are easy to generate.
The primary innovations of our work are summarized as follows:
(1)
We define the PM-based FH waveforms in dual-function MIMO radar-communications configuration and derive the associated MIMO ambiguity function;
(2)
We develop a two-step waveform optimization scheme in the adaptive PM-based dual-function MIMO radar-communications framework;
(3)
We evaluate the performance of the proposed scheme in terms of target response estimation, delay-Doppler resolution and communication symbol error rate (SER).
(4)
We compare the proposed scheme with other radar systems through analysis of the target detection and receiver operating characteristics (ROC) in an interference noise environment.
The organization of this paper is as follows. In Section 2, we describe the dual-function MIMO radar-communications system and the PM-based FH signal model. In Section 3, we derive the MIMO radar ambiguity function of the PM-based FH waveform. In Section 4, we present an adaptive approach to optimizing the proposed information embedding waveform. The transmit waveforms are designed at step 1 of the algorithm and selected based on the criterion presented in step 2. The simulation results demonstrating the proposed scheme are presented in Section 5. Finally, our conclusions and directions for possible future work are drawn in Section 6.
Throughout this paper, the following notations will be used. We use boldface lowercase letters and boldface uppercase letters to denote vectors and matrices, respectively; ( . ) * to denote the complex conjugate operation; ( . ) T to denote the transpose operation; ( . ) H to denotes the Hermitian operation; to denote the Kronecker product; I M N to denote the identity matrix of size M N × M N ; A n g l e ( ) to denote the angle of a complex number.

2. System Configuration and Signal Model

2.1. Phase Modulation (PM)-Based Frequency Hopping (FH) Waveforms

A PM-based method for embedding information into the radar emission was recently proposed in [31]. To deliver a finite number of binary bits per radar pulse, the PM-based method maps the binary data into a phase symbol that belongs to a phase dictionary of an appropriate size. During each radar pulse, the PM-based method embeds one phase symbol into the radar emission toward the intended communication direction. At the communication receiver, a phase detector is used to detect the embedded symbol and, subsequently, deciphers the corresponding binary sequence. Unlike the amplitude modulation (AM) and amplitude shift keying (ASK) methods [30], the PM-based method offers the ability to embed information toward communication receivers, regardless of whether they are located within the sidelobe or the main lobe.
Since target detection is the main task of the dual-function radar-communications system, the transmit waveform should be considered primarily based on the requirements of the radar function. One fundamental requirement of the radar function is the high efficiency of transmitted power. So constant-modulus waveforms are selected. FH waveform enjoys the constant-modulus feature and is easy to generate. Furthermore, FH pulse waveform is immune to multipath channel fading and clutter interference under the circumstances.
The configuration for a joint radar-communications system and the PM-based FH signal model was developed in [44]. In this section, we follow the methods of [44] and further develop a PM-based dual-function MIMO radar-communications system. We express the MIMO FH waveforms as
ϕ m ( t ) = q = 1 Q e j 2 π c m , q Δ f t u ( t Δ t ) , m = 1 , , M
where c m , q , m = 1 , , M ; q = 1 , , Q describes the FH code, M and Q denote the number of transmit antennas and the length of FH code, respectively. Δ f and Δ t respectively denote the frequency step and the hopping interval duration, and
u ( t ) = { 1 , 0 < t < Δ t 0 , o t h e r w i s e
The duration of a pulse is expressed as T 0 = Q Δ t . We assume that the code is c m , q { 1 , , J } , where J denotes a predefined value. As a result, the bandwidth of the radar pulse can be approximately denoted by J Δ f . To realize transmit waveform orthogonality, FH code should be designed to meet the following requirement c m , q = c m , q , q , m = m . Several papers proposed FH waveform optimization for MIMO systems (see [19]; and references therein). Then, we propose that a dual-function MIMO radar-communications system with phase symbols yields the extended virtual data model at the radar receiver. We also present the information-embedding scheme and the associated transmission rate. Let { Ω m , q [ 0 , 2 π ] } , m = 1 , , M , q = 1 , , Q be a set of M Q PM-based symbols. Hence, the set of FH waveforms is expressed as:
x m ( t ) = q = 1 Q e j Ω m , q e j 2 π c m , q Δ f t u ( t Δ t ) , m = 1 , , M
In (3), Δ t should be designed to meet the requirement:
0 Δ t e j 2 π c m , q Δ f t e j 2 π c m , q Δ f t d t = 0 , m = m , q = q
By using (1) and the orthogonality between the FH waveforms stated in (4), it is easy to prove that PM-based FH waveforms x m ( t ) , m = 1 , , M are also orthogonal. The set of transmitted waveforms is presented in Figure 1.
The N × 1 vector of the target scattering signals at the output of the radar receive antenna array can be expressed as:
y ( t , i ) = Η T ( i ) x ( t , i ) + n ( t , i )
where i describes the index. Η ( i ) = [ h m , n i ] M × N indicates the target response matrix during the i-th radar scan. h m , n stands for the channel coefficient between the m-th transmit element and the n-th receive element. x ( t , i ) = [ x 1 ( t , i ) , , x M ( t , i ) ] T represents the M × 1 vector of PM-based waveforms during the i-th radar scan. x m ( t , i ) is the waveform transmitted from the m-th transmit element, and n ( t , i ) is an N × 1 vector of zero-mean white Gaussian noise.
The target-scattering signal components associated with the individual PM-based FH waveforms are achieved by using a matched filter bank. As a result, the target scattering signals observed at the output of the MIMO radar receiver (5) are matched-filtered to the proposed PM-based sequences (3), yielding the M N × 1 vector of the virtual signal. The extended vector can be expressed as:
r ( i ) = v e c ( T 0 y ( t , i ) x H ( t ) d t )   = v e c ( H T ) + n ^ ( i )
In (6), v e c ( ) describes the vectorization operator that stacks the columns of a matrix into one long column vector, and n ^ ( i ) denotes the M N × 1 zero-mean white Gaussian noise with covariance δ z 2 I M N .

2.2. Information-Embedding Scheme

It is assumed that a phase symbol denotes B bits of binary sequence. During the i-th radar pulse, the binary sequence that needs to be embedded is used to select phase symbols Ω m , q ( i ) , m = 1 , , M , q = 1 , , Q from a pre-defined constellation of K = 2 B symbols. It is considered that the constellation is uniformly distributed between [ 0 , 2 π ] , which can be expressed as:
P S K = { 0 , 2 π K , , ( K 1 ) 2 π K }
The set of PM-based FH waveforms at time i can be rewritten as:
x m ( t , i ) = q = 1 Q e j Ω m , q ( i ) e j 2 π c m , q Δ f t u ( t Δ t ) , m = 1 , , M
In (8), Ω m , q ( i ) P S K . To simplify discussion and implementation, it is assumed that a communication receiver equipped with a single antenna is located at a known direction θ c . Therefore, the received signal at the output of the communication receiver can be expressed as:
y ( t , i ) = α c a T ( θ c ) x ( t , i ) + n ( t , i )
where α c denotes the channel attenuation coefficient. It is considered that α c keeps constant during the whole processing interval. a ( θ c ) indicates the steering vector of the transmit antenna array, and n ( t , i ) denotes the zero-mean white Gaussian noise with covariance δ w 2 .
It is assumed that the communication receiver has full knowledge of the code c m , q and the step Δ f . Hence, the communication received signals (5) are matched-filtered to the orthogonal FH sub-pulses yields:
r m , q ( i ) = 0 Δ t y ( t , i ) e j 2 π c m , q Δ f t u ( t Δ t ) d t = α c e j π d m sin θ c e Ω m , q ( τ ) + n m , q ( i )     m = 1 , , M , q = 1 , , Q
In (10), d m denotes the distance between the first and the m-th antennas of the transmit array. The communication receiver has the ability to undo e j π d m sin θ c before it estimates the symbol Ω m , q ( i ) . r m , q ( i ) is viewed as a phase-rotated and noisy version of the m-th entry of a ( θ c ) . Hence, Ω m , q ( i ) that need to be embedded is estimated as:
Ω ^ m , q ( i ) = A n g l e ( r m , q ( i ) ) A n g l e ( α c ) + 2 π d m sin θ c
Then, the communication receiver compares Ω ^ m , q ( i ) to the predefined P S K . Ω m , q ( i ) is restored from r m , q ( i ) at the output of the (m,q)-th matched filter. It allows the communication receiver to determine Ω m , q ( i ) and convert Ω m , q ( i ) into the original sequence. The advantages of the communication information embedding scheme have been presented in [44].

3. Multiple-Input Multiple-Output (MIMO) Radar Ambiguity Function

In this section, we consider a radar target at χ ( τ , v , f ) where τ describes the delay of the target range, v denotes the Doppler frequency, f = 2 π d R λ sin θ indicates the spatial frequency. Here θ and λ are the angle of the target and the wavelength, respectively. d R and d T are the distance between the transmit antennas and between the receiver antennas, respectively. To simplify the discussion, we assume d T = d R in this paper. x m ( t ) , m = 1 , , M denotes the waveform radiated by the m-th transmit antenna. In [17], the MIMO radar ambiguity function can be defined as:
χ ( τ , v , f , f ) = m = 1 M m = 1 M χ m , m ( τ , v ) e j 2 π ( f m f m )
where
χ m , m ( τ , v ) = x m ( t ) x m ( t + τ ) e j 2 π v t d t
In (13), χ m , m ( τ , v ) describes the cross ambiguity function, which implicates two radar waveforms x m ( t ) and x m ( t ) . χ m , m ( τ , v ) is analogous to the SIMO ambiguity function presented in [19]. We now discuss the MIMO radar ambiguity function for the case when x m ( t ) is composed of the shifted forms of a rectangular pulse u m ( t ) .
x m ( t ) = l = 1 L u m ( t T l )
where l denotes the number of the rectangular pulse and T describes pulse repetition interval (PRI). The cross ambiguity function of the waveform x m ( t ) is defined as:
χ m , m ( τ , v ) = l = 1 L l = 1 L χ m , m u ( τ + T l T l , v ) e j 2 π v T l
where χ m , m u ( τ , v ) describes the cross ambiguity function of the rectangular pulses u m ( t ) and u m ( t ) . We assume that the pulse duration T 0 and the Doppler frequency v are small enough so that T 0 v 0 . At a result, the envelope of the Doppler frequency keeps unchanged within an entire pulse period. The cross ambiguity function χ m , m u ( τ , v ) reduces to the cross correlation function r m , m u ( τ ) , which is not a function of v any more, we obtain:
χ m , m u ( τ , v ) r m , m u ( τ )
We further assume that little reflections take place at these second trip ranges. Therefore, the cross ambiguity function of x m ( t ) can be rewritten as:
χ m , m ( τ , v ) r m , m u ( τ ) l = 1 L e j 2 π v T l
The MIMO ambiguity function of x m ( t ) can be expressed as:
χ ( τ , v , f , f ) = m = 1 M m = 1 M r m , m u ( τ ) e j 2 π ( f m f m ) l = 1 L e j 2 π v T l
It is worth noting that the MIMO ambiguity function χ ( τ , v , f , f ) depends on the cross correlation functions r m , m u ( τ ) . Furthermore, the shifted forms of a rectangular pulse { u m ( t ) } just have an impact on the range and spatial resolution. The pulses have no effect on the Doppler resolution. Consequently, to gain a sharp MIMO ambiguity function, these waveforms should be designed such that the function Ω ( τ , f , f ) is denoted as:
Ω ( τ , f , f ) = m = 1 M m = 1 M r m , m u ( τ ) e j 2 π ( f m f m )
Next, we discuss the MIMO radar ambiguity function of the PM-based information embedding FH waveforms. The proposed waveforms can be expressed as (3). We intend to obtain the expression for the function Ω ( τ , f , f ) in terms of { c m , q } and { Ω m , q } . To derive Ω ( τ , f , f ) , we therefore begin with the cross correlation function of the PM-based FH waveform. Making use of (3) and (16), the cross correlation function can be expressed as:
r m , m x ( τ ) = q = 1 Q q = 1 Q χ r e c t ( τ ( q q ) Δ t , ( c m , q c m , q ) Δ f ) e j Ω m , q e j 2 π Δ f ( c m , q c m , q ) q Δ t e j 2 π Δ f c m , q τ
and the function Ω ( τ , f , f ) can be expressed as:
Ω ( τ , f , f ) = m , m = 1 M q , q = 1 Q χ r e c t ( τ ( q q ) Δ t , ( c m , q c m , q ) Δ f ) e j Ω m , q e j 2 π Δ f ( c m , q c m , q ) q Δ t e j 2 π Δ f c m , q τ e j 2 π ( f m f m )
where χ r e c t ( τ , v ) describes the ambiguity function of u ( t ) , which can be denoted by:
χ r e c t ( τ , v ) = { Δ t | τ | Δ t sin c ( v ( Δ t | τ | ) ) e j π v ( τ + Δ t ) ,   if | τ | < Δ t 0 ,   otherwise
For M = 1 , the function is the special case of the SIMO radar. For the general case of the MIMO radar M > 1 , not only the auto-correlation functions but also the cross-correlation functions between the waveforms should be taken into account such that the function (22) is sharp around { ( τ , f , f ) | τ = 0 , f = f } .

4. Waveform Optimization

In this section, we aim to enhance the target detection and feature estimation performance by maximizing the MI between the target response and the target returns in the first step, and then minimizing the MI between successive target scattering signals in the second step. These two stages correspond to the design of the ensemble of excitations and the selection of a suitable signal out of the ensemble, respectively. The two-step cognitive waveform design strategy is based upon continuous learning from the radar scene. The dynamic information about the target feature is utilized to design PM-based FH codes. In this way the transmitter adjusts its probing signals to suit the dynamically changing environment.
Step 1: this step involves the design of PM-based FH waveforms for the dual-function transmit array. The main idea of waveform design is to maximize the MI between the target scattering signal and the target response, subject to the transmit power constraint.
Step 2: once an ensemble of optimal transmit waveforms has been designed, then we select the most reasonable waveform for emission from the ensemble. The key concept of waveform selection is to minimize the MI between the target scattering signals at present and the next target-scattering signals. The step ensures that we continually obtain target scattering signals that are independent of each other in time, in order to achieve more feature information of the target at each time instant of reception. Figure 2 describes the architecture of an adaptive dual-function MIMO radar communication system.
Step 1: during the i-th radar scan, the set of PM-based FH waveforms x m ( t , i ) , m = 1 , , M can be expressed as a matrix X i M × K after discrete sampling, where describes the number domain and K indicates the sample number. N N × K is zero-mean white noise matrix. Therefore, the N × K matric of the target scattering signals is written as:
Y i = H i T X i + N
Here H i Ν ( 0 , R H i ) and N Ν ( 0 , R N ) , and R H i = E { H i H i H } and R N = E { N H N } indicate the covariance matrices of the channel response H i and the zero-mean white noise N , respectively. We intend to maximize the MI between the target-scattering signal and the target response given the transmit waveforms. This involves that the target-scattering signals would be more dependent upon the actual target feature information. According to the definition of MI, we have:
I ( Y i ; H i | X i ) = H ( Y i | X i ) H ( N )
In (24), I ( Y i ; H i | X i ) represents the MI between two random variates Y i and H i given the transmit matrix X i , and H ( Y i | X i ) indicates the conditional entropy that X i conveys about Y i . The main objective of this step is to maximize I ( Y i ; H i | X i ) between Y i and H i given X i . According to the definition of entropy, we have:
H ( Y i | X i ) = p ( Y i | X i ) I n [ p ( Y i | X i ) ] d Y i
Here p ( Y i | X i ) indicates the conditional probability density function (PDF) of the received matrix Y i given transmit matrix X i . The conditional PDF p ( Y i | X i ) can be expressed as follows:
p ( Y i | X i ) = n = 1 N p ( y n , i | X i ) = n = 1 N 1 π K det ( X i H R H i X i + R N ) e y n , i H ( X i H R H i X i + R N ) 1 y n , i   = 1 π N K [ det ( X i H R H i X i + R N ) ] N e t r [ ( X i H R H i X i + R N ) 1 Y i H Y i ]
Solving (25) and (26) gives rise to the following result for the conditional entropy [40]:
H ( Y i | X i ) = N K I n ( π ) + N K + N I n [ det ( X i H R H i X i + R N ) ]
Similarly, the result for the entropy of the noise can be derived as follows:
H ( N ) = N K I n ( π ) + N K + N I n [ det ( R N ) ]
Making use of Equations (24), (27) and (28), the MI between the target-scattering signal and the channel response given transmit waveforms can be rewritten as follows:
I ( Y i ; H i | X i ) = N I n [ det ( X i H R H i X i + R N ) ] N I n [ det ( R N ) ]
Therefore, we can formulate the MI maximization criterion as follows:
max X i { N I n [ det ( X i H R H i X i + R N ) ] N I n [ det ( R N ) ] } n ! r ! ( n r ) ! s . t .   tr [ X i H X i ] P 0           c m , q c m , q ,   for   m = m , q       Δ t Δ f = 1
where P 0 indicates the transmit power. The work of [40] has given a rigorous solution of the above optimization problem (30). Then, we can obtain the ensemble X ^ i out of the whole set of PM-based FH sequences, and the corresponding power allocation vector over diverse dual-function transmit antennas. We start the design procedure with the transmit waveforms from the PM-based FH matrix, and allocate the power on a pulse level as well as across the transmit antennas based on the MI maximization criterion.
Step 2: we then proceed to the waveform selection procedure, in which the successive target scattering signals are different from each other. The step ensures that we achieve more information of the target feature at each time instant of reception. We denote the MI between the successive target scattering signals at time i and at time i + 1 as:
I ( Y i + 1 , Y i ) = H ( Y i + 1 | X i + 1 ) + H ( Y i | X i ) H ( Y i + 1 , Y i | X i + 1 , X i )
In (31), the term H ( Y i | X i ) (or H ( Y i + 1 | X i + 1 ) ) denotes the entropy of Y i (or Y i + 1 ) at time i (or i + 1 ) given the knowledge of X i (or X i + 1 ). The term H ( Y i + 1 , Y i | X i + 1 , X i ) is defined similarly. According to the literature [6,34], the above Equation (31) can be rewritten as follows:
I ( Y i , Y i + 1 ) = N I n ( det ( I ( M × M ) D i , i + 1 2 ) ) = N m = 1 M I n ( 1 ( d i , i + 1 m ) 2 )
where d i , i + 1 m ( d i , i + 1 1 d i , i + 1 2 d i , i + 1 M ) is the diagonal element of the diagonal matrix D i , i + 1 . D i , i + 1 is achieved by singular value decomposition (SVD) of the covariance matrix, which can be denoted as:
R Y i , Y i + 1 = E { Y i H Y i + 1 } = X i H R H i X i + 1
In (33), R Y i , Y i + 1 is the cross-covariance of the expressions for Y i and Y i + 1 . Hence, the MI minimization criterion can be expressed as:
min X i + 1 X ^ i { N m = 1 M I n { 1 ( d i , i + 1 m ) 2 } } s . t .   tr [ X i + 1 H X i + 1 ] P 0 c m , q c m , q ,     for   m = m , q     Δ t Δ f = 1
We can estimate the values for Y i + 1 over all possible values of X ^ i + 1 X ^ i using (23). Thus, we can also form an estimate of all the values of the corresponding D i , i + 1 ( d i , i + 1 m ( d i , i + 1 1 d i , i + 1 2 d i , i + 1 M ) ) and choose the value for X ^ i + 1 that minimizes (34). The above waveform optimization problem (34) is convex. We can obtain the optimal solution directly by using a MATLAB optimization toolbox, such as CVX.
Step 1 designs the optimal ensemble X ^ i based on MI maximization criterion over the spatial domain, and step 2 selects the optimized PM-based FH waveform for each dual-functional transmit antenna from the ensemble X ^ i based on MI minimization criterion over the temporal domain. The proposed adaptive waveform design and selection procedures can be summarized as Algorithm 1.
Algorithm 1. The adaptive waveform design and selection algorithm
Initialize the transmit matrix X 0 and the covariance matrix of the noise R N .
1. At the initial time t = 0 , measure the target scattering signal Y 0 and calculate the covariance of the target-scattering signal R Y 0 . The covariance matrix of the channel response R H 0 can be obtained through successive measurements with uniform power allocation over the transmit antennas by solving (36).
2. The optimized waveform ensemble X ^ 0 can be obtained through the waveform design process by solving (30).
3. At time t = 1 , measure Y 1 and calculate R Y 1 and the cross-covariance matrix R Y 0 , Y 1 . The actual channel response H 1 is obtained by de-convolving the target-scattering signal with the transmit signal by using (23).
4. Calculate the corresponding singular values d 0 , 1 m , m = 1 , , M of R Y 0 , Y 1 and X ^ 1 X ^ 0 can be acquired through the waveform selection process by solving (38).
5. At time t = 1 , emission X ^ 1 and observe the corresponding received signal Y ^ 1 to achieve R H 1 by using (36).
6. If i = I max , the iterative procedure ends; or else, repeat steps 2–5 iteratively.
It is worth noting that adaptation is included in the waveform design and selection procedures through the feedback process and numerous interactions with the radar channel. In summary, the adaptive waveform optimization strategy for target detection is implemented according to the block diagram in Figure 3.

5. Simulation

In this section, numerical results based on Monte Carlo simulations have been provided to validate the effectiveness of the proposed method. Without loss of generality, each entry of the channel matrices follows the standard complex Gaussian distribution. The simulation parameters are based on radar application with a high Pulse Repetition Frequency (PRF), such as in X-band radar. A data rate in the range of dozens of Mbps can be achieved. We provide a comparison between the proposed scheme and the method of [32]. To implement the method in [32], we consider a dual-function MIMO system operating in the X-band with carrier frequency f c = 8.2 GHz and bandwidth B = 500 MHz. The sampling frequency is f s = 10 9 sample/sec, which is taken as the Nyquist rate. The PRI is T 0 = 10   μ s . We assume an arbitrary linear transmit array consisting of M = 16 elements. We further assume that the minimum transmit/receive antenna spacing is sufficiently larger than half wavelength (distributed MIMO configuration). Hence, the correlation introduced by finite antenna element spacing is low enough that the fades associated with two different antenna elements can be considered independent. To implement the radar function, we further assume that the FH step is Δ f = 10   MHz , the length of the FH code is Q = 20 and the FH interval duration is Δ t = 0.1   μ s . We generate a set of 16 FH pulse waveforms. The parameter J = 50 is used. Therefore, the 320 FH code is generated randomly from the set { 1 , 2 , , J } , where J = 50 .
We employ orthogonal sequences of the FH pulse over the transmit antenna elements. The backscatter signals are matched filtered at the receivers and the transmitted signals are later modified by the waveform optimization module as shown in Figure 2. The optimized transmission sequence at one particular transmit antenna after the two-step optimization process. The target response extracted from the received target echoes after matched filtering at the end of 20 iterations of the algorithm, where an excellent performance of the target response estimation can be observed. At each iteration, the scattering coefficients for the target and non-target scatterers in H i vary as described by the Swerling III model [45,46]. This causes the amplitude returns of the backscatter signals to vary at each instance. However, the amplitudes of the echoes from the target are always assumed to be stronger than those from the clutter sources.

5.1. Target Detection Performance

Figure 4a illustrates the detection probability offered by the proposed scheme versus the signal-to-noise ratio (SNR) for different iterations. The iteration process is run 20 times. All optimal waveforms are generated by the proposed two-step algorithm. The value of requested SNR increases as the probability increase for a particular number of iterations. The SNR decreases as the number of iterations increase for a certain detection probability. The detection performance offered by the proposed scheme improves as the number of iterations increases. Simulation results show that, on average, 20 iterations of the waveform optimization algorithm are required in order to achieve convergence of the target response estimation for a wide range of radar scenes.
Figure 4b illustrates the ROC for the following four approaches while the value of received SNR equals to 8 dB. (1) 4 × 4 MIMO system based on conventional maximum a posterior (MAP) approach; (2) 4 × 4 MIMO system using the Kalman filtering [42]; (3) 4 × 4 MIMO system based on MI optimization (step 2) scheme; (4) 4 × 4 dual-function MIMO radar-communications system using the proposed scheme.
The plot is run at the end of 24 iterations. For the probability of false alarm P f a = 0.005 , the of target detection probability generated by the proposed scheme is 0.95 as compared with 0.75 provided by the Kalman filtering method, 0.7 by MI optimization (step 2) criterion and 0.55 by conventional MAP approach. As the two-step scheme can use the temporal correlation of target characteristic during the radar scan interval, the dual-function MIMO system constantly adapts transmit mode to suit the dynamic radar scene. Furthermore, the successive target returns are considered as independent of each other. The property guarantees that information about the radar scene is learned at each instant of reception. In this case, the detection performance offered by the proposed scheme is best.

5.2. Target Response Estimation Performance

Figure 5a displays the MSE achieved by the proposed scheme with regard to the estimation of target response. This plot demonstrates an improved MSE performance for the two-step optimization as compared with the conventional MAP approach, the Kalman filtering method and the MI optimization (step 2) modules, particularly for the first few iterations.
It is evident from Figure 5 that the MSE performance offered by the proposed scheme is superior to the conventional MAP approach and the Kalman filtering method. The MSE performance offered by the proposed scheme is superior to MI optimization (step 2) criterion. Similarly, Figure 5b shows the MSE performance offered by the proposed scheme with respect to the estimate of target response. The MSE of target response estimation provided by the two-step scheme and other several approaches are compared to reveal the benefit of the proposed scheme.

5.3. Delay-Doppler Resolution Performance

Figure 6a displays the 2 × 2 MIMO radar ambiguity function contours of the original waveform. The delay-Doppler resolution deteriorates owing to the presence of environment noise as shown in Figure 6a. This phenomenon will became worse if the interference scatterers are placed in the vicinity along the line linking the radar target and the antenna. The radar target is considered to be located at the origin of the plan. Figure 6b displays the MIMO radar ambiguity function contours of the optimization waveform by the proposed scheme at the end of 20 iterations. As can be seen from Figure 6a, the target discrimination capability becomes significantly enhanced by the increasing iteration number. Specifically, the noise is suppressed by about 2.5 dB.
The delay-Doppler resolution, which is the output of matched filter at the radar receiver, is connected to the MIMO radar ambiguity function of the transmit waveforms. The near-ideal thumbtack response would appear if the statistically independent waveforms are employed for emission. We constantly use the optimization waveforms provided by the proposed scheme to match the estimated target response. It is worth noting that the estimated target response is continuously updated at each iteration. At a result, the matching process ensures that the noise interference can be suppressed over the radar channels and further enhance the SNR of the target-scattering signals. Figure 6b illustrates the enhancement in SNR with multiple iteration number. The phenomenon shows the enhanced ability of the proposed dual-function MIMO system to discriminate the target from radar environment and resolve the target’s range and velocity.

5.4. Symbol Error Rate Performance

We study the SER performance of the communication source embedding scheme using Binary Phase Shift Keying (BPSK), Quadrature Phase Shift Keying (QPSK), 16-Phase Shift Keying (16-PSK) and 256-Phase Shift Keying (256-PSK) constellations and compare the original waveform with the encoded waveform. The original waveform corresponds to data rate of R = 32; 64; 128 and 256 Mbps, respectively. To obtain encoded waveform, a convolutional encoder of rate 2/3 is employed in the original waveform. We use a Viterbi decoder to decode the received encoded waveform at the communication receiver. The communication channel coefficient is taken as | α c | = 1 , and the phase of α c is uniformly distributed within the interval [ 0 , 2 π ] . To test the SER performance, we generate a number of 10 × 10 17 random BPSK, QPSK, 16-PSK and 256-PSK symbols. The SERs versus SNR for all constellation sizes is illustrated in Figure 7.
The results demonstrate that the smaller the constellation size is, the better the SER performance will be. As the constellation size increases, it is more difficult to detect the symbols. It can be explained that the defective cross-correlation between the non-orthogonal transmit sequences influences the detection performance. Obviously, the encoded waveform achieves greater SER performance as compared to the original waveform. Therefore, intersymbol interference is a significant source of detection error resulting in performance degradation. It is expected that the SER performance gets worse if longer FH waveforms are used. It is worth noting that, for all techniques tested, the SER performance offered by the encoded waveform outperforms the performance provided by the original waveform. Therefore, the SER superiority comes at the price of lower data rate.
Figure 8 illustrates the throughput result provided by the proposed optimization waveform versus distance for BPSK, QPSK, 16-PSK, and 256-PSK constellation. 256-PSK waveform provides a data rate of approximately 8 Mbps at a distance of 10 m, which is better than that generated by BPSK, QPSK, 16-PSK constellations. The highest data rate is acquired by the 256-PSK waveform within a distance of 60 m, as the distance between the system nodes increases the data rate decreases.

5.5. Detection Variation Performance

The detection constraint optimization has recently been studied in works such as [38,39], where the authors address the problem of radar code design for target recognition in the presence of colored Gaussian disturbance. The objective function in [38] aims to maximize the weighted average Euclidean distance between the ideal echoes from different target hypotheses. Furthermore, additional practical constraints are considered in [38]. For example, the modulus of the waveform is restricted to be a constant and the detection constraints require that the achievable SNR for each target hypothesis is larger than a desired threshold.
Figure 9 indicates the detection variation offered by the proposed scheme subject to the detection constraint. We assume a radar scene, which has three range-separated targets. The target-scattering signals derived from the radar scene is normalized and the proposed dual-function MIMO system intends to discriminate the scatterers by using a particular detection threshold.
With subsequent iterations of the proposed algorithm, the detection performance of the multiple targets is enhanced. As can be seen from Figure 9, by suppressing noise, the dual-function MIMO system could discriminate three scatterers effectively at the end of 20 iterations. The result proves the performance of target detection presented in Figure 6b as well.
The detection performance is enhanced by providing the waveform design (step 1) procedure in the proposed scheme. The waveform design procedure ensures the maximum of the Euclidean distance between the target echoes from different scatterers based on the MI maximization criterion in (30). The proposed design procedure is similar to the optimization method in [43]. By increasing the number of iterations, we can obtain more accurate estimates of the target responses, which are used to improve detection of the multiple targets.

5.6. Comparison with Other Methods

We have assess the performance of the new approach in terms of the Bit Error Rate (BER) and compare it with the method in [29]. Note that the latter approach uses a single sequence in tandem with 4 sidelobe levels towards the communication direction to deliver 2 bits of information. On the other hand, we use a PM-based FH waveform to deliver two bits per pulse. To test the BER, a sequence of symbols unencoded (two bits each) is used. Furthermore, a convolutional encoder of rate 2/3 is used in the original waveform leading to encoded waveform. Both the unencoded and the encoded waveforms are embedded independently using the approach of [29] and the new approach. The received encoded signal is decoded using a Viterbi decoder. The BERs versus the signal-to-noise ratio (SNR) for the two approaches is presented in Figure 10 for both the unencoded as well as the encoded data sequences. Obviously, the proposed approach achieves better BER performance compared to the method of [29]. Note that the latter approach embeds 25% of the information via each of the four beams. Thus, intersymbol interference is a non-negligible source of detection error leading to capability degradation. The behavior can be expected to be worse if longer pulses are transmitted. Note that, for two approaches tested, the BER with respect to the encoded waveform outperforms the BER with respect to the unencoded waveform. However, this BER superiority comes at the price of slower data transmission rate.
In Figure 11, we compare the performance for adaptive waveform provided by the proposed method to the performance for a static waveform [32] over multiple snapshots. As the proposed method selects specific waveforms, which generate target returns having low correlation over time, the system adapts its transmit waveform better to the fluctuating target Radar Cross-Section (RCS). On the other hand, the static waveform [32] is unable to match the time-varying target response. Therefore, the performance of the static waveform [32] is worse than the proposed adaptive waveform.

6. Conclusions

The problem of dual-functionality system design for joint radar and communications operation was considered. We further presented a two-step waveform optimization algorithm for the proposed system, which combines the waveform optimization and selection processes. The proposed algorithm is based upon constant learning of the radar scene at the receiver and adaptation of the probing codes to suit the dynamic target feature. The adaptive process ensures maximum information extraction from the target of interest. The effectiveness of the proposed technique and its superiority over existing techniques in terms of the BER performance and the target detection performance were demonstrated through extensive simulations. The proposed system would form a joint platform for future intelligent transportation applications for which both environmental perception and establishment of data links are crucial. Future research will look into the tradeoff between the performance improvement offered by the proposed approach and the computational complexity involved.

Author Contributions

Conceptualization, Y.Y. and X.L.; Methodology, Y.Y.; Software, X.L.; Validation, Y.Y., X.L. and L.W.; Formal Analysis, Y.Y.; Investigation, Y.Y.; Resources, X.L.; Data Curation, X.L.; Writing-Original Draft Preparation, Y.Y.; Writing-Review & Editing, X.L.; Visualization, Y.Y.; Supervision, L.W.; Project Administration, L.W.; Funding Acquisition, Y.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (61761019, 61861017, 61861018, 61862024) and the Natural Science Foundation of Jiangxi Province (Jiangxi Province natural Science Fund) (20181BAB211014, 20181BAB211013) and Foundation of Jiangxi Educational Committee of China (GJJ180352).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, J.; Stoica, P. (Eds.) MIMO Radar Signal Processing; Wiley: Hoboken, NJ, USA, 2009. [Google Scholar]
  2. Li, B.; Petropulu, A.P.; Trappe, W. Optimum co-design for spectum sharing between matrix completion based MIMO radars and a MIMO communication system. IEEE Trans. Signal Process. 2016, 64, 4562–4575. [Google Scholar] [CrossRef] [Green Version]
  3. 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]
  4. Liu, S.; Zhang, Y.D.; Shan, T.; Tao, R. Structure-Aware Bayesian Compressive Sensing for Frequency-Hopping Spectrum Estimation with Missing Observations. IEEE Trans. Signal Process. 2018, 66, 2153–2166. [Google Scholar] [CrossRef]
  5. Fan, W.; Liang, J.; Li, J. Constant Modulus MIMO Radar Waveform Design with Minimum Peak Sidelobe Transmit Beampattern. IEEE Trans. Signal Process. 2018, 66, 4207–4222. [Google Scholar] [CrossRef]
  6. Yao, Y.; Zhao, J.; Wu, L. Adaptive Waveform Design for MIMO Radar-Communication Transceiver. Sensors 2018, 18, 1957. [Google Scholar] [CrossRef] [Green Version]
  7. Wu, L.; Babu, P.; Palomar, D.P. Transmit Waveform/Receive Filter Design for MIMO Radar with Multiple Waveform Constraints. IEEE Trans. Signal Process. 2017, 66, 1526–1540. [Google Scholar] [CrossRef]
  8. Stoica, P.; Li, J.; Xie, Y. On Probing Signal Design for MIMO Radar. IEEE Trans. Signal Process. 2007, 55, 4151–4161. [Google Scholar] [CrossRef] [Green Version]
  9. Tang, B.; Li, J. Spectrally Constrained MIMO Radar Waveform Design Based on Mutual Information. IEEE Trans. Signal Process. 2018, 67, 821–834. [Google Scholar] [CrossRef]
  10. Yang, Y.; Blum, R. MIMO radar waveform design based on mutual information and minimum mean-square error estimation. IEEE Trans. Aerosp. Electron. Syst. 2007, 43, 330–343. [Google Scholar] [CrossRef]
  11. Wang, L.; Zhu, W.; Zhang, Y.; Liao, Q.; Tang, J. Multi-Target Detection and Adaptive Waveform Design for Cognitive MIMO Radar. IEEE Sens. J. 2018, 18, 9962–9970. [Google Scholar] [CrossRef]
  12. Fuhrmann, D.R.; San Antonio, G. Transmit beamforming for MIMO radar systems using signal crosscorrelation. IEEE Trans. Aerosp. Electron. Syst. 2009, 44, 171–186. [Google Scholar] [CrossRef]
  13. Song, X.; Zhou, S.; Willett, P. Reducing the Waveform Cross Correlation of MIMO Radar with Space–Time Coding. IEEE Trans. Signal Process. 2010, 58, 4213–4224. [Google Scholar] [CrossRef]
  14. Jajamovich, G.H.; Lops, M.; Wang, X. Space-Time Coding for MIMO Radar Detection and Ranging. IEEE Trans. Signal Process. 2010, 58, 6195–6206. [Google Scholar] [CrossRef]
  15. Cheng, Z.; Liao, B.; He, Z.; Li, Y.; Li, J. Spectrally Compatible Waveform Design for MIMO Radar in the Presence of Multiple Targets. IEEE Trans. Signal Process. 2018, 66, 3543–3555. [Google Scholar] [CrossRef]
  16. Yang, Y.; Blum, R.S.; He, Z.S.; Fuhrmann, D.R. MIMO radar waveform design via alternating projection. IEEE Trans. Signal Process. 2010, 58, 1440–1445. [Google Scholar] [CrossRef]
  17. Antonio, G.S.; Fuhrmann, D.R.; Robey, F.C. MIMO Radar Ambiguity Functions. IEEE J. Sel. Top. Signal Process. 2007, 1, 167–177. [Google Scholar] [CrossRef]
  18. Chen, C.-Y.; Vaidyanathan, P. MIMO Radar Ambiguity Properties and Optimization Using Frequency-Hopping Waveforms. IEEE Trans. Signal Process. 2008, 56, 5926–5936. [Google Scholar] [CrossRef] [Green Version]
  19. Gogineni, S.; Nehorai, A. Frequency-Hopping Code Design for MIMO Radar Estimation Using Sparse Modeling. IEEE Trans. Signal Process. 2012, 60, 3022–3035. [Google Scholar] [CrossRef]
  20. Han, K.; Nehorai, A. Jointly optimal design for mimo radar frequency-hopping waveforms using game theory. IEEE Trans. Aerosp. Electron. Syst. 2016, 52, 809–820. [Google Scholar] [CrossRef]
  21. Griffiths, H.; Cohen, L.; Watts, S.; Mokole, E.; Baker, C.; Wicks, M.; Blunt, S. Radar spectrum engineering and management: Technical and regulatory issues. Proc. IEEE 2015, 103, 85–102. [Google Scholar] [CrossRef]
  22. Chiriyath, A.R.; Paul, B.; Jacyna, G.M.; Bliss, D.W. Inner Bounds on Performance of Radar and Communications Co-Existence. IEEE Trans. Signal Process. 2015, 64, 464–474. [Google Scholar] [CrossRef]
  23. Ahmed, A.; Zhang, Y.D.; Gu, Y. Dual-function radar-communications using QAM-based sidelobe modulation. Digit. Signal Process. 2018, 82, 166–174. [Google Scholar] [CrossRef]
  24. Bliss, D.W. Cooperative radar and communications signaling: The estimation and information theory odd couple. In Proceedings of the 2014 IEEE Radar Conference, Cincinnati, OH, USA, 19–23 May 2014; Institute of Electrical and Electronics Engineers (IEEE): New York, NY, USA, 2014; pp. 0050–0055. [Google Scholar]
  25. Geng, Z.; Deng, H.; Himed, B. Adaptive radar beamforming for interference mitigation in radar-wireless spectrum sharing. IEEE Signal Process. Lett. 2015, 22, 484–488. [Google Scholar] [CrossRef]
  26. Junhui, Z.; Tao, Y.; Yi, G.; Jiao, W.; Lei, F. Power control algorithm of cognitive radio based on non-cooperative game theory. China Commun. 2013, 10, 143–154. [Google Scholar] [CrossRef]
  27. Khawar, A.; Abdelhadi, A.; Clancy, C.; Clancy, T. Target Detection Performance of Spectrum Sharing MIMO Radars. IEEE Sens. J. 2015, 15, 1. [Google Scholar] [CrossRef] [Green Version]
  28. Ahmed, A.; Zhang, Y.D.; Himed, B. Distributed Dual-Function Radar-Communication MIMO System with Optimized Resource Allocation. In Proceedings of the 2019 IEEE Radar Conference (RadarConf), Boston, MA, USA, 22–26 April 2019; Institute of Electrical and Electronics Engineers (IEEE): New York, NY, USA, 2019; pp. 1–5. [Google Scholar]
  29. Hassanien, A.; Amin, M.G.; Zhang, Y.D.; Ahmad, F. Dual-function radar-communications: Information embedding using sidelobe control and waveform diversity. IEEE Trans. Signal Process. 2016, 64, 2168–2181. [Google Scholar] [CrossRef]
  30. Hassanien, A.; Amin, M.G.; Zhang, Y.D.; Ahmad, F. Signaling strategies for dual-function radar communications: An overview. IEEE Aerosp. Electron. Syst. Mag. 2016, 31, 36–45. [Google Scholar] [CrossRef]
  31. Hassanien, A.; Amin, M.G.; Zhang, Y.D.; Ahmad, F. Phase-modulation based dual-function radar-communications. IET Radar, Sonar Navig. 2016, 10, 1411–1421. [Google Scholar] [CrossRef] [Green Version]
  32. Hassanien, A.; Himed, B.; Rigling, B.D. A dual-function MIMO radar-communications system using frequency-hopping waveforms. In Proceedings of the 2017 IEEE Radar Conference (RadarConf), Seattle, WA, USA, 8–12 May 2017; Institute of Electrical and Electronics Engineers (IEEE): New York, NY, USA, 2017; pp. 1721–1725. [Google Scholar]
  33. Ahmed, A.; Gu, Y.; Silage, D.; Zhang, Y.D. Power-Efficient Multi-User Dual-Function Radar-Communications. In Proceedings of the 2018 IEEE 19th International Workshop on Signal Processing Advances in Wireless Communications (SPAWC), Kalamata, Greece, 25–28 June 2018; Institute of Electrical and Electronics Engineers (IEEE): New York, NY, USA, 2018; pp. 1–5. [Google Scholar]
  34. Yao, Y.; Zhao, J.; Wu, L. Cognitive Radar Waveform Optimization Based on Mutual Information and Kalman filtering. Entropy 2018, 20, 653. [Google Scholar] [CrossRef] [Green Version]
  35. Zhao, J.; Guan, X.; Li, X.P. Power allocation based on genetic simulated annealing algorithm in cognitive radio networks. Chin. J. Electron. 2013, 22, 177–180. [Google Scholar]
  36. Bell, M. Information theory and radar waveform design. IEEE Trans. Inf. Theory 1993, 39, 1578–1597. [Google Scholar] [CrossRef] [Green Version]
  37. Romero, R.A.; Bae, J.; Goodman, N.A. Theory and Application of SNR and Mutual Information Matched Illumination Waveforms. IEEE Trans. Aerosp. Electron. Syst. 2011, 47, 912–927. [Google Scholar] [CrossRef]
  38. Tang, B.; Tang, J.; Peng, Y. MIMO Radar Waveform Design in Colored Noise Based on Information Theory. IEEE Trans. Signal Process. 2010, 58, 4684–4697. [Google Scholar] [CrossRef]
  39. Aubry, A.; De Maio, A.; Huang, Y.; Piezzo, M.; Farina, A. A new radar waveform design algorithm with improved feasibility for spectral coexistence. IEEE Trans. Aerosp. Electron. Syst. 2015, 51, 1029–1038. [Google Scholar] [CrossRef]
  40. Chen, Y.; Nijsure, Y.; Yuen, C.; Chew, Y.H.; Ding, Z.; Boussakta, S. Adaptive Distributed MIMO Radar Waveform Optimization Based on Mutual Information. IEEE Trans. Aerosp. Electron. Syst. 2013, 49, 1374–1385. [Google Scholar] [CrossRef]
  41. Naghibi, T.; Behnia, F. MIMO Radar Waveform Design in the Presence of Clutter. IEEE Trans. Aerosp. Electron. Syst. 2011, 47, 770–781. [Google Scholar] [CrossRef]
  42. Yao, Y.; Zhao, J.; Wu, L. Waveform Optimization for Target Estimation by Cognitive Radar with Multiple Antennas. Sensors 2018, 18, 1743. [Google Scholar] [CrossRef] [Green Version]
  43. Wang, X.; Hassanien, A.; Amin, M.G. Dual-Function MIMO Radar Communications System Design Via Sparse Array Optimization. IEEE Trans. Aerosp. Electron. Syst. 2019, 55, 1213–1226. [Google Scholar] [CrossRef] [Green Version]
  44. Yao, Y.; Zhao, J.; Wu, L. Frequency-Hopping Code Design for Target Detection via Optimization Theory. J. Optim. Theory Appl. 2019, 183, 731–756. [Google Scholar] [CrossRef]
  45. Nijsure, Y.; Chen, Y.; Boussakta, S.; Yuen, C.; Chew, Y.H.; Ding, Z. Novel System Architecture and Waveform Design for Cognitive Radar Radio Networks. IEEE Trans. Veh. Technol. 2012, 61, 3630–3642. [Google Scholar] [CrossRef]
  46. Goodman, N.A.; Venkata, P.R.; Neifeld, M.A. Adaptive Waveform Design and Sequential Hypothesis Testing for Target Recognition with Active Sensors. IEEE J. Sel. Top. Signal Process. 2007, 1, 105–113. [Google Scholar] [CrossRef]
Figure 1. Illustrative diagram of a dual-function multiple-input multiple-output (MIMO) transmit platform.
Figure 1. Illustrative diagram of a dual-function multiple-input multiple-output (MIMO) transmit platform.
Sensors 20 00415 g001
Figure 2. The architecture of an adaptive dual-function MIMO radar communication system.
Figure 2. The architecture of an adaptive dual-function MIMO radar communication system.
Sensors 20 00415 g002
Figure 3. The adaptive waveform design and selection scheme for target detection.
Figure 3. The adaptive waveform design and selection scheme for target detection.
Sensors 20 00415 g003
Figure 4. (a) The detection probability for different iterations provided by the proposed scheme; (b) The receiver operating characteristics (ROC) for four approaches.
Figure 4. (a) The detection probability for different iterations provided by the proposed scheme; (b) The receiver operating characteristics (ROC) for four approaches.
Sensors 20 00415 g004
Figure 5. (a) The mean-square error (MSE) of target response estimation subject to power constraint; (b) the MSE of target response estimation subject to power and probability constraint.
Figure 5. (a) The mean-square error (MSE) of target response estimation subject to power constraint; (b) the MSE of target response estimation subject to power and probability constraint.
Sensors 20 00415 g005
Figure 6. Ambiguity Function (AF) contours showing delay-Doppler resolution at (a) iteration 1 and (b) iteration 20, which demonstrates smaller focal area in (b) as compared with (a).
Figure 6. Ambiguity Function (AF) contours showing delay-Doppler resolution at (a) iteration 1 and (b) iteration 20, which demonstrates smaller focal area in (b) as compared with (a).
Sensors 20 00415 g006
Figure 7. The symbol error rates (SERs) versus signal-to-noise ratio (SNR) for Binary Phase Shift Keying (BPSK), Quadrature Phase Shift Keying (QPSK), 16-Phase Shift Keying (16-PSK) and 256-Phase Shift Keying (256-PSK).
Figure 7. The symbol error rates (SERs) versus signal-to-noise ratio (SNR) for Binary Phase Shift Keying (BPSK), Quadrature Phase Shift Keying (QPSK), 16-Phase Shift Keying (16-PSK) and 256-Phase Shift Keying (256-PSK).
Sensors 20 00415 g007
Figure 8. Comparative throughput of BPSK, QPSK, 16-PSK, 256-PSK.
Figure 8. Comparative throughput of BPSK, QPSK, 16-PSK, 256-PSK.
Sensors 20 00415 g008
Figure 9. Target-scattering signal profiles at (a) iteration 1, (b) iteration 5, (c) iteration 10, and (d) iteration 20.
Figure 9. Target-scattering signal profiles at (a) iteration 1, (b) iteration 5, (c) iteration 10, and (d) iteration 20.
Sensors 20 00415 g009
Figure 10. The SERs versus SNR for BPSK.
Figure 10. The SERs versus SNR for BPSK.
Sensors 20 00415 g010
Figure 11. Probability of detection for adaptive waveform and static waveform.
Figure 11. Probability of detection for adaptive waveform and static waveform.
Sensors 20 00415 g011

Share and Cite

MDPI and ACS Style

Yao, Y.; Li, X.; Wu, L. Cognitive Frequency-Hopping Waveform Design for Dual-Function MIMO Radar-Communications System. Sensors 2020, 20, 415. https://doi.org/10.3390/s20020415

AMA Style

Yao Y, Li X, Wu L. Cognitive Frequency-Hopping Waveform Design for Dual-Function MIMO Radar-Communications System. Sensors. 2020; 20(2):415. https://doi.org/10.3390/s20020415

Chicago/Turabian Style

Yao, Yu, Xuan Li, and Lenan Wu. 2020. "Cognitive Frequency-Hopping Waveform Design for Dual-Function MIMO Radar-Communications System" Sensors 20, no. 2: 415. https://doi.org/10.3390/s20020415

APA Style

Yao, Y., Li, X., & Wu, L. (2020). Cognitive Frequency-Hopping Waveform Design for Dual-Function MIMO Radar-Communications System. Sensors, 20(2), 415. https://doi.org/10.3390/s20020415

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