Next Article in Journal
Analysis on BDS-3 Autonomous Navigation Performance Based on the LEO Constellation and Regional Stations
Next Article in Special Issue
Gaussian Process Gaussian Mixture PHD Filter for 3D Multiple Extended Target Tracking
Previous Article in Journal
Comprehensive Ocean Information-Enabled AUV Motion Planning Based on Reinforcement Learning
Previous Article in Special Issue
An Optimization Method for Collaborative Radar Antijamming Based on Multi-Agent Reinforcement Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Anti-Intermittent Sampling Jamming Technique Utilizing the OTSU Algorithm of Random Orthogonal Sub-Pulses

School of Electronics and Communication Engineering, Sun Yat-Sen University, Shenzhen 518000, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(12), 3080; https://doi.org/10.3390/rs15123080
Submission received: 10 April 2023 / Revised: 5 June 2023 / Accepted: 9 June 2023 / Published: 12 June 2023
(This article belongs to the Special Issue Advanced Radar Signal Processing and Applications)

Abstract

:
The utilization of intermittent sampling jamming can engender a lofty verisimilitude false target cluster that exhibits coherence with the transmitted signal. Such an assemblage bears the hallmarks of both suppression jamming and deceitful jamming, capable of inflicting substantial impairment upon the radar, potentially leading to its profound incapacitation. Henceforth, the precise discernment of the target and various forms of intermittent sampling jamming emerges as a novel endeavor. In response to this predicament, this paper posits a pulsed radar waveform featuring intra-pulse random orthogonal frequency modulation (FM) and inter-pulse phase coherence. This innovative approach not only presents formidable challenges for the jammer in acquiring radar waveform parameters, but also bolsters the radar’s low probability of intercept (LPI), while maintaining the phase coherence of sub-pulses between pulses. Furthermore, based on this waveform, the characteristics of the intermittent sampling jamming signal and its differences from the target echo signal are analyzed in the time domain, frequency domain, time-frequency domain, and pulse compression domain. Building upon these findings, this paper proposes the sub-division algorithms for typical types of intermittent sampling jamming under this waveform: the full-pulses multi-level maximum inter-class variance and sub-pulses multi-level maximum inter-class variance anti-intermittent sampling jamming algorithms. Simulation analysis demonstrates that this waveform and the anti-jamming algorithms can accurately identify and effectively counteract different types of intermittent sampling jamming in typical scenarios.

1. Introduction

With the continuous advancement of digital radio frequency memory (DRFM) hardware technology [1], intermittent sampling jamming poses an ever more formidable menace, characterized by amplified transmission power, heightened reception sensitivity, accelerated storage response rates [2], and expanded frequency range coverage [3,4]. Undeniably, this multifaceted threat poses a grave risk to the survival of radar systems on the battlefield and their capacity for target identification. Intermittent sampling jamming intercepts part of the radar transmission signal, which is then carefully modulated and forwarded to the radar to create various interference patterns, including active false targets with high fidelity, controllable quantity, and configurable distribution area [5,6,7]. Traditional coherent pulse radar has limited freedom in waveform parameters and poor interception performance against this type of jamming. Anti-interference signal design also has limited degrees of freedom and is less effective, and severe cases may even disable the radar [8,9]. Nevertheless, as per the existing literature, research on identification and suppression methods for intermittent sampling jamming is still in its nascent stage. The suppression of intermittent sampling jamming is mainly carried out within the domains of polarization, time-frequency, and waveform.
In the time-frequency domain, the discontinuous characteristics of intermittent sampling jamming in comparison to the target signal can be effectively utilized [10,11]. One strategy involves constructing a one-dimensional time-frequency energy map by integrating the frequency domain of the target and jamming signals. Then a bandpass filter can be designed to selectively extract the desired target signal [12]. Additionally, effective feature parameters such as time-frequency separability and third-order Renyi entropy can be extracted to differentiate the interference [13]. These feature extraction methods possess a manageable computational complexity, making them suitable for practical engineering implementation. However, it is worth noting that while these methods are effective in addressing intermittent sampling direct jamming scenarios characterized by a low duty cycle, their performance may deteriorate significantly when faced with intermittent sampling repeater jamming and intermittent sampling cyclic jamming, both of which exhibit a high duty cycle.
In the polarization domain, the optimization criterion is to reduce the energy of intermittent sampling direct jamming, under the simultaneous full-polarization radar system. By considering the constant modulus constraint, a joint design of transmit waveforms and receive filters is conducted in a fully polarized manner to highlight targets and suppress interference [14]. Additionally, target characteristics can be incorporated to accommodate extended targets in wideband radar scenarios [15]. These methods utilize polarization information to distinguish between targets and interference. However, they are specifically designed to address intermittent sampling direct jamming and require prior knowledge of partial interference parameters from the jamming source in order to solve complex waveform optimization problems. Consequently, their applicability in rapidly changing interference scenarios is challenging.
In the waveform domain, conventional frequency agile waveforms exhibit greater effectiveness in countering deceptive jamming over pulse repetition periods [16,17,18], but they are struggling with intermittent sampling jamming. This is because these agile waveforms with varying inter-pulse frequencies make it difficult for deceitful jamming across pulse repetition periods to keep up with the rapid pace of inter-pulse frequency changes [19]. As a result, the frequency band of the jamming signal is isolated from the echo signal of this period, making it relatively easy to separate the signal and interference in the frequency domain using a bandpass filter. On the other hand, intermittent sampling jamming involves slicing and storing the transmitted signal for rapid retransmission, allowing for interference within the current pulse repetition period to disrupt the radar.
Therefore, in recent years, some scholars have presented proposals to improve the aforementioned waveforms. For instance, the utilization of intra-pulse orthogonal linear frequency modulation combined with phase encoding waveforms has been proposed [20]. This method yields waveforms with reduced doppler sensitivity and low sidelobe characteristics. By leveraging the effects of mismatched filtering on intermittent sampling jamming, it becomes possible to differentiate between targets and interference, leading to higher identification rates for intermittent sampling direct jamming. However, obtaining precise analytical expressions for the phase encoding in the waveform proves challenging, and it is also sensitive to the Signal-to-Noise Ratio (SNR), thereby making its application in engineering practice difficult. Moreover, by exploiting the characteristic of intermittent sampling jamming, which fails to capture the complete transmitted pulse, the adoption of intra-pulse stepped linear frequency modulation waveform allows for the extraction of undisturbed pulse echoes [21]. This approach demonstrates effective performance in identifying intermittent sampling direct jamming and intermittent sampling repeater jamming. Nonetheless, the discrimination threshold is typically set based on previous experience, and it proves inadequate in adapting to variations in interference and radar parameters. Furthermore, some scholars have adopted similar concepts to construct inter-pulse and intra-pulse agile waveforms, aiming to identify targets by extracting features such as fuzzy C-means clustering and variance [19,22,23]. However, these waveforms exhibit random frequency hopping between pulses, leading to non-linear phase transitions between pulses. Consequently, the attainment of coherent synthesis of inter-pulse phases in practical engineering becomes challenging [24], thus making them less suitable for application in coherent radar systems such as Pulse Doppler (PD) radar.
Building upon these foundations, this paper proposes a radar waveform that integrates intra-pulse random orthogonal frequency modulation with inter-pulse phase coherence. Initially, within the operating bandwidth, the sub-pulse carrier frequencies are randomly encoded to ensure mutual orthogonality between the frequency bands of the sub-pulses. Subsequently, during the inter-pulse period, these sub-pulses are shuffled and rearranged to ensure phase coherence among sub-pulses with the same encoding. The combination of intra-pulse random orthogonal frequency modulation for sub-pulses effectively restricts the interference system’s ability to extract radar waveform parameters, thereby reducing the probability of interception by the jamming system [19,22]. Moreover, the inter-pulse phase coherence among sub-pulses enables the proposed waveform to be applicable in coherent radar systems, including Moving Target Detection (MTD) radar and PD radar [25].
In the aforementioned literature, it can be observed that these anti-jamming methods are analyzed in four signal domains: time domain, frequency domain, time-frequency domain, and pulse compression domain. These methods mainly focus on distinguishing between the target signal and intermittent sampling direct jamming or intermittent sampling repeater jamming. However, there exists a dearth of comprehensive distinction between the target signal and the three typical types of intermittent sampling jamming, particularly the more challenging intermittent sampling repeater jamming and intermittent sampling cyclic jamming, which can lead to heightened confusion. Hence, building upon this background, this paper presents a comprehensive analysis of the target signal and the three types of intermittent sampling jamming (direct, repeater, and cyclic) within the context of the proposed waveform. The analysis and simulation experiments are conducted in the time domain, frequency domain, time-frequency domain, and pulse compression domain to explore the characteristic differences between the target signal and the interference signals. Moreover, matched filters are constructed based on the frequency encoding of the intra-pulse sub-pulses. After obtaining the variance of each compressed sub-pulses, the paper proposes two novel algorithms: the full-pulses multi-level maximum inter-class variance algorithm and the sub-pulses multi-level maximum inter-class variance algorithm. These algorithms enable precise identification and differentiation of the target signal and the three types of intermittent sampling jamming, facilitating accurate target detection and effective mitigation of different types of intermittent sampling jamming.

2. Models of the Intermittent Sampling Jamming and Interference Hypothesis

From a procedural perspective, intermittent sampling jamming can be classified into three phases: sampling, modulation processing, and transmission. During the sampling phase, the jammer searches and arranges the signal, then samples it. In the processing and modulation stage, the jammer modulates, stores, and delays the sampled signal to generate the jamming signal. During the transmission stage, the jammer amplifies and transmits the jamming signal. By adjusting the sampling timing, duration, delay time, and choosing suitable modulation and processing techniques, the number, strength, and spatial distribution of the jamming false targets can be regulated [26].
Mathematically, the sampling process can be modeled as a signal multiplication to characterize. Let the signal be S ( t ) , and the sampled signal be P ( t ) ; then we have
P ( t ) = rect t τ J n = + δ t n T J
where τ J is the pulse width of the rectangular envelope of the disturbance (and often the sampling pulse width), and T J is the repetition period of the sampled signal. The sampled signal is
X ( t ) = P ( t ) S ( t )
In terms of modulation processing style, intermittent sampling jamming can be divided into three categories: intermittent sampling direct jamming, intermittent sampling repeater jamming, and intermittent sampling cyclic jamming [27,28]. A visual illustration of the intermittent sampling direct jamming is presented in Figure 1.
The mathematical expression of intermittent sampling direct jamming is shown in Equation (3); τ s is the sampling time and τ d is the processing time.
X direct t = X t τ s τ d
To enhance the duty cycle of intermittent sampling jamming, which is hampered by low duty cycles due to direct jamming, the sampled signal can be relayed multiple times within an interference repetition cycle rather than just once. The maximum allowable number of relays can be determined as
H = T J τ s τ d τ J = T w τ J
T w is the actual length of time available for forwarding in an interference cycle T J ; that is, T J = T w + τ d + τ s . This gives the signal form of intermittent sampling repeater jamming as
X repeat = i = 0 H 1 X direct t i τ J
A visual illustration of the intermittent sampling repeater jamming is presented in Figure 2. Building upon this concept, it is feasible to incorporate not only a repeat of the current signal but also the previously sampled signal into the forwarded signal, thereby implementing intermittent sampling cyclic jamming. Figure 3 depicts the typical reconnaissance, interference timing, and modulation pattern of intermittent sampling cyclic jamming. In this instance, the interfering signal can be formulated as
X loop = i = 0 Z 1 X direct t i T J + τ J
Z = min G , H , G = T p T J + 1 , T p is the pulse width of the radar transmit signal, and G is the number of interference cycles within each pulse repetition cycle.

3. Description of the Proposed Method

3.1. Model of the Intra-Pulse Randomly Coded Quadrature FM and Inter-Pulse Phase Coherence Waveform

Let the radar transmit N pulses; each pulse has a pulse width of τ and a repetition period of T s . There are P sub-pulses within each pulse, the operating frequency bandwidth is B, and the carrier frequency difference of P sub-pulses is an integer multiple of Δ f . Then there are B Δ f frequency bands available in the total frequency bandwidth, and for all sub-pulses, their carrier frequencies are selected by random shuffling in the above frequency bands. Then the mathematical expression of the mth sub-pulse within the nth pulse is
S ( t , n , m ) = a ( t ) δ t m + 1 2 τ sub n T s e j 2 π f n m t
a ( t ) can be LFM [29] or non-LFM, pulse, phase coding, and other signals. In this paper, we only discuss LFM signals, in which case a ( t ) can be expressed as
a ( t ) = rect t τ sub e j π k t 2
where τ sub is the sub-pulse pulse width, k is the FM slope, and the sub-pulse bandwidth is B sub . The LFM slope is calculated as k = B sub τ sub ; rect (t) is a rectangular function.
rect t τ sub = 1 , t τ sub 1 2 0 , else
where f n m is the carrier frequency of the mth sub-pulse, f n m = f c + Δ f c ( n , m ) , and f c is the initial carrier frequency. The formation of c ( n , m ) is accomplished in the following manner. Firstly, m elements are randomly chosen from set 0 , 1 , 2 , B Δ f , and these m elements are strictly orthogonal in the frequency domain. Subsequently, within each pulse, the order of these m elements is shuffled and rearranged. c ( n , m ) is also referred to as the frequency modulation coefficient of the mth sub-pulse in the nth pulse.
For the intra-pulse random quadrature FM and inter-pulse phase coherence waveform, the waveform schematic is shown in Figure 4, and the mathematical expression is
S ( t ) = n = 0 N 1 m = 0 P 1 S ( t , n , m )
Thus, the echo signal received by the radar is R ( t ) , shown as Equation (11), where S r ( t ) is the target echo signal, X J ( t ) is the intermittent sampling jamming, and N ( t ) is the noise signal.
R ( t ) = S r ( t ) + X J ( t ) + N ( t )
As intermittent sampling jamming merely divides the transmit signal and fails to store complete pulses, certain sub-pulses are inevitably lost, posing a significant advantage for the detection of intermittent sampling jamming. In the case of the aforementioned waveform, the transmitter possesses knowledge of the frequency modulation number c ( n , m ) of every sub-pulse. Thus, a band-pass filter can be tailored for each sub-pulse, and subsequently, the signal can be down-converted and pulse-compressed to derive the echo signal y ( t , n , m ) of the corresponding sub-pulse, as demonstrated in Equation (12). Consequently, the interference can be identified utilizing the following algorithm.
y ( t , n , m ) = R ( t ) S * ( t , n , m )

3.2. Full-Pulses and Sub-Pulses Multi-Level Maximum Inter-Class Variance Algorithm

Utilizing an intra-pulse randomly coded orthogonal FM waveform to enhance the information entropy of the transmit waveform is akin to constructing a higher-dimensional function space compared to conventional signal. The scattering process of the target and passive environment on the transmit signal is, in turn, equivalent to a linear mapping of the elements within the function space. Nevertheless, due to limitations arising from space-time constraints and detection and binning processes, the signal generated by intermittent sampling jamming corresponds to a nonlinear mapping of the elements within the function space. Upon considering the signal function’s spatial dimension, the interfering signal falls short of the necessary dimensionality, and thus, the disparity in feature dimensions can be leveraged to validate the echoes, enabling the identification and differentiation of interfering echoes.
The echo signal of the randomly coded quadrature FM waveform can be extracted by analyzing its mean, variance, moment kurtosis, and other characteristic elements within various function spaces such as the time domain, frequency domain, time-frequency domain, and fractal domain. Upon conducting a comparative study, it was observed that the variance of the sub-pulses following pulse compression serves as a differentiating factor between the target signal and intermittent sampling jamming. Additionally, the quadrature FM waveform employed in this paper encodes both intra-pulse and inter-pulse sub-pulses, resulting in slight differences in the variance of the three types of intermittent sampling jamming. Thus, the variance of the sub-pulses following pulse compression can be utilized to distinguish between the target signal and the three types of interference with precision.
The OTSU algorithm is a widely used image segmentation algorithm that selects the optimal threshold to maximize inter-class variance and minimize intra-class variance, reducing the probability of misclassification [30,31,32]. The multi-layer OTSU algorithm is an extension of the OTSU algorithm that divides data into multiple classes [33,34,35], rather than just two, with the same aim of minimizing intra-class variance and maximizing inter-class variance.
In this paper, the full-pulses multilevel maximum interclass variance algorithm and the sub-pulses multilevel maximum interclass variance algorithm are proposed based on the multilayer OTSU algorithm, specifically for the waveforms used in the study. The algorithm categorizes the variance of the absolute values of each sub-pulse after pulse compression to maximize the variance between classes. Hence, there are two concepts of “variance” here: one as the input to the algorithm and the other as the objective function of the algorithm. The procedure for the full-pulses multilevel maximum interclass variance algorithm is as follows:
1.
The absolute value of the time domain variance of the n ( n = 1 , 2 , , N ) th sub-pulse of the p ( p = 1 , 2 , , P ) th echo signal of the four echoes after pulse compression is denoted as var ( n , p ) target , var ( n , p ) straight , var ( n , p ) reapet , and var ( n , p ) loop , and these data are arranged from small to large into a one-dimensional array Var of size 4 N P ;
2.
Assume that the first level threshold T 1 is chosen as the ith data Var ( i ) of the array Var , and the second and third level thresholds T 2 , T 3 are the jth data Var ( j ) of Var and the kth data Var ( k ) of Var , and i < j < k < 4 N P , such that the array can be divided into four subintervals, g 1 , g 2 , g 3 , and g 4 , respectively;
3.
Calculate the probability of occurrence of each subinterval, denoted as:
w 1 = i 4 N P
w 2 = j i 4 N P
w 3 = k j 4 N P
w 4 = 4 N P k 4 N P
4.
Calculate the average amplitude of the four intervals and the total interval:
λ 1 = h = 1 i Var ( h ) i
λ 2 = h = i + 1 j Var ( h ) j i
λ 3 = h = j + 1 k Var ( h ) k j
λ 4 = h = k + 1 4 N P Var ( h ) 4 N P k
λ = h = 1 4 N P Var ( h ) 4 N P
5.
Define the interclass variance as:
σ 2 = h = 1 4 w h λ λ h 2
6.
Iterating over the threshold T 1 to T 3 (satisfying T 1  <  T 2  <  T 3 ) such that the maximum inter-class variance σ 2 is the optimal threshold:
T 1 * , T 2 * , T 3 * = arg max T 1 < T 2 < T 3 σ 2
Theflowchart of the full-pulses multilevel maximum interclass variance algorithm is shown in Figure 5.
The full-pulses multilevel maximum interclass variance algorithm considers all sub-pulses for threshold determination, which is effective for distinguishing between target signals and intermittent sampling jamming, but may not perform well when differentiating between different types of intermittent sampling jamming. This is because the quadrature FM waveform used in this paper varies significantly between sub-pulses, and the different types of intermittent sampling jamming also have different impacts on sub-pulses. If the same threshold is applied to all sub-pulses, misclassification of the types of interference is likely to occur. Therefore, to address this issue, differential thresholds can be set for different sub-pulses. Based on this approach, a sub-pulses multilevel maximum interclass variance algorithm is proposed, and the schematic diagram of the algorithm is illustrated in Figure 6.

4. Analysis of the Proposed Method

To validate the proposed method described in the paper, theoretical analysis and simulation experiments are conducted on the proposed waveform and the intermittent sampling jamming under this waveform. Hence, a meticulous analysis is performed in the time domain, time-frequency domain, frequency domain, and pulse compression domain to examine the characteristics of the intra-pulse random orthogonal and inter-pulse coherent frequency-modulated waveform signal, as well as three types of intermittent sampling jamming signals.
Table 1 displays the parameters of the waveform utilized in the experiment. The radar system operates at a carrier frequency of 3 GHz in the S-band, with a pulse width of 41 μ s and a pulse repetition period of 1 ms. The experiment consists of a total of 10 pulses, with each pulse containing 4 sub-pulses. The sub-pulse bandwidth is set at 10 MHz, resulting in a total operating frequency bandwidth of 400 MHz, with a frequency interval of 40 MHz for each hop. The interference pulse width is 1 μ s, with an interference time delay of 0.2 μ s, and an interference repetition period of 10.2 μ s.

4.1. Time Domain Analysis and Simulation

Before conducting the time domain simulation experiments, let us delve into the difference between the signal and the three types of interference in the time domain from a theoretical point of view. A more in-depth study of the time domain expression of the emitted signal waveform is carried out according to Equation (7):
S ( t ) = n = 0 N 1 m = 0 P 1 rect t τ sub e j π k t 2 δ t m + 1 2 τ sub n T s e j 2 π f n m t = n = 0 N 1 m = 0 P 1 + rect u τ sub e j π k u 2 δ t u m + 1 2 τ sub n T s d u e j 2 π f n m t = n = 0 N 1 m = 0 P 1 exp j π k t n T s m + 1 2 τ sub 2 + j 2 π f n m t rect t n T s m + 1 / 2 τ sub τ sub
Hence, the S ( t ) signal is composed of P linearly frequency-modulated sub-pulses with different carrier frequencies within each pulse. Each sub-pulse has a width of τ sub . By utilizing Equations (2) and (3), we can derive the time domain expression of the intermittent sampling direct relay interference as follows:
X ( t ) = rect t τ J n = + δ t n T J S ( t ) = + rect u τ j n = + δ t u n T J d u S ( t ) = n = + rect t n T J τ j S ( t )
Therefore, intermittent sampling direct jamming corresponds to the temporal interception of the transmitted signal S ( t ) with a width of τ j and a period of T J . Intermittent sampling repeater jamming and intermittent sampling cyclic jamming both stem from intermittent sampling direct jamming.The time domain simulation plot is as follows.
From Figure 7, we can observe that within one pulse period, there are four sub-pulses with different carrier frequencies, each marked with a corresponding color. The characteristics of intermittent sampling repeater jamming are particularly evident in the plot, especially in the case of the interference signal for the fourth sub-pulse (indicated by the color purple), which undergoes H repetitions. Additionally, it can be observed that intermittent sampling cyclic jamming not only relays the currently intercepted sub-pulse signal but also relays past intercepted sub-pulse signals. This cyclic relay of past intercepted sub-pulse signals can generate multiple coherent false target echoes in the radar system, which will be analyzed in detail in the subsequent pulse compression analysis.

4.2. Time-Frequency Domain Analysis and Simulation

The analysis in the time-frequency domain builds upon the foundation of frequency domain analysis. According to Equation (24), we can determine the phase of each sub-pulse as follows:
Φ ( t ) = π k t n T s m + 1 2 τ sub 2 + 2 π f n m t
The phase is a quadratic function of time. Taking the derivative of the phase with respect to time, we can obtain the instantaneous frequency as follows:
f ( t ) = 1 2 π d Φ ( t ) d t = k t + f n m
It is evident that the frequency is a linear function of time. Therefore, the frequency range corresponding to each sub-pulse is f n m B sub 2 , f n m + B sub 2 , and the bandwidth for each sub-pulse is B sub . This corresponds to the concept of sub-pulse segmentation filtering discussed in Section 3.2. Therefore, when setting the hopping frequency and total bandwidth, it is necessary to ensure that the frequencies of each sub-pulse do not overlap to extract each sub-pulse without losing information in the spectrum. This lays the foundation for subsequent signal identification and processing.
Regarding the interference, according to Equation (25), we can infer that the phase of the interference signal is approximately the same as the transmitted signal. However, since the interference signal only intercepts a portion of the transmitted signal, its bandwidth is B j = k τ j .
From Figure 8, it can be observed that the intermittent sampling repeater jamming exhibits a relatively continuous pattern in the time-frequency domain. Therefore, distinguishing between the target echo signal and intermittent sampling repeater jamming based on signal continuity in the time-frequency domain is not very effective. As introduced in Section 1, signal and intermittent sampling jamming identification in the time-frequency domain is more effective for low-duty intermittent sampling direct jamming and more challenging for high-duty intermittent sampling repeater jamming.
It is clear that the intermittent sampling cyclic jamming corresponds to a forwarding period of T J + τ J and a forwarding count of G for the interfering sub-signal within a signal pulse width, which aligns with Equation (6).

4.3. Frequency Domain Analysis and Simulation

According to Section 4.1, by applying the Fourier transform to the time domain signals, we can explore the distinctions between the signal and the three types of interference in the frequency domain.
S ( f ) = F n = 0 N 1 m = 0 P 1 exp j π k t n T s m + 1 2 τ sub 2 + j 2 π f n m t rect t n T s m + 1 / 2 τ sub τ sub = n = 0 N 1 m = 0 P 1 F exp j π k t n T s m + 1 2 τ sub 2 + j 2 π f n m t rect t n T s m + 1 / 2 τ sub τ sub = n = 0 N 1 m = 0 P 1 S ( f , n , m )
Let F denote the Fourier transform. According to the nature of the Fourier transform as a linear transform, we can infer that the spectrum of each transmitted pulse signal is the summation of the spectra of all sub-pulses. Let t 0 = m + 1 2 τ sub ; we can obtain
S ( f , n , m ) = F exp j π k t t 0 2 + j 2 π f n m t rect t t 0 τ sub = + exp j π k t t 0 2 + j 2 π f n m t exp ( j 2 π f t ) rect t t 0 τ sub d t = + rect t t 0 τ sub exp j π k t t 0 2 + j 2 π f n m t j 2 π f t d t = + A ( t ) exp j θ ( t ) d t
where A ( t ) = rect t t 0 τ sub , θ ( t ) = π k t t 0 2 + 2 π f n m t 2 π f t .
A ( t ) represents the real envelope, while θ ( t ) represents the signal modulation phase. In comparison to the phase, the envelope is a slowly varying function of time. According to the principle of stationary phase [36,37,38], for a rapidly changing signal in the time domain, apart from the positions where the derivative is zero (stationary points), the positive and negative areas in the remaining regions almost cancel each other out. As a result, the integral of the function is primarily influenced by the stationary point regions. Now, let us solve for the stationary points:
θ ( t ) = d θ ( t ) d t = 2 π k t + 2 π f n m 2 π f
Thus, the stationary point is t k = f f n m k + t 0 , and then the Taylor expansion [39,40] is performed for θ ( t ) in the vicinity of t k :
θ ( t ) = θ ( t k ) + θ ( t k ) ( t t k ) + θ ( t k ) 2 t t k 2 + O ( t t k 2 ) = θ ( t k ) + θ ( t k ) 2 t t k 2 + O ( t t k 2 )
where θ ( t ) = d θ ( t ) d t = 2 π k .
By neglecting higher-order terms beyond the second order and applying the Fresnel integral theorem [41,42], the expression can be simplified as follows:
S ( f , n , m ) = A t k exp [ j θ ( t k ) ] t k δ t k + δ exp j θ ( t k ) 2 t t k 2 d t = 1 k rect f f n m B sub exp j π f f n m 2 k exp j 2 π f f n m t 0 exp j π 4
where δ represents an infinitesimal quantity.
Therefore, from Equation (32), the envelope of S ( f , n , m ) resembles a rectangle with a width distribution of f n m B sub 2 , f n m + B sub 2 . This corresponds to the frequency spectrum bandwidth of the m-th sub-pulse of the n-th pulse, which is consistent with the previous theoretical analysis.
On the basis of Equation (2), for intermittent sampling direct jamming, the spectrum is
X ( f ) = F P ( t ) S ( t ) = F ( P ( t ) ) F S ( t ) = F rect t τ J n = + δ t n T J S ( f ) = τ j f J n = + S a ( π τ j f J n ) δ f n f J S ( f ) = τ j f J n = + S a ( π τ j f J n ) S ( f n f J )
where f J = 1 T J is the sampling frequency of the interference, S a ( x ) = sin ( x ) x .
From Equation (33), it can be observed that the spectrum of intermittent sampling direct jamming is a periodic replication of the transmit signal, weighted by the amplitude and with a weight factor of τ j f J S a ( π τ j f J n ) . The replication period is f J . The spectral analysis of intermittent sampling repeater jamming and cyclic jamming is similar to this, and the following simulations are performed to compare.
From Figure 9, we can see that the carrier frequencies corresponding to the four sub-pulses of the original signal are 90 MHz, 210 MHz, 310 MHz, and 360 MHz, respectively. The bandwidth of each sub-pulse is 10 MHz, which is B sub . It is evident from the figure that the spectral bandwidths of intermittent sampling direct jamming, repeater jamming, and cyclic jamming are similar to the radar’s transmitted signal. This spectral overlap enables highly effective interference to the radar target echoes in the frequency domain. However, due to the fact that the intermittent sampling jamming is a truncation of the transmitted signal, its spectral side lobes are significantly higher and the main lobe is widened. Moreover, in the frequency domain, there is a clear distinction between the target signal and the three types of intermittent sampling jamming, but the differences among the three interference types are relatively small.

4.4. Pulse Compression Domain Analysis and Simulation

In the above experiments, we can validate the simulation experiments through theoretical analysis, revealing that the distinction between the target signals and the three types of intermittent sampling jamming is inadequate in the time domain, time-frequency domain, and frequency domain when employing the waveform incorporating intra-pulse random orthogonal coded frequency modulation and inter-pulse phase coherence.
Therefore, it is crucial to conduct theoretical analysis and experimental verification of the transmitted signal and the three types of intermittent sampling jamming in the pulse compression domain. After pulse compression, the differences between the transmitted signal and the interference can be clearly manifested. Thus, the variance of each sub-pulse after pulse compression can serve as a discriminant. In this study, the proposed Full-Pulse Multilevel Maximum Inter-class Variance algorithm and Sub-Pulse Multilevel Maximum Inter-class Variance algorithm are utilized to analyze their effectiveness in discriminating between targets and interference.
As the frequency bands of different sub-pulses are orthogonal to each other, in signal processing, different filters are constructed for each sub-pulse based on their corresponding modulation index ( c n , m ). Then, the echoes of these sub-pulses are individually pulse-compressed. In this case, we can fix the pulse number (n) and sub-pulse number (m) for analysis. For the sake of demonstration, let us consider n = m = 0. According to Equation (24), we have
S ( t ) = rect t τ sub / 2 τ sub exp j π k t τ sub 2 2 + j 2 π f 0 t
Then the filter corresponding to pulse compression is
h ( t ) = S * ( t ) = rect t τ sub / 2 τ sub exp j π k t + τ sub 2 2 + j 2 π f 0 t
The result of pulse compression is
S 0 ( t ) = h ( t ) S ( t ) = + rect u τ sub / 2 τ sub exp j π k u + τ sub 2 2 + j 2 π f 0 u * rect t u τ sub / 2 τ sub exp j π k t u τ sub 2 2 + j 2 π f 0 ( t u ) d u = sin π k t τ sub 1 t τ sub π k t exp ( j 2 π f 0 t )
The envelope of the pulse compression result is approximately τ sub S a π B sub t . When π B sub t = ± π , we can obtain t = ± 1 B sub . This is the first zero point, that is, the interval between the first zero points is 2 B sub . For intermittent sampling direct jamming, a similar derivation can lead to its pulse compression envelope expression as τ j S a π B j t ; likewise, the interval between its first zero points is 2 B j . Intermittent repeater forwarding interference and intermittently cyclic forwarding interference both involve forwarding the interfering sub-pulses of the intermittent direct interference with different repetitions and time delays. Therefore, the differences between the signal and the three types of interference after pulse compression can be verified through simulation.
From Figure 10, it can be observed that the distance interval between the first zero crossings of the target echo sub-pulses after pulse compression is approximately 30 m, which corresponds to a time interval of 0.2 μ s. This is consistent with the theoretical analysis of interval 2 B sub . On the other hand, the first zero crossing interval of the intermittent sampling direct jamming is approximately 300 m, which matches the theoretical analysis of interval 2 B j .
The spacing between the false targets generated by intermittent sampling repeater jamming is approximately 150 m, which corresponds to a time interval of 1 μ s. This is just the interference pulse width τ j . This correspondence is consistent with Figure 8, as during repeater jamming, each interference sub-pulse is delayed by τ j compared to the previous one. Therefore, the repetition period corresponds to τ j , and due to the close proximity of each dummy target, the sidelobes overlap during pulse compression, resulting in wider sidelobes with higher amplitudes. Similarly, the interval between dummy targets generated by intermittent sampling cyclic jamming is approximately 1680 m, or 11.2 μ s, which is just one interference period plus the interference pulse width, denoted as T J + τ j . Additionally, due to the relatively larger separation between each dummy target, the mutual influence of sidelobes is reduced.
This observation reveals that under this waveform, there exist significant distinctions among the target echo signal, intermittent sampling direct jamming, intermittent sampling repeater jamming, and intermittent sampling cyclic jamming after pulse compression. Consequently, the variance obtained after pulse compression can be utilized for classification, serving as the fundamental basis for the Full-Pulse Multilevel Maximum Inter-class Variance algorithm.
Furthermore, from Figure 11, it is evident that the pulse compression results of intermittent sampling repeater jamming vary among different sub-pulses. In fact, due to the changing relationship between interference and signal in terms of period and pulse width, the effects generated by intermittent sampling jamming differ within each sub-pulse. Therefore, it is possible to perform signal classification separately within each sub-pulse, which forms the fundamental basis of the Sub-Pulses Multilevel Maximum Inter-class Variance algorithm.

5. Echo Certification Simulation Experiment and Analysis

Based on Section 4, we set the number of radar transmitting pulses N to 10, sub-pulse P to 8, and varied the Jamming-to-Signal Ratio (JSR) from 10 dB to 20 dB, as well as the Signal-to-Noise Ratio (SNR) from −20 dB to 20 dB after sub-pulse matched filtering. The simulation analyzed the recognition accuracy of the full-pulses multilevel maximum inter-class variance algorithm and the sub-pulses multilevel maximum inter-class variance algorithm for target signal and three types of interference. The simulation was repeated 300 times using the Monte Carlo method. Figure 10 shows the results of the 300 Monte Carlo simulations when JSR is 10 dB.
Figure 12 demonstrates that with increasing SNR after pulse compression, the recognition accuracy of the target echo signal, intermittent sampling direct jamming, and intermittent sampling cyclic jamming all improve when using the full-pulses shared thresholding algorithm. However, at low SNR levels, the algorithm primarily calculates the variance of noise undulation, leading to lower recognition accuracy.
As the SNR increases, the full-pulses shared threshold algorithm can distinguish the variances of the target and the three types of interference based on their different energy and signal undulation characteristics. The recognition accuracy of the target echo signal is stable at 100% for SNR = −5 dB, while the recognition accuracy of intermittent sampling direct jamming and intermittent sampling cyclic jamming is stable at approximately 85%. The recognition accuracy of all four signals is stable at around 78%.
However, the recognition accuracy of the full-pulses shared threshold algorithm for intermittent sampling repeater jamming is low. From Figure 11, we can observe that the interference effects vary among different sub-pulses. Therefore, if the same threshold is applied to all sub-pulses, it will inevitably lead to misclassifications. Furthermore, the variance of intermittent sampling cyclic jamming is similar to that of intermittent sampling repeater jamming, making it difficult to distinguish the subtle difference in the case of full-pulses shared thresholding. Consequently, intermittent sampling repeater jamming is often recognized as intermittent sampling cyclic jamming, as shown in the fuzzy matrix in Table 2 and Table 3 for SNR = 0 dB.
From the results, it is clear that when 40 sub-pulses are interfered by intermittent sampling repeater jamming, they are misidentified as intermittent sampling cyclic jamming using a common threshold for all pulses. Thus, it is difficult to distinguish intermittent sampling repeater jamming from intermittent sampling cyclic jamming with a common threshold for all pulses. However, the recognition accuracy of the target echo signal is almost perfect, with a recognition accuracy close to 100%.
Based on the information provided in Figure 13 and Figure 14, it is evident that the intermittent sampling repeater jamming has a larger variance than the intermittent sampling cyclic jamming, which has a smaller variance. When applying the full-pulses shared threshold algorithm, it becomes clear that the two types of interference cannot be accurately distinguished at either the first or third sub-pulses, as the variance of the interferences varies across different sub-pulses. Hence, using a global threshold for all sub-pulses would result in a significant reduction in recognition accuracy. In contrast, the sub-pulses adaptive thresholding approach can effectively differentiate between the two types of interference by employing a set threshold.
Figure 15 illustrates that the sub-pulses adaptive thresholding algorithm shows an increasing trend in recognition accuracy as SNR increases for both the signal and the three interferences, with a stabilized recognition accuracy of about 93%. It can be concluded that this algorithm has a significantly higher recognition rate for intermittent sampling repeater jamming compared to the full-pulses shared thresholding algorithm. However, due to the smaller energy and variance fluctuations of the echo signal within each sub-pulse compared to the interference, the recognition accuracy decreases compared to the full-pulses shared threshold algorithm. Overall, the sub-pulses adaptive thresholding algorithm has a higher recognition accuracy for all signals and interferences, but at the cost of a higher SNR for each sub-pulse.
The confusion matrix in Table 4 indicates that the recognition rate of the target signal and the three types of interference is high, and the intermittent sampling repeater jamming and intermittent sampling cyclic jamming can be accurately distinguished.
When the JSR is increased to 20 dB, the results shown in Figure 16 and Figure 17 indicate that the SNR values of both the signal and interference decrease while the recognition accuracy remains stable, as compared to the JSR of 10 dB. This is attributed to the enhanced energy of the interference signal and more distinct variance characteristics, allowing for a stable recognition accuracy to be achieved at a lower SNR.
Based on the theoretical analysis and experimental results, it can be concluded that the proposed algorithms have demonstrated their effectiveness in identifying and distinguishing different types of interference in radar systems. The full-pulses multi-level maximum inter-class variance algorithm is capable of detecting intermittent sampling interference and the target signal with high accuracy, while the sub-pulses multi-level maximum inter-class variance algorithm is more effective in identifying the three types of intermittent sampling interference. By combining these two algorithms in cascade, it is possible to achieve highly accurate identification of the target signal and effective distinction of different types of interference in radar systems.

6. Conclusions

With the purpose of addressing the challenge pertaining to the identification of the target signal and diverse forms of intermittent sampling jamming, this paper puts forth a radar waveform that amalgamates intra-pulse random orthogonal coded frequency modulation with inter-pulse phase coherence. This waveform takes into account the functions of anti-jamming and inter-pulse coherence. Leveraging this proposed waveform, the paper introduces the full-pulses multi-level maximum inter-class variance algorithm, which exhibits a remarkable recognition rate for the target signal. When the Jamming-to-Signal Ratio (JSR) is 10 dB and the Signal-to-Noise Ratio (SNR) is −5 dB, the target recognition rate approximates 100%. Nevertheless, distinguishing between intermittent sampling repeater jamming and intermittent sampling cyclic jamming proves challenging for the algorithm. Building upon this premise, the paper enhances the algorithm and introduces the innovative sub-pulses multilevel maximum interclass variance algorithm. This algorithm further employs adaptive thresholds tailored to individual sub-pulses, thereby improving the accuracy of the recognition process. With respect to typical SNR and JSR conditions, the recognition rate for both the signal and the three types of jamming surpasses 90%. Furthermore, the algorithm’s adaptability to diverse environments enables its practical implementation in engineering. This, in turn, enhances the resistance of coherent phase radar system against intermittent sampling jamming.
In MTI and MTD radar systems, the waveform proposed in this paper can not only obtain better anti-jamming effect, but also be compatible with anti-clutter processing. Therefore, one of the forthcoming objectives of this study involves developing an effective approach to filter out intermittent sampling jamming amidst the presence of significant ground clutter. Furthermore, in this paper, the waveform is designed to vary the sub-pulse carrier frequency within the pulse and the relative order of the sub-pulses between the pulses, without further analysis of the more complex sub-pulse carrier frequency-pulse repetition period joint variation waveform. Hence, the subsequent focus will be on exploring methods to enhance the variation of the pulse repetition interval (PRI) based on the existing waveform. This enhancement aims to expand the range of speed measurement and ensure robust interference recognition capabilities.

Author Contributions

H.Z. performed the theoretical study, conducted the experiments, processed the data, and wrote the manuscript. T.W. helped on the theoretical study, provided research suggestions, and revised the manuscript. T.G. and X.S. helped in performing the experiments and gave suggestions for the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported in part by the Science and Technology Planning Project of Guangdong Province under grants 2019B121203006, in part by the Shenzhen Science and Technology Program under grants 2019ZT08X751.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the fact that it is currently privileged information.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ran, H.; Zhang, H.; Zhang, M. Development Status and Trend of Digital Radio Frequency Memory. Aerodyn. Missile J. 2018, 6, 44–48. [Google Scholar]
  2. Almslmany, A.; Wang, C.; Cao, Q. Advanced Deceptive Jamming Model Based on DRFM Sub-Nyquist Sampling. In Proceedings of the 2016 13th International Bhurban Conference on Applied Sciences and Technology (IBCAST), Islamabad, Pakistan, 12–16 January 2016. [Google Scholar]
  3. Li, Z.; Liao, X.; Shu, D.; Yu, W. A Study on Dynamic Reconfigurable Broadband Digital Channelized RF Storage Technology. J. Abbr. 2020, 42, 42–48. [Google Scholar]
  4. Li, C.; Su, W.; Gu, H.; Ma, C. Improved Interrupted Sampling Repeater Jamming based on DRFM. In Proceedings of the IEEE International Conference on Signal Processing, Guilin, China, 5–8 August 2014. [Google Scholar]
  5. Lan, L.; Liao, G.; Xu, J.; Zhang, Y. Suppression Approach to Main-Beam Deceptive Jamming in FDA-MIMO Radar Using Nonhomogeneous Sample Detection. IEEE Access 2018, 6, 34582–34597. [Google Scholar] [CrossRef]
  6. Feng, D.; Xu, L.; Pan, X.; Wang, S. Jamming Wideband Radar Using Interrupted-Sampling Repeater. IEEE Trans. Aerosp. Electron. Syst. 2017, 53, 1341–1354. [Google Scholar] [CrossRef]
  7. Zhang, J.; Wen, S.; Gao, H.; Sha, M. A New-style Jamming Technology and Implementation Based on Interrupted-sampling. Mod. Radar 2018, 40, 81–85. [Google Scholar]
  8. Wen, C.; Peng, J.; Zhou, Y.; Wu, J. Enhanced Three Dimensional Joint Domain Localized STAP for Airborne FDA MIMO Radar Under Dense False Target Jamming Scenario. IEEE Sens. J. 2018, 18, 4154–4166. [Google Scholar] [CrossRef]
  9. An, T.; Yang, Y.; Liu, J. Research on Self-Defense Jamming Technology of Radar Based on Intermittent Sampling. Ship Electron. Eng. 2020, 40, 79–82. [Google Scholar]
  10. Chen, J.; Wu, W.; Xu, S. A Band Pass Filter Design against Interrupted-Sampling Repeater Jamming based on Time-Frequency Analysis. IET Radar Sonar Navig. 2019, 13, 1646–1654. [Google Scholar] [CrossRef]
  11. Zhang, L.; Wang, G.; Yang, Z.; Zhang, X. Method for Identifying Radar Deception Jamming Using Time-Frequency Correlation Difference of Echo. J. Xi’an Jiaotong Univ. 2021, 55, 136–142. [Google Scholar]
  12. Yang, S.; Tian, B.; Zhou, R. ECCM Against Interrupted Sampling Repeater Jamming Based on Time frequency Analysis. J. Signal Process. 2016, 32, 1244–1251. [Google Scholar]
  13. Yang, S.; Tian, B.; Li, X.; Tan, M. A Recognition Method of LFM Radar Active Deception Jamming Based on SPWVD Figure. J. Air Force Eng. Univ. 2016, 17, 56–59. [Google Scholar]
  14. Wang, F.; Pang, C.; Huang, D.; Li, Y. An anti-jamming method against interrupted sampling repeater jamming based on designing transmitted waveforms and receiving filters for simultaneous polarimetric radar. Sci. Sin. 2022, 52, 1333–1348. [Google Scholar] [CrossRef]
  15. Wang, W.; Li, M.; Wang, F.; Rao, B. An Optimization Method for Transmitting Waveform of Polarimetric Radar Against Interrupted Sampling Repeater Jamming. J. Electron. Inf. Technol. 2023, 45, 1–10. [Google Scholar]
  16. Yu, H.; Li, X.; Li, Q.; Ma, H. Pulse Dopple Fuze Based on Pulse-to-Pulse Random Frequency Agility Technology. Aero Weapon. 2021, 28, 50–54. [Google Scholar]
  17. Quan, Y.; Wu, Y.; Li, Y.; Sun, G. Range–Doppler Reconstruction for Frequency Agile and PRF-Jittering Radar. IET Radar Sonar Navig. 2018, 12, 348–352. [Google Scholar] [CrossRef]
  18. Yang, H.; Liu, Z.; Wei, X. SAR Anti-Jamming Imaging Method Based on Frequency-Modulated Polarity Agility and Inter-Pulse Agility. Digit. Technol. Appl. 2015, 11, 118–120. [Google Scholar]
  19. Dong, S.; Quan, Y.; Sha, M.; Fang, W. Frequency Agile Radar Combined With Intra-Pulse Frequency Coding to Resist Intermittent Sampling Jamming. Syst. Eng. Electron. 2022, 44, 3371–3379. [Google Scholar]
  20. Zhou, C.; Tang, Z.; Yu, F.; Lu, Y. Anti Intermittent Sampling Repeater Jamming Method Based on Intrapulse Orthogonality. Syst. Eng. Electron. 2017, 39, 269–276. [Google Scholar]
  21. Zhang, J.; Mu, H.; Wen, S.; Li, Y. Anti Interrupted-Sampling Repeater Jamming Method Based on Stepped LFM Waveform. Syst. Eng. Electron. 2019, 41, 1013–1020. [Google Scholar]
  22. Dong, S.; Wu, Y.; Fang, W.; Quan, Y. Anti-Interrupted Sampling Repeater Jamming Method Based on Frequency-Agile Radar Joint Fuzzy C-means. J. Radars 2022, 11, 289–300. [Google Scholar]
  23. Liu, Z.; Du, S.; Wu, Y.; Sha, M. Anti-Interrupted Sampling Repeater Jamming Method for Interpulse and Intrapulse Frequency-Agile Radar. J. Radars 2022, 11, 301–312. [Google Scholar]
  24. Quan, Y.; Fang, W.; Sha, M.; Ruan, F. Present Situation and Prospects of Frequency Agility Radar Waveform Countermeasures. Syst. Eng. Electron. 2021, 43, 3126–3136. [Google Scholar]
  25. Ren, L.; Long, T.; Yuan, H. Signal Processing and Parameter Design of HPRF Pulsed-Doppler Stepped Frequency Radar. Acta Electron. Sin. 2007, 09, 1630–1636. [Google Scholar]
  26. Liu, Z.; Wang, X.; Liu, J.; Wang, G. Jamming Technique of Interrupted-Sampling and Periodic Repeater Based on Digital Radio Frequency Memory. Acta Armamentarii 2008, 29, 405–410. [Google Scholar]
  27. Wang, X.; Liu, J.; Zhang, W.; Fu, Q. Mathematical Principle of Interrupted-Sampling Repeater Jamming. Sci. Sin. 2006, 08, 891–901. [Google Scholar]
  28. Liu, Z. Jamming Technique for Countering LFM Pulse Compression Radar Based on Digital Radio Frequency Memory. Ph.D. Thesis, National University of Defense Technology, Changsha, China, 2006. [Google Scholar]
  29. Zhang, J.; Mu, H.; Wen, S.; Li, Y. Anti-Intermittent Sampling Repeater Jamming Method Based on LFM Segmented Pulse Compression. J. Electron. Inf. Technol. 2019, 41, 1712–1720. [Google Scholar]
  30. Otsu, N. A Threshold Selection Method from Gray-Level Histograms. IEEE Trans. Syst. Man Cybern. 1979, 9, 62–66. [Google Scholar] [CrossRef] [Green Version]
  31. Kapur, J.N.; Sahoo, P.K.P.; Wong, A.K.C. A new method for grey-level picture thresholding using the entropy of the histogram. Comput. Vis. Graph. Image Process. 1985, 29, 273–285. [Google Scholar] [CrossRef]
  32. Sezgin, M.; Sankur, B. Survey Over Image Thresholding Techniques and Quantitative Performance Evaluation. J. Electron. Imaging 2004, 13, 146–168. [Google Scholar]
  33. Wang, S.; Haralick, R.M. Automatic multithreshold selection. Comput. Vis. Graph. Image Process. 1984, 25, 46–67. [Google Scholar] [CrossRef]
  34. Yen, J.; Chang, F.; Chang, S. A New Criterion for Automatic Multilevel Thresholding. IEEE Trans. Image Process. 1995, 4, 370–378. [Google Scholar]
  35. Liao, P.; Chen, T.; Chung, P. A Fast Algorithm for Multilevel Thresholding. J. Inf. Sci. Eng. 2001, 17, 713–727. [Google Scholar]
  36. Erdélyi, A. Asymptotic Representations of Fourier Integrals and the Method of Stationary Phase. J. Soc. Ind. Appl. Math. 1955, 3, 17–27. [Google Scholar] [CrossRef]
  37. Mulgrew, B. The Stationary Phase Approximation, Time-Frequency Decomposition and Auditory Processing. IEEE Trans. Signal Process. 2014, 62, 56–68. [Google Scholar] [CrossRef] [Green Version]
  38. Vizitiu, L.; Anton, F.; Popescu, F.; Iubu, G. The synthesis of some NLFM laws using the stationary phase principle. In Proceedings of the 2012 10th International Symposium on Electronics and Telecommunications, Timisoara, Romania, 15–16 November 2012. [Google Scholar]
  39. Vetter, W.J. Matrix calculus operations and Taylor expansions. SIAM Rev. 1973, 15, 352–369. [Google Scholar] [CrossRef]
  40. Shidama, Y. The Taylor expansions. Formaliz. Math. 2004, 12, 195–200. [Google Scholar]
  41. Heald, M.A. Rational approximations for the Fresnel integrals. Math. Comput. 1985, 44, 459–461. [Google Scholar] [CrossRef]
  42. Umul, Y.Z. Equivalent functions for the Fresnel integral. Opt. Express 2005, 13, 8469–8482. [Google Scholar] [CrossRef]
Figure 1. Block diagram of intermittent sampling direct jamming model.
Figure 1. Block diagram of intermittent sampling direct jamming model.
Remotesensing 15 03080 g001
Figure 2. Block diagram of intermittent sampling repeater jamming model.
Figure 2. Block diagram of intermittent sampling repeater jamming model.
Remotesensing 15 03080 g002
Figure 3. Block diagram of intermittent sampling cyclic jamming model.
Figure 3. Block diagram of intermittent sampling cyclic jamming model.
Remotesensing 15 03080 g003
Figure 4. Schematic diagram of intra-pulse random quadrature FM and inter-pulse phase coherence waveform.
Figure 4. Schematic diagram of intra-pulse random quadrature FM and inter-pulse phase coherence waveform.
Remotesensing 15 03080 g004
Figure 5. Full pulses multilevel maximum interclass variance algorithm.
Figure 5. Full pulses multilevel maximum interclass variance algorithm.
Remotesensing 15 03080 g005
Figure 6. Schematic diagram of sub-pulses multi-level maximum inter-class variance algorithm.
Figure 6. Schematic diagram of sub-pulses multi-level maximum inter-class variance algorithm.
Remotesensing 15 03080 g006
Figure 7. Time domain diagram of the transmit signal and three intermittent sampling interferences. (a) Time domain diagram of transmit signal. (b) Time domain diagram of intermittent sampling direct jamming. (c) Time domain diagram of the repeater jamming. (d) Time domain diagram of the cyclic jamming.
Figure 7. Time domain diagram of the transmit signal and three intermittent sampling interferences. (a) Time domain diagram of transmit signal. (b) Time domain diagram of intermittent sampling direct jamming. (c) Time domain diagram of the repeater jamming. (d) Time domain diagram of the cyclic jamming.
Remotesensing 15 03080 g007
Figure 8. Time-frequency diagram of the transmit signal and three intermittently sampled interferences. (a) Time-frequency diagram of the transmit signal. (b) Time-frequency diagram of the direct jamming. (c) Time-frequency diagram of the repeater jamming. (d) Time-frequency diagram of the cyclic jamming.
Figure 8. Time-frequency diagram of the transmit signal and three intermittently sampled interferences. (a) Time-frequency diagram of the transmit signal. (b) Time-frequency diagram of the direct jamming. (c) Time-frequency diagram of the repeater jamming. (d) Time-frequency diagram of the cyclic jamming.
Remotesensing 15 03080 g008aRemotesensing 15 03080 g008b
Figure 9. Frequency domain diagram of the transmit signal and three intermittent sampling jamming. (a) Amplitude–frequency diagram of transmitted signal. (b) Amplitude–frequency diagram of direct jamming. (c) Amplitude–frequency diagram of repeater jamming. (d) Amplitude–frequency diagram of cyclic jamming.
Figure 9. Frequency domain diagram of the transmit signal and three intermittent sampling jamming. (a) Amplitude–frequency diagram of transmitted signal. (b) Amplitude–frequency diagram of direct jamming. (c) Amplitude–frequency diagram of repeater jamming. (d) Amplitude–frequency diagram of cyclic jamming.
Remotesensing 15 03080 g009
Figure 10. Pulse compression domain diagram of the transmit signal and three intermittent sampling jamming. (a) Pulse compression diagram of the transmit s diagram of transmitted signal. (b) Pulse compression diagram of the transmit s diagram of direct jamming. (c) Pulse compression diagram of the transmit s diagram of repeater jamming. (d) Pulse compression diagram of the transmit s diagram of cyclic jamming.
Figure 10. Pulse compression domain diagram of the transmit signal and three intermittent sampling jamming. (a) Pulse compression diagram of the transmit s diagram of transmitted signal. (b) Pulse compression diagram of the transmit s diagram of direct jamming. (c) Pulse compression diagram of the transmit s diagram of repeater jamming. (d) Pulse compression diagram of the transmit s diagram of cyclic jamming.
Remotesensing 15 03080 g010
Figure 11. Schematic diagram of intra-pulse random quadrature FM waveform.
Figure 11. Schematic diagram of intra-pulse random quadrature FM waveform.
Remotesensing 15 03080 g011
Figure 12. Identification results of the full-pulses multilevel maximum interclass variance algorithm at JSR = 10 dB.
Figure 12. Identification results of the full-pulses multilevel maximum interclass variance algorithm at JSR = 10 dB.
Remotesensing 15 03080 g012
Figure 13. Variance and threshold plots for the first sub-pulse using the full-pulses and sub-pulses OTSU. (a) The first sub-pulse uses the full-pulses shared threshold algorithm variance graph. (b) The first sub-pulse uses the sub-pulses adaptive threshold algorithm variance graph.
Figure 13. Variance and threshold plots for the first sub-pulse using the full-pulses and sub-pulses OTSU. (a) The first sub-pulse uses the full-pulses shared threshold algorithm variance graph. (b) The first sub-pulse uses the sub-pulses adaptive threshold algorithm variance graph.
Remotesensing 15 03080 g013
Figure 14. Variance and threshold plots for the third sub-pulse using the full-pulses and sub-pulses OTSU. (a) The third sub-pulse uses the full-pulses shared threshold algorithm variance graph. (b) The third sub-pulse uses the sub-pulses adaptive threshold algorithm variance graph.
Figure 14. Variance and threshold plots for the third sub-pulse using the full-pulses and sub-pulses OTSU. (a) The third sub-pulse uses the full-pulses shared threshold algorithm variance graph. (b) The third sub-pulse uses the sub-pulses adaptive threshold algorithm variance graph.
Remotesensing 15 03080 g014
Figure 15. Identification results of sub-pulses multi-level maximum inter-class variance algorithm at JSR = 10 dB.
Figure 15. Identification results of sub-pulses multi-level maximum inter-class variance algorithm at JSR = 10 dB.
Remotesensing 15 03080 g015
Figure 16. Identification results of full-pulses multi-level maximum inter-class variance algorithm at JSR = 20 dB.
Figure 16. Identification results of full-pulses multi-level maximum inter-class variance algorithm at JSR = 20 dB.
Remotesensing 15 03080 g016
Figure 17. Identification results of sub-pulses multi-level maximum inter-class variance algorithm at JSR = 20 dB.
Figure 17. Identification results of sub-pulses multi-level maximum inter-class variance algorithm at JSR = 20 dB.
Remotesensing 15 03080 g017
Table 1. Theparameters of the waveform and jamming.
Table 1. Theparameters of the waveform and jamming.
ParameterValueParameterValue
T p 41  μ sB400 MHz
τ J μ s B sub 10 MHz
τ d 0.2  μ s Δ f 40 MHz
T J 10.2  μ sP4
Table 2. Confusionmatrix for full-pulses multilevel maximum interclass variance algorithm.
Table 2. Confusionmatrix for full-pulses multilevel maximum interclass variance algorithm.
Mode one78110
Mode two47033
Mode three233540
Mode four43073
Mode oneMode twoMode threeMode four
Table 3. Signal patterns corresponding to the confusion matrix.
Table 3. Signal patterns corresponding to the confusion matrix.
Mode oneTarget echo signal
Mode twoTarget echo signal plus Direct Jamming
Mode threeTarget echo signal plus Repeater Jamming
Mode fourTarget echo signal plus Cyclic Jamming
Table 4. Confusionmatrix for sub-pulses multilevel maximum interclass variance algorithm.
Table 4. Confusionmatrix for sub-pulses multilevel maximum interclass variance algorithm.
Mode one75221
Mode two47402
Mode three23714
Mode four43073
Mode oneMode twoMode threeMode four
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhan, H.; Wang, T.; Guo, T.; Su, X. An Anti-Intermittent Sampling Jamming Technique Utilizing the OTSU Algorithm of Random Orthogonal Sub-Pulses. Remote Sens. 2023, 15, 3080. https://doi.org/10.3390/rs15123080

AMA Style

Zhan H, Wang T, Guo T, Su X. An Anti-Intermittent Sampling Jamming Technique Utilizing the OTSU Algorithm of Random Orthogonal Sub-Pulses. Remote Sensing. 2023; 15(12):3080. https://doi.org/10.3390/rs15123080

Chicago/Turabian Style

Zhan, Haihong, Tao Wang, Tai Guo, and Xingde Su. 2023. "An Anti-Intermittent Sampling Jamming Technique Utilizing the OTSU Algorithm of Random Orthogonal Sub-Pulses" Remote Sensing 15, no. 12: 3080. https://doi.org/10.3390/rs15123080

APA Style

Zhan, H., Wang, T., Guo, T., & Su, X. (2023). An Anti-Intermittent Sampling Jamming Technique Utilizing the OTSU Algorithm of Random Orthogonal Sub-Pulses. Remote Sensing, 15(12), 3080. https://doi.org/10.3390/rs15123080

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop