# Sparse Frequency Waveform Design for Radar-Embedded Communication

Authors: Chaoyun Mai; Jinping Sun; Rui Zhou; Guohua Wang
Journal: Mathematical Problems in Engineering
 (2016)

---

## Abstract

According to the Tag application with function of covert communication, a method for sparse frequency waveform design based on radar-embedded communication is proposed. Firstly, sparse frequency waveforms are designed based on power spectral density fitting and quasi-Newton method. Secondly, the eigenvalue decomposition of the sparse frequency waveform sequence is used to get the dominant space. Finally the communication waveforms are designed through the projection of orthogonal pseudorandom vectors in the vertical subspace. Compared with the linear frequency modulation waveform, the sparse frequency waveform can further improve the bandwidth occupation of communication signals, thus achieving higher communication rate. A certain correlation exists between the reciprocally orthogonal communication signals samples and the sparse frequency waveform, which guarantees the low SER (signal error rate) and LPI (low probability of intercept). The simulation results verify the effectiveness of this method.


---

## Body

## 1. Introduction

The purpose of general radar-communication integration design is to achieve high bandwidth data communications between radar platforms by using radar transmitter/receiver subsystem. Different from former methods, Blunt proposed a method of intrapulse radar-embedded communication applied to the Tag [1], which is a new battlefield IFF (Identification Friend or Foe) equipment. The active communication capability between Tag and radar provides further information about the target area, for example, the exact location. Compared with the general radar-communication integration majoring in the sharing ability between radar waveform and communication signals, radar-embedded communication system focuses more on concealment of communication signals. The best way to avoid being detected by enemy devices near the Tag is to embed the Tag signal into radar echo [2]. The basic principle of radar-embedded communications is shown in Figure 1. The Tag device installed in the radar detection zone receives the radar signals, embeds the communication signals by remodulation, and transmits them as radar scattering echoes to radar receiver and achieves covert communication [3, 4]. The communication signals usually choose orthogonal waveforms, which are widely used in radar waveform design.Figure 1
Radar-embedded communications applied to the Tag.Researches of intrapulse radar-embedded communications in [5–9] focus on linear frequency modulation (LFM) pulses. In [7, 8], Ciuonzo proposed a novel waveform design procedure based on multiobjective optimization. The framework of design strategies for intrapulse radar-embedded communications is summarized in [9]. Due to nearly flat spectral characteristic of LFM, the embedded communication signals can only utilize a small range of spectrum band, which leads to limited available communication waveform samples and low efficiency of communication. In view of this problem, a method for sparse frequency waveform design based on radar-embedded communication is proposed. It can improve the band utilization of the communication signal, increase the number of communication samples, and achieve low signal error rate and low probability of intercept while the detecting performance of the radar will not be degraded.Sparse frequency waveform design is an important research direction of the radar waveform design. For the radar system, sparse frequency waveforms have advantages of suppressing interference and improving detection performance. To design a specific sparse frequency waveform, power spectral density fitting can be used to get the objective function [10–15] by using quasi-Newton method and then the eigenvalue decomposition of sparse frequency waveform sequence can be further utilized to obtain the dominant space. Following that K communication waveforms are designed through the projection of orthogonal pseudorandom vectors in the nondominant subspace. Those K communication waveform samples are orthogonal to each other, with strong independence, but have certain relevance with the sparse frequency waveforms. These features ensure that the information would not be easily intercepted by enemy radar during the communication process. Therefore, the use of sparse frequency waveform stopband can effectively increase the frequency band occupation of communication signals, to gain a higher communication rate.


## 2. Radar-Embedded Communication Model

The Tag device in detection range of pulse radar receives radar signals and embeds communication signals by remodulation. Then the output signals are mixed with echo of surrounding scene and return to the radar receiver. The received signals can be expressed as(1)
y
r
t
=
α
k
c
k
t
+
y
s
t
+
n
t
,where n
(
t
) describes system noise, y
s
(
t
) describes ambient radar echoes of the Tag device, α
k describes the combined effect of the path loss, c
k
(
t
) describes communication signals designed, k
=
1,2
,
…
,
K means the kth communication waveform samples, and K is the total number of samples. Rewrite formula (1) in vector form as(2)
y
r
=
α
k
c
k
+
S
x
+
n
,where c
k is the communication signals, n is the noise vector, x is the range samples of radar scattering around Tag device, and S
x is discrete representation of convolution process between the radar waveform and the scattering response. Given waveform vector s
=
s
0
s
1
⋯
s
N
-
1
T, the Toeplitz matrix S can be expressed as(3)
S
=
s
N
-
1
s
N
-
2
⋯
s
0
0
⋯
0
0
s
N
-
1
s
N
-
2
⋯
s
0
⋯
0
⋯
⋯
⋯
⋯
⋯
⋯
⋯
0
0
⋯
s
N
-
1
s
N
-
2
⋯
s
0
.To ensure the concealment (low probability of intercept) of communication signals, the power ofc
k
(
t
) must be lower than that of the echo signal y
s
(
t
). When there is a correlation between y
s
(
t
) and communication signal c
k
(
t
), it is difficult for standard coherent filter to isolate communication samples from the mixed signals. To solve this problem, a design of correlation receiver was proposed in [1].As communication receiving is relatively simple (embedding only one ofK communication samples) and the relative power of the signal is unknown (decided by α
k), a maximum likelihood receiver (known as decorrelator receiver) can be used to realize the communication receiving function. The length of sequence is set to N. The N
×
(
2
N
+
K
-
1
) matrix can be formed as(4)
C
=
S
c
1
⋯
c
K
.Assuming that the noise is white Gaussian N
(
0
,
σ
v
2
), the probability density function of the received signal y
r parameterized with z
=
x
T
α
1
⋯
α
K can be expressed as(5)
p
y
r
;
z
=
1
2
π
σ
v
2
N
/
2
exp
⁡
-
y
r
-
C
z
H
y
r
-
C
z
2
σ
v
2
.The maximum likelihood estimate of parameter z can therefore be obtained:(6)
z
^
=
C
C
-
1
C
H
y
r
.For embedded communication, only the most likely one of the K communication waveform samples is needed to be found. The decorrelator receiver [4] of the kth sample can be expressed as(7)
w
k
=
C
C
H
-
1
c
k
.After obtaining w
1
,
w
2
,
…
,
w
K, correlate each of them with the received signal y
r, respectively, and the maximum likelihood estimate of k can be calculated by(8)
k
^
=
arg
⁡
max
k
⁡
w
k
H
y
r
.The maximum likelihood receiver of formula (7) is the minimum variance unbiased estimator of parameter z for the linear model in formula (2). And formula (8) is the minimum distance receiver for multiple hypothesis tests.


## 3. Sparse Frequency Waveform Design of Embedded Communications

### 3.1. Sparse Frequency Radar Waveform Design

A sparse frequency waveform with specific passband and stopband could be obtained by power spectrum fitting. Given the distribution of power spectral densityu, the objective function of sparse frequency waveform can be expressed as(9)
m
i
n
Θ
A
s
Θ
⊗
A
s
Θ
∗
-
u
2
,where s represents the designed waveform, whose phase vector Θ
=
[
ϕ
1
⋯
ϕ
N
], A is the discrete Fourier transform matrix.The objective function (9) is a quartic nonconvex optimization problem. It can be solved by computing its local minimum rather than the global minimum with the quasi-Newton method [10]. The quasi-Newton method uses information of the objective function value and its first derivative to construct the curvature approximation of the objective function, so it has the advantages of fast convergence. The gradient of the objective function needs to be calculated in each iteration.The objective function (9) can be rewritten as(10)
f
Θ
=
A
s
Θ
⊗
A
s
Θ
∗
-
u
2
=
a
1
H
s
Θ
s
H
Θ
a
1
-
u
1
⋮
a
N
H
s
Θ
s
H
Θ
a
N
-
u
N
f
2
=
∑
m
=
1
N
s
H
Θ
a
m
a
m
H
s
H
Θ
-
u
m
2
=
∑
m
=
1
N
s
H
Θ
A
m
s
H
Θ
-
u
m
2
.The derivative with respect toΘ of the objective function can be expressed as(11)
g
=
∂
f
Θ
∂
Θ
=
2
∑
m
=
1
N
s
Θ
A
m
s
H
Θ
-
u
m
·
2
Im
⁡
diag
⁡
s
Θ
A
m
H
s
H
Θ
.Specific iterative steps of sparse frequency waveform design are as follows.Step 1.
InitializeΘ
0 and tolerance ε, set k
=
0 and S
0
=
I, and then calculate g
0 by formula (11).Step 2.
Set the search directiond
k
=
-
S
k
g
k, calculate f
Θ
k
+
α
k
d
k, and find the optimal step size α
k by line search method so that f
Θ
k
+
α
k
d
k is minimized; let δ
k
=
α
k
d
k and Θ
k
+
1
=
Θ
k
+
δ
k.Step 3.
Ifδ
k
≤
ε, end the iteration and output Θ
optim
=
Θ
k
+
1; if δ
k
>
ε, go to Step 4.Step 4.
Calculateg
k
+
1, set γ
k
=
g
k
+
1
-
g
k, and then compute(12)
S
k
+
1
=
S
k
+
1
+
γ
k
H
S
k
γ
k
γ
k
H
δ
k
δ
k
δ
k
H
γ
k
H
δ
k
-
δ
k
γ
k
H
S
k
+
S
k
γ
k
δ
k
H
γ
k
H
δ
k
.Set k
=
k
+
1 and repeat from Step 2.Set an appropriate value ofε (e.g., 1
0
-
3); then the sparse frequency radar waveform with a distribution close to that of the ideal power spectral density can be obtained by the quasi-Newton method described above.


### 3.2. Waveform Design of Embedded Communications

First of all, the eigenvalue decomposition of the sparse frequency waveform sequence designed is used to get the dominant space. ThenK communication waveform samples are designed through the projection of orthogonal pseudorandom vectors in the nondominant subspace [4].Set the length of sparse frequency radar waveform sequenceN. Known from formulas (1) and (2), radar scattering echo around the Tag within the radar illumination region can be discretely represented as(13)
y
s
=
S
x
,where S is a N
×
(
2
N
-
1
) Toeplitz matrix. x is a 2
N
-
1 vector, representing range samples of radar scattering.After the eigenvalue decomposition ofS, N eigenvalues λ
0
,
λ
1
,
…
,
λ
N
-
1 and corresponding eigenvectors v
0
,
v
1
,
…
,
v
N
-
1 can be calculated. With eigenvalues in ascending order (i.e., λ
0
≤
λ
1
≤
⋯
≤
λ
N
-
1), constitute the matrix V
=
v
0
v
1
⋯
v
N
-
1 accordingly. The eigenvalue decomposition is shown as(14)
S
S
H
=
V
Λ
V
H
,where Λ is a diagonal matrix which contains N associated eigenvalues.A sequence will be obtained if the eigenvalues are sorted in ascending order and it clearly exhibits a demarcation that the large eigenvalues in the sequence correspond to the passband while small eigenvalues correspond to the stopband. Assuming that the space spanned by eigenvectors corresponding to the firstL eigenvalues is the nondominant space, as s
p
a
n
{
V
o
}
=
s
p
a
n
{
v
0
,
v
1
,
…
,
v
L
-
1
}, the space spanned by eigenvectors corresponding to the remaining N
-
L eigenvalues is the dominant space, as s
p
a
n
{
V
s
}
=
s
p
a
n
{
v
L
,
v
L
+
1
,
…
,
v
N
-
1
}. The frequency spectrum of the stopband should be fully utilized to avoid the effect on the performance when designing communication waveform samples. The design steps are as follows.First of all, generate a set ofK pseudorandom N-dimension vectors denoted by b
k, k
=
1,2
,
…
,
K (both Tag and receiver know the specific values of pseudorandom vectors b
k). Then by projecting the orthogonal pseudorandom vector b
k to its nondominant subspace, according to [4], the first communication sample c
1 can be designed as(15)
c
1
=
I
-
V
s
,
0
V
s
,
0
H
b
1
,where V
s
,
0
=
v
L
v
L
+
1
⋯
v
N
-
1. Matrices S and c
1 are combined to obtain a N
×
2
N matrix S
1
=
S
c
1. Do the eigenvalue decomposition again and N eigenvalues λ
0
′
,
λ
1
′
,
…
,
λ
N
-
1
′ of S
1 will be obtained as well as corresponding eigenvectors v
0
′
,
v
1
′
,
…
,
v
N
-
1
′. The eigenvalue decomposition is shown as(16)
S
1
S
1
H
=
V
1
Λ
1
V
1
H
,where V
1
=
v
0
′
v
1
′
⋯
v
N
-
1
′ and Λ
1 is a diagonal matrix which contains N associated eigenvalues. The nondominant space of V
1 is s
p
a
n
{
V
o
,
1
}
=
s
p
a
n
{
v
0
′
,
v
1
′
,
…
,
v
L
-
2
′
} and the dominant space of V
1 is s
p
a
n
{
V
s
,
1
}
=
s
p
a
n
{
v
L
-
1
′
,
v
L
′
,
…
,
v
N
-
1
′
}.Then the second communication sample can be designed as(17)
c
2
=
I
-
V
s
,
1
V
s
,
1
H
b
2
,where V
s
,
1
=
v
L
-
1
′
v
L
′
⋯
v
N
-
1
′. Thus, matrices S and c
1
,
c
2
,
…
,
c
k
-
1 can be combined to obtain an N
×
(
2
N
+
k
-
1
) matrix S
k
-
1
=
S
c
1
c
2
⋯
c
k
-
1. N eigenvalues λ
0
(
k
-
1
)
,
λ
1
(
k
-
1
)
,
…
,
λ
N
-
1
(
k
-
1
) and corresponding eigenvectors v
0
(
k
-
1
)
,
v
1
(
k
-
1
)
,
…
,
v
N
-
1
(
k
-
1
) of S
k
-
1can be further obtained after the eigenvalue decomposition shown as(18)
S
k
-
1
S
k
-
1
H
=
V
k
-
1
Λ
k
-
1
V
k
-
1
H
,where V
k
-
1
=
v
0
k
-
1
v
1
k
-
1
⋯
v
N
-
1
k
-
1 and Λ
1 is a diagonal matrix which contains N associated eigenvalues. The nondominant space of V
k
-
1 is s
p
a
n
{
V
o
,
k
-
1
}
=
s
p
a
n
{
v
0
(
k
-
1
)
,
v
1
(
k
-
1
)
,
…
,
v
L
-
k
(
k
-
1
)
} and the dominant space of V
k
-
1 is s
p
a
n
{
V
o
,
k
-
1
}
=
s
p
a
n
{
v
0
(
k
-
1
)
,
v
1
(
k
-
1
)
,
…
,
v
L
-
k
(
k
-
1
)
}. Then the kth communication sample can be designed as(19)
c
k
=
I
-
V
s
,
k
-
1
V
s
,
k
-
1
H
b
k
,where V
s
,
k
-
1
=
v
L
-
k
+
1
k
-
1
v
L
-
k
+
2
k
-
1
⋯
v
N
-
1
k
-
1.An important indicator of covert communication is interception [3]. Since the conventional method based on measuring the spectral energy is no longer applicable, Blunt et al. proposed the normalized correlation to measure interception [4]. The intended radar waveform is assumed to be known by the enemy intercept receiver, and the dominant space of waveform is obtained through the eigenvalue decomposition by formula (14). The N
×
l matrix V
~
s
,
l is formed by the eigenvectors corresponding to the largest l eigenvalues. The projection of received signal in the nondominant subspace of the main space can be considered as nonradar echo component. The lth projected value of the received signal y
r is(20)
a
l
=
P
l
y
r
,where P
l
=
I
-
V
~
s
,
l
V
~
s
,
l
H
,
l
∈
1
,
…
,
N is the projection matrix corresponding to different l.Then, normalized correlation of thekth communication waveform is defined as(21)
η
k
,
l
=
c
k
H
a
l
c
k
H
c
k
a
l
H
a
l
.
η
k
,
l, having a range of 0
1, can be used to measure the intercept probability of communications waveforms [4]. The larger value of η
k
,
l indicates higher similarity between the communication signal detected by the intercept receiver and the actual embedded signal, which means higher probability to be intercepted. Communication waveforms are thought to have no covertness when η
k
,
l approaches 1.


## 3.1. Sparse Frequency Radar Waveform Design

A sparse frequency waveform with specific passband and stopband could be obtained by power spectrum fitting. Given the distribution of power spectral densityu, the objective function of sparse frequency waveform can be expressed as(9)
m
i
n
Θ
A
s
Θ
⊗
A
s
Θ
∗
-
u
2
,where s represents the designed waveform, whose phase vector Θ
=
[
ϕ
1
⋯
ϕ
N
], A is the discrete Fourier transform matrix.The objective function (9) is a quartic nonconvex optimization problem. It can be solved by computing its local minimum rather than the global minimum with the quasi-Newton method [10]. The quasi-Newton method uses information of the objective function value and its first derivative to construct the curvature approximation of the objective function, so it has the advantages of fast convergence. The gradient of the objective function needs to be calculated in each iteration.The objective function (9) can be rewritten as(10)
f
Θ
=
A
s
Θ
⊗
A
s
Θ
∗
-
u
2
=
a
1
H
s
Θ
s
H
Θ
a
1
-
u
1
⋮
a
N
H
s
Θ
s
H
Θ
a
N
-
u
N
f
2
=
∑
m
=
1
N
s
H
Θ
a
m
a
m
H
s
H
Θ
-
u
m
2
=
∑
m
=
1
N
s
H
Θ
A
m
s
H
Θ
-
u
m
2
.The derivative with respect toΘ of the objective function can be expressed as(11)
g
=
∂
f
Θ
∂
Θ
=
2
∑
m
=
1
N
s
Θ
A
m
s
H
Θ
-
u
m
·
2
Im
⁡
diag
⁡
s
Θ
A
m
H
s
H
Θ
.Specific iterative steps of sparse frequency waveform design are as follows.Step 1.
InitializeΘ
0 and tolerance ε, set k
=
0 and S
0
=
I, and then calculate g
0 by formula (11).Step 2.
Set the search directiond
k
=
-
S
k
g
k, calculate f
Θ
k
+
α
k
d
k, and find the optimal step size α
k by line search method so that f
Θ
k
+
α
k
d
k is minimized; let δ
k
=
α
k
d
k and Θ
k
+
1
=
Θ
k
+
δ
k.Step 3.
Ifδ
k
≤
ε, end the iteration and output Θ
optim
=
Θ
k
+
1; if δ
k
>
ε, go to Step 4.Step 4.
Calculateg
k
+
1, set γ
k
=
g
k
+
1
-
g
k, and then compute(12)
S
k
+
1
=
S
k
+
1
+
γ
k
H
S
k
γ
k
γ
k
H
δ
k
δ
k
δ
k
H
γ
k
H
δ
k
-
δ
k
γ
k
H
S
k
+
S
k
γ
k
δ
k
H
γ
k
H
δ
k
.Set k
=
k
+
1 and repeat from Step 2.Set an appropriate value ofε (e.g., 1
0
-
3); then the sparse frequency radar waveform with a distribution close to that of the ideal power spectral density can be obtained by the quasi-Newton method described above.


## 3.2. Waveform Design of Embedded Communications

First of all, the eigenvalue decomposition of the sparse frequency waveform sequence designed is used to get the dominant space. ThenK communication waveform samples are designed through the projection of orthogonal pseudorandom vectors in the nondominant subspace [4].Set the length of sparse frequency radar waveform sequenceN. Known from formulas (1) and (2), radar scattering echo around the Tag within the radar illumination region can be discretely represented as(13)
y
s
=
S
x
,where S is a N
×
(
2
N
-
1
) Toeplitz matrix. x is a 2
N
-
1 vector, representing range samples of radar scattering.After the eigenvalue decomposition ofS, N eigenvalues λ
0
,
λ
1
,
…
,
λ
N
-
1 and corresponding eigenvectors v
0
,
v
1
,
…
,
v
N
-
1 can be calculated. With eigenvalues in ascending order (i.e., λ
0
≤
λ
1
≤
⋯
≤
λ
N
-
1), constitute the matrix V
=
v
0
v
1
⋯
v
N
-
1 accordingly. The eigenvalue decomposition is shown as(14)
S
S
H
=
V
Λ
V
H
,where Λ is a diagonal matrix which contains N associated eigenvalues.A sequence will be obtained if the eigenvalues are sorted in ascending order and it clearly exhibits a demarcation that the large eigenvalues in the sequence correspond to the passband while small eigenvalues correspond to the stopband. Assuming that the space spanned by eigenvectors corresponding to the firstL eigenvalues is the nondominant space, as s
p
a
n
{
V
o
}
=
s
p
a
n
{
v
0
,
v
1
,
…
,
v
L
-
1
}, the space spanned by eigenvectors corresponding to the remaining N
-
L eigenvalues is the dominant space, as s
p
a
n
{
V
s
}
=
s
p
a
n
{
v
L
,
v
L
+
1
,
…
,
v
N
-
1
}. The frequency spectrum of the stopband should be fully utilized to avoid the effect on the performance when designing communication waveform samples. The design steps are as follows.First of all, generate a set ofK pseudorandom N-dimension vectors denoted by b
k, k
=
1,2
,
…
,
K (both Tag and receiver know the specific values of pseudorandom vectors b
k). Then by projecting the orthogonal pseudorandom vector b
k to its nondominant subspace, according to [4], the first communication sample c
1 can be designed as(15)
c
1
=
I
-
V
s
,
0
V
s
,
0
H
b
1
,where V
s
,
0
=
v
L
v
L
+
1
⋯
v
N
-
1. Matrices S and c
1 are combined to obtain a N
×
2
N matrix S
1
=
S
c
1. Do the eigenvalue decomposition again and N eigenvalues λ
0
′
,
λ
1
′
,
…
,
λ
N
-
1
′ of S
1 will be obtained as well as corresponding eigenvectors v
0
′
,
v
1
′
,
…
,
v
N
-
1
′. The eigenvalue decomposition is shown as(16)
S
1
S
1
H
=
V
1
Λ
1
V
1
H
,where V
1
=
v
0
′
v
1
′
⋯
v
N
-
1
′ and Λ
1 is a diagonal matrix which contains N associated eigenvalues. The nondominant space of V
1 is s
p
a
n
{
V
o
,
1
}
=
s
p
a
n
{
v
0
′
,
v
1
′
,
…
,
v
L
-
2
′
} and the dominant space of V
1 is s
p
a
n
{
V
s
,
1
}
=
s
p
a
n
{
v
L
-
1
′
,
v
L
′
,
…
,
v
N
-
1
′
}.Then the second communication sample can be designed as(17)
c
2
=
I
-
V
s
,
1
V
s
,
1
H
b
2
,where V
s
,
1
=
v
L
-
1
′
v
L
′
⋯
v
N
-
1
′. Thus, matrices S and c
1
,
c
2
,
…
,
c
k
-
1 can be combined to obtain an N
×
(
2
N
+
k
-
1
) matrix S
k
-
1
=
S
c
1
c
2
⋯
c
k
-
1. N eigenvalues λ
0
(
k
-
1
)
,
λ
1
(
k
-
1
)
,
…
,
λ
N
-
1
(
k
-
1
) and corresponding eigenvectors v
0
(
k
-
1
)
,
v
1
(
k
-
1
)
,
…
,
v
N
-
1
(
k
-
1
) of S
k
-
1can be further obtained after the eigenvalue decomposition shown as(18)
S
k
-
1
S
k
-
1
H
=
V
k
-
1
Λ
k
-
1
V
k
-
1
H
,where V
k
-
1
=
v
0
k
-
1
v
1
k
-
1
⋯
v
N
-
1
k
-
1 and Λ
1 is a diagonal matrix which contains N associated eigenvalues. The nondominant space of V
k
-
1 is s
p
a
n
{
V
o
,
k
-
1
}
=
s
p
a
n
{
v
0
(
k
-
1
)
,
v
1
(
k
-
1
)
,
…
,
v
L
-
k
(
k
-
1
)
} and the dominant space of V
k
-
1 is s
p
a
n
{
V
o
,
k
-
1
}
=
s
p
a
n
{
v
0
(
k
-
1
)
,
v
1
(
k
-
1
)
,
…
,
v
L
-
k
(
k
-
1
)
}. Then the kth communication sample can be designed as(19)
c
k
=
I
-
V
s
,
k
-
1
V
s
,
k
-
1
H
b
k
,where V
s
,
k
-
1
=
v
L
-
k
+
1
k
-
1
v
L
-
k
+
2
k
-
1
⋯
v
N
-
1
k
-
1.An important indicator of covert communication is interception [3]. Since the conventional method based on measuring the spectral energy is no longer applicable, Blunt et al. proposed the normalized correlation to measure interception [4]. The intended radar waveform is assumed to be known by the enemy intercept receiver, and the dominant space of waveform is obtained through the eigenvalue decomposition by formula (14). The N
×
l matrix V
~
s
,
l is formed by the eigenvectors corresponding to the largest l eigenvalues. The projection of received signal in the nondominant subspace of the main space can be considered as nonradar echo component. The lth projected value of the received signal y
r is(20)
a
l
=
P
l
y
r
,where P
l
=
I
-
V
~
s
,
l
V
~
s
,
l
H
,
l
∈
1
,
…
,
N is the projection matrix corresponding to different l.Then, normalized correlation of thekth communication waveform is defined as(21)
η
k
,
l
=
c
k
H
a
l
c
k
H
c
k
a
l
H
a
l
.
η
k
,
l, having a range of 0
1, can be used to measure the intercept probability of communications waveforms [4]. The larger value of η
k
,
l indicates higher similarity between the communication signal detected by the intercept receiver and the actual embedded signal, which means higher probability to be intercepted. Communication waveforms are thought to have no covertness when η
k
,
l approaches 1.


## 4. Experimental Results

Set the length of the sparse frequency radar waveform sequenceN
=
200. The desired power spectral density u has 4 passbands and 3 stopbands. The set of passbands is 5 db and stopbands −30 dB. ε is presumed to be 0.001. According to the iterative steps from Section 3.1, the obtained power spectral density of sparse frequency waveform is shown in Figure 2.Figure 2
Power spectral density of radar waveform after communication signalc
1 embedded.As the sparse frequency waveforms have been obtained, communication waveforms can be designed through the method discussed in Section3.2. Set the total number of communication samples K
=
16 and compute c
1
,
c
2
,
…
,
c
16. Figure 2 shows the power spectral density of radar waveform after the communication signal c
1 is embedded; that leads to the increased power in the stopband obviously. We can also see the change for the radar waveform in time domain when the communication signal c
1 is embedded. The result waveforms in time domain before and after c
1 embedded are shown in Figure 3.Figure 3
Sparse frequency radar waveform in time domain. (a) The real part beforec
1 embedded; (b) the real part after c
1 embedded; (c) the imaginary part before c
1 embedded; (d) the imaginary part after c
1 embedded.
(a)
(b)
(c)
(d)Figures3(a) and 3(c) show the real part and the imaginary part of the original sparse frequency radar waveform, which is a multiphase coded waveform. It is constant modulus, and the modulus value is normalized. Figures 3(b) and 3(d) show the real part and the imaginary part of the waveform after the communication signal c
1 is embedded, and it can be seen that the difference is not obvious in time domain. For the case that c
2
–
c
16 embedded, the results are similar.After obtaining the radar-embedded communication waveforms, we can further analyze the SER with the variation of SNR. This design method containsK
=
16 kinds of communication waveforms and the first one as an example is adopted in simulation (the probability is evaluated when the intended receiver successfully detects the communication signal c
1 transmitted by transmitter). Change the power of noise to vary the value of SNR from −25 dB to 0 dB. Set the number of nondominant eigenvectors used by Tag and radar receiver L
=
80, and the number of large eigenvalues used by intercept receiver j
=
120. The result averaged over 10000 Monte Carlo simulations is shown in Figure 4.Figure 4
Error rate after communication signalc
1 embedded.As is shown, when the SNR is −8 dB, SER can achieve 10−3 by using a decorrelator receiver, far less than that of the intercept receiver. It demonstrates that the communication waveform has good concealment. Figure 4 only shows the result of c
1. The results of samples of c
2
–
c
16 are similar to that of c
1. When SNR ranges from −10 dB to −5 dB, the decorrelator receiver has lower SER.The interception of communication waveforms is further analyzed. Set the SNR to −15 dB. While calculating the normalized correlation, the number of large eigenvalues of the intercept receiverl varies from 0 to 200 (i.e., change from 0–100%, because N
=
200). Each point uses 200 samples to calculate η
k
,
l and get the average. Results are shown in Figure 5.Figure 5
Normalized correlation curve by computation of waveformc
1.In Figure5, the horizontal axis represents the percentage of the number of eigenvalues l; the vertical axis represents the value of η
k
,
l. According to formula (20), the calculation of normalized correlation needs to know one of the communication waveforms. Supposing the communication sample c
1 is known, compute the normalized correlation of c
2
,
c
3
,
…
,
c
16 and also that of c
1, for generality. It can be seen from Figure 5 that the maximum of η
k
,
l corresponding to c
2
–
c
16 is less than 0.2 and the maximum of η
k
,
l corresponding to c
1 is approximately equal to 0.3. Even increasing the percentage of eigenvalues, the waveform that is designed by this method can still maintain a good ability of interception. Simulation in Figure 5 only shows the result provided that the communication sample c
1 is known. In the actual operation, using the communication samples c
2
–
c
16 will get the similar curve and the same conclusion.The impact of embedded communication signals into radar signal by autocorrelation function is further analyzed. Figure6 shows normalized autocorrelation function (ACF) of radar waveform [16] before and after the waveform embedded communication signals, with ordinate in logarithm form. Figure 6(a) shows the ACF under original sparse frequency radar waveform, and Figure 6(b) shows the ACF after adding the communication waveform, which would be c
1. The integral sidelobe level (ISL) is −10.19 dB and peak sidelobe level (PSL) is −15.62 dB before communication waveform is embedded. And the ISL is −10.94 dB and PSL is −16.44 dB when communication waveform is embedded.Figure 6
ACF of sparse frequency radar-embedded communication waveform. (a) Beforec
1 embedded; (b) after c
1 embedded.
(a)
(b)It can be seen that the autocorrelation performances are similar before or after the communication waveformc
1 is embedded, which means the embedded communication waveforms do not affect the detection performance of the radar system. For c
2
–
c
16, the results are basically the same.


## 5. Conclusions

In this paper, we propose the sparse frequency waveform design based on radar-embedded communication to improve the performance of the existing radar Tag system with the function of covert communication. Firstly, sparse frequency waveforms are designed based on power spectral density fitting with the quasi-Newton method and the eigenvalue decomposition is used to get the dominant space of the waveform sequence. Secondly, the communication waveforms are designed through the projection of orthogonal pseudorandom vectors in the nondominant subspace. In the end, we analyze the error rate and the interception of sparse frequency waveform in which communication waveform samples are embedded under different SNR conditions in simulation. The result shows that using sparse frequency waveform can improve the bandwidth occupation of communication signals. Thus, a higher communication rate can be gained. Furthermore, these samples are orthogonal to each other, which guarantees low SER and LPI. It is also verified that the autocorrelation performance of radar signals would not be degraded by embedding communication signals.



---

*Source: 7270301-2016-06-15.xml*