A Segmented Sliding Window Reference Signal Reconstruction Method Based on Fuzzy C-Means

: Reference signal reconstruction serves as a crucial technique for suppressing multipath interference and noise in the reference channel of passive radar. Aiming at the challenge of detecting Low-Slow-Small (LSS) targets using Digital Terrestrial Multimedia Broadcasting (DTMB) signals, this article proposes a novel segmented sliding window reference signal reconstruction method based on Fuzzy C-Means (FCM). By partitioning the reference signals based on the structure of DTMB signal frames, this approach compensates for frequency o ﬀ set and sample rate deviation individually for each segment. Additionally, FCM clustering is utilized for symbol mapping reconstruction. Both simulation and experimental results show that the proposed method signi ﬁ cantly suppresses constellation diagram divergence and phase rotation, increases the adaptive cancellation gain and signal-to-noise ratio (SNR), and in the meantime reduces the computation cost.


Introduction
In recent years, Low-Slow-Small (LSS) targets represented by Unmanned Aerial Vehicles (UAVs) have made full use of low-altitude flight capabilities [1].They posed a threat to government agencies, airports, and important gatherings [2].The major security risks had attracted worldwide attention.At present, the main detection means of LSS targets include radio detection equipment [3], photoelectric equipment [4], and radar [5].The radio detection equipment detects targets by using the image transmission and telemetry signal of the UAV, but it is vulnerable to the interference of a complex electromagnetic environment [6].The photoelectric equipment tracks targets through infrared and visible light; however, it is seriously affected by the weather [7,8].Radar can work in all weather conditions and is the main method for target early warning and detection [9,10].
Due to the strong environment interference and small radar cross section of LSS targets [11], traditional radar has limited detection capabilities.Passive radar utilizes non-cooperative signals as illumination sources to achieve moving target detection and localization through coherent accumulation [12].It is considered to be a potentially effective solution to the above challenges.Exploiting its inherent characteristic of not emitting radio waves, it eliminates the influence of radio interference and effectively avoids being detected and disrupted by enemy [13,14].Compared with traditional radars, its excellent low-altitude coverage capabilities and easy deployment give it obvious advantages in detecting LSS targets [15,16].However, the emission power and waveform parameters of illumination sources are uncontrollable [17].Therefore, appropriate signals must be selected as radiation sources to ensure the detection performance of the system.Currently, research concerning the illumination sources of opportunity signals includes Digital Audio Broadcasting (DAB) [18][19][20], Digital Video Broadcasting-terrestrial (DVB-T) [21][22][23], Digital Terrestrial Multimedia Broadcasting signals (DTMB) [24][25][26], etc.The DTMB, characterized by wide signal bandwidth, high emission power, and extensive coverage, emerges as an ideal illumination source for passive radar [27].
A novel approach was utilized in target detection by simultaneously receiving target echo signals from the surveillance channel and reference signals from the reference channel in passive radar [28,29].After the correlation processing of both channels, crucial target information such as distance, speed, and coordinates were extracted [30].The reference signal serves two important purposes.Firstly, it contributes to clutter cancellation to suppress direct waves and clutter in the echo signals [31].Secondly, it facilitates a coherent accumulation with the echo signals for target detection [32].The effectiveness of interference suppression and the reliability of target detection were significantly influenced by the quality of the reference signal [33,34].In practical scenarios, the reference signal included not only direct waves but also multipath interference and clutter caused by buildings and terrain during the transmission, which posed new challenges for signal reconstruction methods [35,36].
The time-domain reference signal reconstruction method was drawn on the signal demodulation techniques commonly used in communication systems [37].Interference suppression in the reference channel was achieved, and non-coherent interference could be effectively suppressed [38].However, considerable computational resources were consumed in the process of symbol reconstruction by adopting this method [39].Additionally, errors may occur during the estimation of frequency offset and sample rate deviation, resulting in constellation diagram divergence and phase rotation [40].Consequently, the bit error rate (BER) of the reconstructed signal was increased, thus reducing the interference suppression gain.
To address the issues of constellation diagram divergence and phase rotation, and reduce the computation cost of the reconstruction process, a segmented sliding window reference signal reconstruction method based on Fuzzy C-Means (FCM) is proposed in this paper.The proposed method compensates the frequency offset and sampling rate deviation, to solve the problems of constellation divergence and phase rotation caused by error accumulation.The FCM clustering is used in the symbol mapping reconstruction to reduce the computation cost.
The remainder of the paper is organized as follows.The DTMB standard is introduced in Section 2. The traditional reconstruction procedure is overviewed in Section 3, followed by a detailed description of the methods proposed in this article.In Section 4, the effectiveness of the proposed methods is verified by simulated and measured data.A relevant summary and conclusions are given in Section 5.

DTMB Standard 2.1. Signal Structure of DTMB
The DTMB standard, titled 'Framing Structure, Channel Coding and Modulation for Digital Television Terrestrial Broadcasting System', was announced in August 2006 [41].This standard has independent intellectual property rights in China and is a mandatory national standard for digital televisions.The combination of Time-Domain Synchronous Orthogonal Frequency Division Multiplexing (TDS-OFDM) and LDPC (Low-density Paritycheck) coding techniques is used in DTMB signals.In the core of the DTMB system structure, the signal frame is the basic unit.The frame structure of DTMB signals is shown in Figure 1.The frame structure contains two essential components, namely the frame head (FH) and the frame body (FB).The FB is comprised of 36 system information symbols and 3744 data symbols, totaling 3780 symbols.Meanwhile, the FH includes three variations: 420 or 945 symbols in the multi-carrier mode and 595 symbols in the single-carrier mode.
variations: 420 or 945 symbols in the multi-carrier mode and 595 sym carrier mode.

Frame Head
Frame Body Preamble m-Sequence Postamble Taking PN420 mode as an example, the FH is composed of a p order PN sequence, and a postamble.Notably, the preamble and posta of the initial 83 symbols and the concluding 82 symbols of the m-seq A distinct PN sequence with different initial phases is used for each FH identification of the frame number.The linear feedback shift register in Figure 2.For a multi-carrier with frame head PN420, the time-domain expr frame is shown in Equation ( 1), where t N and c N represent the sy and FB, respectively; ( ) PN n is frame head data; k denotes the nu frequency; and ( ) X k represents the FB symbol before modulation.

Transmitting System Scheme
The DTMB transmission system completes the conversion from a to the transmission signal [42,43].Figure 3 shows the transmitter blo consists of the following elements: Taking PN420 mode as an example, the FH is composed of a preamble, an eighth-order PN sequence, and a postamble.Notably, the preamble and postamble are composed of the initial 83 symbols and the concluding 82 symbols of the m-sequence, respectively.A distinct PN sequence with different initial phases is used for each FH signal as a unique identification of the frame number.The linear feedback shift register of PN420 is shown in Figure 2.
Remote Sens. 2024, 16, x FOR PEER REVIEW 3 of 22 variations: 420 or 945 symbols in the multi-carrier mode and 595 symbols in the singlecarrier mode.

Frame Head Frame Body
Preamble m-Sequence Postamble Taking PN420 mode as an example, the FH is composed of a preamble, an eighthorder PN sequence, and a postamble.Notably, the preamble and postamble are composed of the initial 83 symbols and the concluding 82 symbols of the m-sequence, respectively.A distinct PN sequence with different initial phases is used for each FH signal as a unique identification of the frame number.The linear feedback shift register of PN420 is shown in Figure 2.For a multi-carrier with frame head PN420, the time-domain expression of the signal frame is shown in Equation (1), where t N and c N represent the symbol length of FH and FB, respectively; ( ) PN n is frame head data; k denotes the number of the carrier frequency; and ( ) X k represents the FB symbol before modulation.

Transmitting System Scheme
The DTMB transmission system completes the conversion from an input data stream to the transmission signal [42,43].Figure 3 shows the transmitter block diagram, which consists of the following elements:  For a multi-carrier with frame head PN420, the time-domain expression of the signal frame is shown in Equation (1), where N t and N c represent the symbol length of FH and FB, respectively; PN(n) is frame head data; k denotes the number of the carrier frequency; and X(k) represents the FB symbol before modulation.

Transmitting System Scheme
The DTMB transmission system completes the conversion from an input data stream to the transmission signal [42,43].Figure 3 shows the transmitter block diagram, which consists of the following elements:

The Process of Reference Signal Reconstruction
The structure of DTMB signals has the ability to resist multipath interference and noise.Under specific SNR conditions, interference in the reference signal can be effectively suppressed through signal reconstruction.The signal reconstruction process involves various stages including synchronization, channel estimation and equalization, constellation demapping, and decoding.Subsequently, the signal is remapped and shaped according to the standard protocol to obtain the original baseband transmission signal, thus obtaining a purified reference signal.A comprehensive description of the reconstruction scheme at the reception is shown in Figure 4.

Synchronization
In DTMB systems, the accurate synchronization processing of the received data at the receiver is crucial.Without precise synchronization, divergence and phase rotation may occur in the constellation diagram, which will lead to an increased BER of signal reconstruction.The DTMB signal ingeniously utilizes PN sequences as FH, which have exceptional autocorrelation characteristics to facilitate rapid and accurate synchronization.This method has significant advantages in operation simplicity and excellent performance.In the field of DTMB synchronization technology, synchronization is typically divided into frame, carrier frequency, and sample rate synchronization.
The digital television signal transmitter operates on the basis of signal frames as its fundamental unit.To ensure accurate signal recovery at the receiver, frame synchronization must be realized, which requires the identification of the starting position of the signal frame.Considering the fixed PN sequence in the FH of DTMB signals, the identification of the digital television signal FH type relies on the well-correlated characteristics of the PN sequence [44].By conducting cross-correlation between the known PN sequence and the signal, the peak position is served as an estimate of the frame header position.The calculation process is given by

The Process of Reference Signal Reconstruction
The structure of DTMB signals has the ability to resist multipath interference and noise.Under specific SNR conditions, interference in the reference signal can be effectively suppressed through signal reconstruction.The signal reconstruction process involves various stages including synchronization, channel estimation and equalization, constellation demapping, and decoding.Subsequently, the signal is remapped and shaped according to the standard protocol to obtain the original baseband transmission signal, thus obtaining a purified reference signal.A comprehensive description of the reconstruction scheme at the reception is shown in Figure 4.

The Process of Reference Signal Reconstruction
The structure of DTMB signals has the ability to resist multipath interference and noise.Under specific SNR conditions, interference in the reference signal can be effectively suppressed through signal reconstruction.The signal reconstruction process involves various stages including synchronization, channel estimation and equalization, constellation demapping, and decoding.Subsequently, the signal is remapped and shaped according to the standard protocol to obtain the original baseband transmission signal, thus obtaining a purified reference signal.A comprehensive description of the reconstruction scheme at the reception is shown in Figure 4.

Synchronization
In DTMB systems, the accurate synchronization processing of the received data at the receiver is crucial.Without precise synchronization, divergence and phase rotation may occur in the constellation diagram, which will lead to an increased BER of signal reconstruction.The DTMB signal ingeniously utilizes PN sequences as FH, which have exceptional autocorrelation characteristics to facilitate rapid and accurate synchronization.This method has significant advantages in operation simplicity and excellent performance.In the field of DTMB synchronization technology, synchronization is typically divided into frame, carrier frequency, and sample rate synchronization.
The digital television signal transmitter operates on the basis of signal frames as its fundamental unit.To ensure accurate signal recovery at the receiver, frame synchronization must be realized, which requires the identification of the starting position of the signal frame.Considering the fixed PN sequence in the FH of DTMB signals, the identification of the digital television signal FH type relies on the well-correlated characteristics of the PN sequence [44].By conducting cross-correlation between the known PN sequence and the signal, the peak position is served as an estimate of the frame header position.The calculation process is given by

Synchronization
In DTMB systems, the accurate synchronization processing of the received data at the receiver is crucial.Without precise synchronization, divergence and phase rotation may occur in the constellation diagram, which will lead to an increased BER of signal reconstruction.The DTMB signal ingeniously utilizes PN sequences as FH, which have exceptional autocorrelation characteristics to facilitate rapid and accurate synchronization.This method has significant advantages in operation simplicity and excellent performance.In the field of DTMB synchronization technology, synchronization is typically divided into frame, carrier frequency, and sample rate synchronization.
The digital television signal transmitter operates on the basis of signal frames as its fundamental unit.To ensure accurate signal recovery at the receiver, frame synchronization must be realized, which requires the identification of the starting position of the signal frame.Considering the fixed PN sequence in the FH of DTMB signals, the identification of the digital television signal FH type relies on the well-correlated characteristics of the PN sequence [44].By conducting cross-correlation between the known PN sequence and the signal, the peak position is served as an estimate of the frame header position.The calculation process is given by where pn is the known FH, and s * is the complex conjugation of the received signal.Carrier frequency synchronization is essential for ensuring accurate frequency alignment between the receiver and transmitter, as well as among the individual subcarrier frequencies [45].For DTMB signals, only considering the OFDM symbol of the FB, the output of the transmitter after inverse Fourier transform can be expressed as where N c represents number of subcarriers, and S k represents the Fourier transform of s(n) with respect to n.After the signal s(n) passes through the discrete multipath fading channel h(n) with noise w(n), in the case of ideal synchronization, the received signal r(n) is given by where H k is the Fourier transform of h(n) with respect to n.The carrier frequency offset (CFO) originates from the differences in oscillator frequency between the transmitter and receiver, which may arise due to inherent instability in the oscillators or fluctuations of the transmission channel affected by environmental factors [46].Assuming that the CFO is being the symbol rate of OFDM.ε I and ε F are the integer and decimal part of ε, respectively.The received signal r ε (n) can be expressed as The Fourier transform of r ε (n) with respect to n is where W k is the Fourier transform of w(n) with respect to n, and I k is the inter-carrier interference caused by CFO.
When ε F = 0, Equation ( 6) can be simplified as Although the integer partial frequency offset will not destroy the orthogonality between subcarriers, it will cause the sampled signal to shift in the frequency domain, resulting in cyclic shift after OFDM demodulation and demodulation errors.
When ε I = 0, Equations ( 6) and ( 7) can be simplified as where the term sin(πε F )/(N c sin(πε F /N c )) caused by ε F is an amplitude attenuate function on each subcarrier, and the phase is rotated by πε F (N c − 1)/N c on each subcarrier.Therefore, CFO will increase the BER of the reconstructed signal.According to the structure of digital television signal frames, CFO can be estimated by using the phase difference of the PN sequence in adjacent frames on the basis of frame synchronization.The Fast Fourier Transform (FFT) algorithm also serves as a valuable tool for obtaining a rough estimate of the CFO.The estimation accuracy of this method is improved by increasing the number of FFT points N.After the FFT, the synchronization header signal exhibits a characteristic spectrum characterized by obvious peaks.A reliable estimate of the CFO is provided by the displacement of the peak relative to the zero frequency.Due to the inevitable defects in synchronization and clock stability between transmitter and receiver, there exists a difference between the actual sampling rate and the nominal sampling rate at the receiver [47], known as sampling frequency offset (SFO).The transmitted signal can be approximated as where T s is the sample interval, f s = 1/T s is the nominal sampling rate at the transmitter, and the sampling rate at the receiver is Affected by SFO only, the g-th sampling value in the d-th frame of the received DTMB signal can be expressed as where N d is the number of symbols in a frame and N h is the number of symbols in an FH.
The OFDM modulation for the sampled data of the d-th frame is as follows According to Equation ( 13), the SFO will cause the shift in the FFT period, destroy the orthogonality between subcarriers, and lead to inter-carrier interference (ICI).Cumulative clock offset may result in time drift and potential Inter-Symbol Interference (ISI), accompanied by the phase rotation of each subcarrier.These effects are shown as divergence in the constellation diagram of the signal and an increase in the BER of the reconstructed signal.It may affect the performance of interference cancellation, thereby resulting in Doppler sidelobe energy leakage.SFO influences both the time and frequency domains.Consequently, SFO can be estimated in the time domain by using the autocorrelation characteristics of δ-like functions of the DTMB signal's PN sequence.

Channel Estimation and Equalization
The reference signal of digital television is affected by both noise interference in the transmission channel and multipath interference caused by reflections from ground, buildings, and other objects.In order to accurately reconstruct the transmitted signal at the receiver, it is imperative to not only achieve synchronization but also estimate the transmission channel and equalize the reference signal [48].
In the digital television signal transmission system, the original transmitted signal matrix S can be expressed as The reference signal received can be expressed as where the channel response vector The estimated channel can be expressed as The square of the difference between the estimated channel and actual channel is minimized to make the estimated channel close to reality.The objective function J( ∧ d) can be expressed as where the superscript H represents the conjugate transpose.In order to minimize the sum of squared error, taking the first-order partial derivative of the above objective function with respect to The best estimate of the transmission channel is After estimating the transmission channel, equalization methods are typically used to suppress multipath interference.Figure 5 shows the principle diagram of equalization, where S represents the transmitted signal matrix, r represents the received signal passing through the channel with response d and polluted by additive noise, and S ′ represents the equalized signal matrix.After estimating the transmission channel, equalization methods are typically used to suppress multipath interference.Figure 5  The system function ' ( ) d z of the equalization filter and the channel transfer func- tion ( ) Based on the estimated channel response, the system function of the channel equalization filter can be obtained.
Based on the estimated channel response, the system function of the channel equalization filter can be obtained.

Decoding
After the equalization and deinterleaving, there still exist a lot of bit errors.DTMB's coding and decoding mechanism can decrease the errors more.The DTMB standard suggests the use of Bose-Chaudhuri-Hocquenghem (BCH) combined with LDPC coding as the Forward Error Correction (FEC) design.The decoding of LDPC codes uses the minimum sum algorithm of belief propagation (BP).The core idea of the algorithm is to use the received soft information to iterate between the variable node and the check node.
Then, the decoded data stream is mapped, interleaved, and framed with the FH.The signal is converted to the actual sampling rate to obtain the reconstructed signal.For DTMB signals, the received reference signal frame time is as long as seconds [49].The increased divergence and phase rotation are caused by the accumulation of CFO and SFO estimation errors.In the meantime, the computational cost is increased by error correction decoding.This poses significant challenges for engineering applications and real-time processing.Therefore, there is an urgent need to optimize the reconstruction process, aiming to minimize computational cost and facilitate application in engineering scenarios.

Segmented Sliding Window Reference Signal Reconstruction Method Based on Fuzzy C-Means
In order to satisfy the system's real-time processing demands and minimize computational resource consumption, this article proposes a segmented sliding window reconstruction method based on FCM.
The traditional method and the block structure designed for DTMB signals in this article are shown in Figure 6.The received reference signal is taken as one data block for frequency offset and sampling rate deviation estimation in the traditional method.In contrast, the approach proposed effectively works by dividing the signal into different data blocks based on the number of frames, which reduces the accuracy requirement of deviation estimation.In each block, frequency offset and sampling rate deviation are estimated independently, and the data block is compensated based on the estimated values.Notably, there is an overlap of signal frames between different data blocks, which ensures the continuity of the reconstructed signal.Furthermore, FCM clustering is used to map the data to cluster centers instead of LDPC decoding.According to the data membership to the cluster centers, the data are remapped to the center values to complete the reconstruction process.The process of the segmented sliding window reconstruction method base is shown in Figure 8.The process of the segmented sliding window reconstruction method based on FCM is shown in Figure 8.  Figure 7 shows the 4QAM signal frame constellation.The presence of distinct clustering centers in the constellation diagram is readily apparent.The data can be classified into respective cluster centers by utilizing clustering methods.FCM clustering determines the probability of data belonging to a specific class by the size of membership [50,51].This method is suitable for all types of PN headers in DTMB, thus improving the effectiveness of the reconstruction process.The process of the segmented sliding window reconstruction method based on FCM is shown in Figure 8.  Step1: Based on the structural properties of digital television signal frames, autocorrelation analysis is conducted in the reference signal.By combining the relationship between the baseband data sampling rate and 7.56 MHz, the frame length of three distinct PN headers is converted to the frame length under the baseband data sampling rate, which is expressed as where f s is actual baseband data sampling rate, and L PN is the frame length of different types of digital television signal frames at 7.56 MHz.Comparing the peak positions of the signal autocorrelation results and frame length, the frame header pattern and data type of the reference signal can be determined.In the case of PN595 and PN420 with a zero phase, along with PN945, the peak position identified in the cross-correlation is denoted as the precise location of the FH signal.However, for PN headers of the cyclic shift variety, the additional synchronization of the frame header sequence number becomes essential after the determination of the estimated frame header position.
In the example of PN420 in cyclic shift mode, the 225 PN420 sequences exhibit a symmetrical distribution centered around the 112th sequence.Notably, the maximum difference occurs within the interval spanning the 111th and 112th sequences, whereas the minimum difference is observed between the 112th and 113th sequences.The difference between the frame header positions is utilized to estimate frame header position, which is derived from the primary PN sequence correlation.By identifying the positions with the maximum and minimum differences, the frame header sequence position can be obtained.
According to the FH position and data type corresponding to the reference signal, each segment of the reference signal is denoted as x(n), and each segment of the echo signal is denoted as y(n), with N k segments in total.The frame number of each signal segment N m is shown in Equation (22), and the step number of each signal segment N s is shown in Equation (23).
where N PN is the number of signal frames of different types of signals in the acquisition time at the 7.56 MHz.
Step2: The process of N-point FFT on the reference signal x(n) is given by where X(k) represents the FFT results of the reference signal.According to the position corresponding to the peak of spectral line, the frequency offset value can be expressed as The CFO is given by To compensate frequency offset, the reference signal x c (n) and the echo signal y c (n) are given by where n = [1, 2, 3, • • • , i], and i denotes the number of data points in this segment.
Based on the symbol rate B s of the transmitting station, f s represents the nominal value of the sampling frequency, f ′ s represents the actual value, and the theoretical position of the k-th correlation peak corresponding to the k-th frame can be expressed as where N symbol indicates the number of symbols in one signal frame, and the actual position of the correlation peak corresponding to the k-th frame can be expressed as The SFO is given by Step3: The transverse filter, commonly employed for equalization, is widely utilized in practical applications.However, due to the limitations in the filter order, it cannot completely eliminate multipath interference, but it can maximize interference suppression.Figure 9 shows the structure of the 2N + 1 order transversal filter.Step3: The transverse filter, commonly employed for equalization, is widely u in practical applications.However, due to the limitations in the filter order, it canno pletely eliminate multipath interference, but it can maximize interference suppr Figure 9 shows the structure of the 2N + 1 order transversal filter.The Minimum Mean Square Error (MMSE) equalization criterion is used in th cle.MMSE equalization adjusts the coefficients of the transversal filter to minim mean square error of the output signal.Since the digital television signal frame ty be determined at the receiver end, the equalization filter coefficients can be solved the standard PN sequence and the received PN sequence.

The coefficient of 2N + 1 order filter is
), , ( ), ( ) ), , ( ), , ( ) , and the signal of the equalization filter can be expressed as n n q n w p The mean square error is The Minimum Mean Square Error (MMSE) equalization criterion is used in this article.MMSE equalization adjusts the coefficients of the transversal filter to minimize the mean square error of the output signal.Since the digital television signal frame type can be determined at the receiver end, the equalization filter coefficients can be solved using the standard PN sequence and the received PN sequence.
The coefficient of 2N + 1 order filter is and the output signal of the equalization filter can be expressed as The mean square error is where q(n) represents the n-th symbol of a standard PN sequence, q * (n) represents the conjugation of q(n), and re{•} represents the real part of the complex.For stationary signals, the autocorrelation matrix is non-negative, and the mean square error is a quadratic function and only has one minimum value.Then, the MMSE equalization filter coefficient is given by where R xx = E p n p H n represents the autocorrelation matrix of the PN sequence, and r xs = E[q * (n)p n ] represents the cross-correlation vector of PN sequence.Each reference signal is filtered by the equalization filter coefficient, and the equalization process in the frequency domain is expressed as where x e (k) is the frequency domain expression of x e (n), x(k) is the frequency domain expression of x(n), and w k is the frequency domain expression of w opt .⊙ is Hadamard product, which represents corresponding position elements' multiplication of the vector.

Step4:
The x e (n) is resampled to the common multiple of the baseband sampling rate and 7.56 MHz, and it is filtered by using a square root raised cosine (SSRC) filter with a roll-off coefficient of 0.05.The filter frequency response is as follows: where f N is the Nyquist frequency.The data from each segment are sampled down to the baud rate of 7.56 MHz.Meanwhile, SFO is compensated in the down-sampling process.
Then, the FCM method is utilized for data clustering, which obtains cluster centers by minimizing an objective function.The sum of Euclidean distances from each point to each cluster center is essentially represented with this objective function.The clustering process is involved in minimizing this objective function through iterative computations, gradually reducing the error value.Upon convergence of the objective function, the final clustering results are obtained.
The Euclidean distance from each data point to the cluster center can be expressed as where b i represents the i-th data point, c k represents the k-th cluster center, and u ki represents the probability that point i belongs to the k cluster.Set the number of clusters k and threshold ε, initialize membership matrix u, and calculate clustering center of each cluster as The LaGrange multiplier is constructed for optimizing the cluster center, as shown in Equation (39).When L ≥ ε, Equation ( 40) is used to update the membership matrix until the cluster centers are obtained:

Simulation and Experimental Verification 4.1. Simulation Analysis
In order to evaluate the effect of the proposed method in this article on the BER under the conditions characterized by different SNRs, multi-carrier PN420 signals with identical bit rates and four distinct mapping modes are generated through simulation.Four distinct signals are subjected to Gaussian white noise interference at different SNRs.The interfered data are synchronized, equalized, and then converted to the baud rate of 7.56 MHz.
The data are reconstructed utilizing four different methods, including LDPC error correction decoding, K−Means, FCM clustering, and MATLAB's QAM modulation function 'Qamdemod'.We initialized the number of clusters as k = 4, 16, 32, 64, the threshold of Fuzzy C-Means as ε = 0.001, and the membership matrix as u = 0.The results are shown in Figure 10.
slightly higher than the preceding two methods.For the 32QAM and 64QAM modulation methods, the BER of FCM clustering is higher than LDPC decoding in low SNR scenarios, but it is close to LDPC decoding under high SNR conditions.The BER change trend of the K−Means method is similar to that of FCM, but it has higher BER.The Qamdemod method exhibits a significantly higher BER.This observation can be attributed to the fact that simpler constellation mapping has fewer cluster centers and is more easily clustered.In addition, the high SNR constellation is more clustered, which makes it easier to identify the clustering center of each data point.In contrast, in the 32QAM and 64QAM mapping modes, it is difficult to accurately determine cluster centers under low-SNR conditions due to constellation diagram divergence.
Then, for the same data, we compared the running time of FCM clustering and LDPC decoding under the same computer configuration, as shown in Table 1.It can be seen that the operation time consumption of the FCM clustering method is much less than LDPC decoding.The reason is that LDPC decoding involves multiplication with the check matrix, and it needs to be iterated until the result converges to zero.The check matrix properties of It is evident that the BER of the proposed method and LDPC decoding are similar for the 4QAM and 16QAM modulation methods.The BER of K−Means and Qamdemod are slightly higher than the preceding two methods.For the 32QAM and 64QAM modulation methods, the BER of FCM clustering is higher than LDPC decoding in low SNR scenarios, but it is close to LDPC decoding under high SNR conditions.The BER change trend of the K−Means method is similar to that of FCM, but it has higher BER.The Qamdemod method exhibits a significantly higher BER.
This observation can be attributed to the fact that simpler constellation mapping has fewer cluster centers and is more easily clustered.In addition, the high SNR constellation is more clustered, which makes it easier to identify the clustering center of each data point.In contrast, in the 32QAM and 64QAM mapping modes, it is difficult to accurately determine cluster centers under low-SNR conditions due to constellation diagram divergence.
Then, for the same data, we compared the running time of FCM clustering and LDPC decoding under the same computer configuration, as shown in Table 1.It can be seen that the operation time consumption of the FCM clustering method is much less than LDPC decoding.The reason is that LDPC decoding involves multiplication with the check matrix, and it needs to be iterated until the result converges to zero.The check matrix properties of LDPC codes are shown in Table 2. L represents the number of one in the row, M is the line number of the check matrix, and N is the column number of the check matrix.In the process of LDPC decoding using the check matrix, the number of matrix operations involved in one frame is shown in Table 3.The calculation of the addition operation is represented as 2L + N, the multiplication operation is represented as 2L − N, and the comparison operation is represented as L + N − M.These matrix operations may lead to a large computation cost.However, the number of matrix operations with different coding rates in FCM is only affected by mapping modes.It can be represented as k mode •n, where k mode is the number of clustering centers and n is number of iterations.The results show that the FCM clustering method proposed in this paper has great practicability in symbol mapping reconstruction.The segmented sliding window reconstruction method based on FCM proposed in this article involves signal segmentation operation.Different lengths of the segmented sliding block affect the performance of the reconstruction method.To analyze the impact of different block lengths on the performance, one second 16QAM multi-carrier PN420 signal is simulated.The different lengths of the block are selected as 45,90,135,180,225,450,675,900,1125,1350,1575, and 1800 frames.The impact of block length on the performance of the reconstruction method was analyzed by comparing the cancellation gain of the reconstructed signals obtained using different block lengths.The results are shown in Figure 11.
of different block lengths on the performance, one second 16QAM multi-carrier PN420 signal is simulated.The different lengths of the block are selected as 45,90,135,180,225,450,675,900,1125,1350,1575, and 1800 frames.The impact of block length on the performance of the reconstruction method was analyzed by comparing the cancellation gain of the reconstructed signals obtained using different block lengths.The results are shown in Figure 11.When the block length becomes shorter, the Normalized Least Mean Square (NLMS) cancellation gain increases.However, the estimation time of CFO and SFO is increased.When the block length becomes longer, the time required for one data block sampling rate conversion is increased.The comparation of the calculation time between the traditional method and the proposed method with a block length equal to 225 frames is shown in Table 4.When the block length becomes shorter, the Normalized Least Mean Square (NLMS) cancellation gain increases.However, the estimation time of CFO and SFO is increased.When the block length becomes longer, the time required for one data block sampling rate conversion is increased.The comparation of the calculation time between the traditional method and the proposed method with a block length equal to 225 frames is shown in Table 4.During the experimental phase, data were collected from both the reference and the surveillance channels.A one-second segment of received data was selected, known to contain a weak target.Subsequently, two separate datasets were prepared for further analysis: the reference signal using the traditional unpartitioned reconstruction method and using the segmented sliding window reference signal reconstruction method based on FCM proposed in this article.The direct waves and clutter in the surveillance channel were suppressed by clutter cancellation for the two datasets.Cross ambiguity functions were conducted for both datasets, with the result used to obtain Range-Doppler results for the purpose of validating the effectiveness and practicality of the method.
The first set of data used the CCTV Tower in Haidian District, Beijing, as the illuminators.The frame header mode was cyclic-shift PN420, the center frequency was 666 MHz, and the system sampling rate was 10 MHz.
The one-second reference signal was divided into eight segments, and each segment was compensated by frequency offset and sampling rate deviation individually and reconstructed by the FCM clustering method.The 900-th frame constellations are shown in Figure 12.The constellation diagram using the traditional unpartitioned method rotates and diverges.
The first set of data used the CCTV Tower in Haidian District, Beijing, as the illuminators.The frame header mode was cyclic-shift PN420, the center frequency was 666 MHz, and the system sampling rate was 10 MHz.
The one-second reference signal was divided into eight segments, and each segment was compensated by frequency offset and sampling rate deviation individually and reconstructed by the FCM clustering method.The 900-th frame constellations are shown in Figure 12.The constellation diagram using the traditional unpartitioned method rotates and diverges.In response to the reconstructed signals obtained from both methods, NLMS cancellation is adopted to suppress direct waves, multipath interference, and noise in the surveillance channel.Compared to the traditional method, the approach proposed in this article has an increased cancellation gain, as shown by the values presented in Figure 13.In response to the reconstructed signals obtained from both methods, NLMS cancellation is adopted to suppress direct waves, multipath interference, and noise in the surveillance channel.Compared to the traditional method, the approach proposed in this article has an increased cancellation gain, as shown by the values presented in Figure 13.The second set of data is a real digital TV signal of Pinggu TV Tower which is located in Pinggu District, Beijing.The received signal contains a weak target.The frame header mode is PN595, the center frequency is 674 MHz, and the system sampling rate is 10 MHz.
The 864-th frame constellations are shown in Figure 15.The constellation diagram using the traditional unpartitioned method rotates and diverges.NLMS cancellation is adopted to suppress direct and multipath interferences in th surveillance channel.Compared to traditional methods, the approach proposed has a increased cancellation gain, as shown in Figure 16.NLMS cancellation is adopted to suppress direct and multipath interferences in the surveillance channel.Compared to traditional methods, the approach proposed has an increased cancellation gain, as shown in Figure 16.
Figure 17 shows the results of two different reconstruction methods after adaptive cancellation and the cross-ambiguity function.The segmented sliding window reconstruction method based on FCM leads to an enhancement in the detected target amplitude and increases the target SNR from 19.41 dB to 27.89 dB.NLMS cancellation is adopted to suppress direct and multipath interferences in the surveillance channel.Compared to traditional methods, the approach proposed has an increased cancellation gain, as shown in Figure 16.increased cancellation gain, as shown in Figure 16.

Conclusions
In this article, a segmented sliding window reference signal reconstruction method based on FCM was proposed for target detection in passive radar.The frequency offset and sampling rate deviation were estimated for each signal segment, which reduced the accuracy requirement of deviation estimation.This method has effectively suppressed the issues of constellation diagram divergence and phase rotation which arose from error accumulation.The utilization of FCM clustering instead of LDPC error correction decoding for symbol mapping reconstruction contributed to reducing the computation cost.This approach has increased the adaptive cancellation gain and target SNR, which is advantageous for the detection of LSS targets, thus having significant implications for practical engineering applications.

Data Availability Statement:
The data that support the findings of this study are available by request to the corresponding author.The data are not publicly available due to data confidentiality.

Figure 3 .
Figure 3.The digital television transmitter block diagram.

Figure 4 .
Figure 4. Reconstruction scheme at the reception.

Figure 3 .
Figure 3.The digital television transmitter block diagram.

Figure 3 .
Figure 3.The digital television transmitter block diagram.

Figure 4 .
Figure 4. Reconstruction scheme at the reception.

Figure 4 .
Figure 4. Reconstruction scheme at the reception.

Figure 7
Figure 7 shows the 4QAM signal frame constellation.The presence of distinct clustering centers in the constellation diagram is readily apparent.The data can be classified into respective cluster centers by utilizing clustering methods.FCM clustering determines

Figure 6 .
Figure 6.Comparison between this article's blocking and traditional blocking.

Figure 7
Figure7shows the 4QAM signal frame constellation.The presence of distinct clustering centers in the constellation diagram is readily apparent.The data can be classified into respective cluster centers by utilizing clustering methods.FCM clustering determines the probability of data belonging to a specific class by the size of membership[50,51].This method is suitable for all types of PN headers in DTMB, thus improving the effectiveness of the reconstruction process.

Figure 6 .
Figure 6.Comparison between this article's blocking and traditional blocking.

Figure 7
Figure7shows the 4QAM signal frame constellation.The presence of dis tering centers in the constellation diagram is readily apparent.The data can be into respective cluster centers by utilizing clustering methods.FCM clustering d the probability of data belonging to a specific class by the size of membership [5 method is suitable for all types of PN headers in DTMB, thus improving the eff of the reconstruction process.

Figure 8 .
Figure 8.The process of the segmented sliding window reconstruction method based on

Figure 6 .
Figure 6.Comparison between this article's blocking and traditional blocking.

Figure 8 .
Figure 8.The process of the segmented sliding window reconstruction method based on FCM.Figure 8.The process of the segmented sliding window reconstruction method based on FCM.

Figure 8 .
Figure 8.The process of the segmented sliding window reconstruction method based on FCM.Figure 8.The process of the segmented sliding window reconstruction method based on FCM. symbol

Figure 11 .
Figure 11.(a) The cancellation gain versus block length; (b) the algorithm running time changes with the block length.

Figure 11 .
Figure 11.(a) The cancellation gain versus block length; (b) the algorithm running time changes with the block length.

Figure 14
Figure 14 shows the results of two different reconstruction met cancellation and cross-ambiguity function.The segmented sliding win method based on FCM leads to an enhancement in the detected targe creases the target SNR from 30.32 dB to 35.76 dB.

Figure 14
Figure 14 shows the results of two different reconstruction methods after adaptive cancellation and cross-ambiguity function.The segmented sliding window reconstruction method based on FCM leads to an enhancement in the detected target amplitude and increases the target SNR from 30.32 dB to 35.76 dB.

Figure 14 Figure 14 .
Figure14shows the results of two different reconstruction methods after adaptive cancellation and cross-ambiguity function.The segmented sliding window reconstruction method based on FCM leads to an enhancement in the detected target amplitude and increases the target SNR from 30.32 dB to 35.76 dB.

Figure 14 .Figure 15 .
Figure 14.Comparison of LSS target detection effects using PN420 DTMB signal.(a) Traditional unpartitioned reconstruction method; (b) segmented sliding window reconstruction method based on FCM.The second set of data is a real digital TV signal of Pinggu TV Tower which is located in Pinggu District, Beijing.The received signal contains a weak target.The frame header mode is PN595, the center frequency is 674 MHz, and the system sampling rate is 10 MHz.The 864-th frame constellations are shown in Figure15.The constellation diagram using the traditional unpartitioned method rotates and diverges.Remote Sens. 2024, 16, x FOR PEER REVIEW 19 of 2

Figure 17 Figure 16 .
Figure 17 shows the results of two different reconstruction methods after adaptive cancellation and the cross-ambiguity function.The segmented sliding window reconstruction method based on FCM leads to an enhancement in the detected target amplitude and increases the target SNR from 19.41 dB to 27.89 dB.

Figure 17 Figure 17 .
Figure 17 shows the results of two different reconstruction methods after adaptive cancellation and the cross-ambiguity function.The segmented sliding window reconstruction method based on FCM leads to an enhancement in the detected target amplitude and increases the target SNR from 19.41 dB to 27.89 dB.

Author Contributions:
Conceptualization, H.L., Y.F. and Z.W.; methodology, H.L., Y.Z. and X.Q.; software, H.L. and Y.Z.; validation, H.L., T.S. and Y.F.All authors have read and agreed to the published version of the manuscript.Funding: This research was funded by the National Natural Science Foundation of China, Grant No. 62171029 and 61931015.This research was also supported by the foundation of the National Key Laboratory of Electromagnetic Environment (Grant No. 6142403180204).

Table 1 .
Algorithm running time.

Table 1 .
Algorithm running time.

Table 2 .
The check matrix properties of LDPC codes.

Table 3 .
The number of matrix operations with different coding rates in LDPC decoding.

Table 4 .
Calculation times for different reconstruction methods

Table 4 .
Calculation times for different reconstruction methods.