Information Content Based Optimal Radar Waveform Design: LPI’s Purpose

: This paper presents a low probability of interception (LPI) radar waveform design method with a ﬁxed average power constraint based on information theory. The Kullback–Leibler divergence (KLD) between the intercept signal and background noise is presented as a practical metric to evaluate the performance of the adversary intercept receiver in this paper. Through combining it with the radar performance metric, that is, the mutual information (MI), a multi-objective optimization model of LPI waveform design is developed. It is a trade-off between the performance of radar and enemy intercept receiver. After being transformed into a single-objective optimization problem, it can be solved by using an interior point method and a sequential quadratic programming (SQP) method. Simulation results verify the correctness and effectiveness of the proposed LPI radar waveform design method.


Introduction
In the past, short duration pulse train waveforms are widely equipped with most military radars.These pulses have a high ratio of peak to average power without any frequency or phase modulations, and, therefore, can be easily detected by adversary intercept systems, such as electronic support measures (ESMs), electronic attack (EA) systems and radar warning receivers (RWRs), in modern battlefields [1][2][3].For surviving these enemy attacks, complex radar waveforms are designed to reduce the detection performance of passive intercept receiver devices, which can eventually achieve that the radar detection range is greater than the enemy interception range.These radar waveforms are called low probability of interception (LPI) waveforms, which have been developed for several decades [1][2][3][4][5][6][7].
They employ wideband modulation techniques to spread their energy in frequency domain, so that the signal strength in each frequency bin can be effectively reduced.Available wideband modulation techniques include frequency modulation, phase modulation, pseudo-noise modulation and so on.They generate the most common LPI radar waveforms, such as linear frequency modulation (LFM) waveform, Frank code, P1-P4 codes, Costas code and frequency shift keying/phase shift keying (FSK/PSK) signal [7].However, the design of these waveforms does not consider optimal or near-optimal radar performance.Actually, LPI waveform design is a trade-off between radar task performance (such as detection performance, identification performance, location performance, etc.) and LPI performance.One extreme is a waveform with optimal radar performance which completely ignores its LPI performance, while the other extreme gives full consideration of the worst detection performance of enemy intercept receiver without regard to radar performance.The ideal LPI waveform would be to both achieve high radar performance while maintaining its LPI performance.Due to the fact that radar can obtain an additional matched filtering gain with prior information about its transmitted waveform, while the adversary intercept receiver can only capture fewer processing gains by using time-frequency analyses as it has no information about the received signals, the ideal LPI radar waveform can be achievable in practice by making a small reduction of radar performance but a sharply shrinkage of interception performance of adversary intercept receiver.
In this paper, we hope to model and solve the trade-off.To achieve it, the first question is the establishment of separate optimization metrics for radar and the intercept receiver.Many metrics about optimal radar performance have been proposed, such as mutual information (MI), output signal-to-noise ratio (SNR), mean square error and so on [8][9][10][11][12][13][14].Therein, the maximization of MI has been widely employed in radar waveform design [8][9][10].However, the performance metrics of the intercept receiver for radar waveform design are rarely studied.The existing LPI waveform ranking metrics include time-bandwidth product, peak to average power ratio [1,2,15,16] and their integrated function [17].In 2015 and 2016, we presented simpler and more effective LPI ranking metrics [18,19] from the perspective of information distance via ignoring the concrete form of radar waveforms.Based on these metrics, in this paper, we treat the Kullback-Leibler divergence (KLD) between the intercept signal and background noise as a practical performance metric of intercept receiver for LPI radar waveform design.Therefore, based on maximizing MI of radar and minimizing KLD of intercept receiver, an LPI waveform design method is presented and solved under the constraint of a fixed average transmitted power.The simulation results verify its correctness and effectiveness.
The rest of this paper is organized as follows.In Section 2, the MI of radar and KLD of intercept receiver are analyzed, respectively.In Section 3, the optimization model of LPI radar waveform design is developed.In addition, its solution is introduced via using the interior point method and the sequential quadratic programming (SQP) method in combination.In Section 4, the performance of the proposed LPI radar waveform design method is evaluated.Conclusions are drawn in Section 5.

Mutual Information in Radar System
We begin with a brief overview of Bell's MI in the radar system [9], which, in this paper, is used as the radar performance metric in LPI radar waveform design.The task of radar is to obtain information about targets, which will be probed by the transmitted waveforms.It can be depicted by Figure 1.Radar waveform s (t) is emitted by the transmitting antenna with gain G t in the target direction.Then, it can be intercepted and scattered by the target with a target impulse response h (t).The scattered waveform x (t) can be expressed as where, α 1 denotes energy attenuation coefficient, which is a function of emitter wavelength λ, antenna gain G t , propagation distance R, path loss L 1 and target impulse response h (t), which is assumed as a zero-mean Gaussian random process, and the symbol " * " denotes convolution.After propagating in the atmosphere, the scattered waveform can be captured by the radar receiving antenna with gain G r .Due to the existence of additive noise in the propagation path and radar receiver, the received signal y (t) is a noisy signal.It can be expressed as where n (t) is a zero-mean additive white Gaussian noise.Similar to the definition of α 1 , α 2 is also an energy attenuation coefficient related to receiving antenna gain G r , propagation distance R and path loss L 2 , and For notation convenience, S, X, Y and N are denoted as the discrete transmitted waveform, scattered waveform, received waveform and additive white Gaussian noise, respectively.The more amount of information about the target impulse response H the scattered waveform X contains, the more information about the target the radar system will obtain to better detect, locate, identify and track the target.Hence, the waveform optimization problem can be described as According to Equation (1), the scattered waveform X is a Gaussian random process with zero mean and variance σ 2 X , which is independent with zero-mean additive white Gaussian noise N with variance σ 2 N .Hence, the received waveform Y is also a Gaussian random process with zero mean and variance H and N are also assumed to be statistically independent.The optimization problem in Equation ( 3) can be converted to maximize the mutual information I (H; Y|S) between the random target impulse response H and the received waveform Y. Since Bell has proved that the waveforms S that maximize I (X; Y|S) also maximize I (H; Y|S), the waveform optimization problem can be recast as [9] max As described in [9], the transmitted waveform S is considered as a finite energy deterministic waveform with energy E and duration T, and the bandwidth of S is subject to the operation bandwidth W of the radar system.For the convenience of analysis, the frequency interval W is split into a large number of sufficiently small and disjoint sub-bands, whose bandwidth is ∆ f , and S k , X k , Y k , N k denote the components of S, X, Y, N with frequencies in the k-th sub-band.
By transforming the x (t) in Equation (1) from time domain to frequency domain, we have , where X , S and H are frequency domain representations of X, S and H, and in each sub-band can be calculated as where P N ( f ) is the one-sided power spectral density of Gaussian noise N of radar.Let ∆ f → 0, the mutual information I (X; Y|S) is an integral of I (X k ; Y k |S) over the frequency band W,

Kullback-Leibler Divergence in the Intercept Receiver
Suppose the target is equipped with an intercept receiver.When the radar signal is used to detect the target, it can also be intercepted by intercept receiver, which makes the radar face serious threats, such as anti-radiation missiles.Therefore, the transmitted waveform S should be designed to reduce the detection and identification probabilities of the intercept receiver.It can be depicted in Figure 2.After being radiated by the radar emitting antenna and propagated in the atmosphere, the radar waveform can be intercepted by the intercept receiver.The intercept signal can be written as where wherein G I denotes the intercept receiver antenna gain.After time-frequency characteristics are analyzed, the intercept signal is detected and extracted features to recognize its modulation type for achieving the emitter location and recognizing the emitter.Since it is a common view that white noise is the best LPI waveform, the closer the distance between intercept signal and noise signal is, the more difficult it is to detect and recognize the intercept signal.It is logical to compare random data sets by ways of the associated probability density function (PDF).The waveform optimization problem against intercept receiver can be briefly written as where f Z and f N are the PDFs of intercept signal Z and noise signal N, and D denotes a distance measure.There are many PDF distance measures to quantify the disparity between PDFs [20], such as L 1 norm, Kullback-Leibler divergence (KLD), Bhattacharyya distance, f -divergence, Hellinger distance, etc.Among these distance measures, KLD has been confirmed to be a powerful and accurate tool to measure the information of multivariate data with lower complexity and superior performance [21].Thus, we use KLD here as the PDF distance measure, which is denoted as D KL (Z; N).In the frequency domain, the intercept signal z (t) in Equation ( 7) can be transformed as Z ( f ) = γ • S ( f ) + N ( f ), where Z, S and N are frequency domain representations of signals Z, S and N. Like the computation of Bell's MI, we also split the frequency interval W into a large number of sufficiently small and disjoint frequency intervals F k whose bandwidth is ∆ f , and for all f ∈ F k , we have # ( f k ) ≈ # ( f k ), where # can be the waveforms Z, S and N .In addition, Z k denotes the component of Z with frequencies in F k , which is the same for S k and N k .
In each frequency interval F k , we suppose the sample frequency is f s = 2∆ f according to the sampling theory, and, therefore, the sample sizes of S k , Z k , and N k are N = 2∆ f T. Since the distribution of background noise N of the intercept receiver at the frequency f k ∈ F k is a normal distribution with zero mean and variance σ 2 N k , which can be denoted as N ( f k ) ∼ Normal 0, σ 2 N k , the distribution of the intercept signal Z at the frequency f k can be achieved as where P N ( f ) is the one-sided power spectral density of Gaussian noise N of intercept receiver.We can calculate the KLD between Z ( f k ) and N ( f k ), The KLD in Equation ( 11) is a function of signal-to-noise ratio (SNR) P N ( f k ) at the frequency f k .The lower SNR a waveform possesses in each frequency, the harder it is for the intercept receiver to detect it.It meets our common knowledge and experience.The waveform optimization problem against intercept receiver then can be described as the minimization of the KLD In order to facilitate the solving of the optimization problem, the natural exponential function of the KLD D KL (Z ( as the new objective function, which can be denoted as Since the sample size in each frequency interval

Optimization of Radar Waveform for the Purpose of LPI
According to the analysis of Section 2, the mutual information I (X; Y|S) and the KLD D KL (Z; N|S) are tied to the radar and the intercept receiver performance, respectively.We need a transmitted waveform S should not only have good performance of detecting, locating, identifying and tracking target, but also have a good LPI performance against the detection and identification of intercept receiver.Here, the process of LPI radar waveform design can be described as the following optimization problem: where W |S ( f ) | 2 d f is the average transmitted power and P S is fixed.
To maximize the radar performance is just one part of the job in radar waveform design, and the minimization of the performance of the intercept receiver should also be taken into account, as shown in Equation ( 14), which is a multi-objective optimization problem with respect to |S ( f ) | 2 under a fixed average power constraint.However, the primary task of the emitting waveform is to fulfill the target detection.The designed waveform should be considered the LPI capability under the condition of meeting radar performance.Thus, the multi-objective optimization problem in Equation ( 14) can be simplified to a single-objective optimization problem, that is, where I th is the needed mutual information to meet the radar performance, which can be set at various values in different battling scenarios.
The maximization of the mutual information of radar and the minimization of the KLD of intercept receiver are just like lance and shield.To solve the problem, we minimize the KLD in the situation that the mutual information makes some concessions, as shown in Equation (15).Before solving the optimization problem, we first digress to discuss the solutions of two independent optimization problems, which take the mutual information and the KLD as objective functions, respectively.The two optimization problems are expressed as (OP2)

Solution of the OP1
For OP1, Bell has given the concrete solving process [9].For comparing with the solving process of OP2 and determing the constraint value I th easily, we restate the process briefly.
The Lagrange multiplier technique is utilized here, and the objective function can be yielded as Taking the partial derivative of ψ |S ( f ) | 2 , λ 1 with respect to |S ( f ) | 2 and setting it to zero, we have Since the partial derivative of 19) and the constraint W |S ( f ) | 2 d f = P S , which maximizes the mutual information I (X; Y|S), can be written as where , and the value of A can be solved numerically by a search algorithm to satisfy the constraint W max 0, The maximum value of the mutual information I (X; Y|S) can be achieved as Thus, the constraint value I th of MI I (X; Y|S) should satisfy I th ≤ I max (X; Y|S).

Solution of the OP2
Like the solving method of OP1, we reuse the Lagrange multiplier technique, and the objective function can be obtained as This is seen to be equivalent to minimizing Taking the partial derivative of φ |S ( f ) | 2 , λ 2 with respect to |S ( f ) | 2 and setting it to zero, we have Solving for |S ( f ) | 2 , we have where Since the partial derivative of with respect to |S ( f ) | 2 is greater than zero, the solution |S ( f ) | 2 in Equation (25) can minimize the KLD D KL (Z; N|S), which can be written as When the background noise of intercept receiver is a white Gaussian noise, the power spectral density P N ( f ) is a constant.Under this assumption, the optimized waveform |S † ( f ) | 2 that minimizes From the Equation (28), we have a conclusion that it is most difficult to intercept the transmitted waveform whose power is distributed equally in the whole bandwidth W for an intercept receiver under the assumption of white Gaussian noise.
Since the radar has no information about the background noise types of unknown intercept receivers, the assumption of white Gaussian noise is the most reasonable choice.Moreover, the background noise in most of the receivers is white Gaussian noise.Here, it is unnecessary to study the solution of OP2 under the assumption of colored Gaussian noise.In this case, we only need to add another constraint |S ( f ) | 2 ≥ 0, and the interior point method can be used to solve the optimization problem.Here, the solving process about colored Gaussian noise will not be described in detail.The corresponding experiment results will be presented in Section 4.

Solution for the Optimization Problem of LPI Waveform Design
We now give the detailed solving process of the optimization problem in Equation (15).In order to achieve the optimal solution smoothly, the discrete form of the optimization problem first need to be obtained, which can be written as where and N = W b .The classical method for solving the nonlinear optimization problem in Equation ( 29) is the combination of interior point method and sequential quadratic programming (SQP) method [22].Through introducing slack variables c i , i = 1, • • • , 2N + 1, the inequality constraints in Equation ( 29) can be transformed into equality constraints, and correspondingly bounds c i ≥ 0, i = 1, • • • , 2N + 1 are yielded.The transformed optimization problem can be described as According to the interior point method, the barrier function can be constructed as By introducing a penalty factor τ ≥ 0, the optimization problem in Equation (30) can be rewritten as The Lagrangian of the optimization in Equation (32) is where Here, the SQP method can be applied to the optimization problem in Equation (32) and yields the quadratic program min ∇h (z) where Considering the possibility of indefiniteness of M, Richard H. Byrd uses a trust region method to solve such problems [22].The optimization problem can be recast as min ∇h (z) where • denotes a norm, and ∆ > 0 denotes a trust region radius, which will be updated in each iteration.
Here, we uses the algorithm proposed by Richard H. Byrd [22] to solve our optimization problem of Equation (29).The algorithmic framework of LPI radar waveform design is briefly introduced below (See Algorithm 1).Details of the nonlinear programming algorithm can be found in [22].

Simulation Results
In this section, several simulation experiments are provided to verify the effectiveness of the proposed LPI radar waveform design method, which considers not only the maximization of radar performance, but also the minimization of intercept receiver performance.The simulation parameters are listed in Table 1.The background noise in simulations includes white Gaussian noise and colored Gaussian noise.The one-side power spectral density P N ( f ) of white Gaussian noise is 1.9531 × 10 −18 .The one-side power spectral density P N ( f ) of colored Gaussian noise is a sinc function, which is set as P N ( f ) = 1.6531 × 10 −18 × sinc 1.75 × 10 −8 × f − 165.62 + 1.6531 × 10 −18 , for frequencies f ∈ [9.744, 10.252] GHz.For target information, we suppose the target headed for radar is F-16 aircraft, which is modeled in Figure 3a, and the variance σ 2 H ( f ) , f ∈ [9.744, 10.252] GHz of its impulse response with azimuth 0 • and elevation 5 • between the radar and the target has been calculated by electromagnetism simulation software, which is shown in Figure 3b.We also compare the proposed optimized radar waveforms with the common LPI radar waveforms: LFM, Frank and P1-P4 codes.

5
∂L(S,c,r E ,r I ) ∂r I .
6 Loop until: E (S, Set τ ← ατ, ε τ ← αε τ , S ← S k , c ← c k .17 End loop Their transmit power, bandwidth and duration are the same as those of optimized waveforms, and the code numbers of Frank, P1-P4 codes are 64. To test the proposed LPI waveform design method, we run four experiments, including both radar and intercept receiver, are under a white Gaussian noise background, only radar is under a colored Gaussian noise, only the intercept receiver is under a colored Gaussian noise, and both radar and intercept receiver are under a colored Gaussian noise.In each experiment, the SNR losses, which are defined in Equation (36), of radar and the intercepter receiver are simulated, respectively, as a metric of performance degradation.Through comparing the SNR losses of radar and intercepter receiver, the superiority of our proposed LPI waveform design method can be judged easily: where ∆SNR is the difference between the receiving SNR SNR MI of the optimized waveform just considering the maximization of radar performance and the receiving SNR SNR MI−KLD of the proposed optimized waveform.

Both Radar and Intercept Receiver are under White Gaussian Noise Background
Figure 4a shows the optimized waveforms with different mutual information constraints I th used our proposed MI-KLD method, and the optimized waveforms used the methods of maximization of radar performance (see OP1) and minimization of intercept receiver performance (see OP2).As a result of the comparison with Figure 3b, the optimized waveform (see the dotted purple line in Figure 4a) of OP1 has a large value of power spectral density at the frequencies where the variance of the target impulse response also has a larger value to get the maximum matching with the target, so that the radar can achieve the best performance.However, due to the restriction of a limited transmit power P S , the power spectral density of the optimized waveform is equal to zero at most of frequencies where the variance of the target impulse response has a relatively low value, which can be verified by Equation (20).For the solution (see the dotted red line in Figure 4a) of OP2, the power of the optimized waveform is equally distributed into all frequencies f ∈ W to decrease the intercept probability at each frequency, so that the intercept receiver can have a worst performance.
The aim of the proposed LPI waveform design method is to reduce the peak power under the condition of a certain loss of target matching quality.The reduction of peak power is to decrease the interception probability of intercept receivers, and the loss of the target matching quality means the degradation of radar performance.The blue line in Figure 4a shows that the optimized waveform, to a certain extent, not only achieves the target matching, but also realizes the reduction of peak power.With the shrinkage of the constraints I th of mutual information, the target matching quality and peak power of the optimized waveform keeps coming down, and the optimized waveform tends to the solution of OP2.The comparison of the performance degration between radar and intercept receiver with the reduction of I th is shown in Figure 4b.Since the radar has the prior information about its transmitted waveform, it can increase the receiving SNR through the matched filtering technique.Thus, a certain loss of target matching quality may have little impact on the receiving SNR of radar.However, the intercept receiver has no information about the receiving waveform.It can only detect the receiving signal by the signal power in each frequency, so a certain reduction of peak power may have a serious effect on the receiving SNR of intercept receivers.These can be verified by the simulation results in Figure 4b.When the constraints I th of mutual information is 0.0665, which is close to the maximum value 0.0675 of mutual information calculated by Equation ( 21), the SNR reduction for radar is 0.1724 dB, while the SNR reduction is 3.6085 dB for intercept receiver.In addition, the gap between the SNR loss of radar and the SNR loss of intercept receiver keeps increasing with the relaxation of mutual information constraint until I th = 0.0305, where the SNR losses of radar and intercept receiver are 3.4450 dB and 15.4911 dB, respectively.As I th ≤ 0.0305, the optimized waveform is the solution of OP2, which is the red line displayed in Figure 4a, and the SNR losses of both radar and intercept receiver remain unchanged.From Figure 4b, we find that the performance of intercept receiver can be degraded sharply, while the performance of radar has a small reduction.Thus, the optimized waveform has a significant LPI performance.
Compared with the common LPI radar waveforms, our designed waveforms still have a significant LPI superiority, as shown in Table 2.For each compared LPI radar waveform (LFM, Frank, P1-P4), we can find an optimized waveform from different MI constraint values, whose SNR loss of radar is equal to (or less than (since the maximal SNR loss of radar is 3.4450 dB for optimized waveforms)) that of the compared waveform, as the second row shows.However, the SNR losses of intercept receiver are different between each compared waveform and its corresponding optimized waveform.As each column of the third row shows, the SNR loss of intercept receiver of each compared waveform is less than that of its corresponding optimized waveform.Thus, we can draw a conclusion that not only can our proposed optimized waveforms ensure the radar performance, but they can also maximize the performance degradations of intercept receiver.

Only Radar or Intercept Receiver is under a Colored Gaussian Noise Background
In Figure 5, we change the background noise of radar from white Gaussian noise to colored Gaussian noise.Since the power spectral density of the optimized waveform for OP1 is a function of 20)), the power spectral density of the optimized waveform has a large value at the frequencies where σ 2 H ( f ) has a large value and P N ( f ) has a small value.Thus, the optimized waveform (see the dotted purple line in Figure 5a) of OP1 has an extra peak compared with the dotted purple line in Figure 4a.Due to the background noise of intercept receiver is still the white Gaussian noise, the optimized waveform of OP2 remains as a constant P S W , and the optimized waveform of the proposed waveform design method tends to the constant with the reduction of the mutual information constraint I th .From Figure 5b, we find the gap between the SNR loss of radar and the SNR loss of intercept receiver is still huge.As there are such huge gaps, the optimized waveform can achieve a superior LPI performance.
Figure 6a gives the optimized solutions of the proposed optimization problem under the assumption that the background noise of intercept receiver is a colored Gaussian noise.In order to better illustrate the changes of the optimized waveform in this case, we decrease the radar transmit power from 20 KW to 5 KW.Since the solution of OP2 is to make the SNR |S( f )| 2 P N ( f ) in each frequency f ∈ W of intercept receiver reach the minimum, it will match the given background noise of intercept receiver well, as the dotted red line shows in Figure 6a.Like the simulation results above, the target matching quality and peak power of the optimized waveform keeps coming down with the relaxing of the mutual information constraints I th , and the optimized waveform tends to the solution of OP2 eventually.As the Figure 6b shows, such descents have less impact on radar performance, while they hit performance of the intercept receiver hard.The SNR losses of intercepter receiver are above 8 dB and finally reach 17 dB, while the SNR losses of radar keep within 4 dB.Thus, the optimized waveforms in a colored Gaussian noise still have a preferable LPI performance.
Table 3 gives the comparison of SNR losses between the common LPI radar waveforms and our optimized waveforms in the case that the background noise of radar is a colored Gaussian noise; meanwhile, Table 4 gives that in the case that the background noise of intercept receiver is a colored Gaussian noise.Both Tables 3 and 4 show our designed waveforms can further reduce the performance of intercept receiver with almost the same radar performance, compared with common LPI radar waveforms.For example, in Table 3, we find the SNR losses 3.3342 dB of radar for waveforms Frank, P1, P3, P4, are more than the maximal SNR loss 3.3266 dB for optimized waveforms.Although the difference is small, the SNR loss of intercept receiver of optimized waveforms is about 3 dB higher than that of compared waveforms.Therefore, the LPI performance of our designed waveforms is more superior to the common LPI radar waveforms.

Both Radar and Intercept Receiver are under Colored Gaussian Noise
In this subsection, we suppose both the background noises of radar and the intercept receiver are colored Gaussian noise.For better illustrating the changes of the optimized waveforms with the simulation results above, we still reduce the power of transmit waveform from 20 KW to 5 KW.The optimized waveforms should be the combination of the waveforms displayed in Figure 5a and Figure 6a.When the constraints of mutual information are strong, which means I th is close to the maximum mutual information, the optimized waveform should be similar to the corresponding waveforms in Figure 5a.Conversely, as there are almost no constraints of mutual information, that is I th has a smaller value than the maximum mutual information, the optimized waveform is the solution of OP2, as the dotted red lines show in Figures 6a and 7b.For other values of I th , the optimized waveforms should not only have a certain target (σ 2 H ( f )) matching quality, but also be around the solution of OP2, as the purple line shows in Figure 7a.For the comparison between the performance degradations of radar and intercepter displayed in Figure 7b, we still find the SNR losses of intercepter receiver are above 9 dB and finally reach 16 dB, while the SNR losses of radar only keep within 4 dB.Thus, the optimized waveforms on the condition that both radar and intercept receiver are under colored Gaussian noise still have a predominant LPI performance.The comparative results in Table 5 also verify this conclusion, since the optimized waveforms can reduce more performance of the intercept receiver with the same performance loss of radar.

Conclusions
In response to the LPI requirement of modern military radar, waveforms with a fixed average power constraint have been designed in this paper to minimize the interception performance of enemy receivers under the condition of meeting a specified radar performance.In other words, the aim of the proposed LPI radar waveform design method is to reduce the peak power as much as possible under the condition of a certain loss of target matching quality.During the modeling of LPI waveform design, the MI between scattered waveform and radar received waveform is treated as a performance metric of radar, and for enemy receiver, the KLD between intercepted waveform and background noise is presented as its performance metric.Simulation results show that the proposed LPI waveform design method is correct and effective.A certain loss of target matching quality has little impact on radar performance, while a certain reduction of peak power has a serious effect on the performance of enemy receivers.The huge gap of performance loss between radar and intercept receiver is the significance of the proposed LPI waveform design method.

Figure 4 .
Figure 4. (a) optimized waveforms of the proposed mutual information and Kullback-Leibler divergence (MI-KLD) based method with different mutual information constraints I th , and optimized waveforms of OP1 and OP2; (b) performance reductions of radar and intercept receiver versus the reduction of I th (both radar and intercept receiver are under white Gaussian noise background and average transmitted power P S = 20 KW).

Figure 5 .
Figure 5. (a) optimized waveforms of the proposed MI-KLD method with different mutual information constraints I th , and optimized waveforms of OP1 and OP2; (b) performance reductions of radar and intercept receiver versus the reduction of I th (the background noise of radar is a colored Gaussian noise, while the intercept receiver is under the white Gaussian noise background and average transmitted power P S = 20 KW).

Figure 6 .
Figure 6.(a) optimized waveforms of the proposed MI-KLD method with different mutual information constraints I th , and optimized waveforms of OP2; (b) performance reductions of radar and intercept receiver versus the reduction of I th (the background noise of the intercept receiver is a colored Gaussian noise, while the radar is under white Gaussian noise background and average transmitted power P S = 5 KW).

Figure 7 .
Figure 7. (a) optimized waveforms of the proposed MI-KLD method with different mutual information constraints I th , and optimized waveforms of OP2; (b) performance reductions of radar and intercept receiver versus the reduction of I th (both radar and intercept receiver are under a colored Gaussian noise background and average transmitted power P S = 5 KW).
c, 0) ≤ ε 0 Set trust region ∆ k , and calculate r E , r I 8Loop until: E (S k , c k , τ) ≤ ε τ 9 Calculate d by solving Equation (35). 10 If the step d gives sufficient reduction in merit function h (S, c) + ν g E (S) ; g 1 I (S) ; • • • ; g 2N+1 I (S) , 11 set S k+1 ← S k + d S , c k+1 ← c k + d c , 12 Calculate the estimates r E and r I , and increase the trust region, 13 else S k+1 ← S k , c k+1 ← c k , and decrease the trust region.14 Set k ← k + 1. 15 End loop 16

Table 2 .
Signal-to-noise ratio (SNR) losses for the common low probability of interception (LPI) radar waveforms and optimized waveforms (Both radar and intercept receiver are under white Gaussian noise background and average transmitted power P S = 20 KW).

Table 3 .
SNR losses for the common LPI radar waveforms and optimized waveforms (the background noise of radar is a colored Gaussian noise, while the intercept receiver is under white Gaussian noise background and average transmitted power P S = 20 KW).

6468 14.0732 14.0732 13.8503 14.0732 14.0732Table 4 .
SNR losses for the common LPI radar waveforms and optimized waveforms (the background noise of intercept receiver is a colored Gaussian noise, while the radar is under white Gaussian noise background, and average transmitted power P S = 5 KW).

Table 5 .
SNR losses for the common LPI radar waveforms and optimized waveforms (both radar and intercept receiver are under a colored Gaussian noise background, and average transmitted power P S = 5 KW).