Pruning median filtering-based vector sensor pulse interference detection and suppression method

Abstract

The invention discloses a vector sensor pulse interference detection and suppression method based on pruning median filtering, and the method comprises the steps: designing a comparison template, and calculating the correlation degree between the comparison template and a signal in a current sliding region through sliding the comparison template in a search region near pulse interference, and selecting the wavelet coefficient of the most similar part near the pulse interference as a replacement wavelet coefficient of the pulse interference. According to the method, discrete wavelet transform is utilized, signals are decomposed at multiple scales, the calculation efficiency is guaranteed, the threshold function of pulse interference detection is optimized based on pruning median filtering, the influence of strong pulses on threshold selection is reduced, and the detection rate of weak pulses is increased.

Description

Technical Field

[0001] The invention belongs to the technical field of sensors, and in particular relates to a vector sensor pulse interference detection and suppression method based on trimmed median filtering. Background Art

[0002] As an advanced underwater sound signal detection device, the vector sensor can simultaneously and accurately measure the sound pressure and particle velocity of a specific point in the sound field. It has a stable dipole pointing characteristic and is not affected by frequency. The co-oscillation vector sensor indirectly obtains the particle velocity by measuring the fluid particle acceleration. Due to its good low-frequency characteristics, it has been widely used in the field of underwater acoustic measurement and detection. However, in actual applications, the co-oscillation vector sensor will produce a series of vibrations under the instigation of the water flow and the obstruction of the cable, resulting in instantaneous and strong pulse interference. Therefore, detecting and suppressing the pulse interference mixed in the output signal of the vector sensor is an important part of improving its receiving performance.

[0003] Short-time Fourier transform (STFT) in the time-frequency domain is a commonly used signal analysis method. This method obtains the frequency domain characteristics of the local signal by selecting a window function, and then slides the window function along the time axis, so that the signal spectrum at each local moment is presented to achieve joint analysis of the signal in the time domain and frequency domain. The length of the short-time Fourier transform window determines the time-frequency resolution. When processing short pulse signals, too high requirements for time resolution will lead to a decrease in frequency domain resolution. Wavelet analysis can change the size of the analysis "window" according to the decomposition scale, so that the wavelet analysis has a larger time resolution at a small scale and a larger frequency resolution at a large scale, which can better capture the characteristics of the signal. Continuous wavelet transform (CWT) can perform global analysis of the signal, provide a complete time-frequency diagram, and observe the changes of the signal at any scale, but CWT needs to continuously calculate multiple scale transforms in the entire signal domain, which has high computational complexity and too much redundant information. Discrete wavelet transform (DWT) uses a filter array to construct multi-resolution time-frequency to perform the Mallat algorithm commonly used in wavelet decomposition. DWT reduces the amount of calculation by sampling and retains most of the important features of the signal. At the same time, due to the convenience of discrete wavelet transform in convolution, the Mallat algorithm can be easily implemented using digital integrated circuits, so DWT is widely used in practice.

[0004] Accurately detecting pulse interference is the prerequisite for suppressing interference. DWT decomposes the signal at different scales to obtain wavelet coefficients of each scale, which can be used to accurately reconstruct the original signal. The key to detecting pulse interference is to select a suitable threshold function to filter out the wavelet coefficients corresponding to the pulse interference for subsequent suppression processing. It is a commonly used pulse detection method to perform wavelet transform on singular signals such as pulses and perform median filtering on the wavelet coefficients to determine the pulse detection threshold and then detect the location of the pulse. However, when a weak pulse appears near a strong pulse, this method is prone to miss detection of the weak pulse. Therefore, choosing a suitable threshold function is crucial to simultaneously detect strong pulse interference and nearby weak pulses.

[0005] After completing the detection of pulse interference, the wavelet coefficients corresponding to the pulse interference are processed. The common processing methods are hard and soft thresholds and their improved threshold algorithms, that is, adjusting the wavelet coefficients corresponding to the pulse interference to the specified range of the set threshold function. This method is simple and efficient, but it destroys the statistical properties of the original signal and distorts it. Therefore, the focus of the research is to appropriately adjust the wavelet coefficients corresponding to the pulse interference to ensure that the original statistical properties of the signal are not destroyed while suppressing the pulse interference. Summary of the invention

[0006] In order to overcome the shortcomings of the prior art, the present invention provides a vector sensor pulse interference detection and suppression method based on trimmed median filtering, designs a comparison template, calculates the correlation between the comparison template and the signal in the current sliding area by sliding the comparison template in the search area near the pulse interference, and selects the wavelet coefficients of the most similar part near the pulse interference as the replacement wavelet coefficients of the pulse interference. The method of the present invention utilizes discrete wavelet transform to decompose the signal at multiple scales, ensures the calculation efficiency, and optimizes the threshold function of pulse interference detection based on trimmed median filtering, reduces the influence of strong pulses on threshold selection, and improves the detection rate of weak pulses.

[0007] The technical solution adopted by the present invention to solve the technical problem is as follows:

[0008] Step 1: Pulse interference detection;

[0009] Step 1-1: Discrete wavelet transform;

[0010] Assume that the one-dimensional continuous-time signal is , is the set of real numbers The function space composed of square integrable functions is:

[0011] (1)

[0012] In the formula, represents a continuous time variable, The corresponding continuous wavelet transform formula is:

[0013] (2)

[0014] In the formula, By mother wavelet The result obtained through stretching and translation is called wavelet sequence; , is the scale factor, is the translation factor; yes The complex conjugate of

[0015] The scale factor of the continuous wavelet and translation factor Along the binary sequence , and perform discretization sampling to obtain discrete wavelet transform:

[0016] (3)

[0017] in, represents the decomposition scale, represents discrete time points, Represents the wavelet sequence after discretization sampling of scale factor and translation factor;

[0018] The Mallat algorithm is given by equation (4):

[0019] (4)

[0020] In the formula, and are the responses of the low-pass filter and high-pass filter, respectively. represents the total number of scales of wavelet decomposition, represents the index of the filter, is the length of the filter, and Respectively represent Approximation coefficients and detail coefficients of discrete wavelet decomposition under different scales; and Respectively represent Approximation coefficients and detail coefficients of discrete wavelet decomposition under different scales;

[0021] Step 1-2: trimmed median filtering;

[0022] A threshold function is set to detect pulse interference. The threshold function is proportional to the smooth envelope of the absolute value of the wavelet coefficients, but is not affected by the pulse interference. The threshold function is defined as:

[0023] (5)

[0024] In the formula, To smooth the envelope of the signal without being affected by pulse interference, is an adjustable factor set empirically;

[0025] If the absolute value of the wavelet coefficient at a certain position is greater than the set threshold, then this position is judged to be a pulse interference, and the position index of all pulse interferences is obtained, which is recorded as :

[0026] (6)

[0027] Threshold function based on trimmed median filter for:

[0028] (7)

[0029] In the formula is an adjustable factor; Indicates that median filtering is performed on the absolute value of the wavelet coefficients. is the number of sample points of a single pulse interference;

[0030] The threshold function calculated by trimmed median filtering is used to detect pulse interference on the wavelet coefficients of each scale after DWT decomposition, and the detection results of each decomposition scale are obtained. , Indicates that pulse interference is detected. Indicates that no pulse interference is detected;

[0031] Step 1-3: Jointly determine pulse interference at multiple scales;

[0032] Perform differential operation on the detected pulse interference index:

[0033] (8)

[0034] in, is the result of the difference operation, the adjacent and The index is recorded as The scale The starting point of the pulse interference and end point ; If The length of the scaled wavelet coefficient is recorded as , then The length of the scaled wavelet coefficient is recorded as ;

[0035] The wavelet coefficients of each scale are linearly interpolated as shown in formula (9) so that the length of the wavelet coefficients of all scales is The length of the wavelet coefficients of each scale is the same; the detection results after interpolation at each scale are recorded as: , the starting and ending point indices are recorded as: , ;

[0036] (9)

[0037] Combine the detection results of all scales and perform voting statistics, recorded as ,but:

[0038] (10)

[0039] Setting Thresholds If the voting result , then it is determined that a pulse interference occurs at discrete time point n; set It means that at least half of the scales detect that there is a pulse here, which is effective pulse interference. The final detection result of the comprehensive judgment of each scale is:

[0040] (11)

[0041] The start and end indexes of the pulse interference in the comprehensive judgment result are: , ; Downsample the comprehensive judgment result so that the length of each scale wavelet coefficient is the same as the length of the initial scale wavelet coefficient; the updated detection result of each scale is recorded as , and the test results of comprehensive judgment of each scale The relationship between them is:

[0042] (12)

[0043] The starting and ending indexes of the pulse are recorded as: and ;

[0044] Step 2: Pulse interference suppression;

[0045] Step 2-1: pulse interference wavelet coefficient processing;

[0046] Wavelet coefficients in the original signal Subtract the wavelet coefficients at the pulse interference from the wavelet coefficients to obtain the wavelet coefficients without pulse interference components. ,Right now:

[0047] (13)

[0048] Calculate the signal envelope after removing the pulse interference, denoted as :

[0049] (14)

[0050] in:

[0051] (15)

[0052] Similarly, the signal envelope of the pulse interference is calculated and recorded as ;

[0053] The following processing is performed on the wavelet coefficients of the pulse interference and the pulse interference vicinity to obtain the processing template of the wavelet coefficients in the pulse interference area :

[0054] (16)

[0055] The template will be processed Multiply the original signal wavelet coefficients to achieve zero wavelet coefficients at the pulse interference, and reset the wavelet coefficients near the pulse interference in proportion:

[0056] (17)

[0057] Step 2-2: Replacement of pulse interference coefficient;

[0058] The starting and ending points of the following pulse interference no longer distinguish the specific decomposition scale , delete from the symbol , remember Scale The starting point of the pulse interference for ,end for ;

[0059] (a) Design a comparison template for the current pulse interference. The starting range of the comparison template is the same as the length of the processing template. ; Indicates the distance between the start and end points of the pulse interference and is used to define the range of the comparison template and the processing template;

[0060] (b) Define a search area, the starting range of the search area is , The processing template slides in the search area, and the energy in the area is calculated by multiplying it with the current wavelet coefficient; Indicates the distance from the start and end points of the pulse interference to both sides, which is used to define the scope of the search area;

[0061] Through (a) and (b), the energy value in the current pulse interference area calculated by the comparison template is obtained, which is recorded as , is a constant; use the processing template to slide in the search area and calculate the energy value in the current template, recorded as ;Compare and Similarity degree, the similarity calculation formula is:

[0062] (18)

[0063] in:

[0064] (19)

[0065] (20)

[0066] (twenty one)

[0067] In the formula is the energy of the current search window; To compare the energy of the template, for Pulse interference is a constant; is the sliding window Points, Is the first of the comparison templates where the pulse is located Points, is the number of sample points of the comparison template;

[0068] After completing the above similarity calculation, the index of the start of the search window with the maximum similarity is obtained , then the starting position index of the pulse interference to be replaced and for:

[0069] (twenty two)

[0070] in, Indicates the starting point of the optimal search window;

[0071] Extended on both sides of the replacement part Sample points and , then replace part of the new area with :

[0072] (twenty three)

[0073] For the processing template in equation (16), find the minimum value and record the index of the first point as , the index of the last point is , then the replacement smoothing template is:

[0074] (twenty four)

[0075] In the formula

[0076] (25)

[0077] The wavelet coefficients of each scale after the final pulse interference suppression are:

[0078] (26).

[0079] Preferably, the This is achieved using a median filter.

[0080] Preferably, the median filter calculation is defined as:

[0081] ;

[0082] The length of the median filter Need to meet:

[0083] .

[0084] A computer program enables a computer to execute the above-mentioned sensor pulse self-noise detection and suppression method.

[0085] An electronic device comprises: a processor and a memory; the memory is used to store a computer program, and the processor is used to execute the computer program stored in the memory, so that the electronic device executes the above-mentioned sensor pulse self-noise detection and suppression method.

[0086] A computer-readable storage medium stores a computer program, which, when executed by a processor, implements the above-mentioned sensor pulse self-noise detection and suppression method.

[0087] A chip includes: a processor for calling and running a computer program from a memory, so that a device equipped with the chip executes the above-mentioned sensor pulse self-noise detection and suppression method.

[0088] A computer program product, the computer program product comprising a computer storage medium, the computer storage medium storing a computer program, the computer program comprising instructions executable by at least one processor, and when the instructions are executed by the at least one processor, the above-mentioned sensor pulse self-noise detection and suppression method is implemented.

[0089] The beneficial effects of the present invention are as follows:

[0090] The method of the present invention utilizes discrete wavelet transform to decompose the signal at multiple scales, thereby ensuring computational efficiency, and optimizes the threshold function of pulse interference detection based on trimmed median filtering, thereby reducing the influence of strong pulses on threshold selection and improving the detection rate of weak pulses. BRIEF DESCRIPTION OF THE DRAWINGS

[0091] Figure 1 Schematic diagram of Mallat algorithm.

[0092] Figure 2 A diagram showing the relationship between processing templates, replacement templates, and search areas.

[0093] FIG3( a ) is a time domain waveform diagram of a Y channel signal according to an embodiment of the present invention.

[0094] FIG3( b ) is a time-frequency diagram of a Y channel signal according to an embodiment of the present invention.

[0095] FIG. 4( a ) is a diagram showing a median filter cA4 threshold curve and a detection result diagram in an embodiment of the present invention.

[0096] FIG4( b ) is a diagram showing a trimmed median filter cA4 threshold curve and a detection result diagram according to an embodiment of the present invention.

[0097] FIG. 5( a ) is a diagram showing a trimmed median filter cA4 threshold curve and a detection result diagram according to an embodiment of the present invention.

[0098] FIG5( b ) is a diagram showing a trimmed median filter cD4 threshold curve and a detection result diagram according to an embodiment of the present invention.

[0099] FIG. 5( c ) is a diagram showing a cD3 threshold curve and a detection result of a trimmed median filter according to an embodiment of the present invention.

[0100] FIG5( d ) is a diagram showing a trimmed median filter cD2 threshold curve and a detection result diagram according to an embodiment of the present invention.

[0101] Figure 6 Detection results at each scale and multi-scale superposition results of the embodiment of the present invention.

[0102] Figure 7 A diagram of multi-scale comprehensive decision detection results according to an embodiment of the present invention.

[0103] Figure 8 The signal envelope and pulse interference envelope of the embodiment of the present invention.

[0104] Fig. 9 Schematic diagram of a processing template according to an embodiment of the present invention.

[0105] Fig.10 Schematic diagram of replacing a template according to an embodiment of the present invention.

[0106] Fig.11The time domain diagram of the original signal in the embodiment of the present invention and the time domain diagram after median filtering and processing by the method.

[0107] Fig.12 The time-frequency diagram of the original signal in the embodiment of the present invention and the time-frequency diagram after being processed by the method. DETAILED DESCRIPTION

[0108] The present invention is further described below in conjunction with the accompanying drawings and embodiments.

[0109] In order to solve the problem of instantaneous and strong pulse interference caused by the sudden vibration of the vector sensor under the instigation of water flow and the obstruction of the cable, an improved median filtering method - trimmed median filtering is used to detect strong and weak pulses, improve the detection rate of weak pulses, and solve the problem of missed detection of weak pulses near strong pulses.

[0110] In order to improve the signal quality after pulse interference suppression, the present invention designs a comparison template, which calculates the correlation between the comparison template and the signal in the current sliding area by sliding the comparison template in the search area near the pulse interference, and selects the wavelet coefficients of the most similar part near the pulse interference as the replacement wavelet coefficients of the pulse interference.

[0111] The technical solution of the present invention is composed of two processing modules, namely a pulse interference detection module and a pulse interference suppression module.

[0112] 1) Pulse interference detection;

[0113] Pulse interference detection can be divided into the following two steps:

[0114] 1.1 Discrete Wavelet Transform;

[0115] Assume a one-dimensional continuous-time signal , Real number set The function space composed of square integrable functions is:

[0116]

[0117] The corresponding continuous wavelet transform formula is:

[0118]

[0119] The scale factor of the continuous wavelet and translation factor Along the binary sequence , and perform discretization sampling to obtain discrete wavelet transform:

[0120]

[0121] It can be seen that the discrete wavelet transform has an effect on the translation factor A secondary extraction is performed. The Mallat algorithm is given by equation (4):

[0122]

[0123] Mallat algorithm Figure 1 shown.

[0124] 1.2 Pruning median filtering;

[0125] A threshold function is set to detect pulse interference. The threshold function is proportional to the smooth envelope of the absolute value of the wavelet coefficients, but is not affected by the pulse interference. The threshold function is defined as:

[0126]

[0127] If the absolute value of the wavelet coefficient at a certain position is greater than the set threshold, it is considered that this position is a pulse interference, and the position index of all pulses is obtained, which is recorded as :

[0128]

[0129] The threshold function used in formula (5) is It is usually implemented using a median filter because it preserves step function type signals and is robust to pulse interference. The traditional median filter calculation is defined as:

[0130]

[0131] The length of the median filter Need to meet:

[0132]

[0133] Traditional median filtering can remove strong pulse interference with large amplitude. However, when there is a weak pulse near the strong pulse, this method is prone to miss the weak pulse. The present invention proposes a threshold selection method based on trimmed median filtering. First, a certain proportion of extreme value data in the wavelet coefficients is removed to reduce the influence of strong pulses on threshold selection. Then, median filtering is performed on the wavelet coefficients to determine the detection threshold, thereby improving the detection rate of weak pulses. Threshold function based on trimmed median filtering for:

[0134]

[0135] In the formula Is an adjustable factor. It is used to sort the tail of the median filter in ascending order. The wavelet coefficients corresponding to the strong pulse interference of the proportion are removed to reduce the influence of the strong pulse on the threshold selection, so that the weak pulse can be detected and the detection rate of the weak pulse can be improved. The threshold function calculated by trimming median filtering is used to realize pulse interference detection for the wavelet coefficients of each scale after DWT decomposition, and the detection results of each decomposition scale are obtained. , Indicates that pulse interference is detected. Indicates that no pulse interference was detected.

[0136] 1.3 Multi-scale joint determination of pulse interference;

[0137] In order to reduce the missed detection rate and effectively suppress weak pulses, the following joint processing is performed:

[0138] Perform differential operation on the detected pulse interference index:

[0139]

[0140] Due to the characteristics of the Mallat algorithm, the length of the wavelet coefficients is halved each time the decomposition is performed. The length of the scaled wavelet coefficient is recorded as , then The length of the scaled wavelet coefficient is recorded as .

[0141] Since the lengths of wavelet coefficients at different scales are different, the wavelet coefficients at each scale are linearly interpolated as shown in formula (9) so that the lengths of wavelet coefficients at all scales are equal to The length of the wavelet coefficients of each scale is the same. The detection results after interpolation at each scale are recorded as: , the starting and ending point indices are recorded as: , .

[0142]

[0143] Combine the detection results of all scales and perform voting statistics, recorded as ,but:

[0144]

[0145] Setting Thresholds If the voting result , then it is determined that a pulse interference occurs at discrete time point n; set It means that at least half of the scales detect the presence of pulses here to be considered effective pulse interference. Therefore, the final detection result of the comprehensive judgment of each scale is:

[0146]

[0147] The start and end indexes of the pulse interference in the comprehensive judgment result are: , The comprehensive judgment result is downsampled so that the length of each scale wavelet coefficient is the same as the length of the initial scale wavelet coefficient. The updated detection results of each scale are recorded as , and the test results of comprehensive judgment of each scale The relationship between them is:

[0148]

[0149] The starting and ending indexes of the pulse are recorded as and .

[0150] 2) Pulse interference suppression;

[0151] Pulse interference suppression can be divided into the following two steps:

[0152] 2.1 Pulse interference wavelet coefficient processing;

[0153] Complete the detection of pulse interference in 1) and obtain the detection results of pulse interference of various scales In order to avoid the artifacts (mutation signals or noise introduced by improper processing) generated near the pulse interference by simply setting the wavelet coefficients at the location of the pulse interference to zero, which will affect the effect of pulse interference suppression, the wavelet coefficients near the pulse interference also need to be replaced. Subtract the wavelet coefficients at the pulse interference from the wavelet coefficients to obtain the wavelet coefficients without pulse interference components. ,Right now:

[0154]

[0155] Calculate the signal envelope after removing the pulse interference, denoted as :

[0156]

[0157] in

[0158]

[0159] Similarly, the signal envelope of the pulse interference is calculated and recorded as The following processing is performed on the wavelet coefficients of the pulse interference and the pulse interference vicinity to obtain the processing template of the wavelet coefficients in the pulse interference area :

[0160]

[0161] Multiply the processing template with the original signal wavelet coefficient to achieve zero wavelet coefficients at the pulse interference and reset the wavelet coefficients near the pulse interference in proportion:

[0162]

[0163] 2.2 Replacement of pulse interference coefficient;

[0164] In order to find the most suitable wavelet coefficients to replace the wavelet coefficients of the pulse interference suppressed in 2.1, the following processing steps are designed. For the sake of symbol simplicity, the starting point and end point of the pulse interference below no longer distinguish the specific decomposition scale , delete from the symbol , remember Scale The starting point of the pulse interference for ,end for .

[0165] (a) Design a comparison template for the current pulse interference. The starting range of the comparison template is the same as the length of the processing template above. This step aims to calculate The energy of the region without pulse interference is used as a reference for subsequent calculation of the energy of the replacement region.

[0166] (b) Define a search area, the starting range of the search area is , The processing template slides in the search area, and the energy in the area is calculated by multiplying it with the current wavelet coefficient.

[0167] Through the above two steps, the energy value of the current pulse interference area calculated by the comparison template is obtained, which is recorded as , is a constant; use the processing template to slide in the search area and calculate the energy value in the current template, recorded as .Compare and Similarity degree, the similarity calculation formula is:

[0168]

[0169] in

[0170]

[0171]

[0172]

[0173] After completing the above similarity calculation, the index of the start of the search window with the maximum similarity is obtained , then the starting position index of the pulse interference to be replaced and for:

[0174]

[0175] However, to minimize artifacts, the replaced part needs to blend smoothly at the boundary. Therefore, the replaced part needs to be extended on both sides to overlap in order to achieve smooth blending. Sample points and , then replace part of the new area with :

[0176]

[0177] In order to smoothly replace the wavelet coefficients of the pulse interference area, a replacement smoothing template needs to be designed to achieve arrive Smooth replacement. For the processing template in equation (16), find the minimum value. The point with the minimum value is usually not unique, and the index of the first point is , the index of the last point is , then the replacement smoothing template is:

[0178]

[0179] In the formula

[0180]

[0181] The relationship diagram of processing template, replacement template and search area is as follows Figure 2 As shown. Then the wavelet coefficients of each scale after the final pulse interference suppression are:

[0182]

[0183] Example:

[0184] In order to verify the effective performance of the present invention, the embodiments of the present invention verify the actual effect of the above-mentioned pulse interference detection and suppression method on actual data.

[0185] The data comes from the Y channel of a domestic isotropic vector sensor. The target signal in the data is CW signal, the original sampling rate of the signal , the signal duration is , downsample the signal, and the sampling rate after downsampling is The duration of the pulse interference is , the number of sampling points corresponding to the time domain is .

[0186] The signal time domain waveform is shown in Figure 3 (a), and the time-frequency diagram is shown in Figure 3 (b). From the time-frequency diagram in Figure 3 (b), it can be seen that the pulse interference continues from the low frequency band to the high frequency band. Scale discrete wavelet transform, according to the discrete wavelet decomposition characteristics, the length of each scale pulse interference after decomposition is , The signal frequency range corresponding to the wavelet coefficients decomposed at each scale is , where the approximation coefficient The corresponding frequency range is , detail coefficient The corresponding frequency range is .by Scaling wavelet approximation coefficients Take the median filter and pruning coefficients as an example to calculate them respectively. The threshold function of the trimmed median filter is used to obtain the detection effect as shown in Figure 4 (a) and Figure 4 (b). The upper figure of Figure 4 (a) is the threshold function calculated by the median filter, and the lower figure of Figure 4 (a) is the detection result. The upper figure of Figure 4 (b) is the threshold function calculated by the trimmed median filter, and the lower figure of Figure 4 (b) is the detection result. By comparison, it can be seen that the threshold function calculated by the trimmed median filter can detect more weak pulse interference and improve the detection rate of weak pulse interference.

[0187] Figure 5(a), Figure 5(b), Figure 5(c), and Figure 5(d) respectively show the threshold function curves calculated by the trimmed median filtering method for cA4, cD4, cD3, and cD2 and the pulse interference detection results. Figure 6 The results of each scale detection and the superposition of all scale detections are shown. In the multi-scale joint determination of pulse interference step, according to the scale of wavelet decomposition , we can know the threshold of multi-scale joint , that is, it is considered that voting statistics are performed on each position at four scales, and the statistical value reaches If the number is less than two, it is considered as a false alarm and is discarded. Figure 7 The figure above shows the result after superposition of all scale detections, and the figure below shows the pulse interference detection result after multi-scale joint determination.

[0188] Suppress and replace the detected pulse interference. Draw the signal envelope and the pulse interference envelope together Figure 8 In the above example, the processing template of (16) is designed according to the envelope of the signal and pulse interference. Taking the first pulse as an example, the processing template of the wavelet coefficients in the pulse interference region is plotted. See Fig. 9, processing template extension length The length of the pulse search area on both sides . Pulse interference replacement template such as Fig.10 As shown, the extended length of the replacement template The signal time domain diagram and the original signal waveform diagram after the pulse interference is detected and suppressed by two threshold functions are shown in the figure below. Fig.11 As shown, Fig.11 The above figure is the time domain diagram of the original signal. Fig.11 The middle figure is the time domain figure after the threshold function detection and suppression calculated by the median filter. Fig.11 The following figure shows the time domain diagram of pulse interference detection and suppression by this method. Compared with the median filtering method, this method has a better pruning coefficient. , 10 additional weak pulses were suppressed, marked by dashed arrows. This proves that this method reduces the missed detection rate and effectively suppresses weak pulses. The time-frequency diagram of the signal before and after processing by this method is compared. Fig.12 As shown, Fig.12 The left picture is the time-frequency diagram of the original signal. Fig.12 The right figure is the time-frequency diagram after pulse suppression using this method.

[0189] According to the embodiment, it can be considered that the pulse interference detection and suppression method proposed in the present invention is feasible. The discrete wavelet transform is used to reduce the amount of calculation, which can meet the requirements for pulse interference detection and optimize the threshold function of pulse interference detection based on the trimmed median filter, reduce the influence of strong pulses on the threshold selection, and improve the detection rate of weak pulses. Actual data processing shows that under the condition that the trimming ratio is between 10% and 20%, the method of the present invention significantly improves the effect of pulse interference detection.

Claims

1. A vector sensor pulse interference detection and suppression method based on trimmed median filtering, characterized in that: The steps include: Step 1: Pulse interference detection; Step 1-1: Discrete wavelet transform; Assume that the one-dimensional continuous-time signal is , is the set of real numbers The function space composed of square integrable functions is: ； In the formula, represents a continuous time variable, The corresponding continuous wavelet transform formula is: ； In the formula, By mother wavelet The result obtained through stretching and translation is called wavelet sequence; , is the scale factor, is the translation factor; yes The complex conjugate of The scale factor of the continuous wavelet and translation factor Along the binary sequence , and perform discretization sampling to obtain discrete wavelet transform: ； in, represents the decomposition scale, represents discrete time points, Represents the wavelet sequence after discretization sampling of scale factor and translation factor; The Mallat algorithm is given by equation (4): ； In the formula, and are the responses of the low-pass filter and high-pass filter, respectively. represents the total number of scales of wavelet decomposition, represents the index of the filter, is the length of the filter, and Respectively represent Approximation coefficients and detail coefficients of discrete wavelet decomposition under different scales; and Respectively represent Approximation coefficients and detail coefficients of discrete wavelet decomposition under different scales; Step 1-2: trimmed median filtering; A threshold function is set to detect pulse interference. The threshold function is proportional to the smooth envelope of the absolute value of the wavelet coefficients, but is not affected by the pulse interference. The threshold function is defined as: ； In the formula, To smooth the envelope of the signal without being affected by pulse interference, is an adjustable factor set empirically; If the absolute value of the wavelet coefficient at a certain position is greater than the set threshold, then this position is judged to be a pulse interference, and the position index of all pulse interferences is obtained, which is recorded as : ； Threshold function based on trimmed median filter for: ； In the formula is an adjustable factor; Indicates that median filtering is performed on the absolute value of the wavelet coefficients. is the number of sample points of a single pulse interference; The threshold function calculated by trimmed median filtering is used to detect pulse interference on the wavelet coefficients of each scale after DWT decomposition, and the detection results of each decomposition scale are obtained. , Indicates that pulse interference is detected. Indicates that no pulse interference is detected; Step 1-3: Jointly determine pulse interference at multiple scales; Perform differential operation on the detected pulse interference index: ； in, is the result of the difference operation, the adjacent and The index of The scale The starting point of the pulse interference and end point ; If The length of the scaled wavelet coefficient is recorded as , then The length of the scaled wavelet coefficient is recorded as ; The wavelet coefficients of each scale are linearly interpolated as shown in formula (9) so that the length of the wavelet coefficients of all scales is The length of the wavelet coefficients of each scale is the same; the detection results after interpolation at each scale are recorded as: , the starting and ending point indices are recorded as: , ; ； Combine the detection results of all scales and perform voting statistics, recorded as ,but: ； Setting Thresholds If the voting result , then determine the discrete time point There is pulse interference; set It means that at least half of the scales detect that there is a pulse here, which is effective pulse interference. The final detection result of the comprehensive judgment of each scale is: ； The start and end indexes of the pulse interference in the comprehensive judgment result are: , ; Downsample the comprehensive judgment result so that the length of each scale wavelet coefficient is the same as the length of the initial scale wavelet coefficient; the updated detection result of each scale is recorded as , and the test results of comprehensive judgment of each scale The relationship between them is: ； The starting and ending indexes of the pulse are recorded as: and ; Step 2: Pulse interference suppression; Step 2-1: Pulse interference wavelet coefficient processing; Wavelet coefficients in the original signal Subtract the wavelet coefficients at the pulse interference from the wavelet coefficients to obtain the wavelet coefficients without pulse interference components. ,Right now: ； Calculate the signal envelope after removing the pulse interference, denoted as : ； in: ； Similarly, the signal envelope of the pulse interference is calculated and recorded as ; The following processing is performed on the wavelet coefficients of the pulse interference and the pulse interference vicinity to obtain the processing template of the wavelet coefficients in the pulse interference area : ； The template will be processed Multiply the original signal wavelet coefficients to achieve zero wavelet coefficients at the pulse interference, and reset the wavelet coefficients near the pulse interference in proportion: ； Step 2-2: Replacement of pulse interference coefficients; The starting and ending points of the following pulse interference no longer distinguish the specific decomposition scale , delete from the symbol , remember Scale The starting point of the pulse interference for ,end for ; (a) Design a comparison template for the current pulse interference. The starting range of the comparison template is the same as the length of the processing template. ; Indicates the distance between the start and end points of the pulse interference and is used to define the range of the comparison template and the processing template; (b) Define a search area, the starting range of the search area is , The processing template slides in the search area, and the energy in the area is calculated by multiplying it with the current wavelet coefficient; Indicates the distance from the start and end points of the pulse interference to both sides, which is used to define the scope of the search area; Through (a) and (b), the energy value in the current pulse interference area calculated by the comparison template is obtained, which is recorded as , is a constant; use the processing template to slide in the search area and calculate the energy value in the current template, recorded as ;Compare and Similarity degree, the similarity calculation formula is: （18）； in: （19）； （20）； （21）； In the formula is the energy of the current search window; To compare the energy of the template, for Pulse interference is a constant; is the sliding window Points, Is the first of the comparison templates where the pulse is located Points, is the number of sample points of the comparison template; After completing the above similarity calculation, the index of the start of the search window with the maximum similarity is obtained , then the starting position index of the pulse interference to be replaced and for: （22）； in, Indicates the starting point of the optimal search window; Extended on both sides of the replacement part Sample points and , then replace part of the new area with : （23）； For the processing template in equation (16), find the minimum value and record the index of the first point as , the index of the last point is , then the replacement smoothing template is: （24）； In the formula （25）； The wavelet coefficients of each scale after the final pulse interference suppression are: （26）。 2. The method for detecting and suppressing pulse interference of a vector sensor based on trimmed median filtering according to claim 1 is characterized in that: Said This is achieved using a median filter.

3. The method for detecting and suppressing pulse interference of a vector sensor based on trimmed median filtering according to claim 2 is characterized in that: The median filter calculation is defined as: ； The length of the median filter Need to meet: 。 4. A computer program, characterized in that The computer program enables a computer to execute the method according to any one of claims 1 to 3.

5. An electronic device, characterized in that: include: Processor and memory; The memory is used to store a computer program, and the processor is used to execute the computer program stored in the memory, so that the electronic device executes the method as claimed in any one of claims 1 to 3.

6. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method according to any one of claims 1 to 3 is implemented.

7. A chip, characterized in that: include: A processor, configured to call and run a computer program from a memory, so that a device equipped with the chip executes a method as claimed in any one of claims 1 to 3.

8. A computer program product, characterized in that The computer program product comprises a computer storage medium storing a computer program, wherein the computer program comprises instructions executable by at least one processor, and when the instructions are executed by the at least one processor, the method according to any one of claims 1 to 3 is implemented.

Images