# System Optimization for Temporal Correlated Cognitive Radar with EBPSK-Based MCPC Signal

Authors: Peng Chen; Lenan Wu
Journal: Mathematical Problems in Engineering
 (2015)

---

## 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%.


---

## Body

## 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 [2–8].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 [17–20]. 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 [21–23]. 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 1 from that of code 0, 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 Section2, 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.(
·
)
H, diag
·, I
L, N
0
,
R, ·, ·
2, det
⁡
·, E
{
·
} and Tr
⁡
· denote the conjugate transpose (Hermitian), the diagonal matrix, the identity matrix of size L, the Gaussian distribution with zero mean and covariance being R, the absolute value, the l
2 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 Figure1, where the transmitted waveform is EBPSK-based MCPC signal s
k
t and the TIR of ET is h
k
t during the kth pulse. Then the echo signal y
k
t can be described as(1)
y
k
t
=
s
k
t
*
h
k
t
+
n
k
t
,
k
=
1,2
,
3
,
…
,where * denotes the convolution operation, and n
k
t denotes the AWGN. The EBPSK-based MCPC signal can be expressed as(2)
s
k
t
=
∑
l
=
1
L
‍
a
k
,
l
∑
r
=
1
R
‍
r
k
,
l
r
t
,where r
k
,
l
r
t
≜
C
k
,
l
r
s
1
t
-
T
s
r
,
f
l
+
1
-
C
k
,
l
r
s
0
t
-
T
s
r
,
f
l, C
k
,
l
r
∈
0,1 denotes the entry at lth row and rth column of the coding matrix C
k, a
k
,
l denotes the amplitude of the lth subcarrier, and f
l
≜
f
c
+
l
-
1
Δ
f denotes the carrier frequency of the lth subcarrier, f
c denotes the carrier frequency of the first subcarrier, and Δ
f denotes the frequency interval between subcarriers. s
1
t
,
f and s
0
t
,
f are the waveform of code 0 and 1 with EBPSK modulation, respectively,(3)
s
0
t
,
f
=
cos
⁡
2
π
f
t
,
0
≤
t
≤
T
s
s
1
t
,
f
=
-
cos
⁡
2
π
f
t
,
0
≤
t
≤
τ
cos
⁡
2
π
f
t
,
τ
≤
t
≤
T
s
,where T
s
≜
N
/
f
c
N
=
1,2
,
3
,
… denotes the time interval of the waveform for one code and τ
≜
K
/
f
c
K
=
1,2
,
3
,
…
,
K
≤
N denotes the jump interval of the EBPSK modulation. Then we can get the spectrum of transmitted waveform by the Fourier transform(4)
S
k
f
=
∑
l
=
1
L
‍
a
k
,
l
∑
r
=
1
R
‍
R
k
,
l
r
f
exp
⁡
-
j
2
π
f
T
s
r
,where R
k
,
l
r
f
≜
C
k
,
l
r
S
1
f
,
f
l
+
1
-
C
k
,
l
r
S
0
f
,
f
l and S
i
f
,
f
l
i
∈
0,1 is the Fourier transform of s
i
t
,
f
l; then the spectrum of echo signal Y
f is(5)
Y
k
f
=
H
k
f
S
k
f
+
N
k
f
.The discrete form in the frequency domain can be obtained as(6)
y
k
=
S
k
g
k
+
n
k
,where the length of y
k is M, n
k
~
N
0
,
σ
n
2
I
M denotes the discrete Gaussian noise, σ
n
2 denotes the variance of noise, y
k denotes the discrete form of Y
k
f, g
k denotes the TSC, that is, the discrete form of H
k
f, and the diagonal matrix S
k
≜
diag
s
k, where s
k is the discrete form of S
k
f. According to [14], the exponential correlation model of the TSC is(7)
g
k
=
e
-
T
/
τ
g
k
-
1
+
u
k
-
1
,where g
k denotes the TSC during the kth pulse, T denotes the pulse repetition interval (PRI) in the radar system, τ denotes the temporal decay constant, and u
k
-
1
~
N
0
,
1
-
e
-
T
/
τ
Σ
g follows the Gaussian distribution. When g
0
=
0,(8)
g
k
=
∑
m
=
1
k
-
1
‍
e
-
m
T
/
τ
u
k
-
1
-
m
,and g
k
~
N
0
,
R
k
≜
1
-
e
-
2
T
/
τ
Σ
g. To simplify the analysis in this work, the correlation between the individual scatters is small enough to assume that Σ
g is a diagonal matrix; then σ
k
,
G
m
2 denotes the variance of g
k
,
m, where g
k
,
m is the mth entry of g
k.Figure 1
The system model of cognitive radar with EBPSK-based MCPC signal.This work proposes a new radar working scheme for the CRS with EBPSK-based MCPC signal. As shown in Figure2, 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±
1 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.Figure 2
The proposed working scheme for the CRS with EBPSK-based MCPC signal.


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

### 3.1. The Initial Target Detection Based on GLRT

Assuming that the target is presentH
1 and absent H
0, the echo waveform is(9)
H
0
:
y
k
=
n
k
,
H
1
:
y
k
=
S
k
g
k
+
n
k
.With the Gaussian assuming of g
k and n
k, the distribution of echo signal is also Gaussian, and the logarithmic likelihood ratio is(10)
l
y
k
∣
g
k
≜
ln
⁡
p
y
k
∣
g
k
,
H
1
p
y
k
∣
H
0
iiiiiiiiiiiiiii
=
1
σ
n
2
2
y
k
H
-
g
k
H
S
k
H
S
k
g
k
,where p
y
k
∣
H
0 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 as(11)
l
y
k
≜
l
y
k
∣
g
¯
k
⋛
H
1
H
0
λ
,where(12)
g
¯
k
=
arg
max
⁡
g
k
l
y
k
∣
g
k
g
¯
k
=
arg
max
⁡
g
k
y
k
-
S
k
g
k
2
=
S
k
-
1
y
k
.Then we can obtain the GLRT detection(13)
l
y
k
=
1
σ
n
2
y
k
H
y
k
⋛
H
1
H
0
λ
.Therefore, without the knowledge of TSC, the GLRT detection is equivalent to the energy detecting. If the energy y
k
2
2 is greater than the threshold λ
σ
n
2, the target exists, otherwise not.In the CRS, CFAR is adopted as the target detection criterion, and the echo signal meets theχ
2 distribution with M degrees of freedom when the target is absent, so the false alarm rate is(14)
P
F
A
=
P
l
y
k
∣
H
0
≥
λ
=
1
-
1
Γ
M
/
2
γ
M
2
,
λ
2
,where Γ
x is the gamma function and γ
s
,
t is the lower incomplete gamma function. When λ
>
M, the Chernoff bound of false alarm rate can be expressed as(15)
P
F
A
≤
λ
M
exp
⁡
1
-
λ
M
M
/
2
.Since the right side of (15) is the Lambert W function, the solution is not easy to obtain, so we use the Taylor series to approximate it, and we can get the detection threshold(16)
λ
P
F
A
≈
M
P
F
A
2
/
M
+
M
1
-
P
F
A
2
/
M
-
3
M
2
1
-
P
F
A
2
/
M
iiiiiiiiiiiiii
+
8
M
3
1
-
P
F
A
2
/
M
4
-
125
M
24
1
-
P
F
A
2
/
M
5
.When the target is present, the process of target detection is(17)
l
y
k
∣
H
1
=
1
σ
n
2
S
k
g
k
+
n
k
H
S
k
g
k
+
n
k
⋛
H
1
H
0
λ
P
F
A
.Since s
k
,
m
g
k
,
m
+
n
k
,
m
~
N
0
,
s
k
,
m
2
σ
G
m
2
+
σ
n
2, where s
k
,
m, g
k
,
m, and n
k
,
m, respectively, denote the mth entry of s
k, g
k, and n
k, x
≜
l
y
k
H
1 meets generalized chi-squared distribution and the probability density function is [32](18)
f
x
=
∑
m
=
1
M
‍
exp
⁡
-
x
/
ξ
m
2
ξ
m
2
∏
i
=
1
,
i
≠
m
M
‍
1
-
ξ
i
2
/
ξ
m
2
,where ξ
m
2
=
1
+
s
k
,
m
2
σ
G
m
2
/
σ
n
2 and ξ
m
2
≠
ξ
n
2
∀
m
≠
n. According to the Neyman_Pearson lemma, we can attain the probability of GLRT detection based on CFAR(19)
P
D
=
∫
λ
P
F
A
∞
‍
f
x
d
x
P
D
=
∑
m
=
1
M
‍
exp
⁡
-
λ
P
F
A
/
ξ
m
2
∏
i
=
1
,
i
≠
m
M
‍
1
-
ξ
i
2
/
ξ
m
2
.


### 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 have(20)
g
^
k
=
arg
max
⁡
g
k
p
g
k
∣
y
k
.According to Bayes’ theorem, we have(21)
p
g
k
∣
y
k
=
p
g
k
,
y
k
p
y
k
=
p
y
k
∣
g
k
p
g
k
p
y
k
,where p
y
k
∣
g
k and p
g
k are both the function of Gaussian distribution and the distribution of y
k is(22)
p
y
k
=
∫
‍
p
y
k
∣
g
k
p
g
k
d
g
k
.We can obtain the posterior probability by substituting (22) into (21) as(23)
p
g
k
y
k
=
p
y
k
∣
g
k
p
g
k
p
y
k
=
exp
⁡
g
k
H
1
/
σ
n
2
S
k
H
y
k
-
1
/
2
S
k
g
k
-
1
/
2
R
k
-
1
g
k
exp
⁡
1
/
2
σ
n
4
y
k
H
S
k
B
k
-
H
S
k
H
y
k
det
⁡
π
B
k
-
1
,where B
k
≜
1
/
σ
n
2
S
k
H
S
k
+
R
k
-
1 and B
k
-
H
≜
B
k
H
-
1. We can get the optimization objective function by substituting (23) into (20) as(24)
g
^
k
=
arg
max
⁡
g
k
-
1
2
g
k
H
B
k
g
k
+
1
σ
n
2
S
k
g
k
H
y
k
.Therefore, the estimated TSC by MAP estimation is(25)
g
^
k
=
1
σ
n
2
B
k
-
1
S
k
H
y
k
.The receiver filter can be written in the matrix form(26)
Q
k
≜
1
σ
n
2
B
k
-
1
S
k
H
.Then the MSE matrix of the estimated TSC is(27)
P
=
E
g
^
k
-
g
k
g
^
k
-
g
k
H
P
=
Q
k
S
k
R
k
S
k
H
+
σ
n
2
I
M
Q
k
H
iiiiiii
-
1
σ
n
2
B
k
-
1
S
k
H
S
k
R
k
-
1
σ
n
2
R
k
S
k
H
S
k
B
k
-
H
+
R
k
.


### 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 Algorithm1.Algorithm 1:The TSC estimation based on KF.
(1)  Set the pulse indexk = 1, and the initial TSC estimation is based on MAP, where the estimated
TSCg
^
k
∣
k can be calculated from (25) and the MSE matrix P
k
∣
k can be obtained from (27).
Set the maximum number of pulsesK
max
⁡.
(2)  Letk = k + 1;
(3)while  k ≤ K
max
⁡ and k > 1 do
(4) According to the temporal correlation of TSC (4), the prediction of TSC is
g
^
k
∣
k
-
1
=
e
-
T
/
τ
g
^
k
-
1
∣
k
-
1,     (i)
where g
^
k
-
1
∣
k
-
1 is the estimated TSC during the (k  − 1)th pulse;
(5) The MSE matrix of the predicted TSC is
P
k
∣
k
-
1
=
e
-
2
T
/
τ
P
k
-
1
∣
k
-
1
+
(
1
-
e
-
2
T
/
τ
)
Σ
g; (ii)
(6) The KF gain matrix is defined as
Φ
k
≜
P
k
∣
k
-
1
S
k
H
(
σ
n
2
Q
k
+
Q
k
S
k
P
k
∣
k
-
1
S
k
H
)
-
1;  (iii)
(7)  Update the estimated TSC
g
^
k
∣
k
=
g
^
k
∣
k
-
1
+
Φ
k
(
g
^
k
-
Q
k
S
k
g
^
k
∣
k
-
1
),    (iv)
where g
^
k
=
Q
k
y
k;
(8) Update the MSE matrix
P
k
∣
k
=
P
k
∣
k
-
1
-
Φ
k
Q
k
S
k
P
k
∣
k
-
1;   (v)
(9) Let k = k + 1.
(10)end whileIn order to improve the performance of KF estimation, the trace of the MSE matrixP
k
∣
k is utilized to measure this performance,(28)
f
s
k
=
Tr
⁡
P
k
∣
k
.The transmitted waveform should be optimized to minimize the MSE at each iteration of KF, which is equivalent to the following optimization problem:(29)
s
k
*
=
arg
min
⁡
s
k
f
s
k
s
.
t
.
Tr
⁡
S
k
H
S
k
=
E
s
i
i
i
I
I
i
i
i
S
k
=
diag
s
k
,where E
s is the power of radar signal. The objective function can be simplified as(30)
f
s
k
=
Tr
⁡
P
k
∣
k
f
s
k
=
∑
m
=
1
M
‍
s
m
2
P
k
∣
k
-
1
,
m
2
σ
n
2
+
P
k
∣
k
-
1
,
m
s
m
2
,where P
k
∣
k
-
1
,
m denotes the entry at mth row and mth column of P
k
∣
k. Then we can get the Lagrangian function of this optimization problem(31)
g
s
k
,
ζ
=
f
s
k
-
ζ
∑
m
=
1
M
‍
s
m
2
-
E
s
.Let the derivation of (31) equal 0;(32)
∂
g
s
k
,
ζ
∂
s
m
2
=
0
.Therefore, the optimal spectrum of transmitted waveform should satisfy the following condition:(33)
s
m
*
2
=
max
⁡
0
,
α
σ
n
2
-
σ
n
2
P
k
∣
k
-
1
,
m
,where α is the constraint of transmitted power Tr
⁡
S
k
H
S
k
=
E
s and the optimal signal obtained in (34) is similar to thewater-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 Section3.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 s
k
*; 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 problem(34)
s
k
′
*
=
arg
max
⁡
s
k
′
=
s
k
′
2
-
s
k
*
2
2
2
,where s
k
′ 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:(35)
a
m
*
a
n
*
=
s
k
,
m
*
s
k
,
n
*
,where s
k
,
m
* donates the mth entry of s
k
*. 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:(36)
C
k
*
=
arg
min
⁡
C
k
s
~
k
H
x
s
~
k
2
x
2
s
.
t
.
iiii
C
l
r
∈
0,1
,where s
~
k denotes the discrete form the transmitted waveform S
~
k
f
=
∑
l
=
1
L
a
k
,
l
∑
r
=
1
R
C
k
,
l
r
S
1
f
,
f
l
+
1
-
C
k
,
l
r
S
0
f
,
f
l
exp
⁡
-
j
2
π
f
T
s
r and x 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 isg
^
k
-
1
∣
k
-
1, the predicted value is(37)
g
^
k
k
-
1
=
e
-
T
/
τ
g
^
k
-
1
∣
k
-
1
.However, the real one is(38)
g
k
=
e
-
T
/
τ
g
k
-
1
+
u
k
-
1
.We should design the transmitted signal x, which minimizes the estimation MSE with the assumption of  g
^
k
-
1
≈
g
k
-
1. 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. Let(39)
x
=
g
^
k
∣
k
-
1
+
1
+
i
1
2
1
-
e
-
2
T
/
τ
Σ
g
,where i is the imaginary unit. A simple iteration algorithm is proposed to calculate the optimal coding matrix, which is described in Algorithm 2.Algorithm 2:The algorithm to optimize the coding matrix.
(1) Initialize the coding matrixC
* with a random ±
1 matrix,
and set the iteration indexi
=
1 and the maximum number of iterations I
max
⁡;
(2)while
i
≤
I
max
⁡ and the coding matrix changes do
(3)  Letr
=
i
mod
⁡
M
L, C
*
r
/
M
+
1
,
r
-
M
r
/
M
=
1
-
C
*
r
/
M
+
1
,
r
-
M
r
/
M,
whereC
*
m
,
n denotes the entry at the mth row and nth column of C
*;
(4)  Generate the EBPSK-based MCPC waveforms
~
i with coding matrix C
* and calculate
z
i
=
s
~
i
2
-
x
2
2
2;             (vi)
(5)  Ifz
i
>
z
i
-
1, i
=
i
+
1, and go to Step  3,
otherwise, letC
*
r
/
M
+
1
,
r
-
M
r
/
M
=
1
-
C
*
r
/
M
+
1
,
r
-
M
r
/
M, i
=
i
+
1 and go to Step  3;
(6)  Leti
=
i
+
1.
(7)end while


## 3.1. The Initial Target Detection Based on GLRT

Assuming that the target is presentH
1 and absent H
0, the echo waveform is(9)
H
0
:
y
k
=
n
k
,
H
1
:
y
k
=
S
k
g
k
+
n
k
.With the Gaussian assuming of g
k and n
k, the distribution of echo signal is also Gaussian, and the logarithmic likelihood ratio is(10)
l
y
k
∣
g
k
≜
ln
⁡
p
y
k
∣
g
k
,
H
1
p
y
k
∣
H
0
iiiiiiiiiiiiiii
=
1
σ
n
2
2
y
k
H
-
g
k
H
S
k
H
S
k
g
k
,where p
y
k
∣
H
0 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 as(11)
l
y
k
≜
l
y
k
∣
g
¯
k
⋛
H
1
H
0
λ
,where(12)
g
¯
k
=
arg
max
⁡
g
k
l
y
k
∣
g
k
g
¯
k
=
arg
max
⁡
g
k
y
k
-
S
k
g
k
2
=
S
k
-
1
y
k
.Then we can obtain the GLRT detection(13)
l
y
k
=
1
σ
n
2
y
k
H
y
k
⋛
H
1
H
0
λ
.Therefore, without the knowledge of TSC, the GLRT detection is equivalent to the energy detecting. If the energy y
k
2
2 is greater than the threshold λ
σ
n
2, the target exists, otherwise not.In the CRS, CFAR is adopted as the target detection criterion, and the echo signal meets theχ
2 distribution with M degrees of freedom when the target is absent, so the false alarm rate is(14)
P
F
A
=
P
l
y
k
∣
H
0
≥
λ
=
1
-
1
Γ
M
/
2
γ
M
2
,
λ
2
,where Γ
x is the gamma function and γ
s
,
t is the lower incomplete gamma function. When λ
>
M, the Chernoff bound of false alarm rate can be expressed as(15)
P
F
A
≤
λ
M
exp
⁡
1
-
λ
M
M
/
2
.Since the right side of (15) is the Lambert W function, the solution is not easy to obtain, so we use the Taylor series to approximate it, and we can get the detection threshold(16)
λ
P
F
A
≈
M
P
F
A
2
/
M
+
M
1
-
P
F
A
2
/
M
-
3
M
2
1
-
P
F
A
2
/
M
iiiiiiiiiiiiii
+
8
M
3
1
-
P
F
A
2
/
M
4
-
125
M
24
1
-
P
F
A
2
/
M
5
.When the target is present, the process of target detection is(17)
l
y
k
∣
H
1
=
1
σ
n
2
S
k
g
k
+
n
k
H
S
k
g
k
+
n
k
⋛
H
1
H
0
λ
P
F
A
.Since s
k
,
m
g
k
,
m
+
n
k
,
m
~
N
0
,
s
k
,
m
2
σ
G
m
2
+
σ
n
2, where s
k
,
m, g
k
,
m, and n
k
,
m, respectively, denote the mth entry of s
k, g
k, and n
k, x
≜
l
y
k
H
1 meets generalized chi-squared distribution and the probability density function is [32](18)
f
x
=
∑
m
=
1
M
‍
exp
⁡
-
x
/
ξ
m
2
ξ
m
2
∏
i
=
1
,
i
≠
m
M
‍
1
-
ξ
i
2
/
ξ
m
2
,where ξ
m
2
=
1
+
s
k
,
m
2
σ
G
m
2
/
σ
n
2 and ξ
m
2
≠
ξ
n
2
∀
m
≠
n. According to the Neyman_Pearson lemma, we can attain the probability of GLRT detection based on CFAR(19)
P
D
=
∫
λ
P
F
A
∞
‍
f
x
d
x
P
D
=
∑
m
=
1
M
‍
exp
⁡
-
λ
P
F
A
/
ξ
m
2
∏
i
=
1
,
i
≠
m
M
‍
1
-
ξ
i
2
/
ξ
m
2
.


## 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 have(20)
g
^
k
=
arg
max
⁡
g
k
p
g
k
∣
y
k
.According to Bayes’ theorem, we have(21)
p
g
k
∣
y
k
=
p
g
k
,
y
k
p
y
k
=
p
y
k
∣
g
k
p
g
k
p
y
k
,where p
y
k
∣
g
k and p
g
k are both the function of Gaussian distribution and the distribution of y
k is(22)
p
y
k
=
∫
‍
p
y
k
∣
g
k
p
g
k
d
g
k
.We can obtain the posterior probability by substituting (22) into (21) as(23)
p
g
k
y
k
=
p
y
k
∣
g
k
p
g
k
p
y
k
=
exp
⁡
g
k
H
1
/
σ
n
2
S
k
H
y
k
-
1
/
2
S
k
g
k
-
1
/
2
R
k
-
1
g
k
exp
⁡
1
/
2
σ
n
4
y
k
H
S
k
B
k
-
H
S
k
H
y
k
det
⁡
π
B
k
-
1
,where B
k
≜
1
/
σ
n
2
S
k
H
S
k
+
R
k
-
1 and B
k
-
H
≜
B
k
H
-
1. We can get the optimization objective function by substituting (23) into (20) as(24)
g
^
k
=
arg
max
⁡
g
k
-
1
2
g
k
H
B
k
g
k
+
1
σ
n
2
S
k
g
k
H
y
k
.Therefore, the estimated TSC by MAP estimation is(25)
g
^
k
=
1
σ
n
2
B
k
-
1
S
k
H
y
k
.The receiver filter can be written in the matrix form(26)
Q
k
≜
1
σ
n
2
B
k
-
1
S
k
H
.Then the MSE matrix of the estimated TSC is(27)
P
=
E
g
^
k
-
g
k
g
^
k
-
g
k
H
P
=
Q
k
S
k
R
k
S
k
H
+
σ
n
2
I
M
Q
k
H
iiiiiii
-
1
σ
n
2
B
k
-
1
S
k
H
S
k
R
k
-
1
σ
n
2
R
k
S
k
H
S
k
B
k
-
H
+
R
k
.


## 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 Algorithm1.Algorithm 1:The TSC estimation based on KF.
(1)  Set the pulse indexk = 1, and the initial TSC estimation is based on MAP, where the estimated
TSCg
^
k
∣
k can be calculated from (25) and the MSE matrix P
k
∣
k can be obtained from (27).
Set the maximum number of pulsesK
max
⁡.
(2)  Letk = k + 1;
(3)while  k ≤ K
max
⁡ and k > 1 do
(4) According to the temporal correlation of TSC (4), the prediction of TSC is
g
^
k
∣
k
-
1
=
e
-
T
/
τ
g
^
k
-
1
∣
k
-
1,     (i)
where g
^
k
-
1
∣
k
-
1 is the estimated TSC during the (k  − 1)th pulse;
(5) The MSE matrix of the predicted TSC is
P
k
∣
k
-
1
=
e
-
2
T
/
τ
P
k
-
1
∣
k
-
1
+
(
1
-
e
-
2
T
/
τ
)
Σ
g; (ii)
(6) The KF gain matrix is defined as
Φ
k
≜
P
k
∣
k
-
1
S
k
H
(
σ
n
2
Q
k
+
Q
k
S
k
P
k
∣
k
-
1
S
k
H
)
-
1;  (iii)
(7)  Update the estimated TSC
g
^
k
∣
k
=
g
^
k
∣
k
-
1
+
Φ
k
(
g
^
k
-
Q
k
S
k
g
^
k
∣
k
-
1
),    (iv)
where g
^
k
=
Q
k
y
k;
(8) Update the MSE matrix
P
k
∣
k
=
P
k
∣
k
-
1
-
Φ
k
Q
k
S
k
P
k
∣
k
-
1;   (v)
(9) Let k = k + 1.
(10)end whileIn order to improve the performance of KF estimation, the trace of the MSE matrixP
k
∣
k is utilized to measure this performance,(28)
f
s
k
=
Tr
⁡
P
k
∣
k
.The transmitted waveform should be optimized to minimize the MSE at each iteration of KF, which is equivalent to the following optimization problem:(29)
s
k
*
=
arg
min
⁡
s
k
f
s
k
s
.
t
.
Tr
⁡
S
k
H
S
k
=
E
s
i
i
i
I
I
i
i
i
S
k
=
diag
s
k
,where E
s is the power of radar signal. The objective function can be simplified as(30)
f
s
k
=
Tr
⁡
P
k
∣
k
f
s
k
=
∑
m
=
1
M
‍
s
m
2
P
k
∣
k
-
1
,
m
2
σ
n
2
+
P
k
∣
k
-
1
,
m
s
m
2
,where P
k
∣
k
-
1
,
m denotes the entry at mth row and mth column of P
k
∣
k. Then we can get the Lagrangian function of this optimization problem(31)
g
s
k
,
ζ
=
f
s
k
-
ζ
∑
m
=
1
M
‍
s
m
2
-
E
s
.Let the derivation of (31) equal 0;(32)
∂
g
s
k
,
ζ
∂
s
m
2
=
0
.Therefore, the optimal spectrum of transmitted waveform should satisfy the following condition:(33)
s
m
*
2
=
max
⁡
0
,
α
σ
n
2
-
σ
n
2
P
k
∣
k
-
1
,
m
,where α is the constraint of transmitted power Tr
⁡
S
k
H
S
k
=
E
s and the optimal signal obtained in (34) is similar to thewater-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 Section3.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 s
k
*; 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 problem(34)
s
k
′
*
=
arg
max
⁡
s
k
′
=
s
k
′
2
-
s
k
*
2
2
2
,where s
k
′ 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:(35)
a
m
*
a
n
*
=
s
k
,
m
*
s
k
,
n
*
,where s
k
,
m
* donates the mth entry of s
k
*. 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:(36)
C
k
*
=
arg
min
⁡
C
k
s
~
k
H
x
s
~
k
2
x
2
s
.
t
.
iiii
C
l
r
∈
0,1
,where s
~
k denotes the discrete form the transmitted waveform S
~
k
f
=
∑
l
=
1
L
a
k
,
l
∑
r
=
1
R
C
k
,
l
r
S
1
f
,
f
l
+
1
-
C
k
,
l
r
S
0
f
,
f
l
exp
⁡
-
j
2
π
f
T
s
r and x 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 isg
^
k
-
1
∣
k
-
1, the predicted value is(37)
g
^
k
k
-
1
=
e
-
T
/
τ
g
^
k
-
1
∣
k
-
1
.However, the real one is(38)
g
k
=
e
-
T
/
τ
g
k
-
1
+
u
k
-
1
.We should design the transmitted signal x, which minimizes the estimation MSE with the assumption of  g
^
k
-
1
≈
g
k
-
1. 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. Let(39)
x
=
g
^
k
∣
k
-
1
+
1
+
i
1
2
1
-
e
-
2
T
/
τ
Σ
g
,where i is the imaginary unit. A simple iteration algorithm is proposed to calculate the optimal coding matrix, which is described in Algorithm 2.Algorithm 2:The algorithm to optimize the coding matrix.
(1) Initialize the coding matrixC
* with a random ±
1 matrix,
and set the iteration indexi
=
1 and the maximum number of iterations I
max
⁡;
(2)while
i
≤
I
max
⁡ and the coding matrix changes do
(3)  Letr
=
i
mod
⁡
M
L, C
*
r
/
M
+
1
,
r
-
M
r
/
M
=
1
-
C
*
r
/
M
+
1
,
r
-
M
r
/
M,
whereC
*
m
,
n denotes the entry at the mth row and nth column of C
*;
(4)  Generate the EBPSK-based MCPC waveforms
~
i with coding matrix C
* and calculate
z
i
=
s
~
i
2
-
x
2
2
2;             (vi)
(5)  Ifz
i
>
z
i
-
1, i
=
i
+
1, and go to Step  3,
otherwise, letC
*
r
/
M
+
1
,
r
-
M
r
/
M
=
1
-
C
*
r
/
M
+
1
,
r
-
M
r
/
M, i
=
i
+
1 and go to Step  3;
(6)  Leti
=
i
+
1.
(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 Table1, 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, 27–29], 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.Table 1
Simulation parameters.
Parameter
Value
Carrier frequencyf
c
100 MHz
Sample frequencyf
s
1 GHz
The parameter of EBPSK modulationK
2
The parameter of EBPSK modulationN
20
The number of subcarriersL
16
The length of codes of each subcarrierN
C
32
Frequency intervalΔ
f
20 MHz
The SNR of echo signal
[−20, 20] dBFigure3 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 as(40)
e
g
^
k
≜
g
^
k
-
g
k
2
2
g
k
2
2
.We 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.Figure 3
The estimation performance based on MAP.


### 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 Section4.1 except SNR = 5dB, and the index of iteration is from 1 to 50. Figure 4 demonstrates the normalized MSE of estimated TSC, including the approach based on KF and MAP. Apparently, the estimation performance is improved about 25
% 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 g
k
=
e
-
T
/
τ
g
k
-
1
+
u
k
-
1.Figure 4
The estimation of TSC based on KF.Figure5 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 5
The optimized spectrum for the KF estimation.Figure6 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 90
%, and the performance of KF estimation is also improved by 90
%. Moreover, the KF estimation with optimized waveform is 18
% better than the MAP one with optimized waveform.Figure 6
The comparison of estimation performance, including (a) the MAP approach with nonoptimized waveform; (b) the MAP approach with optimized waveform; (c) the KF approach with nonoptimized waveform; (d) the KF approach with optimized waveform.


### 4.3.  The Optimization of EBPSK-Based MCPC Signal

Figure7 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 7
% 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 7
The spectrum of EBPSK-based MCPC signal with optimized power allocation.Figure 8
The comparison of estimation performance by KF approach, including (a) the EBPSK-based MCPC signal with uniform distributed power; (b) the EBPSK-based MCPC signal with optimized power allocation; (c) the optimal signal.Figure9 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 8
% by simultaneously optimizing these two parameters. Furthermore, the algorithm proposed in this work to optimize the coding matrix is easy to realize.Figure 9
Optimize the power and coding matrix of the EBPSK-based MCPC signal. The comparison of estimation performance by KF approach, including (a) the EBPSK-based MCPC signal with optimized power allocation and coding matrix; (b) the EBPSK-based MCPC signal with optimized power allocation; (c) the optimal signal.


## 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 Table1, 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, 27–29], 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.Table 1
Simulation parameters.
Parameter
Value
Carrier frequencyf
c
100 MHz
Sample frequencyf
s
1 GHz
The parameter of EBPSK modulationK
2
The parameter of EBPSK modulationN
20
The number of subcarriersL
16
The length of codes of each subcarrierN
C
32
Frequency intervalΔ
f
20 MHz
The SNR of echo signal
[−20, 20] dBFigure3 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 as(40)
e
g
^
k
≜
g
^
k
-
g
k
2
2
g
k
2
2
.We 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.Figure 3
The estimation performance based on MAP.


## 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 Section4.1 except SNR = 5dB, and the index of iteration is from 1 to 50. Figure 4 demonstrates the normalized MSE of estimated TSC, including the approach based on KF and MAP. Apparently, the estimation performance is improved about 25
% 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 g
k
=
e
-
T
/
τ
g
k
-
1
+
u
k
-
1.Figure 4
The estimation of TSC based on KF.Figure5 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 5
The optimized spectrum for the KF estimation.Figure6 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 90
%, and the performance of KF estimation is also improved by 90
%. Moreover, the KF estimation with optimized waveform is 18
% better than the MAP one with optimized waveform.Figure 6
The comparison of estimation performance, including (a) the MAP approach with nonoptimized waveform; (b) the MAP approach with optimized waveform; (c) the KF approach with nonoptimized waveform; (d) the KF approach with optimized waveform.


## 4.3.  The Optimization of EBPSK-Based MCPC Signal

Figure7 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 7
% 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 7
The spectrum of EBPSK-based MCPC signal with optimized power allocation.Figure 8
The comparison of estimation performance by KF approach, including (a) the EBPSK-based MCPC signal with uniform distributed power; (b) the EBPSK-based MCPC signal with optimized power allocation; (c) the optimal signal.Figure9 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 8
% by simultaneously optimizing these two parameters. Furthermore, the algorithm proposed in this work to optimize the coding matrix is easy to realize.Figure 9
Optimize the power and coding matrix of the EBPSK-based MCPC signal. The comparison of estimation performance by KF approach, including (a) the EBPSK-based MCPC signal with optimized power allocation and coding matrix; (b) the EBPSK-based MCPC signal with optimized power allocation; (c) the optimal signal.


## 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 is25
% better than that of MAP. In addition, the transmitted waveform with optimal spectrum is obtained, which can improve the estimation performance of KF by 90
%. 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 7
% and 8
% 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.



---

*Source: 302083-2015-10-04.xml*