Abstract
This study deals with the problem of radar waveform design based on the weighted mutual information (MI) and the difference of two mutual information metrics (DMI) in signal-dependent interference. Since the target and clutter information are included in the received signal at the beginning of the design, DMI-based waveform is designed according to the following criterion: maximizing the MI between the received signal and target impulse response while minimizing the MI between the received signal and the clutter impulse response. This criterion is equivalent to maximizing the difference between the first MI and the second MI. Then maximizing the difference of two types of MI is used as the objective function, and the optimization model with the transmitted waveform energy constraint is established. In order to solve it, we resort to maximum marginal allocation (MMA) method to find the DMI-based waveform. Since DMI-based waveform does not allocate energy to the frequency band where the clutter power spectral density (PSD) is greater than the target PSD, we propose to weight the MI-based waveform and DMI-based waveform to synthesize the final optimal waveform. It could provide different trade-offs between two types of MI. Simulation results show the proposed algorithm is valid.
1. Introduction
Radar performance is not only dependent on the signal processing at receiving end, but also closely related to the transmitted waveform [1]. The transmitter of conventional radar transmits the same fixed waveforms every time regardless of environment factors. Different from the conventional radar, cognitive radar (CR), which was first proposed by Simon Haykin in 2006, has a closed loop feedback structure from the receiver to the transmitter [2]. It could improve the radar performance by sensing the environment and transmitting adaptive waveforms. Waveform optimization, as one of the important technologies of CR, has attracted attention by lots of scholars in the past years.
There are lots of metrics for adaptive waveform design, such as maximum signal-to-noise ratio (SNR) or maximum signal-to-interference-plus-noise ratio (SINR), mean square error, and information theory [3–24]. Pillai proposed to maximize the output SINR as criterion to optimize the waveform by jointly transmitted waveform and the receiver [4]. In [5], the authors used maximum SINR as the objective function to design waveform to improve target detection performance in the case of unknown target impulse response (TIR). The maximum SINR metric is also applied to the robust waveform design, such as [6, 7]. Some scholars extended maximum SNR/SINR criterion to MIMO radar waveform design with different signal models or different constraints [8–11].
Information theory is one of the main methods. In [3], Bell first proposed to maximize the received signal and target impulse response (TIR) to optimize the transmitted waveform. In [12], the authors used a linear combination of the MI between target signals and the corresponding received beams to optimize transmitted waveform for multiple targets. Goodman et al. conducted in-depth research on the optimization of transmission waveform by using the MI method [13–19]. Their studies mainly aimed to design optimal waveform at the target identification and improved the waveform optimization method based on maximizing SNR criterion and maximizing MI criterion. In [13], the authors used sequential probability ratio test to weight each target PSD to optimize the waveform at each transmission. In [17], the authors summarized SNR-based waveform design and MI-based waveform design for known target and stochastic target, respectively; they also derived the relationship of MI and SNR for stochastic signals. In [18], the authors derived an enhanced waveform design algorithm, which derived a strict lower bound of MI based on energy spectral variance and optimized an enhanced waveform by maximizing the lower bound. In [19], the authors proposed a waveform design algorithm based on maximizing MI between the measurements and the target signals for Gaussian mixture target range profiling. Yang and Blum [20] extended the maximizing MI criterion to MIMO radar waveform design for target identification and classification. They also proposed the waveform optimization based on minimum mean square error (MMSE) criterion and pointed out that optimal waveforms based on MI criterion and MMSE criterion have the same solution under the same energy constraint. In [21], the authors aimed to maximize the MI between target echoes and TIR to design waveform for the MIMO radar with STAP to detect target both considering the range cell including target and influence of neighboring range cells with different spatial distributions. Maximizing MI criterion has also been extended to robust waveform design. In the actual radar scene, the target prior information is uncertain, which leads to the waveform designed based on the target prior information and cannot satisfy the requirements of parameter estimation; thus, Wang et al. proposed a robust waveform design method for parameter estimation. They also maximized MI between the radar target echo and the random target spectrum response under the hierarchical game model of radar and jammer [22]. Due to uncorrelation of the successive target echoes in some actual scene, a novel waveform design method is presented, which minimized the MI between successive target echoes as criterion to optimize the transmitted waveform [23]. In [24], the authors minimized the MSE of the target scattering coefficients estimation and then designed the transmitted waveform by using minimization of MI between the radar successive target echoes. Tang et al. proposed an algorithm to design MIMO radar waveforms via relative entropy in the clutter; the proposed algorithm is devised based on minorization-maximization approach, which could reduce complexity. In [25], the authors recommended to maximize Kullback–Leibler divergence locally most powerful (KLD-LMP) metric to optimize the transmitted waveform and indicated the relationship between KLD, MI, and SNR. In [26], the authors developed a close-loop adaptive OFDM waveform design approach based on maximizing MI criterion in the presence of signal-dependent clutter. They updated the complex weights of the transmitted OFDM signal to enhance target detection. In [27], they illustrated the relationship of SINR, MI, and probability of detection criterion. With the same constraints, three criteria could optimize the similar waveform. Numerical results part of [27] indicates the mutual information may not be the largest when the SINR reaches the maximum, which is consistent with the simulation results in our paper.
According to the above analysis, the waveform design methods based on mutual information are mainly to maximize the target information in the received signal or minimize the uncorrelated parameters as the optimization criterion. In this paper, we proposed a new CR waveform design scheme, which combine DMI-based waveform algorithm and MI-based waveform algorithm to obtain an optimization waveform; the proposed scheme is divided into two steps as follows: Step 1: The step is to design transmitted waveform based on DMI criterion with the transmitted energy constraint. The main objective is to maximize MI between the received signal and TIR and minimize MI between the received signal and clutter impulse response (CIR). Step 2: The final optimal waveform is obtained by weighting the DMI-based waveform and MI-based waveform. Let the sum of the weights of DMI-based waveform and MI-based waveform be equal to 1, so as to ensure the transmitted energy constraint. The main innovations of our work can be described as follows:(1)We first present a waveform design method based on DMI criterion, which ensures that the receive signals could contain more target information and less clutter information compared with conventional radar.(2)Although the DMI-based waveform can improve the radar performance, it has a disadvantage. DMI-based waveform will not allocate energy to the frequency band where the clutter PSD is greater than the target PSD, which results in received signals containing little target information within these bands. In order solve it, we propose to weigh the DMI-based waveform and MI-based waveform to obtain the final optimized waveform. By adjusting the weight, the final optimized waveform can allocate the main energy to all frequency bands with strong target PSD. Even if the clutter PSD is greater than the target PSD, more target information can be obtained from the received signal while the clutter information can be reduced. Compared with MI-based waveform, the final optimized waveform could reduce more clutter information when the target information is little loss.
The structure of this paper is as follows. The signal model is given in Section 2. In Section 3, we derive the waveform design based on DMI criterion and use MMA to solve the optimization model for waveform. In Section 4, we present an improved method to solve the problem of DIM-based waveform. In Section 5, we present the simulation results of the proposed method and discuss them. Finally, we conclude the paper in Section 6.
2. Signal Model
The target channel model in clutter environment is shown in Figure 1. The transmitted waveform has the energy within time interval . denotes the target impulse response, which is a Gaussian random process with zero mean and variance . is additive Gaussian noise with zero mean and one-sided PSD . Denotes the clutter impulse response, which is still a Gaussian random process with zero mean and variance . , , and are statistically independent. The received signal could be written bywhere , , and the symbol denote convolution.

3. DMI-Based Waveform
3.1. Waveform Design Based on DMI Criterion
The classical waveform optimization method based on MI criterion is to maximize the mutual information between the received signal and TIR by optimizing the transmitted waveform, that is, , where denotes mutual information.
In (1), the received signal contains two components: the target echo and the interference . The radar system hopes that the received signal contains more target information and less clutter information. Hence, the MI between the received signal and TIR, , should be maximum. The MI between the received signal and CIR, , should be minimum. The optimized waveform based on DMI criterion should satisfy the following:
Expression (2) could be equivalent towhere is the Fourier transform of , is the transmitted waveform bandwidth, and is the duration of the received signal. The derivation of (3) is shown in Appendix.
Equation (3) is the optimization criterion based on DMI. Considering the total energy constraint of transmitter, the optimization model for DMI-based waveform is
If the clutter does not exist, the optimization model (4) is consistent with the optimization model in literature [3].
3.2. MMA for DMI-based Waveform Solution
The solution of optimization model (4) is solved by MMA method. First, we discretize model (4):
Let
Substituting (6)–(8) into (5), we have
The waveform optimization problem can be transformed to seek the maximum of with the constraint , whereand could be divided equally into parts, that is, , and take values in the set for all , and is defined as the minimum energy unit. The energy unit is allocated at each step of MMA until all energy units are allocated. If , then . Repeat the same procedure, and only in the procedure we need to select the value of to make (11) reach the maximum or the to make the maximum marginal increase:
For example, let , , , and . We wish to maximize
Either 0, 1, 2, 3, or 4 units could be allocated to , , and with the energy constraint. The possible values of corresponding to the values of for are shown in Table 1. When the energies are initially allocated units, the values of are 0.2877, 0.2231, and 0.1823, respectively. The maximum value 0.2877 is selected. Thus, units are allocated to in the first step 1.
After units are allocated to , the new marginal energies are obtained at as shown in Table 2. The marginal energies are 0.0488, 0.2231, and 0.1823. The maximum value 0.2231 is selected. Thus, units are allocated to in the second step.
Similarly, the new energies are obtained at in the third step as shown in Table 3. The marginal energies are 0.0488, 0.0646, and 0.1823. the maximum value 0.1823 is selected. Thus, units are allocated to . Table 4 shows the last energies allocation. The maximum value 0.0646 is selected. Table 5 shows the final energies allocation. units are allocated to , units are allocated to , and units are allocated to . More details on MMA method could refer to literature [28].
4. Weighted MI- and DMI-Based Waveform
4.1. Problem Formulation
In the previous section, DMI-based waveform could allocate the energy reasonably according to the PSDs of target and clutter. However, DMI-based waveform has two disadvantages.(1)DMI-based waveform only allocates the energy to the frequency bands where target PSD is greater than that of clutter (see Figure 2 in Section 4).(2)There may be no solution to DMI-based waveform in the frequency domain.
Proof. The optimization model for DMI-based waveform is (4). If using Lagrange multiplier to model (4), we have Equation (13) can be equivalently maximized asThe first derivatives of with respect to areThe second derivatives of with respect to arewhereThe denominator of (16) is greater than 0, and . If DMI-based waveform has solution for optimization model (4), (16) should be less than 0. Thus, must be greater than 0. That is, the reason for DMI-based waveform allocates the energy to the frequency band where target PSD is greater than that of the clutter PSD; see Figure 2 in Section 5. If , the DMI-based waveform has no solution in all frequency bands.

4.2. Final Scheme
In order to solve the above disadvantages of DMI-based waveform, the final scheme for waveform design based on weighted MI and DMI criterion is proposed.
The optimization model for waveform based on MI criterion with energy constraint is [29].
Assuming that is the optimal solution to model (18) (it can be solved by Lagrange multiplier, MMA, or other methods).
Assuming that is the optimal solution to model (4) (if ), then, the final optimized waveform is obtained by weighting MI-based waveform and DMI-based waveform ; that is,where and are weights and
Equation (20) could ensure that final optimal waveform still has the total energy .
When , the optimal waveform has no solution. In this situation, let and . The final optimal waveform is MI-based waveform .
The final optimization scheme for waveform design is as follows: Step 1: The radar system estimates the PSDs of the target and clutter. Step 2: MI-based waveform e is derived according to model (18). Step 3: Determine whether is greater than 0, if , DMI-based waveform is derived according to model (28), and Step 4 is executed. Otherwise, step 5 is executed. Step 4: According to (19) and (20), the final optimal waveform is derived by weighting and . Step 5: If , the final optimal waveform is MI-based waveform .
The algorithm flow is shown in Figure 3.

5. Simulation and Results
In this section, we demonstrate the effectiveness of the proposed method. The signal energy is (energy unit). The noise PSD is . The clutter-to-noise ratio (CNR) is and the target-to-noise ratio (TNR) is . is set to 1000. Figure 4 shows the spectrum signatures of the target and clutter.

Figure 2 shows the energy spectral densities (ESD) of MI-based waveform and DMI-based waveform, respectively. The main characteristics of DMI-based waveform are as follows:(1)DMI-based waveform could allocate the most energy to the frequency bands with the stronger target PSD. Although the target PSD is symmetric about the 0 frequency point, the clutter PSD is stronger in the frequency band than that in the frequency band ; the DMI-based waveform allocates more energy to the sub two frequency bands with the stronger target PSD in the range of frequency band .(2)Since DMI-based waveform takes into account the minimum mutual information between the received signal and clutter, the energy of DMI-based waveform is only allocated to the frequency bands whose target PSD is stronger than the clutter PSD. When the clutter PSD is greater than target PSD, DMI-based waveform will not allocate energy even if target PSD is still strong.
The performance of DMI-based waveform is verified as shown in Figure 5, and the linear frequency modulated (LFM) waveform is used as a benchmark. Figure 5 shows detection probability curves corresponding to the DMI-based waveform, MI-based waveform, and LFM waveform, respectively. The performance of DMI-based waveform and MI-based waveform are better than that of LFM waveform. However, the detection probability of DMI-based waveform is slightly lower than that of MI-based waveform. This is due to the fact that, in order to minimize the mutual information between the received signal and CIR, DMI-based waveform does not allocate energy to the frequency band, wherein the clutter PSD is greater than that of target PSD. To solve this problem, weighted MI- and DMI-based waveforms is proposed.

MMA algorithm is applied to solve DMI-based waveform. Since , the larger is, the smaller is. Figure 6 shows the DMI-based waveform based on MMA algorithm with different values. is set to 10, 100, and 1000, respectively. If is growing larger, will become smaller, and the optimal waveform will be more accurate. Figure 7 shows the performance of different DMI-based waveforms corresponding to different values. The performance of DMI-based waveform with is the worst. The DMI-based waveform with performs the best. It demonstrates that is great significance in MMA algorithm; however, when increases to a certain value, the performance improvement is not obvious.


Figure 8 shows weighted MI- and DMI-based waveforms corresponding to , , and , respectively. The main characteristics of weighted MI- and DMI-based waveforms are as follows:(1)The ESD of weighted MI- and DMI-based waveforms is between ESDs of MI-based waveform and DMI-based waveform at each frequency point. Weighted MI- and DMI-based waveforms allocate the energy to all frequency bands with strong PSD of target.(2)When both target PSD and clutter PSD are strong and clutter PSD is greater than target PSD, compared with MI-based waveform and DMI-based waveform, weighted MI- and DMI-based waveforms will allocate less energy in these frequency bands rather than no energy, which can not only ensure the acquisition of target information but also reduce clutter information.

The performance of the weighted MI- and DMI-based waveforms is shown in Figure 9; we plot the curves of detection probabilities of several weighted MI- and DMI-based waveforms corresponding to , , and , respectively. As can be seen from Figure 9, all weighted MI- and DMI-based waveforms with different weights perform better than DMI-based waveform. The performance of the weighted MI- and DMI-based waveforms with different weights is similar to that of the MI-based waveform in our simulation scheme. Depending on different weights selected, some weighted MI- and DMI-based waveforms have better performance than MI-based waveform, while some are slightly inferior to MI-based waveform such that when and , the performance of the weighted MI- and DMI-based waveforms is very approximate to that of MI-based waveform, and when and , the performance of the weighted MI- and DMI-based waveforms is better than that of MI-based waveform, which indicates that different weights affect the performance of optimal waveforms. How to select optimal weights to achieve the best performance of optimal waveform is our future research.

Table 6 shows the SINR, MI between received signal and TIR and MI between received signal and CIR corresponding to different optimal waveforms in our simulation scheme. DMI-based waveform could obtain the minimum value of MI between the clutter and received signal, while the SINR is also the worst, since, at the same time of reducing the clutter information, the target information also lost the most. Thus, the DMI-based waveform has the worst performance. Weighted waveform (, ) has the best SINR.
6. Conclusions
In this paper, we have developed radar waveform design based on weighted MI and DMI criterion in the presence of clutter, which takes into account the fact that the received signal contains more target information and as little clutter information as possible. Firstly, the DMI criterion, the maximizing MI between the received signal and TIR and the minimizing MI between the received signal and CIR, is used to establish optimization model with transmitted energy constraint, which is solved by MMA. Since DMI-based waveform is solved only when the target PSD is stronger than the clutter PSD in the same frequency band, we propose weighted MI and DMI-based waveform; it is a trade-off between maximum target information and minimum clutter information. If the clutter PSD is stronger than the target PSD in the entire frequency band, the weighted MI and DMI-based waveform is equal to the MI-based waveform. The performance of weighted MI and DMI-based waveform is dependent on selection of weights. How to optimize the weights adaptively is our future research.
Appendix
Proof for Equation (3):
Considering a random variable with probability density function (PDF) , the differential entropy is defined by
Assume there exists another random variable . and are jointly distributed. The conditional differential entropy is defined bywhere is the PDF of conditioned on . The mutual information between two jointly distributed random variables and is defined by [3].
Assuming the random variable is Gaussian variable with zero mean and variance , its PDF is
Substituting (A.4) into (A.1), the differential entropy of is
Equation (A.5) demonstrates the entropy of the Gaussian variable with zero mean is dependent on the variance of the Gaussian random variable.
Figure 10 shows the channel of additive Gaussian noise and signal-dependent clutter. is the Gaussian random variable with zero mean and variance and is the signal-dependent clutter, which is Gaussian random process with zero mean and variance . is the additive Gaussian noise with zero mean and variance . The output is the sum of , , and .

The mutual information between and is
Assuming that , , and are statistically independent, hence is still Gaussian random variable and its variance is
According to (A.5), the entropy is
According to literature [3], the conditional differential entropy is
Substituting (A.8) and (A.9) into (A.6), the authors have
Discussing Figure 1 again, the signals , , , and with frequency components confined to the frequency interval , according to sampling theorem; each signal can be replaced by a series of samples obtained from uniform sampling rate of . Assuming that the spectra , , , and are smooth and have a constant value for all , the samples of the Gaussian process sampled are statistically independent.
The samples are independent and identically distributed random variables with zero mean and variance . Over the time interval , the energy in is
The energy is evenly spread among statistically independent samples. Hence, the variance of each sample, , is
Similarly, the energy in is
The energy is evenly spread among statistically independent samples. Hence, the variance of each sample, , is
The energy of noise is
The energy is evenly spread among statistically independent samples. Hence, the variance of each sample, , is
Substituting (A.13), (A.14), and (A.16) into (A.10), for each sample of and corresponding sample of , the mutual information between and is
Similarly, for each sample of and corresponding sample of , the mutual information between and is
In the time interval , and contain statistically independent samples. Thus,
The bandwidth is divided into many continuous disjoint intervals, and the interval bandwidth is . According to literature [3], when , the mutual information between and can be obtained as
The mutual information between and can be obtained as
Substituting (A.20) and (A.21) into , the authors have
Data Availability
All data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported by the Fundamental Research Funds for the Central Universities (no. N182304022), the National Natural Science Foundation of China (nos. 61601109 and 61473066), and the Natural Science Foundation of Hebei Province (no. F2018501051).