An anti-interference method in the detection of radiated noise line spectrum

Abstract

The present invention discloses an anti-interference method in radiation noise line spectrum detection, comprising the following steps: (1) initializing data processing parameters, a radiation noise line spectrum information matrix and a determination threshold; (2) reading the line spectrum information obtained from a frame of line spectrum detection and updating the radiation noise line spectrum information matrix, where the line spectrum information in the radiation noise line spectrum information matrix includes line spectrum frequency, number of detections, frame number of first occurrence, frame number of last occurrence, and interference confidence level; (3) determining the value of the interference confidence level of the current line spectrum according to the number of detections of the current line spectrum and its proportion; (4) judging whether the line spectrum in the radiation noise line spectrum information matrix and the current line spectrum form a multiple frequency relationship and updating the line spectrum interference confidence level in the radiation noise line spectrum information matrix; (5) when the number of processed frames reaches a certain amount, judging whether the current line spectrum is an interference line spectrum according to the interference confidence level in the radiation noise line spectrum information matrix.

Description

Technical Field

[0001] The present invention relates to an anti-interference technology in the detection of radiated noise line spectra, and belongs to the technical field of feature extraction and recognition of underwater acoustic target incident noise signals. Background Art

[0002] Underwater acoustic systems such as sonars use hydrophone arrays to receive underwater acoustic signals and achieve functions such as underwater target detection, positioning, and recognition. Among them, target recognition and positioning are particularly important links, which are related to the correct judgment and decision-making of commanders. Due to many interferences existing in the equipment itself and the environment, the anti-interference technology in the detection of radiated noise line spectra plays an important supporting role in distinguishing targets and interferences.

[0003] Underwater acoustic target recognition mainly distinguishes target types and species information by extracting target feature quantities. Target feature information is the information that can accurately and simply indicate the target state and identity and is contained in or can be extracted from the target original data. Underwater acoustic targets mainly include features such as noise, motion, wake, and geometric structure, and the feature information of different underwater acoustic targets is different. Submerged vehicles such as submarines and surface ships inevitably generate radiated noise during navigation. Therefore, the radiated noise of underwater targets is the main information source for the current passive sonar work, and the characteristics of radiated noise are the main basis for passive sonar target recognition.

[0004] However, due to interference factors such as structural vibration and circuit noise, the signals received by hydrophones may have interferences with harmonic relationships or continuous fixed interferences. For the above two types of interferences, the present invention performs statistical analysis on the line spectrum information matrix output by the passive sonar processing, calculates the interference confidence level and the duration of interference of the line spectra in the line spectrum information matrix, and determines whether the line spectrum is caused by interference by judging whether the interference confidence level and the interference duration are greater than the threshold values. Summary of the Invention

[0005] Technical Problem: The purpose of the present invention is to provide an anti-interference method in the detection of radiated noise line spectra for the anti-interference problem in the detection of radiated noise line spectra. This method first records all detected line spectrum information (line spectrum frequency, number of detections, first appearance frame number, last appearance frame number) in the line spectrum information matrix, and sets the interference confidence value according to the number of detections, the proportion of the number of detections, and the harmonic relationship. When the line spectrum information matrix accumulates to a certain duration, it is judged whether the line spectrum source is an interference signal according to the confidence level and the duration of the current line spectrum. While ensuring that the line spectrum detection probability remains unchanged, reduce the interference with harmonic relationships and the constant interference in the environment, so as to reduce the false alarm probability and meet the anti-interference requirements.

[0006] Technical Solution: An anti-interference method in the detection of radiated noise line spectra adopted by the present invention includes the following steps:

[0007] Step 1, initialize the data processing parameters, the radiation noise line spectrum detection information matrix, the maximum confidence value, and the decision threshold;

[0008] Step 2, read a frame of radiation noise line spectrum detection information, including the line spectrum frequency, the number of detections, the frame number of the first occurrence, and the frame number of the last occurrence, update it to the line spectrum information matrix, and determine the interference confidence value of each line spectrum according to the number of detections of the line spectrum and its proportion;

[0009] Step 3, traverse each line spectrum in the current frame of radiation noise line spectrum detection information. If the number of line spectra in the radiation noise line spectrum information matrix that form a multiple-frequency relationship with the line spectra in the current frame of radiation noise line spectrum detection information exceeds the threshold, set the interference confidence value of the corresponding line spectrum and its multiple-frequency line spectra in the radiation noise line spectrum information matrix to the maximum confidence value, otherwise the interference confidence value remains the initial value;

[0010] Step 4, determine whether the number of processed frames reaches the processed frame threshold. If it reaches, go to Step 5, otherwise go to Step 6;

[0011] Step 5, traverse each line spectrum in the current frame of radiation noise line spectrum detection information. According to the line spectrum interference confidence value corresponding to the line spectrum in the current frame of radiation noise line spectrum detection information in the radiation noise line spectrum information matrix, determine whether the line spectrum in the current frame of radiation noise line spectrum detection information is an interference line spectrum. If the interference confidence value is greater than the confidence threshold and the duration is greater than the duration threshold, then determine that the line spectrum is interference, otherwise determine that the line spectrum is a suspected target line spectrum;

[0012] Step 6, if it is selected to end the interference processing, then end the processing, otherwise return to Step 2 to continue reading the line spectrum detection information of the next frame after determination.

[0013] Among them,

[0014] The following parameters are initialized in Step 1:

[0015] Step 1.1, set the sampling rate f s 、the frequency resolution Δf and the detection duration T of each frame, the lower limit f l of the radiation noise line spectrum detection frequency and the upper limit f h of the detection frequency;

[0016] Step 1.2, set the matrix G for storing the radiation noise line spectrum information. G is a matrix with N rows and 5 columns, where Among them is the symbol for rounding up; the first column of G is the line spectrum frequency f n , the second column is the number of line spectrum detections M n , the third column is the frame number F fn, the fourth column is the frame number L where the line spectrum appears for the last time fn , the fifth column is the line spectrum interference confidence level B n ; the initial value of each element in the first column of the line spectrum information matrix is f l +(n - 1)×Δf, n = 1, 2, 3, …, N, and the initial values of the other four columns are all assigned 0;

[0017] Step 1.3, set the maximum value K of the number of line spectra detected in each frame m , the frame number threshold M λ , the current frame number is r, and let r = 1;

[0018] Step 1.4, set the detection times threshold γ and the detection times ratio threshold α of the initialization confidence level of the radiation noise line spectrum i , β i , i = 1, 2, 3, 4, 5;

[0019] Step 1.5, set the confidence level threshold B of the radiation noise line spectrum λ , the maximum confidence level B max , the number of multiple frequency relationship thresholds N λ , the upper limit value Q of the number of multiple frequencies;

[0020] Step 1.6, set the duration threshold T1 and the duration threshold T2 of the radiation noise line spectrum;

[0021] Step 1.7, set the allowable error f when judging the equality or harmonic relationship of the radiation noise line spectrum ε .

[0022] The specific steps of the said Step 2 include the following steps:

[0023] Step 3.1, read the radiation noise line spectrum detection information of the r-th frame, that is, the line spectrum frequency f k , k = 1, 2, 3, …, K, K represents the total number of line spectra detected in the current frame, and K ≤ K m ;

[0024] Step 3.2, for the K line spectra detected in the r-th frame, if there exists a positive integer p, and 1 ≤ p ≤ N, such that the elements in the p-th row of the matrix G satisfy |f k -f p |<f ε , then let M p = M p + 1, L fp = r; at this time, if F fp = 0, then let F fp = r, otherwise F fp remains unchanged; if α i <M p / L fp <βi and M p > γ, i = 1, 2, 3, 4, 5, then let B p = i.

[0025] The specific steps of step 3 are as follows:

[0026] Step 4.1, for each line spectrum f in matrix G n , use N n to represent the number of line spectra in matrix G that form a multiple frequency relationship with f n ; if L fn > 0, let N n = 0, variable j = 1, and the Q elements of array H are initialized to 0;

[0027] Step 4.2, if there exists a positive integer q, 2 ≤ q ≤ Q + 1, such that the elements in the j-th row of matrix G satisfy |q × f n - f j | < f ε and M j > 0, then let N n = N n + 1, and let the q-th element of H, H(q) = j;

[0028] Step 4.3, if j < N, then let j = j + 1 and return to step 4.2;

[0029] Step 4.4, if N n ≥ N λ , then let B n = B max , and let all B H(q) corresponding to H(q) ≠ 0 be B max .

[0030] The specific steps of step 4 are as follows: Determine whether the number of detected frames reaches the threshold M λ , that is, if r ≥ M λ then proceed to step 5, otherwise enter step 6.

[0031] The specific steps of step 5 are as follows: For the K line spectra detected in the r-th frame, in the matrix G of radiation noise line spectrum information, the line spectrum frequency is f p , and p is the row number corresponding to this line spectrum in matrix G:

[0032] Step 6.1, for the K line spectra f k detected in the r-th frame, k = 1, 2, 3,..., K, if there exists a positive integer l, and 1 ≤ l ≤ N, such that the elements in the l-th row of matrix G satisfy |f k - f l | < f ε , then proceed to step 6.2;

[0033] Step 6.2, if B l = B max and (L fl - F fl ) × T > T1, then it is determined that the radiated noise line spectrum with the line spectrum frequency of f k is an interference line spectrum; otherwise, the next judgment is carried out;

[0034] Step 6.3, if B l > B λ and (L fl - F fl ) × T > T2, then it is determined that the radiated noise line spectrum with the line spectrum frequency of f k is an interference line spectrum; otherwise, it is determined that the radiated noise line spectrum is a suspected target line spectrum.

[0035] The specific steps of Step 6 are as follows: If it is selected to end the interference processing, the anti-interference processing is terminated, and the radiated noise line spectrum information matrix no longer changes; otherwise, let r = r + 1, return to Step 2, and continue to update the radiated noise line spectrum information matrix and perform anti-interference processing.

[0036] Beneficial effects: Compared with the prior art, the method disclosed in the present invention has the following advantages: 1. Without changing the line spectrum detection probability, the interference with harmonic relationship and the constant interference in the environment are reduced to reduce the false alarm probability and meet the anti-interference requirements; 2. The setting of the line spectrum information matrix can not only record the radiated noise line spectrum situation for a period of time, but also facilitate the suppression of the line spectrum with a specific frequency in some environments. Brief Description of the Drawings

[0037] Figure 1 is the implementation flowchart of the method of the present invention.

[0038] Figure 2 is the interference confidence of the radiated noise line spectrum in the embodiment.

[0039] Figure 3 is the single-frame detection result of the radiated noise line spectrum in the embodiment.

[0040] Figure 4 is the multi-frame detection result of the radiated noise line spectrum after anti-interference processing in the embodiment.

[0041] Figure 5 is the detection accuracy rate and false alarm rate of the radiated noise line spectrum after anti-interference processing in the embodiment.

[0042] Figure 6 is the multi-frame detection result of the radiated noise line spectrum without anti-interference processing.

[0043] Figure 7The detection accuracy rate and false alarm rate of the radiation noise line spectrum without anti-interference processing. Detailed implementation mode

[0044] The present invention will be further described in detail below in conjunction with the accompanying drawings and specific implementation modes:

[0045] The present invention analyzes the long sequence line spectrum information output by the passive sonar processing, counts the line spectrum information in the sequence to obtain a line spectrum information matrix, and finally combines whether the confidence level, detection ratio at this time, and duration in the line spectrum information matrix reach the set threshold rules to determine whether the line spectrum is interference or a target.

[0046] Embodiment 1:

[0047] The present invention will be further clarified below in conjunction with the accompanying drawings and specific implementation modes.

[0048] As Figure 1 shown, the anti-interference method in the radiation noise line spectrum detection of the present invention includes the following steps:

[0049] Step 1, initialize the data processing parameters, the radiation noise line spectrum information matrix, and the determination threshold, specifically:

[0050] (1) Set the sampling rate fs = 2 kHz, the frequency resolution Δf = 0.06 Hz, and the detection duration T of each frame = 32 s. The lower limit f of the radiation noise line spectrum detection frequency l = 50 Hz, and the upper limit f of the detection frequency h = 500 Hz. During the lake experiment, the target frequency values are 81 Hz, 165 Hz, 204 Hz, 311 Hz, and 500 Hz.

[0051] (2) Set the radiation noise line spectrum information matrix to be saved as an 825×5 matrix G. Among them represents the total capacity of the line spectrum. The line spectrum information matrix has a total of 5 columns. The first column of G is the line spectrum frequency f n , the second column is the line spectrum detection times M n , the third column is the frame number F when the line spectrum first appears fn , the fourth column is the frame number L when the line spectrum last appears fn , and the fifth column is the line spectrum interference confidence level B n . The first column of the line spectrum information matrix represents all the line spectrum frequencies within the detection frequency range. Therefore, the initial value of this column is an array from 50 Hz to 500 Hz with an interval of 0.03 Hz. The other four columns are real-time updated quantities, so the initial values are assigned 0;

[0052] (3) Set the maximum value K of the number of line spectra detected in each frame m = 100, and the initial accumulated frame number threshold M λ= 20 frames, i.e., 640 s, the current frame number is r, and let r = 1;

[0053] (4) Set the detection times threshold γ of the initial confidence level of the radiated noise line spectrum to 20 times, the detection times ratio threshold α i , β i , i = 1, 2, 3, 4, 5, and their values are shown in Table 1 below;

[0054] Table 1 Confidence level values of the detection times ratio under different thresholds

[0055]

[0056] (5) Set the confidence level threshold B of the radiated noise line spectrum λ = 3, the maximum confidence level B max = 6, the number of octave relationship thresholds N λ = 6, the upper limit value Q of the number of octaves = 8;

[0057] (6) Set the duration threshold T1 of the radiated noise line spectrum to 320 s and the duration threshold T2 to 1000 s;

[0058] (7) Set the allowable error f when judging the equality or harmonic relationship of the radiated noise line spectrum ε = 0.01 Hz;

[0059] Step 2 is specifically as follows:

[0060] (1) Read the radiated noise line spectrum detection information of the r-th frame, i.e., the line spectrum frequency f k , k = 1, 2, 3,..., K, where K represents the total number of line spectra detected in the current frame, and K ≤ 100;

[0061] (2) For the K line spectra detected in the r-th frame, if there exists a positive integer p, and 1 ≤ p ≤ 825, such that the elements in the p-th row of the matrix G satisfy |f k - f p | < 0.01 Hz, then let M p = M p + 1, L fp = r; At this time, if F fp = 0, then let F fp = r, otherwise F fp remains unchanged; If α i < M p / L fp < β i , and M p > 20, i = 1, 2, 3, 4, 5, then let B p = i;

[0062] Step 3 is specifically as follows:

[0063] (1) For each line spectrum f in the matrix G n , use N n Indicates the matrix G with f n The number of line spectra that form the frequency relationship; if L fn >0, let N n =0, variable j=1, and the 9 elements of array H are initialized to 0;

[0064] (2) If there exists a positive integer q, 2≤q≤9, such that the elements in the jth row of the matrix G satisfy |q×f n -f j |<0.01Hz and M j > 0, then let N n =N n +1, and let the qth element of H H(q) = j;

[0065] (3) If j<825, set j=j+1 and return to the previous step;

[0066] (4) If N n ≥N λ , then let B n =B max , and let all B corresponding to H(q)≠0 H(q) =B max .

[0067] Step 4 is as follows:

[0068] Determine whether the number of detection frames reaches the threshold of 20 frames, that is, if r ≥ 20, proceed to step 5, otherwise proceed to step 6; so far, the method completes the update of the line spectrum interference confidence in the line spectrum information matrix, such as Figure 2 and Figure 3 As shown, it can be seen that the frequencies with harmonic relationship (62.5Hz double frequency line spectrum) and constant interference frequency (for example, around 30Hz) are effectively shielded.

[0069] Step 5 is as follows:

[0070] For the K line spectra detected in the rth frame, the line spectrum frequency corresponding to the matrix G of the radiation noise line spectrum information is f p , p is the row coordinate corresponding to the line spectrum in the matrix G:

[0071] (1) The K line spectra f detected in the rth frame k , k=1,2,3,…,K, if there exists a positive integer l, and 1≤l≤N, such that the elements of the lth row of the matrix G satisfy |f k -f l |<0.01Hz, proceed to the next step;

[0072] (2) If Bl = B max and (L fl - F fl ) × 32s > 320s, then it is determined that the line spectrum frequency of f k is the interference line spectrum of the radiated noise line spectrum, otherwise the next judgment is carried out;

[0073] (3) If B l > 3 and (L fl - F fl ) × 32s > 1000s, then it is determined that the line spectrum frequency of f k is the interference line spectrum of the radiated noise line spectrum, otherwise it is determined that the radiated noise line spectrum is the suspected target line spectrum.

[0074] Step 6 is specifically as follows:

[0075] If it is selected to end the interference processing, the processing ends; otherwise, let r = r + 1 and return to Step 2.

[0076] First, according to Figure 2 as shown, when encountering interference frequencies with harmonic relationships with a confidence level of B max , and constant interference frequencies with a confidence level greater than 3, these interference frequencies will be blocked after detection, achieving a good anti-interference effect. And the anti-interference effect will gradually stabilize over time. In this example, only a time interval after 1000s (the 32nd frame) of anti-interference processing is selected as an example.

[0077] Through Figure 4 is the line spectrum detection result after anti-interference processing, Figure 6 is the line spectrum detection result without anti-interference processing. The solid line represents the radiated noise line spectrum emitted during the lake experiment, and the hollow origin represents the detected line spectrum frequency points. It can be seen that without anti-interference processing ( Figure 6 ), there are many harmonic interferences and constant interferences on the time-frequency diagram. After using the anti-interference method of the present invention ( Figure 4 ), the interference frequencies with harmonic relationships and constant interference frequencies are blocked. The detection accuracy rate and false alarm rate of the radiated noise line spectrum after anti-interference processing are as Figure 5 shown, and the detection accuracy rate and false alarm rate of the radiated noise line spectrum without anti-interference processing are as Figure 7 shown. On the premise that the line spectrum detection accuracy rate remains unchanged, the false alarm probability drops significantly. Through calculation, in the time region where the target appears, the average false alarm probability per frame drops from the original 82.33% to 58.15%, and the average number of false alarms per frame drops from the original 25.92 to 7.27.

[0078] The above embodiments show that the detection points calculated and screened by using the anti-interference method of the present invention can well remove the interference with harmonic relationship and the continuous interference in the environment, greatly reduce the false alarm probability, and achieve a good anti-interference effect.

[0079] The above is only a preferred embodiment of the present invention, and it is not a limitation of the present invention in any other form. Any modification or equivalent change made according to the technical essence of the present invention still belongs to the scope of protection required by the present invention.

Claims

1. An anti-interference method in the detection of radiated noise line spectrum, characterized in that, The method includes the following steps: Step 1, initialize data processing parameters, a radiation noise line spectrum detection information matrix, the maximum confidence level, and a determination threshold; Step 2, read a frame of radiation noise line spectrum detection information, including line spectrum frequency, detection times, first occurrence frame number, and last occurrence frame number, update it to the line spectrum information matrix, and determine the interference confidence level value of each line spectrum according to the detection times of the line spectrum and its proportion; Step 3, traverse each line spectrum in the current frame of radiation noise line spectrum detection information. If the number of line spectra in the radiation noise line spectrum information matrix that form a multiple frequency relationship with the line spectra in the current frame of radiation noise line spectrum detection information exceeds the threshold value, set the interference confidence level values of the corresponding line spectra and their multiple frequency line spectra in the radiation noise line spectrum information matrix to the maximum confidence level; otherwise, the interference confidence level value remains the initial value; Step 4, determine whether the number of processed frames reaches the processed frame number threshold. If it reaches, proceed to Step 5; otherwise, enter Step 6; Step 5, traverse each line spectrum in the current frame of radiation noise line spectrum detection information. According to the line spectrum interference confidence level corresponding to the line spectrum in the current frame of radiation noise line spectrum detection information in the radiation noise line spectrum information matrix, determine whether the line spectrum in the current frame of radiation noise line spectrum detection information is an interference line spectrum. If the interference confidence level is greater than the confidence level threshold and the duration is greater than the duration threshold, it is determined that the line spectrum is interference; otherwise, it is determined that the line spectrum is a suspected target line spectrum; Step 6, if it is selected to end the interference processing, end the processing; otherwise, after determination, return to Step 2 to continue reading the line spectrum detection information of the next frame.

2. The anti-interference method in the detection of radiated noise line spectrum according to claim 1, wherein The parameters initialized in Step 1 are as follows: Step 1.1, set the sampling rate f s , the frequency resolution Δf, and the detection duration T of each frame, the lower limit f l of the radiation noise line spectrum detection frequency and the upper limit f h ; Step 1.2, set up a matrix G for storing the line spectrum information of the radiation noise. G is a matrix with N rows and 5 columns, where where is the symbol for rounding up; the first column of G is the line spectrum frequency f n , the second column is the number of times the line spectrum is detected M n , the third column is the frame number F when the line spectrum first appears fn , the fourth column is the frame number L when the line spectrum last appears fn , and the fifth column is the interference confidence level B of the line spectrum n ; the initial value of each element in the first column of the line spectrum information matrix is f l +(n - 1)×Δf, n = 1, 2, 3, …, N, and the initial values of the other four columns are all assigned 0; Step 1.3, set the maximum value K of the number of line spectra detected in each frame m , process the frame number threshold M λ , the current frame number is r, and let r = 1; Step 1.4, set the detection times threshold γ for initializing the confidence of the radiated noise line spectrum, and the detection times ratio threshold α i , β i , i = 1, 2, 3, 4, 5; Step 1.5, set the confidence threshold B of the radiated noise line spectrum λ , the maximum confidence value B max , the threshold N for the number of octave relationships λ , the upper limit Q of the number of octaves; Step 1.6, set the radiation noise line spectrum duration threshold T1 and the duration threshold T2; Step 1.7, set the allowable error f when judging the equality or harmonic relationship of the radiation noise line spectra ε .

3. The anti-interference method in the detection of radiated noise line spectrum according to claim 1, characterized in that, The specific steps of Step 2 include the following steps: Step 3.1, read the detection information of the radiated noise line spectrum of the r-th frame, i.e., the line spectrum frequency f k , k = 1, 2, 3, …, K, where K represents the total number of line spectra detected in the current frame, and K ≤ K m ; Step 3.2, for the K line spectra detected in the r-th frame, if there exists a positive integer p, and 1 ≤ p ≤ N, such that the elements of the p-th row of matrix G satisfy |f k -f p |<f ε , then let M p = M p +1, L fp = r; at this time, if F fp = 0, then let F fp = r, otherwise F fp remains unchanged; if α i < M p / L fp < β i , and M p > γ, i = 1, 2, 3, 4, 5, then let B p = i.

4. A method for anti-interference in the detection of radiated noise line spectrum according to claim 1, characterized in that The specific steps of Step 3 include the following steps: Step 4.1, for each line spectrum f in matrix G n , use N n to represent the number of line spectra in matrix G that form a multiple-frequency relationship with f n ; if L fn >0, set N n = 0, variable j = 1, and initialize the Q elements of array H to 0; Step 4.2, if there exists a positive integer q, 2 ≤ q ≤ Q + 1, such that the elements in the j-th row of matrix G satisfy |q × f n - f j | < f ε and M j > 0, then let N n = N n + 1, and let the q-th element of H, H(q) = j; Step 4.3, if j < N, then let j = j + 1, and return to Step 4.2; Step 4.4, if N n ≥ N λ then set B n = B max , and set all B corresponding to H(q) ≠ 0 H(q) = B max .

5. A method for anti-interference in the detection of radiation noise line spectrum according to claim 1, characterized in that Step 4 specifically includes the following steps: Determine whether the number of detected frames reaches the threshold M λ , that is, if r ≥ M λ then proceed to Step 5, otherwise enter Step 6.

6. The anti-interference method in the detection of radiated noise line spectrum according to claim 1, characterized in that The specific steps of step 5 are as follows: For the K line spectra detected in the r-th frame, at the line spectrum frequencies f corresponding to the matrix G of radiation noise line spectrum information p , where p is the row number corresponding to this line spectrum in matrix G: Step 6.1, for the K line spectra f detected in the r-th frame k , k = 1, 2, 3, …, K, if there exists a positive integer l, and 1 ≤ l ≤ N, such that the elements of the l-th row of matrix G satisfy |f k - f l | < f ε , then proceed to Step 6.2; Step 6.2, if B l = B max and (L fl - F fl ) × T > T1, then it is determined that the line spectrum of the radiated noise with a line spectrum frequency of f k is an interference line spectrum; otherwise, the next judgment is carried out; Step 6.3, if B l > B λ and (L fl - F fl ) × T > T2, then it is determined that the line spectrum of the radiated noise with frequency f k is an interference line spectrum, otherwise it is determined that the radiated noise line spectrum is a suspected target line spectrum.

7. The anti-interference method in the detection of radiation noise line spectrum according to claim 1, characterized in that The specific steps of Step 6 include the following steps: If it is selected to end the interference processing, terminate the anti-interference processing, and the radiation noise line spectrum information matrix no longer changes; otherwise, let r = r + 1, return to Step 2, continue to update the radiation noise line spectrum information matrix and perform anti-interference processing.

Images