1. Introduction
With the increasing types and quantities of targets near space, including aircraft, ballistic missiles, decoys, satellites, near space target, and space debris, inverse synthetic aperture radar (ISAR) plays an increasingly important role in target imaging and identification [
1,
2]. At the same time, the electronic countermeasures (ECM) technology against ISAR is developing so rapidly with the development of digital radio frequency memory (DRFM). In particular, the interrupted sampling repeater jamming (ISRJ) method based on DRFM can realize deceiving the victim radar by generating amounts of realistic false targets [
3,
4]. The response speed of ISRJ is fast, and the jamming signal is correlated with the target signal [
5,
6,
7]. Therefore, extensive research of ISRJ has been carried out due to its significant application in electronic warfare.
ISRJ is a novel jamming method against ISAR by interrupting sampling the target echo signal and retransmitting the intercepted segments. The improved ISRJ inherits the traditional “storage-repeater” model and adjusts the sampling interval and retransmission period to conquer the deficiencies of the amplitude decaying rapidly and main false target hysteresis [
8,
9,
10]. Besides, the sub-Nyquist sampling technology is another kind of improved ISRJ method whose sampling rate of the intercepted target echo signal is less than the Nyquist sampling rate [
11]. The false targets are induced along with the different orientations, respectively, by interrupted sampling the echo signal along with the corresponding different time domain at sub-Nyquist sampling rate in [
12,
13,
14]. In [
14], multiple false targets will be induced along with the range orientation while the CS sampling number should satisfy inequation condition with sparsity. By introducing the sub-Nyquist sampling and transforming the convolution into a parallel formation, the structure for generating the jamming signal is simplified, and the computation complexity is reduced in [
7]. The jamming methods mentioned above are mainly focused on one-dimensional jamming. In the field of radar target feature extraction, ISAR imaging using a broadband signal will still be interfered by the ISRJ signal. However, there is a lack of work to study the characteristics of ISRJ jamming signals from the two dimensions of ISAR imaging. Further, this paper investigated the 2D jamming method that the ISRJ imposed both in the fast-time domain and slow-time domain and derived the explicit signal form.
Meanwhile, many anti-jamming methods against ISRJ have attracted much attention. In [
15], the jamming signal parameter is estimated firstly, and then an adaptive parameter-adjusted intra-pulse frequency coded signal is transmitted to realize ISRJ suppression. In [
16], the jamming signal parameters of intercepted slice number, forwarding times, and slice width was estimated firstly, and then the parameters were used for iterative cancellation of suppression. However, these methods’ performances depend on the precision of parameter estimating. Based on the discontinuous characteristics of the ISRJ signal, [
17] proposed a bandpass filter based on the discontinuous pulse duration. A novel method was proposed based on the radar signal energy function, and the results of the function were compared with the threshold to obtain the non-disturbed segments in [
18]. According to the energy distribution of the de-chirped signal, a function was constructed to extract the jamming-free signal segments, and the bandpass filter based on these segments was generated in [
19]. The filter results are suitable for constant false alarm rate (CFAR) detectors to improve the suppression performance. Nevertheless, the filter method is constructed complicatedly, and it is easily affected by the low signal-to-noise ratios (SNR) or high jam to signal ratio (JSR).
CS theory can reconstruct the sparse data under limited measurements [
20,
21]. The paper [
22] proposed a phase-aided distributed compressive sensing (DCS) method to preserve the target and suppress randomly distributed jamming signals. In [
23], a CS-based method for high-resolution ISAR images is proposed under limited measurements. These methods in the aforementioned five papers solve the reconstruction problem based on Compressed Sensing. However, the anti-jamming performance of low SNR and high JSR has not been analyzed.
Since the radar echo signal and ISRJ signal are sparse, the CS method can suppress ISRJ signals [
24,
25,
26]. An imaging method with a non-periodic interrupted sampling-linear frequency modulated signal (NIS-LFM) is proposed in [
24]. The NIS-LFM is piecewise sparse, and the target echo can be extracted based on CS theory, while parameters analyses on the anti-jamming performance are limited to restricted isometry property (RIP) conditions. The high-resolution range profile (HRRP) is reconstructed, and this method’s minimum sub-Nyquist sampling rate countering jamming is addressed in [
25]. According to the variable sparsity influence, the reconstruction effect, a Bayesian compress sensing-based method is proposed to reconstruct the multiple targets signal by optimizing the sparse model in [
26]. This model is based on the unjammed discontinuous segments extracted by the ISRJ strategy already. However, only the range-dimension anti-jamming method is discussed in the above methods, and parameter analysis of CS-based reconstruction is incomplete.
Most of the methods against ISRJ focus on parameter estimation or filter principle, while the anti-jamming method of parameter estimation requires a higher real-time response of hardware. The filter-based anti-jamming method mainly handles the segments without jamming. There is a lack of research on processing directly the radar echo containing jamming signal by means of an algorithm. At present, the available CS-based methods against ISRJ are to filter out the signal without jamming first and then reconstruct the signal with CS theory only using the compressed sensing method to reconstruct the signal. Inspired by the promising CS theory, this paper proposed an adaptive anti-jamming method to suppress the ISRJ signal by adjusting vector sparsity. According to the restricted isometry property (RIP) condition of CS theory, the signal reconstruction effect is relative to the number of the measurements, which varies with the sparsity. In addition, the ISRJ signal is piecewise sparse in the frequency-domain, and the reconstruction sparsity of the target signal is different from that of the ISRJ signal. Hence, the ISRJ signal can be suppressed primarily by searching the number of optimal measurements. This paper first derivates the ISRJ signal model along with range orientation and azimuth orientation and then investigates the 2D anti-jamming method. This method can realize anti-jamming with only a few signal measurements at the receiver and does not need complicated parameter estimation and filter design, which is applicable in various scenarios.
The structure of this paper is organized as follows. In
Section 2, a 2D ISRJ signal model is built. In
Section 3, the procedure of the 2D anti-jamming method is derived in detail. In
Section 4, numerous simulations are carried out. Finally, in
Section 5, conclusions are drawn.
2. Signal Model
The geometry of the jamming configuration is illustrated in
Figure 1. The jammer is assumed on the same platform as the radar, and this jamming belongs to main lobe interference. The process of radar receiving jamming signals can be divided into two stages: Firstly, the radar transmits the linear frequency modulation (LFM) wideband waveform, some of which are intermittently intercepted by the jammer and then transmitted to the target. Secondly, the jamming signals are received by the radar receiver, and then multiple false targets are induced in ISAR image.
In
Figure 1, the coordinate
XOY is established on the point
, which denotes the reference point of the target. The x-axis is the line of the radar transmitting signal.
,
and
represents the distance between jammer and reference point, radar and reference point, radar, and jammer, respectively.
is the incident angle of the jammer. Assume that the scatter
of the rigid object rotates with speed
and the instantaneous rotation angle of the scatter
during the observation period is
, which can be expressed as
with slow time
. Hence, the distance history of ISRJ signal is accumulated after sampling by jammer and reflecting on the scatter
of target, which can be expressed below
where
represents the displacement distance of scatter along with the range orientation, while
represents the displacement distance along with the azimuth orientation. The formula transformation of the ISRJ signal during the two stages is described below. Assume that the signal transmitted by radar is as below
where
is the width of the pulse,
is the carrier frequency and
is the chirp rate.
represents the fast time while
is the slow time and
is the full time. The function
can be described as
At the jammer, the sampling function in the fast time and in the slow time are, respectively, defined as
where
represents sampling pulse duration and
is the sampling interval in the fast time. Then the equivalent
represents the duty ratio, where
is the sampling frequency.
is the sampling interval in the slow time domain and
is the sampling frequency. This sampling function makes the jamming signal equivalent to the target signal weighted with different Doppler shifts.
The ISRJ signal in the range dimension can be generated by interrupted transmitting the target signal in the fast time domain. The generation processes of ISRJ signal in azimuth orientation are the same as those in range dimension, except that intermittent sampling signal is carried out in the slow time domain with
. In this section, we take the range dimension as an example. The jamming signal can be expressed as
where
is the scattering coefficient of the
point and
is the electromagnetic propagation velocity. Assume the history distance of the reference point (i.e.,
point) is
. Then the jamming signal of reference point can be expressed as below
Cast
into Equation (5) and de-chirp it with the reference signal, which reduces the bandwidth. Then the jamming signal can be expressed as
where
When
is a slight angle, there is an approximation of Equation (8) as
Omitting the constant term and calculating the Fourier transform of
in Equation (7), then the ISRJ signal after range compression can be expressed as
Based on the frequency-distance relationship of de-chirping,
, the jamming signal is equivalent to
where
represents the range unit and
is the bandwidth of LFM signal. Actually, intermittent sampling in the fast time domain is equivalent to applying ISRJ jamming along the range orientation, and the ISAR image of false target can be obtained along the range orientation with ISRJ. Equation (11) illustrates that the jamming signal consists of a few sinc functions, which represent the strong scattering centers. This equation describes multiple peaks of the sinc function along the range orientation at
, which represents the results of the ISRJ signal with the function
. Assuming the radar receives the target echo with the same parameters, the target echo signal can also be processed along the cross-range direction with an interrupted sampling repeater.
The ISRJ signal in the azimuth dimension is similar with that in the range dimension that the jamming signal can be generated by interrupted transmitting the radar signal in the slow time domain. Multiplying
by Equation (11), the jamming signal along two dimensions can be expressed as
Casting Equations (9)–(12) and applying the Fourier transform in slow time, the jamming signal can be obtained. The distance variable along range orientation can be omitted when the variation of the Doppler Domain is discussed. The jamming signal can be expressed as
where
is the number of pulse strings and
is the pulse repetition interval. Equation (13) shows that peaks of the
sinc function along the azimuth orientation are also located at
after pulse compression.
According to ISAR principle, the ISAR image of the target is obtained along with the range orientation and azimuth orientation [
27]. Thus, the ISRJ signal will induce a set of false-targets symmetrically centered on the real target along two dimensions of the RD image. The range interval of the false target is
, and the azimuth interval of the false target is
. The number and amplitude of false-target are related to jamming sampling interval and duty cycle. Even though the number of targets on the RD image are increased after interrupted sampling repeater jamming, the false target on the RD can also be deemed to be sparse.
3. The Theory of Adaptive CS-Based Method against ISRJ
According to the basic idea of the CS, the sparse signal can be reconstructed by a few measurements. A few scattered points of the target contained the most energy of target echo, thus the target echo signal is discrete and sparse. In that case, both the target echo and ISRJ signal can be reconstructed based on CS theory. In the following part, the ISAR imaging method based on CS theory is briefly introduced firstly. Then it is the CS reconstruction of the ISRJ signal. The last subsection describes and analyzes the anti-jamming procedure of the CS-based method.
3.1. ISAR Imaging via CS
The CS imaging of ISAR includes two steps. The first step is the measured process. The ISAR measured data
, which is accepted by the radar receiver, can be expressed as
where
is the measurement matrix with
and
denotes the measured noise. Assume
represents the interested original signal. Generally,
can be represented as
where
is the orthonormal base matrix and
is a sparse weighting coefficients vector. If there are only
nonzero coefficients for
at most with
, then the signal
is called owning
sparsity in the base matrix
. Actually,
and
are two representations of the same signal that
is the complete signal representation, while
is the sparse representation of
in the sparse domain
. Substituting Equations (15) to (14) and assuming that
is the dictionary matrix, which is represented as
, we have
The second step is signal recovery. The sparse weighting coefficients vector can be obtained from limited measurement by calculating the convex optimization problem below
where
stands for
norm and
represents the minimization,
represents the coefficient vector to calculate, and
is the noise level. Equation (16) will have a solution in the condition of matrix
obeying the RIP condition [
20]. Then the ISAR image can be obtained by reconstructing pulse signal firstly and range unit signal secondly based on CS.
3.2. ISRJ Signal Imaging via CS
The ISRJ signal can be deemed sparse because the multiple false-targets consist of a few strong scattering centers in ISAR imaging. Thus, the ISRJ jamming signal can be collected and recovered with CS theory. This section analyzes the CS-based anti-jamming model of 2D ISRJ jamming signal. The anti-jamming process is divided into two the 2D CS reconstruction.
The first step is countering the range orientation jamming by reconstructing the jamming signal with the CS method. The down-range orientation of the ISRJ signal is applied with de-chirping and Fourier transform, then the ISRJ signal model with 2D interrupted sampling can be simplified from Equation (15) as
where
. Assuming that the observation region along range orientation is
, the
represents the range resolution, and the number of bins along range orientation is
. In the ISAR image domain, the total number of scatters is
. Then the observation matrix can be constructed as
Suppose the number of bins along azimuth orientation is
, then the dimension of the matrix
is
. The scattering echo coefficient vector of one pulse is
where
represents the n-th backscattering coefficient along the range orientation. Moreover,
means there is no target in the range units. If the false targets occupy only a small part of the observation area, the number of nonzero elements in
is much smaller than the number of distance units
. The de-chirping echo signal of the jamming signal can be expressed as
where the
is the additive noise and
is
-sparse vector, which can be estimated by calculating the convex optimization problem below
where
represents one pulse echo signal.
The second step is the ISAR azimuth processing of the ISRJ signal based on CS. After range alignment and phase compensation, the ISRJ signal in one range unit after compression can be expressed as
where
represents the i-th backscattering coefficient along with the azimuth orientation. The reflection coefficient of the background is zero. Substituting Equation (9) into Equation (23) and ignoring the RVP, the jamming signal model of ISAR image in one range unit after compression can be represented as
Assuming that the observation region along azimuth orientation is
, the
represents the azimuth resolution, and the number of bins along range orientation is
. Then the observation matrix of jamming signal in one range unit can be described as
where
is the observation matrix with
, which is different from the matrix
. Hence, the vectors of scatters along azimuth dimension can be obtained by calculating the optimization problem below.
3.3. Anti-Jamming Procedure on CS Theory
According to the theory of CS-based ISAR imaging, the sparsity actually represents the scattered number of targets. Hence, the target can be reconstructed with the optimum sparsity on the condition of RIP, which is
where
is the number of the measurements,
is the sparsity of vector, and
is a small constant.
represents the dimension of the coefficient vector, which equals the number of distance units
and
respectively in
Section 3.2. Assuming the minimum number of measurements satisfies the equation
, the signal can be reconstructed as long as the measurements are bigger than
. The relationship between sparsity and measurements number is proved and simulated in [
22]
The probability of reconstruction increases with adding the measurements under invariant sparsity;
With the increasing sparsity, the reconstruction is guaranteed only by more measurements.
This conclusion is significant for suppressing the ISRJ signal below.
Section 2 describes that the ISRJ signal is generated by periodically sampling and forwarding the original target signal. Therefore, the ISRJ signal can induce multiple false targets with high similarity in the frequency domain.
Figure 2 shows the CS-based anti-jamming method against ISRJ, and those characteristics of the ISRJ signal are revealed in it. The ISAR imaging process of the ISRJ signal is divided into two steps.
Firstly,
points are randomly sampled along with the range orientation for each echo signal, that is, along the direction of the yellow arrow in
Figure 2a. Those red points represent the value that is sampled. Each echo’s high-resolution range profile (HRRP) can be obtained by reconstructing the vector based on the proposed CS method.
Next, different reconstructed echo values in the range unit are randomly sampled
points along the azimuth direction, which is the direction of the green arrow in
Figure 2b. Those green blocks indicate the point being sampled. Then every column vector is reconstructed by the proposed CS method, and the imaging results of the ISRJ signal can be obtained.
Figure 2c shows the dense false targets induced by the ISRJ signal, in which the blue aircraft represents the true target, and the red aircraft represents the false targets. Multiple false targets with high similarity are distributed as an equal interval. In [
4], the main false target was in the middle of all false targets, and the amplitude of the main false target was stronger than that of other false targets. Besides, the main false target overlaps with the true target. These characteristics of ISRJ can be utilized for anti-jamming. The proposed CS method obtains K-sparse optimal solutions, which can express the most energy of the observation vector by solving optimization problems. Since the ISRJ signal is sparse in ISAR image domain, it can be reconstructed by CS method only with a few random measurements.
Figure 2 exhibits the scattering number of the ISRJ signal is more than those of the target signal without jamming. Hence, the sparsity of jamming signal
is bigger than that of target signal
. Then the condition of reconstructing jamming signal can be represented from Equation (20)
Let
indicate the minimum measurements needed to reconstruct the jamming signal, then
. There is a minimum anti-sampling number
, which satisfies the equality
. The amplitude of multiple false targets induced by the ISRJ signal is different. The energy of the main false target in the middle is the largest, and the energy of other false targets decreases in turn. Furthermore, the main false target overlaps with the real target. Therefore, the jamming signal can be suppressed with the number of the optimal measurements
, by which the target can be reconstructed while the jamming signal cannot be reconstructed. Moreover, the false middle target and real target are reconstructed firstly with the largest energy because the energy is the largest. The one-dimensional procedure of the algorithm is shown in
Figure 3. The 2-D imagery repeats this algorithm along with range and azimuth orientation.
In this flowchart, parameters of system are initialized firstly, and a threshold is set up. is an empirical constant calculated from a distance between the target and the radar and the transmitting power of the jammer. Secondly, the number of measurements is calculated from , which can ensure complete reconstruction of the vector . Thirdly, the vector is reconstructed by the proposed CS method. It can be judged whether the ISRJ signal has been removed by the inequality . The overlapping of the main false targets and the true target will be reconstructed, which is mentioned earlier. Therefore, the energy of the jamming signal can be filtered by adjusting the number of measurements . The energy of the reconstructed signal also includes the energy of the ISRJ signal when . Thus, the energy of sampled signal should be reduced by reducing the number of measurements (i.e., ). Until , the overlapping of true target and main false target will be reconstructed.