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; to denote the transpose operation; to denotes the Hermitian operation; to denote the Kronecker product; to denote the identity matrix of size ; to denote the angle of a complex number.
3. Multiple-Input Multiple-Output (MIMO) Radar Ambiguity Function
In this section, we consider a radar target at
where
describes the delay of the target range,
denotes the Doppler frequency,
indicates the spatial frequency. Here
and
are the angle of the target and the wavelength, respectively.
and
are the distance between the transmit antennas and between the receiver antennas, respectively. To simplify the discussion, we assume
in this paper.
denotes the waveform radiated by the
m-th transmit antenna. In [
17], the MIMO radar ambiguity function can be defined as:
where
In (13),
describes the cross ambiguity function, which implicates two radar waveforms
and
.
is analogous to the SIMO ambiguity function presented in [
19]. We now discuss the MIMO radar ambiguity function for the case when
is composed of the shifted forms of a rectangular pulse
.
where
denotes the number of the rectangular pulse and
describes pulse repetition interval (PRI). The cross ambiguity function of the waveform
is defined as:
where
describes the cross ambiguity function of the rectangular pulses
and
. We assume that the pulse duration
and the Doppler frequency
are small enough so that
. At a result, the envelope of the Doppler frequency keeps unchanged within an entire pulse period. The cross ambiguity function
reduces to the cross correlation function
, which is not a function of
any more, we obtain:
We further assume that little reflections take place at these second trip ranges. Therefore, the cross ambiguity function of
can be rewritten as:
The MIMO ambiguity function of
can be expressed as:
It is worth noting that the MIMO ambiguity function
depends on the cross correlation functions
. Furthermore, the shifted forms of a rectangular pulse
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
is denoted as:
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
in terms of
and
. To derive
, 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:
and the function
can be expressed as:
where
describes the ambiguity function of
, which can be denoted by:
For , the function is the special case of the SIMO radar. For the general case of the MIMO radar , 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 .
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
can be expressed as a matrix
after discrete sampling, where
describes the number domain and
indicates the sample number.
is zero-mean white noise matrix. Therefore, the
matric of the target scattering signals is written as:
Here
and
, and
and
indicate the covariance matrices of the channel response
and the zero-mean white noise
, 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:
In (24),
represents the MI between two random variates
and
given the transmit matrix
, and
indicates the conditional entropy that
conveys about
. The main objective of this step is to maximize
between
and
given
. According to the definition of entropy, we have:
Here
indicates the conditional probability density function (PDF) of the received matrix
given transmit matrix
. The conditional PDF
can be expressed as follows:
Solving (25) and (26) gives rise to the following result for the conditional entropy [
40]:
Similarly, the result for the entropy of the noise can be derived as follows:
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:
Therefore, we can formulate the MI maximization criterion as follows:
where
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
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
and at time
as:
In (31), the term
(or
) denotes the entropy of
(or
) at time
(or
) given the knowledge of
(or
). The term
is defined similarly. According to the literature [
6,
34], the above Equation (31) can be rewritten as follows:
where
is the diagonal element of the diagonal matrix
.
is achieved by singular value decomposition (SVD) of the covariance matrix, which can be denoted as:
In (33),
is the cross-covariance of the expressions for
and
. Hence, the MI minimization criterion can be expressed as:
We can estimate the values for over all possible values of using (23). Thus, we can also form an estimate of all the values of the corresponding and choose the value for 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
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
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 and the covariance matrix of the noise . 1. At the initial time , measure the target scattering signal and calculate the covariance of the target-scattering signal . The covariance matrix of the channel response can be obtained through successive measurements with uniform power allocation over the transmit antennas by solving (36). 2. The optimized waveform ensemble can be obtained through the waveform design process by solving (30). 3. At time , measure and calculate and the cross-covariance matrix . The actual channel response is obtained by de-convolving the target-scattering signal with the transmit signal by using (23). 4. Calculate the corresponding singular values of and can be acquired through the waveform selection process by solving (38). 5. At time , emission and observe the corresponding received signal to achieve by using (36). 6. If , 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
GHz and bandwidth
MHz. The sampling frequency is
sample/sec, which is taken as the Nyquist rate. The PRI is
. We assume an arbitrary linear transmit array consisting of
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
, the length of the FH code is
and the FH interval duration is
. We generate a set of 16 FH pulse waveforms. The parameter
is used. Therefore, the 320 FH code is generated randomly from the set
, where
.
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
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)
MIMO system based on conventional maximum a posterior (MAP) approach; (2)
MIMO system using the Kalman filtering [
42]; (3)
MIMO system based on MI optimization (step 2) scheme; (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 , 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
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
= 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
, and the phase of
is uniformly distributed within the interval
. To test the SER performance, we generate a number of
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.