Radar Reconnaissance Pulse-Splitting Modeling and Detection Method

: When radar receivers adopt digital channelization, it is prone to generating a cross-channel split signal, the rabbit ear effect, and a transition band repeated signal, leading to errors in radar signal sorting or identification. The pulse-splitting model and detection method proposed in this paper can model split pulses and identify them in radar pulse streams, facilitating the merging of split pulses to enhance sorting and identification performance. Firstly, the mechanism of splitting pulse generation is deeply analyzed, and the splitting site theory is proposed. Then, the split pulse signal model and the split pulse number statistical model based on geometric distribution are constructed, which are used to guide the construction of simulation data of split pulse flow with different characteristics. Furthermore, a time-domain convergence degree (TCD) index is proposed to characterize the pulse split phenomenon. At the same time, in order to avoid a large number of threshold searching problems in pulse-splitting detection, an empirical formula for the pulse-splitting detection threshold based on the TCD is given to quickly determine whether there is a pulse train split problem. The selected measured radar pulse stream is verified to follow a geometric distribution at a significance level of 0.05. The proposed method achieved a detection accuracy of at least 99.55% on the simulation dataset and at least 95.68% on the experimental dataset, validating the rationality of pulse-splitting modeling and the effectiveness of the detection method.


Introduction
The reliable and accurate detection of non-cooperative radar pulses is one of the key tasks of radar reconnaissance.The detection results are crucial for estimating radar frequency, pulse width (PW), and other parameters.The quality of parameter estimation directly impacts subsequent tasks such as sorting recognition and jamming decision making [1,2].
For the detection of radar pulse signals in reconnaissance signals, the basic principle is to utilize the differences in statistical characteristics between random noise and signals in various transform domains.The most classic algorithm is the time-domain energy detection method, whose performance is significantly affected by noise [3].Xi Yimin [4] and Zou Simin [5] achieved pulse detection under low signal-to-noise ratio (SNR) conditions by computing the delayed autocorrelation of the received signal.Li Cheng [2] utilizes the Difference of Boxes (DOB) filter for pulse edge detection, and its detection performance surpasses that of an energy detection method [3].On the basis of Li Cheng [2], Xinqun Liu [6] utilizes the fine-tuned Short-Time Fourier Transform (STFT) to refine the rough pulse edge, and then the Constant False Alarm Rate (CFAR) method is used to implement adaptive threshold detection in the frequency domain, which is better than the DOB pulse detection algorithm [2] and the traditional energy detection method [3].Wenmao Luo [7] utilizes the disparity in state variables between two generalized time-scale transformed Duffing oscillators to detect pulses, thereby addressing the issue of weak pulse signal detection.Existing radar pulse detection algorithms [2][3][4][5][6][7] focus on detecting individual pulses, but prior to this paper, no one has studied the detection of split radar pulses.
There are several factors that affect pulse detection performance, and CAO Chunyan [8] has conducted a comprehensive analysis focusing on SNR, pulse overlap, and processing channels.However, the complex and variable modern electromagnetic environment results in a significant difference in correlation power, various modulation forms, and a serious time-domain overlap in the receiver [9].The received pulse signal is prone to splitting due to multipath effects, equipment impedance mismatch, and colored noise [10].The use of ultra-wideband radar [11] with high-resolution detection capability will cause a serious split phenomenon in both the time-domain and frequency-domain detection methods based on digital channalization, which is not conducive to reliable detection and analysis of ultra-wideband radar.In particular, Guangliang Zou [12], Ying Liu [13] and Haiqing Jiang [14] point out that wideband pulse amplitude (PA) is susceptible to noise, fluctuates greatly in its interior, and the intra-pulse data segment is prone to be lower than the minimum detection sensitivity, thus causing pulse splitting.
In view of the pulse split phenomenon, Guangliang Zou [12] uses the correlation of adjacent pulse description words (PDW) in the time-frequency domain and the consistency of intra-pulse modulation types to solve the merging of split PDWs.Ying Liu [13] considers the misjudgment of intra-pulse modulation type caused by insufficient linearity when there are fewer intra-pulse sample points.It also removes the constraint condition of the intra-pulse modulation type mentioned in reference [12].Ying Liu [13] designs the merging conditions for split pulses under different modulation types.The corresponding parameters after merging are obtained by averaging the PA and frequency-modulated slope of the split pulses that satisfy the merging condition.Then, the corresponding expression for the combination of multiple split pulses in the time-frequency domain parameter is constructed.Different from Ying Liu [13], Haiqing Jiang [14] considers the actual measurement error and introduces the error coefficient into the LFM parameters and frequency domain merging conditions.Liming Tang [15] proposes a serial combination judgment method that judges the signal type, the time of arrival (TOA) difference of adjacent pulses, the difference of pulse slope, and the center frequency in turn.When the combination condition is met step by step, the pulse combination process is carried out with the weight of the split PW.Chao Wang [16] provides a theoretical basis for constructing the split pulse signal model through mathematical modeling of the split linear frequency modulation signal.
The pulse-splitting model proposed in this paper analyzes the mechanisms behind three pulse-splitting phenomena: the cross-channel split signal, the rabbit ear effect, and the transition band repeated signal.Furthermore, this paper establishes correspondences between the three pulse-splitting phenomena and the temporal relationships among the split pulses: the temporal relationship of cross-channel split signal pulses is overlapping, the temporal relationship of rabbit ear effect split pulses is independent, and the temporal relationship of transition band repeated signal split pulses is inclusive.Subsequently, a split pulse signal model covering the above three pulse-splitting phenomena is proposed.Finally, this paper also discovers and verifies that the statistical distribution of split pulse counts follows a geometric distribution.Wu Zhenyu [17] and Gui Youlin [18] conducted a mechanistic study on the three scenarios of pulse splitting described above.Guangliang Zou [12], Ying Liu [13], Haiqing Jiang [14], and Liming Tang [15] outlined the merging conditions and methods of split pulses.However, none of them, including Guangliang Zou [12], Ying Liu [13], Haiqing Jiang [14], Liming Tang [15], Wu Zhenyu [17], and Gui Youlin [18], has mathematically modeled the phenomenon of pulse splitting.Wang Chao [16]'s proposed pulse-splitting model only addresses independently timed split pulses and does not consider other splitting scenarios.Additionally, none of these researchers have investigated the statistical regularity of split pulse numbers.
The advantages and disadvantages of existing radar pulse-splitting methods are shown in Table 1.Aiming to address the lack of research on the pulse-splitting model and detection method, the main contributions of the proposed pulse-splitting modeling and detection method are as follows: 1.
The split-site theory is proposed to describe the mechanism of split pulses in different ways.Then, based on the split pulse law of the experimental data, the split pulse signal model and the statistical model of the number of split pulses using the geometric distribution are constructed, respectively.Subsequently, simulation data of split pulses with various characteristics are generated using the proposed model.

2.
A convergence degree in the time domain is proposed to represent pulse splitting.Additionally, based on the TCD, an empirical formula is provided for the threshold of pulse-splitting detection.This formula helps to avoid the low detection efficiency that can be caused by a large number of threshold searches in pulse splitting detection, enabling the reliable assessment of pulse splitting in a radar pulse stream.
They have not studied radar pulse-splitting detection algorithms.[17,18] They have conducted mechanistic studies on the three scenarios of pulse splitting.
They have not mathematically modeled the phenomenon of pulse splitting.
They have not studied the statistical regularities of split pulse numbers. [16] The proposed pulse splitting model only considers mutually independent split pulses.
The model does not consider other splitting scenarios, and there has been no investigation into the statistical regularities of split pulse numbers.
This paper presents the concept of split pulse fusion processing and signal parameter estimation.It can improve the sorting and identification performance, supports the efficient and reliable implementation of radar reconnaissance systems, and has significant engineering application value.

Pulse-Splitting Model
When receiving the non-cooperative radar pulses, the pulse-splitting phenomenon will be caused by the environment and the detection equipment, which makes the reconnaissance process difficult.By analyzing the split law of the experimental data, the generation mechanism of the split pulse is expounded according to the problems brought by digital channelization, and the split site theory is proposed to describe the split pulse of different split modes.Then, a statistical model is constructed for the number of split pulses based on the geometric distribution.The split site theory is then incorporated to guide the construction of split pulse flows with different characteristics.

Split Pulse Generation Mechanism
In this paper, the phenomenon of pulse splitting caused by the receiver device is primarily considered.The digital channelization receiver device divides the receiver bandwidth into multiple subchannels and outputs signals of different frequencies through a set of filters that can cover the entire frequency band.Digital channelization has many advantages, such as the ability to handle wideband signals and signals arriving at different frequencies simultaneously.However, there are still problems such as the cross-channel signal split phenomenon, the rabbit ear effect, and the repetition of signals in the transition band, which are collectively referred to as the pulse split phenomenon [17,18].In the following, the split-site theory is proposed, and the mechanism of the above pulse split phenomenon is described, respectively.

• Split Site Theory
As shown in Figure 1, the original pulse is located on the upper time axis, while the split pulse generated by the original pulse is located on the lower time axis.The split sites 1 to 10 represent the leading edge and the trailing edge of the split pulse, respectively, for the original pulse; the number of split pulses is 5.The relationship between the number of split sites and the number of split pulses of a single original pulse can be determined as follows: According to the experimental data of split pulses, the following provisions are made: 1.
The split sites in a complete pulse are not allowed to be reused; that is, a split site can only be used as the leading edge or trailing edge of a split pulse.

2.
The current split site can only generate a split pulse with the next split site.

3.
Split pulses are only allowed to be generated in a single complete pulse.4.
The first split site must be located at the leading edge of the complete pulse, and the last split site must be located at the trailing edge of the complete pulse.

5.
All split sites within a complete pulse must be used.
Let any adjacent split pulse split by the same original pulse with earlier arrival time be the n number subpulse and the later one be the n + 1 number subpulse.The n number subpulse arrival time front is TOA n , the front site is a, the back edge of the arrival time is TOE n , and the back edge site is a + x.The n + 1 subpulse arrival time front is TOA n+1 , the front site is b, the arrival time back edge is TOE n+1 , and the back edge site is b + y, a, b, x, y ∈ N.

• Cross-channel splitting signal
When channelizing a wideband signal, such as a linear frequency modulation signal, its frequency band often covers multiple subchannels, potentially resulting in the detection of split pulses with different TOAs in adjacent channels.The split pulses detected by multiple channels are plotted on the same graph, revealing the phenomenon known as the cross-channel signal split with the split sites located in the time domain corresponding to the transition bands between subchannels.The schematic diagram and the time domain diagram of the experimental data are shown in Figures 2 and 3.  1 ∼ M in the upper timeline of Figure 4 represent the split sites of the original pulse.The split sites are denoted as a, b, a + x and b + y, which correspond to the overlapping split sites of the leading edge TOA n , TOA n+1 , and the trailing edge TOE n , TOE n+1 of the nth and n + 1th split pulse in the lower timeline.The timing relationship between the leading and trailing edges of adjacent split pulses and the size of split sites must satisfy the following condition:

• Rabbit Ear Effect
Due to the abrupt change in PA and phase near the rising and falling edges of the pulse signal, channelization will spread the energy of the abrupt change to several subchannels near the actual output subchannel of the signal, resulting in a waveform similar to the transient output of a rabbit ear.The rabbit ear effect appears near the beginning or end of the main signal.The PW is much smaller than the main signal PW, and the PA is slightly smaller as well.Additionally, the adjacent split pulses are independent in the time domain.The schematic diagram and the time-domain diagram of the experimental data are shown in Figures 5 and 6.Similar to Figure 4, the upper timeline of Figure 7 shows independent split sites labeled as a, b, a + x, and b + y.These split sites correspond to the leading edge TOA n , TOA n+1 , and trailing edge TOE n , TOE n+1 in the lower timeline.The timing relationship of adjacent split pulses and the size of split sites satisfy the following:

• Transition Band Repeated Signal
When the signal frequency is located at or near the intersection of adjacent subchannels, it may appear almost simultaneously in two adjacent subchannels, leading to repeated signals in the transition band.Its PA is slightly smaller than that of the real signal, and its leading edge and trailing edge arrival times are approximately the same as those of the real signal.In the time domain, it is mostly encompassed by the real signal.The schematic diagram and time-domain plot of experimental data are shown in Figures 8 and 9.  Similar to Figure 4, the upper timeline of Figure 10 shows mutually contained split sites, which are represented by a, b, a + x, and b + y.These split sites correspond to the leading edge TOA n , TOA n+1 , and trailing edge TOE n , TOE n+1 in the lower timeline.The timing relationship of adjacent split pulses and the size of split sites must satisfy the following:

• Split PDW Value Range
The range of PDW values corresponding to the split original pulse is given below.
(1) Carrier Frequency The carrier frequency at the split site has a value range of generating M random values f site that follow a uniform distribution. (

Split Pulse Signal Model
When the original pulse splits in the manner described above, this section provides a mathematical description of the pulse-splitting phenomenon and develops a model for the split pulse signal by constructing the original pulse signal model.
Let the pulse number of the observed original pulse stream be 1, • • • , e, • • • , E, and the pulse envelope shape is assumed to be ideal, resembling a rectangle.The type of modulation within the pulse is not limited.Therefore, the signal model of the e th original pulse can be expressed as shown below: where rect(•) represents the unit rectangular window function, A e denotes the PA, TOA e refers to the TOA, τ e represents the PW, ω 0 represents the angular frequency of the signal, and ϕ(t) represents the modulation phase of the signal pulse.
A splitting pulse signal model is then constructed based on the original pulse signal model.The sequence number of N splitting pulses generated by the eth original pulse is 1, • • • , n, • • • , N, the intra-pulse modulation type is consistent with the original pulse.The signal model of the nth split pulse is expressed as: where A e,n represents the PA of the splitting pulse, TOA e,n denotes the TOA of the splitting pulse, and τ e,n refers to the width of the splitting pulse.

Statistical Model of Split Pulse Number Based on Geometric Distribution
Through the analysis of the experimental data, it is known that pulses transmitted by the same radar and pulses transmitted by different radars will produce a different number of split pulses when passing through the receiver.Count the number of original pulses corresponding to different split pulse numbers from the same radar: The value range of the number of split pulses corresponding to the original pulse of this radar is represented by 1, • • • , l + 1. F 1 represents the number of original pulses without splitting, and F l+1 represents the number of original pulses with the maximum number of split pulses l + 1.
It can be inferred that the number of split pulses follows a geometric distribution based on statistical analysis and practical significance.Therefore, a statistical model of the number of split pulses based on the geometric distribution is constructed.First of all, the probabilities of success and failure of events in the context of the pulse split problem, which are based on the n-weight Bernoulli experiment, need to be explained.The random variable X is defined to follow a Bernoulli distribution with parameter p ∈ (0, 1].X can take on the values 1 and 0, representing whether the split pulse can form a complete pulse or not.p ∈ (0, 1] indicates the probability that each split pulse can form a complete pulse.The distribution law is as follows: Here, p ∈ (0, 1), P(X = 1) represents the probability that the split pulse can form a complete pulse, and P(X = 0) represents the probability that the split pulse cannot form a complete pulse.In addition, when p = 1, the split pulse is the complete pulse.Based on the above setting, let L random variables X 1 , • • • , X l , X l+1 , • • • , X L be independent and identically distributed according to the Bernoulli distribution with parameter p ∈ (0, 1].The random variable Y is defined as follows: for any l ∈ F, Y = l is true if and only if indicates that the 1 ∼ l-th split pulse cannot form a complete pulse, and X l+1 = 1 indicates that the l + 1-th split pulse can form a complete pulse.Then, the event Y = l means that its neighboring l split pulses cannot form a complete pulse until the l + 1 th split pulse arrives.When the l + 1 th split pulse appears, it successfully forms a complete pulse with the first l split pulses exactly.The distribution law of the random variable Y is as follows: • Effect of PW and PA on Split Pulse Number As shown in Figure 11, the p values of different radars vary due to differences in parameters such as PW and PA.The higher the probability p that the split pulse constitutes a complete pulse, the larger the average slope |k| of the curve, and the faster the curve declines simultaneously.The minimum P in f = 1.5 × 10 −4 close to 0 is set as the discriminative threshold for determining the maximum number of split pulses needed to successfully construct a complete pulse.P(Y = l − 1) > P in f if and only if P(Y = l) < P in f .In this case, the abscissa l + 1 represents the maximum number of split pulses required for the radar to form a complete pulse.The calculated results are shown in Table 2.It is clear that the larger the value of p, the fewer split pulses are needed to form a complete pulse.Then, the influence of PW on the splitting of pulses is discussed.The experimental data are used to calculate the average number of splitting pulses corresponding to different PWs.The results are shown in Figure 12, which indicates that the average number of splitting pulse increases as the PW increases.According to Table 2, the smaller the probability p of the split pulse forming a complete pulse, the flatter the corresponding curve of the geometric distribution law becomes.Finally, the influence of PA on pulse splitting is analyzed.The experimental data of pulse splitting from the same radar are selected to calculate the average number of pulse splitting with the same PW but different PA.As shown in Figure 13, the average number of pulse splitting gradually decreases with the increase in PA until it remains unchanged after reaching −82 dB.Assuming that the ambient noise of the receiving device does not change over time, the radar pulse SNR is primarily influenced by the PA.The larger the PA, the less impact pulse splitting has on the environment, and the smaller the maximum number of pulse splitting required to form a complete pulse.Combined with Table 2, it can be observed that as the probability p of the split pulse forming a complete pulse increases, the curve of the geometric distribution law becomes steeper.

Building the Split Pulse Flow
According to the generation mechanism of split pulse, the split pulse signal model, and the statistical model of split pulse number, a split pulse flow is constructed to verify the performance of the split pulse detection method described below.The pulse flow density, pulse splitting degree, pulse staggering degree, and pulse recurrence interval (PRI) modulation complexity are defined as characteristics of pulse flow to describe the changing pattern of split pulse flow in a complex environment.

Pulse Flow Density
Pulse flow density ρ is defined as the number of pulses received by the receiver per unit time and can be expressed as shown below: where N p is the number of pulses received within a given timeframe, and T f represents the duration of that timeframe.The pulse flow density can be reflected by the length of the pulse repetition interval.U radars are assumed to transmit pulse streams with densities The general formula for calculating pulse flow densities is as follows: The PRI of the radar u pulse stream is PRI u , and the pulse flow density ρ u is the reciprocal of PRI u .The split PDW datasets with different pulse flow densities are generated by combining the pulse-splitting model with other pulse flow characteristics unchanged.

Pulse-Splitting Degree
The degree of pulse splitting P(X = 0) = q, q ∈ [0, 1) is defined as the probability that split pulses cannot form a complete pulse.As can be seen from Figure 11 and Table 2, the larger q is, the higher the maximum number of split pulses required to form a complete pulse.Suppose that when receiving the pulse streams arriving successively from the same radar, the corresponding pulse splitting degree is The split PDW datasets with varying degrees of pulse splitting are generated by combining the pulse-splitting model with other unchanged pulse flow characteristics.

Pulse Staggering Degree
The pulse staggering degree RN is defined as the number of radars receiving the pulse stream within a given timeframe.The higher the number of radars, the greater the pulse staggering degree.Set the staggering degree corresponding to the arriving pulse streams as Keep the characteristics of other pulse streams unchanged and combine the pulse-splitting model to generate split PDW datasets with different pulse staggering degrees.

PRI Modulation Complexity
The PRI modulation complexity c is defined as the modulation type of the PRI corresponding to the radar pulse stream within a basic timeframe.The higher the complexity of the modulation type, the higher the PRI modulation complexity.Set the modulation complexity of the PRI corresponding to the incoming pulse streams as Keep the characteristics of other pulse streams unchanged and combine the pulse-splitting model to generate split PDW datasets with varying PRI modulation complexities.

Integrated Pulse Flow Characteristics
Let the pulse flow density, pulse-splitting degree, pulse-staggering degree, and PRI modulation complexity, which change simultaneously in a time sequence, be: Combined with the pulse-splitting model, a split PDW dataset is generated, which exhibits multiple pulse flow characteristics changing simultaneously.

Reference Label Generation Method
The generation method for timeframes with reference labels is described as follows: First, the original PDW dataset is generated, and then different numbers of split pulses corresponding to the original PDW are generated according to the pulse-splitting model.The original PDW data with split pulse numbers is processed in frames, and the timeframes are labeled according to Table 3.

Timeframe Feature Label
There exists at least one original PDW with a split pulse number greater than 1. 1 The split pulse number corresponding to the original PDW is at most 1. 0 No pulse is received.−1

Split Pulse Detection
A method for detecting splitting pulses based on a split pulse model is proposed.Firstly, the pulse stream is divided into frames of equal length, and then the degree of convergence in the time domain is introduced to characterize the splitting of pulses.An empirical formula for determining the detection threshold based on the convergence degree in the time domain is provided to obtain an appropriate threshold for each timeframe.Formed by calculating the derivatives of the convergence degree in the time domain at adjacent thresholds, the convergence derivative matrix is used to efficiently identify any potential splitting problems in the pulse flow.The flow of the algorithm is shown in Figure 14.

Select the split pulse detection index
Calculate the time domain convergence degree for different thresholds, and construct the convergence derivative matrix Determine the presence of split pulses in the pulse stream according to the convergence derivative matrix Frame processing pulse stream Obtain the optimal detection threshold by using the empirical formula for pulse splitting detection threshold Introduce the time-domain convergence degree

Framing the Pulse Stream
The pulse stream with a timeframe length T f is the fundamental processing unit of the split pulse detection algorithm and is calculated as follows: where T PDW is the duration of interception time for PDW and N f is the number of timeframes.The processing of pulse stream framing can reduce the algorithm's processing time, enabling the localization of split pulse detection within each timeframe and facilitating the subsequent processing of timeframes containing split pulses.

TCD
In order to determine whether there is a split pulse in each frame, this paper introduces the full-window standard deviation (FSTD) based on the standard deviations.It then improves the single-box moving step of FSTD to a multi-box moving step of independent full-window standard deviations (IFSTD).Finally, the concept of TCD is introduced, which is a statistical distribution.
The standard deviation S can reflect the dispersion of a dataset and is expressed as follows: where represents each sample of the dataset, X denotes the sample mean, and V represents the sample size.The smaller the standard deviation, the more stable the dataset.As shown in Figure 15, the standard deviation of the window length for a fixed value w (referred to as the FSTD) is calculated sequentially for each box in the histogram, starting from the left and moving toward the right along the horizontal coordinate.When each FSTD is calculated, the window will shift to the right by the length of a single box.Here is how to calculate each FSTD: where the pulse cumulant per box (histogram height) is Then, by improving the single-move step length of FSTD to the window length, TCD is introduced into the calculation to represent the convergence of the pulse stream parameter distribution histogram.As shown in Figure 16, the FSTD is calculated from left to right along the abscissa with a fixed window length and independent coverage of the histogram, which is defined as IFSTD.When calculating each IFSTD, the window will shift w elements to the right.When calculating the last IFSTD, the window will automatically truncate at the end of the abscission value.IFSTD is calculated as follows: IFSTD reflects the standard deviation of each window length in the parametric histogram, and the results obtained are more independent and significant compared to FSTD.Finally, IFSTD is averaged to obtain the TCD: The histogram window number is N w = ⌈ N box w ⌉, and the meaning of the remaining variables is consistent with the description provided above.

Split Pulse Detection Indicator
In order to perform split pulse detection on the received PDWs, the following time domain indicators are selected to draw the histogram and calculate the TCD.

• PW
The PW (τ) is the duration for which the pulse can reach its maximum value.

• Duty Cycle (DC)
DC (D) is the ratio of PW to the adjacent pulse time interval (DTOA): For split pulses, the calculation of the DC can be divided into two cases: pulse independence and pulse non-independence.For independent splitting pulses, as shown in Figure 17, the sequential relationship between adjacent splitting pulses is TOA d+1 > TOA d + τ d , and the DC calculation method is consistent with Equation ( 19): For splitting pulses that are not independent of each other, as shown in Figure 18, the sequential relationship between adjacent splitting pulses is TOA n+1 < TOA n + τ n , n = d, z, and the DC D d = D z = 1 of pulses d and z is specified.In particular, as shown in Figure 19, if the current timeframe has only one pulse, the DC of the pulse is specified as the ratio of the PW to the current timeframe length The degree of convergence of PW and DC in the time domain over a period of time can indicate the splitting of the pulse train.Therefore, these parameters are used as indicators to measure whether the pulse train is split.
TCD can reflect the extent of overlap and splitting of the pulse flow.Take PW as an example.As shown in Figure 20, for a complete pulse flow, the TCD calculated from the PW distribution histogram is lower.As shown in Figures 21 and 22, the TCD calculated from the PW distribution histogram is larger for overlapping or splitting pulse flows.

Empirical Formula of Splitting Pulse Detection Threshold
The difficulty of the algorithm lies in selecting the appropriate detection threshold for each PDW frame to achieve high accuracy in split pulse detection.In order to clarify the measurement object of the detection threshold ths, the missing pulse length T g is defined as the time difference between the arrival time of the back edge of the first arrival pulse and the arrival time of the front edge of the second arrival pulse [16].It can be expressed as follows: where TOA n+1 is the arrival time of the front edge of the second arrival pulse, and TOE n is the arrival time of the back edge of the first arrival pulse.As shown in Figure 23, it can be seen that within a given timeframe, the number of missing pulse lengths is one less than the number of splitting pulses.In order to facilitate matrix calculation, the first splitting pulse and the second splitting pulse share the same missing pulse length Tg 1 = Tg 2 to ensure that the number of missing pulse lengths is consistent with the number of splitting pulses.For split pulses, Tg is calculated differently when the pulses are independent of each other and when they are not independent of each other.For independent splitting pulses, as shown in Figure 17, Tg is calculated in the same way as Equation (20).For split pulses that are not independent of each other, as shown in Figure 18, the missing pulse length of pulses t and z is specified as Tg t = Tg z = 0.In particular, as shown in Figure 19, if the current timeframe has only one pulse, the missing pulse length Tg i is equal to 0.
Calculate the absolute difference between the mean missing pulse length and the mean PW of each PDW frame, multiply this absolute value by the mean duty ratio, and round the result to obtain the empirical formula for the detection threshold of split pulses.The mean values of the missing pulse length (Tg), PW (τ), DC (D) and the detection threshold for split pulses (ths max ) are calculated as follows: where round(x) is the rounding function used to ensure that the reasonable detection threshold ths max of each frame PDW is an integer.The reasonable detection threshold of each frame data is ths max .Tg i , τ i , DTOA i and D i represent the missing pulse length, PW, adjacent pulse time interval, and DC per frame, respectively.Additionally, n b denotes the number of pulses per frame.When there are no split pulses within the timeframe, Tg, τ and |Tg − τ| are relatively large.As the number of split pulses increases, both Tg and τ gradually decrease, and |Tg − τ| decreases compared to the timeframe with only complete pulses.Coefficient D ensures that ths max falls within the range (τ max , DTOA o min ) in the majority of cases, ensuring accurate and effective detection of split pulses.The number of ths max falling within the (τ max , DTOA o min ) range as a percentage of the total number of detection thresholds is consistent with the split pulse detection accuracy, where τ max is the maximum pulse width within the frame, and DTOA o min is the minimum adjacent pulse time interval between complete pulses.

Radar Split Pulse Detection Method
The attempts to merge split pulses with different threshold values (ths) are performed on the PDW frame sequentially.If Tg < ths, the adjacent pulses should be merged, and the pulses are considered split pulses.If Tg ≥ ths, they should not be merged, and the pulses are considered complete pulses.The detection threshold range ths ∈ [ths min , ths max ] is set to adjust the pulse detection threshold incrementally from small to large with a fixed step size.Usually, the starting threshold ths min = 0 is set to indicate that pulse merging is not performed.
As shown in Figure 24, the number of pulses in each timeframe can be set as n b , and the pulse merging processing pointer is positioned to the current pulse i, i ∈ [1, n b − 1] that needs to be processed.If the current pulse, denoted by the pulse merging processing pointer, and the next pulse i + 1 satisfy the merging condition, the pulses will be merged.Otherwise, the original state of the pulse will remain unchanged.Then, the position of the pulse merging processing pointer is updated, and the search for the next pulse that satisfies the pulse merging condition continues until the pointer points to the pulse with sequence number n b − 1, indicating that the merging process is complete.Then, the position of the pulse merging processing pointer is updated to locate the next pulse that satisfies the pulse merging condition until the pointer points to the pulse with sequence number n b − 1.According to the given rules, the PDW sequences generated by different detection thresholds are merged.The PW dimension of the merged PDW sequence is extracted, and the DC dimension is calculated using Equation (19).The histogram of the PW and DC dimensions is plotted, and the time domain convergence degree is calculated to create the time domain convergence degree matrix: where i = 1, • • • , I is the sequence number of the detection threshold.Finally, the convergence degree of adjacent detection thresholds in the time domain is derived to reflect the change rate k i,i+1 of a fixed step size m: The results of the derivation are then formed into the time domain convergence degree derivative matrix: where i = 2, • • • , I, k(1, 2) and k(I − 1, I) represent the derivative results of the time domain convergence degree of the first two and the last two detection thresholds, respectively.According to the above calculation process, the convergence degree derivative matrixes of the PW and DC dimensions are calculated as k τ and k DC .If the time domain convergence degree TCD i remains unchanged, the change rate k(i, i + 1) is equal to 0 and the convergence degree derivative matrices k τ and k DC are all zero matrices with the increase in the detection threshold, indicating that the pulse timing position does not change, and no new merging pulse is generated.If the time domain convergence degree changes, the convergence derivative matrices k τ and k DC are not all-zero matrices, indicating that the new merging pulses are generated.The discriminative labels are generated according to the above algorithm shown in Table 4.

The Convergence Derivative Matrices k τ and k DC Label
All zero matrices 0 Not all zero matrices 1 No pulse is received −1

Parameter Setting
The simulation data in Table 5 and experimental data in Table 6 are derived from multiple radars with fixed PRI, staggered PRI, and group-staggered PRI modulation types.These data are numerically translated based on the original data.Then, the methods for setting data under various conditions are described.Pulse flow density: Set the PRI of radars A, B and C to generate different PRI settings, resulting in pulse flow densities of magnitudes in the ranges of [10 2 , 10 3 ), [10 3 , 10 4 ) and [10 4 , 10 5 ) pulses per second.Only a fixed PRI signal is transmitted.The pulse-splitting degree q is set to 1/3 for each radar, while the pulse-staggering degree RN is kept constant at 1.
Pulse-splitting degree: For radar D, set q = 1/4, 1/3 and 1/2 to generate different pulse splitting degrees.The pulse flow density ρ remains in the same order of magnitude, with [10 3 , 10 4 ) pulse numbers per picosecond.Only a fixed PRI signal is transmitted, and the pulse staggering degree RN is set to 1 remain unchanged.
Pulse-staggering degree: Combine radars E, F and G for interleaving to generate different pulse-staggering degrees within a single frame with RN = 1, 2, 3, respectively.For each radar, the signal with the same pulse-splitting degree q = 1/3 and fixed PRI modulation type is set.The data frame is taken with a pulse flow density in the same order of magnitude [10 3 , 10 4 ) pulse numbers per picosecond.
PRI modulation cmplexity: Set different PRI modulation complexities for radars H, I and J: fixed PRI, staggered PRI and grouped staggered PRI.Each radar is set with the same pulse flow density order of magnitude [10 3 , 10 4 ) pulse numbers per picosecond.Additionally, a pulse-splitting degree of q = 1/3 and a pulse-staggering degree of RN = 1 are set for each radar.
Integrated pulse flow characteristics: Set different radar parameters for radars K, L and M with varying pulse flow density, pulse-splitting degree, pulse-staggering degree and PRI modulation complexity to reasonably replicate experimental data in a complex electromagnetic environment.
Experimental data: The labeled experimental data are used to verify the validity of the proposed model and algorithm.

Verification of the Statistical Distribution of the Split Pulse Number
The number of split pulses corresponding to the original pulse is counted using the experimental data from the same radar.The percentage of the original pulse number f i corresponding to the different split pulse number i in the total number E of the original pulse is then calculated.The statistical results are shown in Table 7.The histogram is drawn based on the statistical distribution of the number of split pulses and is compared to the geometric distribution law, as shown in Figure 25.In order to validate the rationality of the statistical model for split pulse number, the K-S (Kolmogorov-Smirnov) test [19] is employed to determine if the split pulse number of the experimental pulse flow follows a geometric distribution.Let the null hypothesis H 0 to be tested be: the distribution of the split pulse number Y of the experimental pulse flow belongs to the geometric distribution family, i.e., Y ∼ Geo(p).It is verified that by accepting H 0 at the significance level α = 0.05, the distribution of population Y is considered to belong to the geometric distribution family, i.e., Y ∼ Geo(p).

Performance Analysis of Algorithms
This section analyzes whether the variation of window length and frame length affects the split pulse detection result and verifies the validity of the empirical formula used for determining the split pulse detection threshold.

Window Length
As one of the parameters that determines the level of convergence in the time domain, the window length can be adjusted in the split pulse detection algorithm.Therefore, it is necessary to test whether the degree of convergence in the time domain, as calculated using different window lengths, affects the result of split pulse detection.Keep other parameters unchanged, set the window length as w = 1, 2, 3, 4, and conduct comparative experiments using simulated data.When w = 1, as shown in Figures 26 and 27, regardless of whether split pulses exist in the frame, the time domain convergence degree of PW and DC is 0 and does not change with the detection threshold.When w > 1, as shown in Figure 26, the convergence degree of PW and DC in the time domain will change with the increase in the detection threshold.The derivative matrices of convergence are not all-zero matrices, and split pulses exist within the frame.As shown in Figure 27, the convergence degree of PW and DC in the time domain remains unchanged with the increase in the detection threshold.The derivative matrices of convergence are all-zero matrices, indicating that no split pulse exists in the frame.Therefore, when the window length is larger than 1, consistent split pulse detection results will be obtained regardless of any changes in the window length.

Frame Length
The frame length determines the length of the data sequence processed by the split pulse detection algorithm in a single iteration.Below, the test will be conducted to determine whether the frame length affects the result of the split pulse detection.The range of frame length variation for different pulse flow characteristics is set as [0.1 : 0.1 : 0.9, 1 : 1 : 10] ms.The results of comparative experiments using simulation data are shown in Figures 28-32.The accuracy of split pulse detection is calculated as the percentage of correctly detected timeframes out of the total number of timeframes.Because the split pulse is randomly generated, it is necessary to set up 100 Monte Carlo simulation experiments, and the average of the detection results generated many times under the same experimental condition is used to ensure the stability of the detection results.Since the split pulse is generated randomly, it is necessary to set 100 Monte Carlo simulation experiments and calculate the mean value of the detection results generated repeatedly under the same experimental conditions to ensure the stability of the detection results.
Figure 28 shows the detection effect with only pulse flow density changes, and 'All' means that pulse flow density is not distinguished.When the pulse flow density and frame length are both relatively small, it becomes difficult to detect whether the pulse is split.When the pulse flow density is [10 2 , 10 3 ) pulse numbers per picosecond, and the frame length T f ≤ 1 ms, the detection accuracy is 0. When the pulse flow density is [10 3 , 10 4 ) pulse numbers per picosecond and the frame length T f = 0.1 ms, the detection accuracy is 0. It can be seen that the smaller the pulse flow density, the more sensitive the frame length selection.The reason for the above phenomenon is that once the pulse flow density becomes smaller, the pulse flow will become sparse when it arrives.When the frame length is relatively small, it is easy to have an extreme situation in which only one split pulse exists in the frame, so that a single split pulse is misjudged as a complete pulse.In addition, regardless of the pulse flow density in the range, with the increase in frame length, the detection accuracy will gradually increase and exceed 99.9% and then remain unchanged.Comparing the change trend of detection effect of different pulse flow density, it can be seen that the higher the pulse flow density, the higher the detection accuracy.When the frame length is large, the algorithm has an accurate detection effect on different pulse flow densities.
Figure 29 shows the detection effect of only pulse splitting degree changes, and 'All' indicates that the pulse-splitting degree is not distinguished.If the frame length is relatively small, it is impossible to detect whether the pulse is split.When T f = 0.1 ms, the detection accuracy is 0. This is because when the frame length is relatively small, regardless of the degree of pulse splitting, the single split pulse remaining in the frame may be misjudged as a complete pulse, leading to detection errors.It can be seen that the selection of frame length has similar sensitivity to different degrees of pulse splitting.In addition, as the frame length increases, the detection accuracy gradually increases to 100% and remains constant.Comparing the trend of detection accuracy for different pulse-splitting degrees, it can be observed that the detection effect becomes more significant as the pulse-splitting degree decreases.When the frame length is large, the algorithm has a reliable detection effect for various degrees of pulse splitting.Figure 30 shows the detection effect with only pulse staggering degree changes, and 'All' indicates that the pulse-staggering degree is not distinguished.If the frame length is small, the split pulse cannot be detected.When RN = 1 and T f = 0.1 ms, RN = 2 and T f ≤ 0.2 ms, and RN = 3 and T f ≤ 4 ms, the detection accuracy is 0. It can be observed that as the pulse staggering degree increases, the requirement for frame length selection becomes more stringent.The above phenomenon can be attributed to two reasons.One reason is that when the frame length is reduced to the point where it no longer meets the agreed degree of pulse interleaving within a single timeframe, the calculation detection accuracy rate becomes meaningless.Therefore, the detection accuracy rate is set to 0. The second reason is that a frame length that is too small will cause the algorithm to mistake the only split pulse in a single timeframe for a complete pulse, resulting in detection errors.In addition, regardless of the degree of pulse staggering, the detection accuracy rate gradually increases to 100% as the frame length increases.It can be seen that when the frame length selection meets the detection requirements for pulse staggering degree within a single timeframe, the algorithm demonstrates a remarkable detection effect across various pulse staggering degrees.
Figure 31 shows the detection effect of only PRI modulation complexity changes, and 'All' indicates that PRI modulation complexity is not distinguished.When T f = 0.1 ms, the detection accuracy of all PRI modulation types is 0, and the reason for this phenomenon is consistent with Figure 29.It can be seen that selecting a frame length that is too small for different PRI modulation complexities will result in detection errors.In addition, as the frame length increases, the detection accuracy gradually increases and exceeds 99.8%.By comparing the change trend of detection accuracy for different PRI modulation complexities, it can be observed that the detection accuracy is higher for simpler PRI modulation types.Therefore, when the frame length is large, the algorithm has a reliable detection effect on different PRI modulation types.Figure 32 shows the detection effect of the comprehensive change in pulse flow characteristics.According to the trend of detection effect, as the frame length decreases, there is a higher likelihood of misjudging a single split pulse within the frame as a complete pulse, leading to a worse detection effect.Conversely, a longer frame length results in a better detection effect.
Figure 33 shows the changing trend of the detection effect for the selected six groups of experimental datasets.It can be observed that the change pattern of the experimental dataset is similar to that of the simulation dataset.The detection accuracy rate gradually improves and consistently exceeds 96% as the frame length increases, which confirms that the algorithm can reliably produce accurate results on the experimental dataset.To summarize, as long as the frame length is reasonable, the split pulse detection effect remains reliable regardless of any changes in the pulse flow characteristics.In engineering applications, the frame length is generally not less than 5 ms.

Split Pulse Detection Threshold
Table 8 presents the simulation data results for the mean p mean , maximum p max , minimum p min and the difference between the maximum and minimum values ∆p of the split pulse detection accuracy for different pulse flow characteristics.The results are obtained from 100 Monte Carlo simulation experiments.It can be seen that 99.55% of the minimum detection accuracy of p mean comes from the fixed PRI modulation type, while 99.96% of the minimum detection accuracy of p max , 95.55% of the minimum detection accuracy of p min , and 4.41% of the peak value of ∆p all derive from the comprehensive pulse flow characteristics.The results show that when the split pulse detection algorithm is applied to the changing complex electromagnetic environment, compared with the changing single pulse flow characteristics, the detection accuracy and the stability of the algorithm are slightly worse.
Then, the detection algorithm is applied to the experimental data, as shown in Table 9.The detection accuracy p is maintained above 95.68%, which verifies the accuracy and stability of the empirical formula for the split pulse detection threshold.

Conclusions
This paper solves the problem of split pulse detection in a radar pulse stream, and it proposes a method of split pulse modeling and detection.Firstly, the generation mechanism of split pulse is analyzed, and the theory of split site is proposed to describe different splitting modes.Then, the split pulse signal model and the statistical model of split pulse number based on geometric distribution are constructed, respectively, to construct the simulation data of a split pulse under different characteristics.Furthermore, in order to avoid a large number of threshold searches during pulse splitting detection, an empirical formula of the pulse-splitting detection threshold based on the convergence degree in a time domain is given to quickly determine whether there is a split pulse in the radar pulse stream.

Figure 1 .
Figure 1.The schematic of pulse split site.

Figure 3 .
Figure 3.The simulation diagram of cross-channel signal split.

Figure 4 .
Figure 4.The overlapping split sites of the split pulses.

Figure 5 .
Figure 5.The schematic of rabbit ear effect.

Figure 6 .
Figure 6.The simulation diagram of rabbit ear effect.

Figure 7 .
Figure 7.The independent split sites of the split pulses.

Figure 8 .
Figure 8.The schematic of transition band repeated signal.

Figure 9 .
Figure 9.The simulation diagram of transition band repeated signal.

Figure 10 .
Figure 10.The mutually contained split sites of the split pulses.
) PW The value range of split PW is τ b ∈ [τ min , τ o ], and τ min is the minimum value of split PW, generating N random values τ b that obey a uniform distribution.(3) Arrival Time Front and Back Edge Let the original pulse arrival time front edge be TOA o , the back edge be TOE o , and the split pulse front site range be denoted as TOA b ∈ [TOA o , TOE o − τ min ], where TOE o − τ min represents the maximum arrival time front that the last split pulse can take.The range of the split pulse back edge site TOE b ∈ [TOA o + τ min , TOE o ], where TOA o + τ min represents the minimum arrival time back edge that the first split pulse can take.Then, N random values values for TOA b and TOE b are generated, which follow a uniform distribution.(4) PA The original PA is PA o , the corresponding split PA range is denoted by PA b ∈ [PA min , PA o ], and PA min represents the minimum value that the split PA can reach, generating N random integers PA b that follow a uniform distribution.

Figure 11 .
Figure 11.Geometric distribution law of splitting pulse number.

Figure 12 .
Figure 12.Average number of split pulses with different PWs.

Figure 13 .
Figure 13.Average number of split pulses of different PA.

Figure 17 .
Figure 17.The DC of independent split pulses.

Figure 18 .
Figure 18.The DC of dependent split pulses.

Figure 19 .Figure 20 .
Figure 19.The DC of only one pulse in a timeframe.

Figure 21 .
Figure 21.Pulse stream with a higher TCD (presence of split pulses).

Figure 22 .
Figure 22.Pulse stream with a higher TCD (presence of overlapping pulses).

Figure 25 .
Figure 25.Comparison of experimental data distribution law and geometric distribution law.

Figure 26 .
Figure 26.The existence of split pulses.

Table 1 .
Advantages and disadvantages of existing radar pulse-splitting methods.

Table 2 .
The maximum number of split pulses required to form a complete pulse.

Table 4 .
Discriminative labels generated by the algorithm.

Table 7 .
Statistical table of split pulse numbers corresponding to the original pulses.

Table 8 .
Under different pulse flow characteristics, the accuracy of split pulse detection is used to verify the validity of the empirical formula for determining the split pulse detection threshold.

Table 9 .
Under the experimental dataset, the accuracy of split pulse detection is used to verify the rationality of the empirical formula of split pulse detection threshold.