阅读说明：本技术 一种无线电信号频谱的平滑滤波方法及系统 (Smoothing filtering method and system for radio signal frequency spectrum ) 是由 张吉楠 向俊 王萌 孙恩元 于 2021-07-21 设计创作，主要内容包括：本发明公开了一种无线电信号频谱的平滑滤波方法及系统,方法包括以下步骤：获取无线电信号的频谱X(n),对频谱X(n)进行左右翻转得到翻转后的频谱X-(R)(n)；根据预设的平滑滤波算法对频谱X(n)进行平滑滤波得到第一频谱Y-(F)(n),且对翻转后的频谱X-(R)(n)进行平滑滤波得到第二频谱Y-(R0)(n),对第二频谱Y-(R0)(n)进行左右翻转得到第三频谱Y-(R)(n)；将第一频谱Y-(F)(n)与第三频谱Y-(R)(n)的数据求和取平均值得到第四频谱Y(n),并将第四频谱Y(n)作为平滑滤波后的频谱输出。本发明的方法调节简单且计算量小,可以抵消平滑滤波的延迟特性。(The invention discloses a smoothing filtering method and a system for a radio signal frequency spectrum, wherein the method comprises the following steps: obtaining a frequency spectrum X (n) of a radio signal, and turning the frequency spectrum X (n) left and right to obtain a turned frequency spectrum X R (n); according to a preset smoothing filtering algorithm, smoothing filtering is carried out on the frequency spectrum X (n) to obtain a first frequency spectrum Y F (n) and for the inverted spectrum X R (n) smoothing the filtered spectrum to obtain a second spectrum Y R0 (n) for the second spectrum Y R0 (n) left-right turning to obtain a third frequency spectrum Y R (n); the first frequency spectrum Y F (n) and a third spectrum Y R And (n) summing and averaging the data to obtain a fourth frequency spectrum Y (n), and outputting the fourth frequency spectrum Y (n) as the smooth filtered frequency spectrum. The method of the invention has simple adjustment and small calculation amount, and can counteract smooth filteringA delay characteristic.)

1. A method of smoothing the spectrum of a radio signal, comprising the steps of:

s1, obtaining the frequency spectrum X (n) of the radio signal, and turning the frequency spectrum X (n) to obtain the turned frequency spectrum X_R(n)；

S2) smoothing and filtering the frequency spectrum X (n) to obtain a first frequency spectrum Y_F(n) and for the flipped spectrum X_R(n) smoothing the filtered spectrum to obtain a second spectrum Y_R0(n) for said second spectrum Y_R0(n) turning over again to obtain a third spectrum Y_R(n)；

S3) dividing the first spectrum Y_F(n) and the third spectrum Y_R(n) combining to obtain a fourth spectrum y (n), and outputting the fourth spectrum y (n) as a smooth filtered spectrum.

2. The method for smoothing the spectrum of a radio signal according to claim 1, wherein in step S1) and step S2), the inverting specifically comprises: and reversing the data sequence of the frequency spectrum to be processed to obtain a reverse data sequence, and taking the reverse data sequence as the data sequence of the processed frequency spectrum.

3. The method for smoothing the spectrum of a radio signal according to claim 1, wherein the step S2) of smoothing filtering includes the steps of:

s201) calculating the calculation results of the first exponential smoothing filter of the current point in the frequency spectrum to be processed respectivelyQuadratic exponential smoothing filter calculation resultAnd cubic exponential smoothingResults of filtering calculations

S202) calculating a result according to the first exponential smoothing filterQuadratic exponential smoothing filter calculation resultAnd cubic exponential smoothing filter calculation resultsRespectively calculating intermediate calculation results A of current points in the frequency spectrum to be processed_i、B_i、C_i；

S203) calculating the result A according to the intermediate_i、B_i、C_iCalculating to obtain the value Y after filtering of the current point in the frequency spectrum to be processed_i。

4. The method for smoothing filter of spectrum of radio signal as claimed in claim 3, wherein in step S201), the result of one exponential smoothing filter calculation for the current point in the spectrum to be processedThe calculation result of the first exponential smoothing filter of the current point value and the previous point isWeighted sum of, quadratic exponential smoothing filter calculation of the current point in the spectrum to be processedComputing results for a first exponential smoothing filter of a current pointSecond order exponential smoothing filter calculation result with last pointThe result of the cubic exponential smoothing filter calculation of the current point in the spectrum to be processedSecond order exponential smoothing filter calculation result for current pointCubic exponential smoothing filter calculation result with last pointIs calculated as a weighted sum of.

5. The method of claim 4, wherein the result of the first exponential smoothing filter calculation at the current point in the spectrum to be processed is the result of the smoothing filter calculationQuadratic exponential smoothing filter calculation resultAnd cubic exponential smoothing filter calculation resultsThe function of (a) is expressed as follows:

in the above formula, α is a smoothing coefficient, i is the number of the current point, and x_iIs the value of the current point in the spectrum to be processed,the result of the first exponential smoothing filtering for a point in the spectrum to be processed,the result of the second exponential smoothing filtering for a point in the spectrum to be processed,and smoothing the filtering result for the third exponential of the upper point in the frequency spectrum to be processed.

6. The method according to claim 3, wherein the intermediate calculation result A of the current point in step S202) is_i、B_i、C_iThe function of (a) is expressed as follows:

in the above formula, α is a smoothing coefficient, i is the number of the current point,for a first exponential smoothing filtering result of the current point in the spectrum to be processed,is the result of the quadratic exponential smoothing filter of the current point in the spectrum to be processed,and smoothing the filtering result for the cubic exponent of the current point in the frequency spectrum to be processed.

7. The method for smoothing the spectrum of a radio signal according to claim 3, wherein in step S203), the current point in the spectrum to be processed is filtered to obtain the value Y_iThe function is expressed as follows:

Y_i＝A_i+B_i+C_i

in the above formula, A_i、B_i、C_iRespectively, the intermediate calculation results of the current point in the spectrum.

8. Method for smoothing the spectrum of a radio signal according to claim 1, characterized in that the first spectrum Y is filtered_F(n) and the third spectrum Y_R(n) combining specifically said first spectrum Y_F(n) and the third spectrum Y_RAnd (n) summing the corresponding data to obtain an average value.

9. The method of smoothing a spectrum of a radio signal according to claim 5 or 6, wherein the smoothing coefficient α is 0.01 to 0.05.

10. A system for smoothing the spectrum of a radio signal, comprising a computer device programmed or configured to perform the method for smoothing the spectrum of a radio signal according to any one of claims 1 to 9.

Technical Field

The present invention relates to the field of wireless communication technologies, and in particular, to a method and a system for smoothing a spectrum of a radio signal.

Background

With the gradual advance of the national information strategy, the wireless communication technology is rapidly developed, and has played an important basic supporting role in a plurality of fields of resident life, military industry and the like. In the fields of radio monitoring and electronic investigation, signal detection is one of the most basic and important functions, and is a prerequisite for tasks such as subsequent signal identification, signal main parameter estimation, signal extraction, signal characteristic statistics and the like. The deviation of the detection result will directly affect the accuracy of these important monitoring indexes and investigation parameters.

In the signal transmission process, the quality of the signal is often affected by the interference of noise, and in order to improve the situation, the signal is often required to be subjected to noise elimination processing. Fig. 1 shows the frequency spectrum of the real actual acquisition signal of a certain satellite, and it can be seen that the number of the satellite signals is large, the bandwidth range span is large, and the satellite signals are concentrated on the same frequency band. Meanwhile, under the condition of limited non-cooperative receiving conditions, the fluctuation of some signal envelopes is large, the background noise change is large, and false alarm or missing judgment of signal detection is easily caused. The frequency point signals indicated by the marks (i) and (ii) in fig. 1 are respectively composed of a plurality of signals, and can be easily considered as the same signal if observed from the full frequency band. Therefore, a certain method is needed for smoothing filtering, otherwise, the accuracy of signal identification and parameter detection is very low, and the subsequent radio signal parameter estimation and demodulation calculation processes are seriously affected.

For smoothing filtering of spectral data of radio signals, a method of smoothing an image in digital image processing is generally used for reference. Generally, the most widely used, most representative smoothing operation is mainly mean smoothing filtering and gaussian smoothing filtering in linear smoothing filtering, and median filtering in nonlinear smoothing filtering.

For mean smoothing filtering, the mean filter used to smooth the signal is often referred to as a moving average filter. It can be seen from the name of the filter that this is an operation of averaging the input signal to generate the input signal, and for the mean value smoothing filtering, the filter has a delay due to the correlation between the smoothed filter output value of the current point i and the M/2 points before and after the current point i. If the smoothing is calculated from the left side to the right side of the spectrogram, the smoothed spectrum has a phenomenon of "frequency shift left" due to the delay effect of the filter, which can be counteracted by compensation as shown by the dotted line in fig. 2, but the magnitude of the delay effect is related to the parameter M of the average smoothing filter, and the signal condition under the actual environment is complicated and variable, and it is difficult to counteract the delay effect by simple compensation. It is difficult to determine reasonable filter parameters and compensation methods in practical situations. More importantly, for the detection of radio signals, the carrier frequency of the signals is an important parameter, whether the carrier frequency estimation is accurate or not directly affects the demodulation process of the subsequent processing of the signals, and the delay effect seriously affects the carrier frequency estimation of the signals.

For gaussian smoothing filtering, also referred to as gaussian blurring, it refers to an operation of blurring a signal using a gaussian function. As shown in fig. 2, the frequency spectrum smoothed by the gaussian smoothing filter does not have the phenomenon of "frequency shift left", but has a low-pass filtering characteristic, so that the contour and edge information of the signal are damaged while high-frequency noise is filtered. In addition, in the Gaussian smoothing filter calculation process, convolution calculation content is adopted. This greatly increases the amount of computation relative to other smoothing filtering methods.

For median filtering, median refers to a critical point value that can segment the larger half and the smaller half of a value in a certain data set or probability distribution. For example, the median value of the data set {1,3,3,6,7,8,9} is 6. The concept of median is introduced into signal processing, and the filter is used as a nonlinear filter and is mainly used for eliminating isolated noise points such as impulse noise, salt and pepper noise and the like. As shown in fig. 2, the median filtering method protects the detail features, contours and edges of the signal, but for the U-shaped data or concave data portions, the median filtered boundary contour is blurred, which is not favorable for distinguishing the middle contour of two adjacent signals.

In the prior art, chinese patent CN201610913626 proposes that after a frequency spectrum after smooth filtering is left-right flipped, smooth filtering is performed again and left-right flipped, and a frequency spectrum of a signal is even-numbered smooth filtered and flipped to obtain a final signal after smooth filtering, thereby preventing the frequency spectrum from shifting due to a delay characteristic of the smooth filtering. In addition, a section of frequency spectrum is subjected to multiple mean filtering, so that data is more and more blurred, the outline and edge information of a signal are greatly damaged, and meanwhile, the calculated amount is also remarkably increased.

Disclosure of Invention

The technical problem to be solved by the invention is as follows: aiming at the technical problems in the prior art, the invention provides a smoothing filtering method and a smoothing filtering system for a radio signal frequency spectrum, which are simple to adjust and small in calculation amount and can offset the delay characteristic of smoothing filtering.

In order to solve the technical problems, the technical scheme provided by the invention is as follows:

a method of smoothing the spectrum of a radio signal, comprising the steps of:

s1, obtaining the frequency spectrum X (n) of the radio signal, and turning the frequency spectrum X (n) to obtain the turned frequency spectrum X_R(n)；

S2) smoothing and filtering the frequency spectrum X (n) to obtain a first frequency spectrum Y_F(n) and for the flipped spectrum X_R(n) smoothing the filtered spectrum to obtain a second spectrum Y_R0(n) for said second spectrum Y_R0(n) turning over again to obtain a third spectrum Y_R(n)；

S3) dividing the first spectrum Y_F(n) and the third spectrum Y_R(n) combining to obtain a fourth spectrum y (n), and outputting the fourth spectrum y (n) as a smooth filtered spectrum.

Further, in step S1) and step S2), the turning specifically includes: and reversing the data sequence of the frequency spectrum to be processed to obtain a reverse data sequence, and taking the reverse data sequence as the data sequence of the processed frequency spectrum.

Further, in step S2), the performing of the smoothing filtering includes the following steps:

s201) calculating the calculation results of the first exponential smoothing filter of the current point in the frequency spectrum to be processed respectivelyQuadratic exponential smoothing filter calculation resultAnd cubic exponential smoothing filter calculation results

S202) calculating a result according to the first exponential smoothing filterQuadratic exponential smoothing filter calculation resultAnd cubic exponential smoothing filter calculation resultsRespectively calculating intermediate calculation results A of current points in the frequency spectrum to be processed_i、B_i、C_i；

S203) calculating the result A according to the intermediate_i、B_i、C_iCalculating to obtain the value Y of the frequency spectrum to be processed after the current point is filtered_i。

Further, in step S201), a calculation result of a first exponential smoothing filter of a current point in the spectrum to be processedThe calculation result of the first exponential smoothing filter of the current point value and the previous point isWeighted sum of, quadratic exponential smoothing filter calculation of the current point in the spectrum to be processedComputing results for a first exponential smoothing filter of a current pointSecond order exponential smoothing filter calculation result with last pointThe result of the cubic exponential smoothing filter calculation of the current point in the spectrum to be processedSecond order exponential smoothing filter calculation result for current pointCubic exponential smoothing filter calculation result with last pointIs calculated as a weighted sum of.

Further, a calculation result of a primary exponential smoothing filter of a current point in a spectrum to be processedQuadratic exponential smoothing filter calculation resultAnd cubic exponential smoothing filter calculation resultsThe function of (a) is expressed as follows:

in the above formula, α is a smoothing coefficient, i is the number of the current point, and x_iIs the value of the current point in the spectrum to be processed,the result of the first exponential smoothing filtering for a point in the spectrum to be processed,the result of the second exponential smoothing filtering for a point in the spectrum to be processed,and smoothing the filtering result for the third exponential of the upper point in the frequency spectrum to be processed.

Further, the intermediate calculation result a of the current point in step S202)_i、B_i、C_iThe function of (a) is expressed as follows:

in the above formula, α is a smoothing coefficient, i is the number of the current point,for a first exponential smoothing filtering result of the current point in the spectrum to be processed,is the result of the quadratic exponential smoothing filter of the current point in the spectrum to be processed,and smoothing the filtering result for the cubic exponent of the current point in the frequency spectrum to be processed.

Further, in step S203), the current point filtered value Y of the spectrum to be processed_iThe function is expressed as follows:

Y_i＝A_i+B_i+C_i

in the above formula, A_i、B_i、C_iRespectively, the intermediate calculation results of the current point in the spectrum.

Further, the smoothing coefficient α is 0.01 to 0.05.

The invention also proposes a system for smoothing the spectrum of a radio signal, comprising a computer device programmed or configured to carry out said method for smoothing the spectrum of a radio signal.

Compared with the prior art, the invention has the advantages that:

the invention turns the frequency spectrum to obtain the turned frequency spectrum, the data sequences in the frequency spectrum and the turned frequency spectrum are opposite to each other, the forward data and the reverse data of the frequency spectrum are respectively filtered smoothly to obtain the first frequency spectrum and the second frequency spectrum, the second frequency spectrum is turned again to obtain the third frequency spectrum, at the moment, the moving positions of the first frequency spectrum and the third frequency spectrum on the frequency spectrum graph relative to the original frequency spectrum are opposite and the distances of the first frequency spectrum and the third frequency spectrum are the same, the data in the first frequency spectrum and the third frequency spectrum correspond to each other one by one, then, the data in the first frequency spectrum and the third frequency spectrum are summed and averaged to obtain the fourth frequency spectrum, and the fourth frequency spectrum just offsets the delay influence of the smooth filtering.

Drawings

Fig. 1 is a full bandwidth spectrum diagram of a real signal from a certain satellite.

Fig. 2 is a comparison graph of the effect of the conventional smoothing filtering method.

Fig. 3 is a flowchart illustrating a smoothing filtering method according to an embodiment of the present invention.

Fig. 4 is a diagram illustrating an effect of applying the smoothing filtering method according to the embodiment of the present invention to a satellite real signal.

Fig. 5 is a diagram illustrating the effect of applying the mean smoothing filtering method to the actual satellite signal.

Fig. 6 is a diagram illustrating the effect of applying the gaussian smoothing filtering method to the actual satellite signal.

Fig. 7 is a diagram illustrating the effect of applying the median filtering method to the satellite real signal.

FIG. 8 is a graph illustrating the effect of the smoothing filtering method according to the embodiment of the present invention compared with the effect of the conventional smoothing filtering method.

Detailed Description

The invention is further described below with reference to the drawings and specific preferred embodiments of the description, without thereby limiting the scope of protection of the invention.

The idea is that the predicted value is the weighted sum of the previous observed values, different weights are given to different data, the new data is given a larger weight, and the old data is given a smaller weight. Compared with the mean value smoothing filtering, the method has smaller influence on the original signal and better filtering effect, and although the method also has a delay characteristic, the time sequence data does not generally consider the influence of slight delay.

In order to achieve better effect based on the exponential smoothing filtering and offset the effect caused by the delay of the exponential smoothing filtering, the invention provides a smoothing filtering method for the spectrum of a radio signal, as shown in fig. 3, comprising the following steps:

s1, obtaining the frequency spectrum X (n) of the radio signal, and turning the frequency spectrum X (n) left and right to obtain the turned frequency spectrum X_R(n) so that the inverted spectrum X_RThe data sequence in (n) and the data sequence in the spectrum x (n) are opposite to each other, for example, the original data sequence in the spectrum is {1,2,3,4,5,6,7}, and the reversed data sequence is {7,6,5,4,3,2,1}, that is, the reversed data sequence is used as the data sequence of the reversed spectrum, so that the delay effect caused by the smoothing filtering is mainly counteracted, and the equalization compensation can be performed for counteracting the delay effect, but the design and the calculation process are very complex, and the data sequence is very sensitive to the change of the smoothing coefficient;

s2) smoothing and filtering the frequency spectrum X (n) to obtain a first frequency spectrum Y_F(n) and for the inverted spectrum X_R(n) smoothing the filtered spectrum to obtain a second spectrum Y_R0(n) for the second spectrum Y_R0(n) left-right turning to obtain a third frequency spectrum Y_R(n) at this time, the first spectrum Y_F(n) and a third spectrum Y_R(n) the moving directions are opposite to each other on the spectrogram relative to the spectrum X (n), and the moving distances are the same;

s3) the first spectrum Y_F(n) and a third spectrum Y_R(n) summing and averaging the data to obtain a fourth frequency spectrum Y (n), and outputting the fourth frequency spectrum Y (n) as the frequency spectrum after smooth filtering, wherein the frequency spectrum X after overturning_RThe data sequence in (n) and the data sequence in the spectrum X (n) are in inverse one-to-one correspondence, so that the first spectrum Y is obtained by smooth filtering_F(n) and a second frequency spectrum Y_R0The data sequences in (n) should also be in inverse one-to-one correspondence with each other, and the first spectrum Y_F(n) and a third frequency spectrum Y_RAnd (n) the data sequences are in one-to-one correspondence, and the fourth frequency spectrum obtained by summing and averaging the corresponding data in the first frequency spectrum and the third frequency spectrum just offsets the delay influence of the smoothing filtering.

In step S2) of this embodiment, an exponential smoothing filtering method is used for smoothing, and the exponential smoothing filtering method is classified into a first-order exponential smoothing method, a second-order exponential smoothing method, a third-order exponential smoothing method, and the like according to different smoothing times, and we find that the third-order exponential smoothing method can achieve an optimal effect on smoothing filtering of a radio signal spectrum through an experiment of three kinds of exponential smoothing and through effect analysis after the experiment, where the third-order exponential smoothing method in this embodiment includes the following steps:

s201) calculating the calculation results of the first exponential smoothing filter of the current point in the frequency spectrum to be processed respectivelyQuadratic exponential smoothing filter calculation resultAnd cubic exponential smoothing filter calculation resultsFirst exponential smoothing filter calculation result of current point in frequency spectrum to be processedThe calculation result of the first exponential smoothing filter of the current point value and the previous point isWeighted sum of, quadratic exponential smoothing filter calculation of the current point in the spectrum to be processedComputing results for a first exponential smoothing filter of a current pointSecond order exponential smoothing filter calculation result with last pointThe result of the cubic exponential smoothing filter calculation of the current point in the spectrum to be processedSecond order exponential smoothing filter calculation result for current pointCubic exponential smoothing filter calculation result with last pointWeighted sum of, a one-time exponential smoothing filter calculation of the current point in the spectrum to be processedQuadratic exponential smoothing filter calculation resultAnd cubic exponential smoothing filter calculation resultsThe function of (a) is expressed as follows:

in the above formula, α is a smoothing coefficient, i is the number of the current point, and x_iIs the value of the current point in the spectrum to be processed,the result of the first exponential smoothing filtering for a point in the spectrum to be processed,the result of the second exponential smoothing filtering for a point in the spectrum to be processed,smoothing a filtering result for a cubic exponential of a previous point in a frequency spectrum to be processed;

s202) calculating a result according to the first exponential smoothing filterQuadratic exponential smoothing filter calculation resultAnd cubic exponential smoothing filter calculation resultsRespectively calculating intermediate calculation results A of current points in the frequency spectrum to be processed_i、B_i、C_i，A_i、B_i、C_iThe function of (a) is expressed as follows:

in the above formula, α is a smoothing coefficient, i is the number of the current point,for a first exponential smoothing filtering result of the current point in the spectrum to be processed,is the result of the quadratic exponential smoothing filter of the current point in the spectrum to be processed,the three-time exponential smoothing filtering result of the current point in the frequency spectrum to be processed is obtained;

s203) calculating the result A according to the intermediate_i、B_i、C_iAnd calculating to obtain a value after filtering of the current point in the frequency spectrum to be processed, wherein the function expression is as follows:

Y_i＝A_i+B_i+C_i (3)

in the above formula, A_i、B_i、C_iRespectively, the intermediate calculation results of the current point in the spectrum to be processed.

Each point on the frequency spectrum X (n) can be calculated to obtain a filtered value through processing of a cubic exponential smoothing method, and a data sequence formed by the filtered values is a first frequency spectrum Y on a frequency spectrogram_F(n) corresponding, inverted spectrum X_R(n) the second spectrum Y is also obtained by processing with cubic exponential smoothing_R0(n), first frequency spectrum Y_F(n) and a second frequency spectrum Y_R0(n) with respect to frequency spectrum X (n) and inverted frequency spectrum X, respectively_R(n) there is a slight delay due to the third spectrum Y_R(n) is the second frequency spectrum Y_R0(n) is left-right flipped so that there is a slight advance with respect to the spectrum X (n), on the spectrogram, the third spectrum Y_R(n) distance to right shift and first spectrum Y_F(n) the distance to the left is the same, thus shifting the first spectrum Y_F(n) and a third spectrum Y_RAnd (n) summing and averaging the data to obtain a fourth frequency spectrum Y (n), namely the result after smooth filtering of the frequency spectrum X (n) after the delay is cancelled.

In this embodiment, a section of actual data of a certain satellite is used as original data, 262144-point fast fourier transform is performed on the data to obtain a frequency spectrum x (n), and the effect of performing smoothing filtering on the frequency spectrum x (n) by using the method in this embodiment is shown in fig. 4, where in this embodiment, a smoothing coefficient α is 0.01 to 0.05, and a specific value is 0.018. In order to offset the delay characteristic of the exponential smoothing filter, the filter coefficients adopted by the forward and backward data of the frequency spectrum x (n) are the same. As can be seen from fig. 4, the smoothing filtering performed on the frequency spectrum x (n) by the method in this embodiment can remove the high-frequency noise to the maximum extent, and also ideally protect the detailed features, contours, and edges of the signal. Meanwhile, the method is also ideal for the processing result of U-shaped data between two adjacent signals, and can clearly distinguish the boundary between the signals, which is very favorable for the extraction and the resolution of the subsequent signals and the parameter estimation.

For comparison with the effect of the method in this embodiment, the conventional smoothing filtering method is used to perform smoothing filtering on the spectrum x (n), and the following contents are provided:

adopting a mean smoothing filtering method, wherein the function expression of the mean smoothing filtering method is as follows:

in the above formula, y (i) is the filtering output value of the current point i of the spectrum x (n), M is the average point number calculated each time as the filtering coefficient, and x (n) is the value from the previous M/2 to the next M/2 of the current point i. In this embodiment, M is 101 points, and the filtering effect is as shown in fig. 5, it can be found that the spectrum smoothing curve has an obvious offset phenomenon, and the offset phenomenon of reducing M points is weakened, but a lot of burrs are added at the same time, and the smoothing effect is very undesirable. For radio signals, the central frequency point or carrier frequency of the signal is a key parameter, and if the data after smooth filtering has deviation, the subsequent carrier frequency estimation has deviation, which affects the overall signal monitoring and detection results.

Adopting a Gaussian smoothing filtering method, wherein the function expression of the Gaussian smoothing filtering method is as follows:

in the above equation, y (i) is a filter output value of the current point i of the spectrum x (n), and σ and u represent a standard deviation and a mean of a gaussian function, respectively, as filter coefficients. In this embodiment, the number of data points involved in each calculation is 1025, x (τ) is the value of the τ -th point of the spectrum x (n), and σ can be 15 by referring to the convolution calculation formula. The filtering effect is shown in fig. 6, and it can be found that the high-frequency noise is filtered, and meanwhile, the contour and edge information of the signal are damaged. The overall filtering effect tends to be fuzzy, and the damage to the signal frequency spectrum details and the contour is large. Adjusting the number of data points and the value of sigma involved in each calculation can obtain different calculation amounts and smoothing effects, but tends to blur details and outlines of the data in general.

Adopting a median filtering method, and the function expression of the median filtering method is as follows:

in the above formula, y (i) is the filtering output value of the current point i of the spectrum x (n), M is the number of points calculated each time, and is used as the filtering coefficient, which is generally odd,represents the median of the values corresponding to points i-M/2 to i + M/2 in spectrum x (n), such as med (3,7,6,1,3,2) ═ 3. In this embodiment, M is 101 points, and the filtering result is shown in fig. 7. It has been found that this method preserves the detail features, contours and edges of the signal. However, between two adjacent signals, the smooth results show that the spectrum baseline is lifted by processing the U-shaped data or the concave data part. This would be detrimental to distinguishing the mid-profile of the two signals and to calculating the bandwidth of the two signals.

In order to more intuitively compare the effect of the method of the embodiment with the effect of the existing smoothing filtering method, as shown in fig. 8, after summarizing and comparing result curves of smoothing filtering of several methods, it is found by marking the lowest point of U-shaped data between two signals that the average smoothing filtering reaches-88.48 dB, the gaussian smoothing filtering reaches-89.51 dB, the median filtering reaches-87.88 dB, and the smoothing filtering method of the embodiment reaches-91.25 dB. 1.74dB is optimized for gaussian smoothing filtering.

The calculation time for each filtering method is statistically calculated, and it is found that the mean smoothing filter time is 0.106125 seconds, the gaussian smoothing filter time is 0.521634 seconds, the median filtering filter time is 0.385380 seconds, and the filtering time of the method of the present embodiment is 0.072542 seconds. The calculation process of the gaussian smooth filtering adopts convolution calculation, so that the calculation amount is large, the calculation amount is also a certain amount when the average value is obtained by the secondary summation of the average filtering, and the calculation resources are consumed when the median value is obtained by the sorting every time the median value is obtained by the median filtering. According to the scheme, the calculation amount is fixed every time and does not change along with the change of the filter coefficient alpha, only the data corresponding to the current point and the previous point in the frequency spectrum are calculated by adopting the iterative calculation process, and the whole calculation amount is greatly reduced.

In summary, the method provided by the embodiment is applied to the smoothing filtering of the radio signal spectrum, the performance is generally superior to that of the existing smoothing filtering method, and the method has the characteristics of simple adjustment of the filter coefficient and small calculation amount. Meanwhile, the triple exponential smoothing filtering is used for acting on the frequency spectrum data in the positive direction and the reverse direction, so that the delay characteristic of the exponential filtering is skillfully counteracted. The test verifies that the high-frequency noise can be removed to the maximum extent, and the detail characteristics, the outline and the edge of the signal are protected very ideally, so that the early data preparation is prepared for subsequent signal identification and main parameter estimation and calculation.

The invention also proposes a system for smoothing the spectrum of a radio signal, comprising a computer device programmed or configured to carry out the method for smoothing the spectrum of a radio signal of the present embodiment.

The present invention also proposes a computer-readable storage medium storing a computer program programmed or configured to execute the smoothing filtering method of a radio signal spectrum of the present embodiment.

The invention is also suitable for smoothing filtering processing of image signals, voice signals and other digital signals.

The foregoing is considered as illustrative of the preferred embodiments of the invention and is not to be construed as limiting the invention in any way. Although the present invention has been described with reference to the preferred embodiments, it is not intended to be limited thereto. Therefore, any simple modification, equivalent change and modification made to the above embodiments according to the technical spirit of the present invention should fall within the protection scope of the technical scheme of the present invention, unless the technical spirit of the present invention departs from the content of the technical scheme of the present invention.