Abstract

The system optimization is considered in cognitive radar system (CRS) with extended binary phase shift keying- (EBPSK-) based multicarrier phase-coded (MCPC) signal. A novel radar working scheme is proposed to consider both target detection and estimation. At the detection stage, the generalized likelihood ratio test (GLRT) threshold is deduced, and the GLRT detection probability is given. At the estimation stage, an approach based on Kalman filtering (KF) is proposed to estimate target scattering coefficients (TSC), and the estimation performance is improved significantly by exploiting the TSC temporal correlation. Additionally, the optimal waveform is obtained to minimize the mean square error (MSE) of KF estimation. For the practical consideration, iteration algorithms are proposed to optimize the EBPSK-based MCPC signal in terms of power allocation and coding matrix. Simulation results demonstrate that the KF estimation approach can improve the estimation performance by 25% compared with maximum a posteriori MAP (MAP) method, and the KF estimation performance can be further improved by 90% by optimizing the transmitted waveform spectrum. Moreover, by optimizing the power allocation and coding matrix of the EBPSK-based MCPC signal, the KF estimation performances are, respectively, improved by 7% and 8%.

1. Introduction

The cognitive radar system (CRS) as the future trend of the radar systems, compared with the traditional ones, mainly includes three different aspects [1]: (1) CRS can sense the targets and environment; (2) the transmitted waveform is adaptively optimized to improve the detection and estimation performance; (3) the transmitter, environment, and receiver form a closed loop feedback system. Therefore, optimizing the transmitted waveform based on the working environment becomes a popular research direction in the CRS [28].

For the modeling of CRS, a point is utilized to model the target [9]. The echo waveform from the point target is only the delay or Doppler frequency shift of the original one. However, when the size of target is large enough to occupy many resolution cells, the echo signals from different scattering cells are superimposed, and this type of target is modeled as an extended target (ET) [10]. In addition, the target impulse response (TIR), which is often referred to as the high resolution range profile (HRRP) in the automatic target recognition (ATR) problem [11], can be used to describe ET based on the assumption of the linear time-invariant target [12, 13]. However, the view angle between the target and radar changes, and the ET does not satisfy this assumption. Hence, an exponential correlation model is proposed to describe this time dynamic characteristic, where the TIR is stationary with time and uncorrelated among different resolution cells, namely, wide sense stationary-uncorrelated scattering (WSSUS). Therefore, the TIR of this type target can be estimated by Kalman filtering (KF) [5, 14], where only the estimation is taken into consideration and under the present assumption of the target. However, the practical radar systems should detect the target before estimating.

The performance of the target estimation and detection can be significantly improved by optimizing the transmitted waveform in the CRS. During the stage to estimate the target scattering coefficients (TSC), which is essential the Fourier transform of TIR, the performance can be improved by optimizing the transmitted waveform to maximize the mutual information (MI) between the echo signal and TSC [12]. Furthermore, when the precise priori knowledge of the target cannot be obtained, an eigensubspace projection-based method is proposed to enlarge the separation between the echo signals from different targets in [15, 16]. During the detection stage, in order to maximize the probability of target detection, more power is concentrated on frequencies with relatively large TSC in the additive white Gaussian noise (AWGN) systems or relatively high ratio between the target and clutter in the clutter interference systems [1720]. However, these optimized waveforms have arbitrary nature and are not conveniently generated in the practical radar systems.

The multicarrier phase-coded (MCPC) signal is first applied in the wideband radar systems by Levanon [2123]. Furthermore, the MCPC waveform can achieve the optimal target detection performance with higher spectrum efficiency than the linear frequency modulation (LFM) signal [24, 25]. Unlike the P3 or P4 signal, the MCPC signal has a thumbtack-shaped ambiguity function and is easily generated in the practical radar systems. However, the nonconstant amplitude of MCPC waveform cannot take full advantage of the nonlinear amplifier [26]. On the modulation of MCPC waveform, an extended binary phase shift keying (EBPSK) modulation proposed by Wu et al. is adopted in this work [27, 28]. EBPSK is more flexible than BPSK, where a small angle phase and jump time are utilized to distinguish the modulation waveform of code from that of code , which tightens the spectrum of the transmitted waveform [25, 29, 30]. In the wireless communication system, reference [31] has confirmed that EBPSK can achieve the same theoretical bit error rate (BER) performance as BPSK and achieve higher spectral efficiency with the same bit rate by tuning modulation parameters. Therefore, we first utilize the EBPSK-based MCPC signal in the CRS, since the tight spectrum is easy to optimize.

In this work, the problem of system optimization is considered in the CRS with EBPSK-based MCPC signal, and we propose a new radar working scheme, where both the target detection and estimation are taken into account. During the initial stage of the target detection, the generalized likelihood ratio test (GLRT) algorithm based on constant false alarm rate (CFAR) is utilized. In the presence of ET, the approach based on KF is proposed to estimate TSC by exploiting the temporal correlation, and the optimal spectrum is obtained to minimize the trace of the mean square error (MSE) matrix of KF estimation. In addition, the iteration algorithms are proposed to optimize the amplitudes and coding matrix of EBPSK-based MCPC signal, respectively.

The remainder of this work is organized as follows. In Section 2, the system model of the CRS with EBPSK-based MCPC signal and the new radar working scheme are given. In Section 3, we optimize this radar system, including the GLRT target detection based on CFAR, the maximum a posteriori probability (MAP) receive filter, the TSC estimation based on KF, the optimal power allocation, and coding matrix. Section 4 gives the simulation results. Finally, Section 5 concludes this work.

The notations used in this work are defined as follows. Symbols for matrices (upper case) and vectors (lower case) are in boldface. , , , , , , , and denote the conjugate transpose (Hermitian), the diagonal matrix, the identity matrix of size , the Gaussian distribution with zero mean and covariance being , the absolute value, the norm, the determinant of a matrix, the expectation, and the trace of a matrix, respectively.

2. The System Model of Cognitive Radar with EBPSK-Based MCPC Signal

The system model of cognitive radar considered in this work to detect and estimate the ET is shown in Figure 1, where the transmitted waveform is EBPSK-based MCPC signal and the TIR of ET is during the th pulse. Then the echo signal can be described aswhere denotes the convolution operation, and denotes the AWGN. The EBPSK-based MCPC signal can be expressed aswhere , denotes the entry at th row and th column of the coding matrix , denotes the amplitude of the th subcarrier, and denotes the carrier frequency of the th subcarrier, denotes the carrier frequency of the first subcarrier, and denotes the frequency interval between subcarriers. and are the waveform of code and with EBPSK modulation, respectively,where denotes the time interval of the waveform for one code and denotes the jump interval of the EBPSK modulation. Then we can get the spectrum of transmitted waveform by the Fourier transformwhere and is the Fourier transform of ; then the spectrum of echo signal isThe discrete form in the frequency domain can be obtained aswhere the length of is , denotes the discrete Gaussian noise, denotes the variance of noise, denotes the discrete form of , denotes the TSC, that is, the discrete form of , and the diagonal matrix , where is the discrete form of . According to [14], the exponential correlation model of the TSC iswhere denotes the TSC during the th pulse, denotes the pulse repetition interval (PRI) in the radar system, denotes the temporal decay constant, and follows the Gaussian distribution. When ,and . To simplify the analysis in this work, the correlation between the individual scatters is small enough to assume that is a diagonal matrix; then denotes the variance of , where is the th entry of .

This work proposes a new radar working scheme for the CRS with EBPSK-based MCPC signal. As shown in Figure 2, the target detection and estimation stages in this radar system are described as follows.(1)Initially, the nonoptimized EBPSK-based MCPC signal with power evenly distributed across all subcarriers is transmitted, and the coding matrix is a random matrix.(2)The CRS detects the target via the method based on GLRT, and the detection threshold is obtained by the theory of CFAR. If the power of echo signal is greater than threshold, then the CRS detects the presence of target, otherwise not.(3)When the GLRT detects the presence of target, we propose an approach to estimate the TSC from the echo signal based on KF.(4)In order to improve the performance of KF estimation, this work optimizes the EBPSK-based MCPC signal in terms of power allocation among subcarriers and coding matrix.(5)With the knowledge of the estimated TSC, the signal-to-noise ratio (SNR) of the echo signal can be further improved, which attains a better detection probability for the particular target. More details to describe this working scheme will be given in the following sections.

3. The Optimization of CRS with EBPSK-Based MCPC Signal

3.1. The Initial Target Detection Based on GLRT

Assuming that the target is present and absent , the echo waveform isWith the Gaussian assuming of and , the distribution of echo signal is also Gaussian, and the logarithmic likelihood ratio iswhere is the conditional probability distribution, and the constant term is eliminated in the second equation of (10). According to the GLRT theory, without the knowledge of TIR, the detection process can be expressed aswhereThen we can obtain the GLRT detectionTherefore, without the knowledge of TSC, the GLRT detection is equivalent to the energy detecting. If the energy is greater than the threshold , the target exists, otherwise not.

In the CRS, CFAR is adopted as the target detection criterion, and the echo signal meets the distribution with degrees of freedom when the target is absent, so the false alarm rate iswhere is the gamma function and is the lower incomplete gamma function. When , the Chernoff bound of false alarm rate can be expressed asSince the right side of (15) is the Lambert function, the solution is not easy to obtain, so we use the Taylor series to approximate it, and we can get the detection thresholdWhen the target is present, the process of target detection isSince , where , , and , respectively, denote the th entry of , , and , meets generalized chi-squared distribution and the probability density function is [32]where and . According to the Neyman_Pearson lemma, we can attain the probability of GLRT detection based on CFAR

3.2. The TSC Estimation Based on MAP

After detecting the presence of target, the TSC can be estimated from the echo waveform. When MAP is adopted as the receiver filter, we haveAccording to Bayes’ theorem, we havewhere and are both the function of Gaussian distribution and the distribution of isWe can obtain the posterior probability by substituting (22) into (21) aswhere and . We can get the optimization objective function by substituting (23) into (20) asTherefore, the estimated TSC by MAP estimation isThe receiver filter can be written in the matrix formThen the MSE matrix of the estimated TSC is

3.3. The TSC Estimation Based on KF

We proposed an approach based on KF to estimate TSC according to the temporal correlation of TSC, when the GLRT detects the presence of target. Furthermore, we optimize the transmitted waveform to further improve the estimation performance by minimizing the trace of MSE matrix. The estimation approach based on KF is described in Algorithm 1.

(1)  Set the pulse index = 1, and the initial TSC estimation is based on MAP, where the estimated
   TSC can be calculated from (25) and the MSE matrix can be obtained from (27).
   Set the maximum number of pulses .
(2)  Let = + 1;
(3)  while   and > 1 do
(4)   According to the temporal correlation of TSC (4), the prediction of TSC is
        ,     (i)
    where is the estimated TSC during the (  − 1)th pulse;
(5)   The MSE matrix of the predicted TSC is
    ; (ii)
(6)   The KF gain matrix is defined as
    ;  (iii)
(7)   Update the estimated TSC
    ,    (iv)
    where ;
(8)   Update the MSE matrix
        ;   (v)
(9)   Let = + 1.
(10) end while

In order to improve the performance of KF estimation, the trace of the MSE matrix is utilized to measure this performance,The transmitted waveform should be optimized to minimize the MSE at each iteration of KF, which is equivalent to the following optimization problem:where is the power of radar signal. The objective function can be simplified aswhere denotes the entry at th row and th column of . Then we can get the Lagrangian function of this optimization problemLet the derivation of (31) equal ;Therefore, the optimal spectrum of transmitted waveform should satisfy the following condition:where is the constraint of transmitted power and the optimal signal obtained in (34) is similar to the water-filling method in the communication systems [12].

3.4. The Optimization of EBPSK-Based MCPC Signal

In the practical radar systems, the waveforms are generated by setting the modulation parameters or choosing from the preset waveforms, so it is difficult to attain the optimal waveform calculated in Section 3.3. Therefore, we propose a suboptimal signal by tuning the parameters of EBPSK-based MCPC waveform in terms of power allocation among subcarriers and coding matrix. In addition, the process of waveform optimization is simple while maintaining the estimation performance of TSC.

The optimal waveform obtained in (33) is referred as ; then the best power allocation scheme for EBPSK-based MCPC signal should approximate the spectrum to this optimal one. Therefore, the MSE approximation can be expressed as an optimization problemwhere denotes EBPSK-based MCPC signal in the frequency domain. Because of the sparsity of the EBPSK spectrum [30], we only need to consider the power allocation in carrier frequencies of subcarriers. The amplitudes should satisfy the following condition:where donates the th entry of . The optimized power allocation (35) of EBPSK-based MCPC signal can improve the estimation performance of TSC via KF, and the simulation results will be given in Section 5.

We will also optimize the coding matrix in the following, which can further reduce the MSE of the estimated TSC based on KF. The optimal coding matrix should be the solution of the following optimization problem:where denotes the discrete form the transmitted waveform and denotes the optimal reference signal, which will be obtained by the estimated TSC and the target scattering model in the following.

When the estimation value of scattering coefficients is , the predicted value isHowever, the real one isWe should design the transmitted signal , which minimizes the estimation MSE with the assumption of  . Since we have proposed the optimal power allocation scheme by the water-filling method, the estimation scattering vector is chosen as the optimal reference waveform, and (37) is the optimization objective function. Letwhere is the imaginary unit. A simple iteration algorithm is proposed to calculate the optimal coding matrix, which is described in Algorithm 2.

(1) Initialize the coding matrix with a random matrix,
  and set the iteration index and the maximum number of iterations ;
(2) while   and the coding matrix changes do
(3)  Let , ,
   where denotes the entry at the th row and th column of ;
(4)  Generate the EBPSK-based MCPC waveform with coding matrix and calculate
       ;             (vi)
(5)  If , , and go to Step  3,
   otherwise, let , and go to Step  3;
(6)  Let .
(7) end while

4. Simulation Results

4.1. The TSC Estimation Based on MAP

First, we simulate the MSE of estimated TSC based on the MAP receiver filter and compare the simulation results with the theoretical ones. The parameters are set in Table 1, where the transmitted waveform is the EBPSK-based MCPC signal. In this work, all the simulation parameter are given based on the published literatures. For example, the parameters of EBPSK modulation are according to [20, 21, 2729], and the parameters for the characteristics of target and radar are based on [5, 14]. Although the realistic data is not adopted in our work, the simulation results are reliable and practicable.

Figure 3 depicts a comparison between theoretical and simulation results of the normalized MSE of estimated TSC, where the theoretical results are calculated by the trace of MSE matrix in (28), and the normalized MSE is defined asWe can observe from this figure that the simulation and theoretical results have a relatively high degree of agreement, which verifies the efficiency to measure the estimation performance of MAP by the trace of MSE matrix.

4.2. The TSC Estimation Based on KF

This subsection gives the simulation results of the estimated TSC based on KF, in which the simulation parameters are the same with Section 4.1 except SNR = dB, and the index of iteration is from to . Figure 4 demonstrates the normalized MSE of estimated TSC, including the approach based on KF and MAP. Apparently, the estimation performance is improved about by using the KF approach instead of MAP, since the temporal correlation is exploited by KF. In addition, the theoretical results given in (28) have high degree of agreement with the simulation results. Therefore, the trace of MSE matrix can be utilized as an objective function, which accurately describes the estimation performance of KF, in the waveform optimization problem. However, during the simulation, the decreasing normalized MSE of MAP estimation is caused by the exponential correlation model .

Figure 5 illustrates the spectrum of transmitted waveform, including the power uniform distributed EBPSK-based MCPC signal and the optimal spectrum. The optimal spectrum, which can minimize the trace of MSE matrix and improve the estimation performance of KF approach, is obtained by the water-filling method, where the parameter is determined by the dichotomy algorithm. In addition, as shown in this figure, the variance of TSC is also given, which is similar to the optimal spectrum. Therefore, more power should be concentrated on the frequencies with relatively large variance of TSC, which provides guidance on the optimization of EBPSK-based MCPC signal.

Figure 6 shows the estimation performance of the methods based on MAP and KF, where the waveform with optimal spectrum and the random one is compared. We can observe that the simulated MSE and theoretical one have a relatively high degree of agreement. Moreover, the accuracy of estimation based on KF is more precise than the one based on MAP. In addition, optimizing the spectrum of transmitted waveform can further improve the estimation performance both by the approaches based on MAP and KF though the criterion of waveform optimization is the MSE matrix of the KF estimation. The performance of MAP estimation is improved by , and the performance of KF estimation is also improved by . Moreover, the KF estimation with optimized waveform is better than the MAP one with optimized waveform.

4.3. The Optimization of EBPSK-Based MCPC Signal

Figure 7 depicts the spectrum of the optimal power allocation for EBPSK-based MCPC signal, where the KF estimation is adopted. The optimized power allocation shows that more power should be concentrated on the frequencies with relatively large variance of TSC. Therefore, we should design the EBPSK-based MCPC signal to approximate the signal with optimal spectrum, which can be calculated in (33). Figure 8 compares the normalized MSE of EBPSK-based MCPC signal with optimized power allocation with the one with uniform distributed power. We can observe that the estimation performance in the CRS can be improved significantly by optimizing the transmission power of each subcarrier, where the KF estimation with optimized power allocation of EBPSK-based MCPC signal achieves better than the one with uniform power allocation. Although the estimation performance of the optimized EBPSK-based MCPC signal is worse than the waveform with optimal spectrum (33), we can generate the optimal EBPSK-based MCPC signal by tuning the amplitudes of each subcarrier, which is more simple and convenient.

Figure 9 compares the estimation performance of two types of optimized EBPSK-based MCPC signal, where the first one optimizes the power allocation and coding matrix simultaneously, and the second one only optimizes the power allocation. As the figure shows, the estimation performance is improved by simultaneously optimizing these two parameters. Furthermore, the algorithm proposed in this work to optimize the coding matrix is easy to realize.

5. Conclusions

In this work, we have investigated the problem of system optimization in the CRS with EBPSK-based MCPC signal, and a new radar working scheme has been proposed to take both the target detection and estimation into consideration. At the detection stage, we have deduced the theoretical threshold and the detection probability based on GLRT. At the estimation stage, the TSC is estimated by the approach based on KF. By exploiting the temporal correlation of TSC, the estimation performance of KF is better than that of MAP. In addition, the transmitted waveform with optimal spectrum is obtained, which can improve the estimation performance of KF by . For the practical consideration, the iteration algorithms are proposed to optimize the EBPSK-based MCPC signal in terms of power allocation and coding matrix, which can achieve and performance improvement, respectively. In the future work, we will put our force on the system optimization of the temporal correlated CRS in the presence of clutter and jam.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publishing of this paper.

Acknowledgment

This work was supported in part by the National Natural Science Foundation of China (Grant no. 61271204) and the National Key Technology R&D Program (Grant no. 2012BAH15B00).