Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Projected Increase in Heatwaves under 1.5 and 2.0 °C Warming Levels Will Increase the Socio-Economic Exposure across China by the Late 21st Century
Next Article in Special Issue
Evaluation of GNSS-TEC Data-Driven IRI-2016 Model for Electron Density
Previous Article in Journal
Application of Shannon Entropy in Assessing Changes in Precipitation Conditions and Temperature Based on Long-Term Sequences Using the Bootstrap Method
Previous Article in Special Issue
Investigation and Validation of Short-Wave Scattering in the Anisotropic Ionosphere under a Geomagnetic Field
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Polynomial Fitting-Based Noise Reduction for Correlation Functions in Medium-Frequency Radar

1
National Key Laboratory of Electromagnetic Environment, China Research Institute of Radiowave Propagation, Qingdao 266107, China
2
Kunming Electro-Magnetic Environment Observation and Research Station, Qujing 655500, China
3
Qujing Electro-Magnetic Environment Observation and Research Station, Qujing 655500, China
4
College of Electronics and Information, Guangxi Minzu University, Nanning 530006, China
5
College of Information and Commnunication, National University of Defense Technology, Changsha 410073, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Atmosphere 2024, 15(8), 899; https://doi.org/10.3390/atmos15080899
Submission received: 9 June 2024 / Revised: 18 July 2024 / Accepted: 24 July 2024 / Published: 27 July 2024

Abstract

:
In the theoretical calculation of atmospheric wind fields using the cross-correlation analysis method of Medium-Frequency radar, it is necessary to compute a series of correlation parameters from the received echo signals, such as autocorrelation and cross-correlation functions, within the main lobe range of the antenna array to retrieve atmospheric parameters. However, both theoretical analysis and practical applications have shown that the shape of correlation functions can be affected by atmospheric conditions and receiver noise, leading to significant biases in the estimated correlation parameters within the main lobe range. In this study, we theoretically analyze the influence of noise on the amplitude of autocorrelation and cross-correlation functions. We propose a noise reduction method based on the characteristics of correlation functions at the zero-delay point to calculate the noise factor and process the correlation functions within the main lobe range. Furthermore, we conduct simulation analysis to evaluate the performance of this noise reduction method and summarize the effects of the number of fitting points and fitting methods on the noise reduction performance.

1. Introduction

The mesosphere and lower thermosphere (MLT) region between 60 and 110 km forms the boundary between the lower atmosphere and space. This region is dominated by atmospheric dynamics, including planetary waves, tides and gravity waves [1,2,3,4]. Ground-based radio radar like Medium-Frequency (MF) radar and VHF meteor radar systems have become the most widely used instrument to routinely observe MLT winds among ground-based techniques because they have the advantages of being low cost, easy to install, and that they operate automatically and continuously under all weather conditions [5,6]. One of the primary advantages of radar systems, in comparison to other remote sensing techniques for the MLT, is their ability to provide continuous measurements regardless of weather conditions [7,8,9]. A crucial aspect of measuring winds is the reliability of each technique. This requires a thorough understanding of the underlying scattering processes, potential instrumental effects, analysis-related simplifications, and assumptions that could introduce biases or systematic errors in the derived winds [10].
MF Radar achieves continuous 24 h observations of atmospheric wind fields at altitudes of 60 to 100 km [11,12,13,14]. Due to its advantages of being unattended and having low observation costs, MF Radar has become a primary tool for observing the wind fields in the mesosphere [15]. Additionally, MF Radar observations of atmospheric wind fields furnish extensive statistical data for atmospheric dynamics and climate research. Studies based on MF Radar wind field observations have significantly deepened our understanding of mesosphere environmental characteristics. The technique for wind field retrieval typically employs the Full Correlation Analysis (FCA) method based on the spaced antenna (SA) theory, which is widely used domestically and internationally due to its high accuracy and stability. Research consistently shows that the FCA method is an effective means of retrieving wind fields in the middle and upper atmosphere using MF Radar [16,17,18,19]. P. T. May used the FCA method to determine the maximum value in the correlation function and the error at the correlation lag value [20], thereby estimating wind fields. Guifu Zhang et al. employed various methods to estimate wind fields [21,22], thereby confirming the efficiency of the FCA method.
However, when using the FCA method for retrieval, autocorrelation and cross-correlation functions are first calculated from the time series of echoes obtained by spaced antennas to determine correlation parameters, followed by computing the vector wind field. In practical applications, noise generated by atmospheric conditions, or the radar itself, can distort the waveform of the autocorrelation and cross-correlation functions, causing significant errors in calculating correlation parameters. These errors can propagate through the iterative FCA retrieval process, severely affecting the accuracy of the calculations. Therefore, special analysis and noise reduction treatment in the correlation function are crucial for the accurate retrieval of atmospheric wind fields using the FCA method with MF Radar. Richard J. Doviak et al. employed the least squares method to estimate parameters of the correlation function [21,22]. David A. Holdsworth et al. provided fundamental principles for addressing noise reduction in correlation functions [14]. This paper examines the performance of this method through numerical simulations, providing important references for practical applications.
This study primarily addresses the issue of amplitude correction by calculating noise factors for the correlation function. The significance of the correlation function in FCA is introduced in Section 2. The theoretical noise issues encountered during data analysis are also discussed in Section 2. Detailed data analysis methods are presented in Section 3. The performance of the PFNR algorithm is evaluated in Section 4. Finally, conclusions are drawn in Section 5.

2. The Important Role of Correlation Functions in the FCA Method and Its Noise Theory

The computation of correlation functions is a crucial step in the FCA method. The principle of the FCA method can be summarized as follows:
Assuming the spatiotemporal correlation function corresponding to the ground diffraction pattern generated by radar echoes is a family of ellipsoids [23]:
ρ ( ξ , η , τ ) = ρ ( A ( ξ V x τ ) 2 + B ( η V y τ ) 2 + K τ 2 + 3 H ( ξ V x τ ) ( η V y τ ) )
where ξ and η , respectively, represent the distances along the latitude and longitude directions in the two-dimensional plane reference system, τ stands for time, V x and V y denote the atmospheric component of the wind speed along the latitude and longitude directions and A , B , K , H represent coefficients. Once these coefficients are determined, the ground diffraction pattern can be described using the parameters of this ellipsoid [24,25,26,27]. Expanding Equation (1):
ρ ( ξ , η , τ ) = ρ ( α ξ 2 + b η 2 + 2 f ξ τ + 2 g η τ + 2 h ξ η )
By comparison, the relationship between coefficients α , b , f , g , h and ( V x , V y ) can be derived as follows α V x 2 + b V y 2 + 2 h V x V y + k = 1 , α V x h V y = f , h V x b V y = g .
Suppose there is a set of antennas spaced apart by ( ξ 0 , η 0 ) . According to Equation (1), the corresponding autocorrelation function is denoted by ρ ( 0 , 0 , τ ) , and the cross-correlation function is denoted by ρ ( ξ 0 , η 0 , τ ) , where ρ ( 0 , 0 , τ ) and ρ ( ξ 0 , η 0 , τ ) are collectively referred to as the correlation functions of this antenna array. Using these correlation functions, a set of correlation parameters, τ , τ 0.5 , τ , and τ , can be determined. Here, τ represents the time delay corresponding to when the autocorrelation function decreases to ρ ( 0 , 0 , τ ) and the cross-correlation function has a value of ρ ( ξ 0 , η 0 , τ ) at zero delay; τ 0.5 represents the time delay corresponding to when the autocorrelation function ρ ( 0 , 0 , τ ) decreases to half of its peak value; τ represents the time delay corresponding to the peak of the cross-correlation function ρ ( ξ 0 , η 0 , τ ) ; τ represents the time delay corresponding to when the autocorrelation function ρ ( 0 , 0 , τ ) increases to the same value as the cross-correlation function ρ ( ξ 0 , η 0 , τ ) at its peak. Specifically, as shown in Figure 1, using these correlation functions, the coefficients in Equation (2) can be determined.
Therefore, the key to the FCA method lies in determining the correlation functions. The calculation methods for the autocorrelation function and the cross-correlation function are as follows: assuming the received echo signal at antenna α = ( x i , y i ) is f α ( t ) , the autocorrelation function f α ( t ) is defined as [14]:
ρ α ( τ k ) = [ f α ( t ) f α ( t + τ k ) ] [ | f α ( t ) | 2 ]
where the symbol [ · ] denotes statistical averaging. By expanding the continuous signal into a time series and performing statistical averaging, the autocorrelation function can be expressed as:
ρ α ( τ k ) = n = 1 N r | k | f α ( t ) f α ( t + τ k ) ( N r | k | ) n = 1 N r | f α ( t ) | 2
Similarly, assuming the received echo signals at antenna α = ( x i , y i ) are f α ( t ) , and the signals at antenna b = ( x l , y m ) are f b ( t ) , the cross-correlation function between f α ( t ) and f b ( t ) is defined as:
ρ α b ( τ k ) = [ f α ( t ) f b ( t + τ k ) ] [ | f α ( t ) | 2 ] [ | f b ( t ) | 2 ]
The cross-correlation function, after statistical averaging and expansion, is expressed as:
ρ α b ( τ k ) = n = 1 N r | k | f α ( t ) f b ( t + τ k ) ( N r | k | ) n = 1 N r | f α ( t ) | 2 n = 1 N r | f b ( t ) | 2
The accuracy of computing the correlation function directly determines the precision of the final calculated values. The influence of noise on the waveform of the correlation function can cause deviations in the computed correlation function from the true values. Below, we illustrate this with a simple numerical simulation.
First, two sinusoidal sequences, s ( t ) = sin ( t ) , s ( t ) = sin ( t + 5 ) , ( t = 0 , 0.1 , , 20 ) , each with 200 points, are generated. Then, additive Gaussian white noise with a signal-to-noise ratio (SNR) of 0 dB, 2 dB, and 4 dB is added to s ( t ) and s ( t ) to generate noisy sinusoidal sequences, s n ( t ) and s n ( t ) . Using Equation (4), the amplitudes of the autocorrelation functions of s n ( t ) at different SNRs are calculated, as shown in the left graph of Figure 2. Similarly, using Equation (5), the amplitudes of the cross-correlation functions between s n ( t ) and s ( t ) at different SNRs are calculated, as shown in the right graph of Figure 2.
From Figure 1 and Figure 2, it can be observed that noise significantly impacts the amplitudes of both types of correlation functions, the autocorrelation function and the cross-correlation function, within the main lobe range. The amplitudes decrease and exhibit distortion as the SNR decreases. This distortion in the correlation functions can introduce errors in the FCA retrieval process. David’s theory addresses this issue by correcting the amplitudes of the correlation functions to determine the attenuation coefficient, referred to as the noise factor. This noise factor aims to minimize the deviation of the autocorrelation function values under noisy conditions from those under noise-free conditions [24].
Assuming antenna α i receives an echo signal composed of both the signal s i ( t ) and noise n i ( t ) , the received signal can be represented as f i ( t ) :
f i ( t ) = s i ( t ) + n i ( t )
Since the signal and noise are uncorrelated, their statistical mean is zero [14]. Where [ s i ( t ) n j ( τ k ) ] = 0   i , j , k .
Based on specific statistical characteristics, noise can be further categorized into two types, correlated noise and uncorrelated noise, as follows:
(1) Noise that has certain correlations between different antennas due to coupling after reception n c i ( t ) , where [ n c i n c i ] 0   i j .
This type of noise includes atmospheric background noise, power supply noise, and artificial noise, among others.
(2) This type of noise, which includes receiver noise and galactic background noise, lacks correlation noise n u i ( t ) between different antennas after reception. The subdivided noise can be expressed as n i ( t ) = n c i ( t ) + n u i ( t ) .
According to the definitions of n c i ( t ) and n u i ( t ) , and utilizing their respective statistical properties, the following statistical laws can be derived, where i = j represents the same antenna:
[ n i ( t ) n j ( t + τ k ) ] = 0   τ k 0   or   τ k = 0   and   i j [ n c i ( t ) n c i ( t + τ ) ] = 0   τ k 0 [ n c i ( t ) n u i ( t + τ k ) ] = 0   i , j , τ k [ n c i ( t ) n u i ( t + τ k ) ] = 0   τ k 0   or   τ k = 0   and   i j
Utilizing these statistical laws of noise, we can further derive the expanded forms of the autocorrelation and cross-correlation functions.
Since any signal is correlated with itself, expanding f i ( t ) in Equation (7) as a combination of signal and noise, the autocorrelation function can be rewritten as Equation (9) (where τ k is abbreviated for noise τ ).
ρ ( τ ) = [ ( s i ( t ) + n i ( t ) ) ( s i ( t + τ ) + n i ( t + τ ) ) ] [ | s i ( t ) + n i ( t ) | 2 ]
Expanding the autocorrelation function yields:
ρ ( τ ) = 1 , τ = 0 [ ( s i ( t ) + n i ( t ) ) ( s i ( t + τ ) + n i ( t + τ ) ) ] [ | s i ( t ) + n i ( t ) | 2 ] , τ 0
Comparing Equations (4) and (17), when τ 0 .
ρ ( τ ) = [ | s i ( t ) | 2 ] [ | s i ( t ) | 2 ] + [ | n i ( t ) | 2 ] [ s i ( t ) s i ( t + τ ) ] [ | s i ( t ) | 2 ] = [ | s i ( t ) | 2 ] [ | s i ( t ) | 2 ] + [ | n i ( t ) | 2 ] ρ i ( τ k )
We define:
F i = [ | s i ( t ) | 2 ] + [ | n i ( t ) | 2 ] [ | s i ( t ) | 2 ]
Here, F i is called the noise factor of s i ( t ) . Multiplying all non-zero values of the expression for the autocorrelation function under noise in Equation (11) by F i yields the estimate of the autocorrelation function without noise:
ρ ( τ ) = 1 , τ = 0 ρ i ( τ k ) F i , τ 0
The value of the autocorrelation function at zero delay must be 1. As τ approaches 0, the autocorrelation function under noise tends towards 1 / F i , so F i can be determined from this.
Similarly, the cross-correlation function f i ( t ) and f j ( t ) in Equation (6), when expanded as a combination of signal and noise, can be rewritten as:
ρ i j ( τ ) = [ ( s i ( t ) + n i ( t ) ) ( s j ( t + τ ) + n j ( t + τ ) ) ] [ | s i ( t ) + n i ( t ) | 2 ] [ | s j ( t ) + n j ( t ) | 2 ]
Finally, the cross-correlation function can be expressed as:
ρ i j ( τ ) = [ s j ( t ) s j ( t ) ] + [ n i ( t ) n j ( t ) ] [ | s i ( t ) | 2 ] + [ | n i ( t ) | 2 ] [ | s j ( t ) | 2 ] + [ | n j ( t ) | 2 ] , τ = 0 s i ( t ) s j ( t + τ ) [ | s i ( t ) | 2 ] + [ | n i ( t ) | 2 ] [ | s j ( t ) | 2 ] + [ | n j ( t ) | 2 ] , τ 0
Analyzing Equation (15), when τ 0 .
ρ i j ( τ ) = s i ( t ) s j ( t + τ ) [ | s i ( t ) | 2 ] + [ | n i ( t ) | 2 ] [ | s j ( t ) | 2 ] + [ | n j ( t ) | 2 ] = [ | s i ( t ) | 2 ] [ | s i ( t ) | 2 ] + [ | n i ( t ) | 2 ] [ | s j ( t ) | 2 ] [ | s j ( t ) | 2 ] + [ | n j ( t ) | 2 ] s i ( t ) s j ( t + τ ) [ | s i ( t ) | 2 ] [ | s j ( t ) | 2 ]
According to the definition of the noise factor in Equation (13), Equation (15) can be rewritten as:
ρ i j ( τ ) = s i ( t ) s j ( t + τ ) F i F j [ | s i ( t ) | 2 ] [ | s j ( t ) | 2 ] = ρ a b ( τ ) F i F j
If F i F j is known, radar parameters ρ i j ( 0 ) can be calculated at τ = 0 using Equation (17). Here, is the noise factor s i ( t ) for antenna F i , and is the noise factor s j ( t ) for antenna F j , with each noise factor originating from the noise of the signals received by the respective antennas.
The key to estimating the autocorrelation and cross-correlation functions lies in computing the noise factor. The definitions and computational formulas for the correlation functions are presented in Table 1.

3. Correction Steps for the Amplitude of Correlation Functions

The autocorrelation function curve can be fitted using an even-degree polynomial to determine the noise factor. Assuming the signal to be tested is a standard sine wave y ( x ) = α sin ( ω 0 x + ϕ ) , for simplicity, the sampling frequency of x is denoted as f , ω = ω 0 f , and y ( m ) = α sin ( ω m + ϕ ) , where m is an integer. If the measurement time T is long enough, the expression for its autocorrelation function is:
R x x ( k ) = 1 T 0 T α sin ( ω m + ϕ ) α sin [ ω ( m + k ) + ϕ ] d x = α 2 2 cos ( ω k )
To avoid interference from the coefficient values α , the autocorrelation value at τ = 0 can be normalized, yielding ρ ( k ) = cos ( ω k ) .
Similarly, if the signal under test is the superposition of n sine waves (this is a commonly used assumption for narrow band signals as per [28]), y ( x ) = i = 1 n α i sin ( ω o i x + ϕ i ) , where f is the sampling frequency of x, ω i = ω o i f , so that y ( m ) = i = 1 n α i sin ( ω i m + ϕ i ) , m is an integer. If the measurement time T is long enough, its autocorrelation function expression is as follows:
R x x ( k ) = 1 T 0 T y ( m ) y ( m + k ) d x = i = 1 n α i 2 2 cos ( ω i k )
After normalization, the expression becomes as follows:
ρ ( k ) = i = 1 n α i 2 2 cos ( ω i k ) i = 1 n α i 2 2
From the results, it is observed that:
(1) When n is a finite integer value, the expression for ρ ( k ) is a sum of finite cosine functions. Using the Taylor series expansion of cosine functions, it can be determined that there exists a unique constant array C satisfying:
ρ ( k ) = c 0 + c 2 x 2 + c 4 x 4 + o ( x 4 )
Then, based on Equation (14), F i can be estimated using a fitting method. From Equations (20) and (21), it can be seen that the fitting function can be selected from the following three functions: α cos ( b x ) , c 1 x 2 , c 1 x 2 + c 2 x 4 (cosine fitting, quadratic fitting, and quartic fitting). In the fitting process, the fitting points are evenly distributed around τ = 0 but do not include the point 0.
(2) The cosine functions ρ ( k ) in Equation (20) have periods related to the angular frequency of the original signal and the sampling frequency. Therefore, the accuracy of the fitting results can be improved by suppressing noise through an appropriate choice of sampling frequency. In cosine fitting, if the angular frequencies of the parameters in ρ ( k ) are known, the fitting process can be simplified. Otherwise, it is necessary to estimate the angular frequencies of each component, complicating the fitting function and making the process more time-consuming. However, radar data typically involve a large amount of data, and complex fitting calculations are not suitable for real-time radar data processing.
Taking two 200-point sine wave sequences s ( t ) = sin ( t ) and s ( t ) = sin ( t + 5 ) ( t = 0 , 0.1 , , 20 ) as an example, additive Gaussian white noise with SNRs of 0 dB, 2 dB, and 4 dB is added to s ( t ) and s ( t ) to generate noisy sine wave sequences s n ( t ) and s n ( t ) .
The autocorrelation functions, obtained after adding additive Gaussian white noise, are then fitted with a second-order polynomial using 10 points around the zero-delay point to calculate the noise factor. Finally, the curves of the autocorrelation and cross-correlation functions in Figure 2 are denoised and corrected, as shown in Figure 3.
Comparing Figure 2 with Figure 3, it is observed that the noise factor fitted by the polynomial significantly corrects the amplitude of the autocorrelation function. According to the definition of SNR, an estimation of its value can be made.
Regarding the application analysis of fitting methods, although the radar can receive thousands of raw data points per second, it requires hundreds of original data to calculate a single data point. Therefore, in this study of cos ( 2 π f k ) , the range of f is generally from 2 < f < 10 corresponding to a large SNR error. However, for the FCA method, the academic community is more concerned with the calculation of the autocorrelation function. The amplitude correction of the distorted correlation function is performed using the proposed polynomial fitting-based noise reduction (PFNR) method, as shown in Algorithm 1.
Algorithm 1. PFNR Algorithm
Input: Receive signal; sampling frequency; fitting order
Output: Correlation function amplitude correction
1. Calculate correlation function.
 According to Equations (4) and (6) to calculate correlation function.
2. Polynomial fitting.
 2.1 Remove the point where the autocorrelation function is at 0 lags.
 2.2 Set an appropriate fitting order.
 2.3 Take an appropriate number signal data points around 0 lags for fitting.
 2.4 Use fitting to obtain the value of autocorrelation function at 0 lags.
3. Calculate the noise factor F i .
 According to Equations (11) and (12) to calculate F i .
4. Correct the amplitude of the correlation function.
 Multiply F i by autocorrelation function without the value of 0 lags to obtain amplitude correction.
Taking τ 0.5 as an example, its equivalent point is the delay point corresponding to the half-wave width point of the autocorrelation function curve. Figure 4 illustrates the comparison results between the original autocorrelation function (left) and the autocorrelation function after noise reduction using the quadraticfitting method (right).
It can be observed that, although, theoretically, fitting methods result in larger errors in low-frequency sampled signals, they still exhibit significant effectiveness in estimating the radar autocorrelation function. Moreover, these methods are straightforward and do not significantly increase radar data processing time, making them suitable for practical applications. In the application of the proposed PFNR algorithm, the fitting function’s order n and frequency can influence fitting performance, and hence the estimation performance of τ 0.5 . In the following section, we will investigate this problem.

4. Performance of the PFNR Algorithm

To evaluate the performance of various methods for estimating atmospheric wind speed, certain quantitative indicators are necessary. The Mean Squared Error (MSE) measures the difference between an estimated value and an actual value. It calculates the average of the squares of the differences between the estimated and actual values, which can be used to evaluate the prediction accuracy of the model. The smaller the MSE, the higher the accuracy of the model’s predictions. Therefore, MSE is used to evaluate performance. Currently, existing methods for estimating wind speed include direct estimation and indirect estimation through relevant parameters [21,22]. In the process of the proposed PFNR algorithm, F i , or, equivalently, ρ ^ ( k | k = 0 ) , is a byproduct in the estimation of delay parameters and hence the wind speed. We find that:
SNR e s t i m a t e = 10 log 10 [ | s i ( t ) | 2 ] [ | n i ( t ) | 2 ] = 10 log 10 [ | s i ( t ) | 2 ] [ | s i ( t ) | 2 ] + [ | n i ( t ) | 2 ] [ | s i ( t ) | 2 ] + [ | n i ( t ) | 2 ] [ | n i ( t ) | 2 ] = 10 log 10 ρ ^ ( k | k = 0 ) 1 ρ ^ ( k | k = 0 )
in which, ρ ^ ( k | k = 0 ) indicates that it is an estimated value using the polynomial fitting method according to the proposed PFNR algorithm. Thus, the value of S N R e s t i m a t e can be a measure of the proposed PFNR algorithm. Using Equation (22), the estimated value of SNR can be calculated from the value at the zero-delay point of the autocorrelation function ρ ^ ( k | k = 0 ) . The difference between the estimated S N R e s t i m a t e and the true value S N R t r u e is S N R e r r o r = S N R e s t i m a t e S N R t r u e , and it can serve as a metric to evaluate the performance of the noise reduction process. The closer the estimation error S N R e r r o r is to zero, the better the performance of the noise reduction. In this work, the polynomial fitting method is used to correct the amplitude of the correlation function and evaluate it by the MSE of SNR. The factors affecting the evaluation mainly include the order of fitting and the sampling frequency. These two influencing factors are analyzed in the following sections.
For a fair comparison, the range of data used is fixed from −20 to 20 lags. Simulations are conducted for two cases of additive white noise signals:
(1) First, a sinusoidal sequence s ( t ) = sin ( t )   ( t = 0 , 0.1 , , 20 ) is generated. Gaussian white noise with a SNR of 0 dB is added to s ( t ) . Finally, based on Equation (22), the mean square value of the error data for 10,000 cases of s n ( t ) under SNR of 0 dB is calculated as the number of fitting points increases. Figure 5 corresponds to s ( t ) . From Figure 5, it can be observed that, under the same conditions, the S N R e s t i m a t e obtained by the quadratic fitting is the lowest. Simultaneously, Table 2 presents the fitting time complexity of different methods under the same conditions and a fitting point of 10, indicating that the quadratic fitting exhibits the lowest time complexity.
(2) To better simulate real-world scenarios, the signal y = sin ( 2 x + 5 ) + 2 sin ( x )   ( t = 0 , 0.1 , , 20 ) is chosen. In view of the results of simulation (1), quadratic fittings are adopted, and MSE is calculated, as shown in Figure 6. From Figure 6, it is observed that, under the same conditions, the range of MSE decreases with the increase in sampling frequency. This is mainly attributed to the cosine factor cos ( 2 π f k ) in the theoretical expression of the autocorrelation function ρ ^ ( k ) , which becomes finer “interpolation” for cos ( k ) with increasing sampling frequency, thus better suppressing noise. According to the simulation results, when f = 5 , the minimum value of MSE is around 1 (Figure 6); when f = 10 , the minimum value of MSE is around 0.7; and when f = 100 , the minimum value of MSE is around 0.2. Since the error curve of quadratic fitting approximates a parabola, quadratic fitting should be preferred when the sampling frequency is low. Otherwise, when fitting with more points, quartic fitting should be preferred.

5. Conclusions

This work proposes an amplitude correction method for the noise factor of sinusoidal sequences using the correlation function, which can be applied to data processing of the autocorrelation and cross-correlation functions of MF Radar. The MSE of noise reduction with different fitting methods is analyzed as the number of fitting points changes. In the scenario set in this study, since the error curve of quadratic fitting approximates a parabola, quadratic fitting performs better than other fitting methods. For a deeper analysis, the impact of different sampling frequencies on polynomial fitting is further examined through numerical simulation under the same SNR. Simulation results indicate that a higher data sampling frequency leads to better amplitude correction of the correlation function, while a smaller number of fitting points is suitable for quadratic fitting.

Author Contributions

Conceptualization, J.C.; methodology, L.W. and J.C.; formal analysis, N.L. and G.K.; resources, J.C.; writing—original draft preparation, J.C. and Y.Z.; writing—review and editing, J.W. and J.C.; supervision, J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research work was supported by the JKW173 Fund (A072201405), the Academician Studio Fund (A392301036), Guangxi Minzu University Research Fund (2023KJQD37), Guangxi University Young and Middle-aged Teachers’ Basic Scientific Research Ability Improvement Project (2024KY0171), Xiangsi Lake Youth Scholar Innovation Team of Guangxi Minzu University (2023GXUNXSHQN06), Graduate Education Innovation Program of Guangxi University for Nationalities (GXUN-CHXS2024122).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Forbes, J.M.; Garrett, H.B. Theoretical studies of atmospheric tides. Rev. Geophys. 1979, 17, 1951–1981. [Google Scholar] [CrossRef]
  2. Fritts, D.C.; Alexander, M.J. Gravity wave dynamics and effects in the middle atmosphere. Rev. Geophys. 2003, 41. [Google Scholar] [CrossRef]
  3. Wang, J.; Yi, W.; Wu, J.; Chen, T.; Xue, X.; Zeng, J.; Vincent, R.; Reid, I.; Batista, P.; Buriti, R.; et al. Coordinated observations of migrating tides by multiple meteor radars in the equatorial mesosphere and lower thermosphere. J. Geophys. Res. Space Phys. 2022, 127, e2022JA030678. [Google Scholar] [CrossRef]
  4. Yang, C.; Lai, D.; Yi, W.; Wu, J.; Xue, X.; Li, T.; Chen, T.; Dou, X. Observed Quasi 16-Day Wave by Meteor Radar over 9 Years at Mengcheng (33.4° N, 116.5° E) and Comparison with the Whole Atmosphere Community Climate Model Simulation. Remote Sens. 2023, 15, 830. [Google Scholar] [CrossRef]
  5. Yi, W.; Xue, X.; Lu, M.; Zeng, J.; Ye, H.; Wu, J.; Wang, C.; Chen, T. Mesopause temperatures and relative densities at midlatitudes observed by the Mengcheng meteor radar. Earth Planet. Phys. 2023, 7, 665–674. [Google Scholar] [CrossRef]
  6. Zeng, J.; Stober, G.; Yi, W.; Xue, X.; Zhong, W.; Reid, I.; Adami, C.; Ning, B.; Li, G.; Dou, X. Mesosphere/lower thermosphere 3-dimensional spatially resolved winds observed by Chinese multistatic meteor radar network using the newly developed VVP method. J. Geophys. Res. Atmos. 2024, 129, e2023JD040642. [Google Scholar] [CrossRef]
  7. Wilhelm, S.; Stober, G.; Chau, J.L. A comparison of 11-year mesospheric and lower thermospheric winds determined by meteor and MF radar at 69 N. Annales Geophysicae. Copernic. GmbH 2017, 35, 893–906. [Google Scholar]
  8. Thayaparan, T.; Hocking, W.K. A long-term comparison of winds and tides measured at London, Canada (43° N, 81° W) by co-located MF and meteor radars during 1994–1999. J. Atmos. Sol.-Terr. Phys. 2002, 64, 931–946. [Google Scholar] [CrossRef]
  9. Tsutsumi, M.; Aso, T. MF radar observations of meteors and meteor-derived winds at Syowa (69° S, 39° E), Antarctica: A comparison with simultaneous spaced antenna winds. J. Geophys. Res. Atmos. 2005, 110. [Google Scholar] [CrossRef]
  10. Zeng, J.; Yi, W.; Xue, X.; Reid, I.; Hao, X.; Li, N.; Chen, J.; Chen, T.; Dou, X. Comparison between the mesospheric winds observed by two collocated meteor radars at low latitudes. Remote Sens. 2022, 14, 2354. [Google Scholar] [CrossRef]
  11. Gaikwad, H.P.; Fadnavis, S.; Gurav, O.B.; Bhosale, J.L.; Sharma, A.K.; Patil, P.T.; Ghodpage, R.N.; Sutar, M.M.; Chavan, G.A.; Birajdar, A.T.; et al. The day-to-day variability in the mesosphere and lower thermosphere in low latitudes: A study using MF radar. Adv. Space Res. 2023, 71, 199–215. [Google Scholar] [CrossRef]
  12. Liu, L.; Wan, W. Recent ionospheric investigations in China (2018–2019). Earth Planet. Phys. 2020, 4, 179–205. [Google Scholar] [CrossRef]
  13. Tang, Z.; Li, N.; Wang, J.; Ding, Z.; Dai, L.; Zhao, L.; Chen, J. Variation of Electron Density in the D-Region Using Kunming MF Radar under Low Solar Activity. Atmosphere 2023, 14, 1764. [Google Scholar] [CrossRef]
  14. Wan, W.; Xu, J. Recent investigation on the coupling between the ionosphere and upper atmosphere. Sci. China Earth Sci. 2014, 57, 1995–2012. [Google Scholar] [CrossRef]
  15. Zhao, J. Research on General Design and Key Technologies of Wuhan MF Radar. Ph.D. Thesis, Wuhan University, Wuhan, China, 2009. [Google Scholar]
  16. Hu, X.; Zeng, Z.; Zhang, D. Observations of the mesosphere and lower thermosphere angular spectrum with an MF radar at Wuhan. Chin. J. Space Sci. 2003, 23, 256–261. [Google Scholar]
  17. Li, F.; Hu, X.; Zhang, D.; Huang, Z.; Xiong, J.; Ceng, Z. Mesospheric MF radar in Wuhan and its preliminary observation results. Chin. J. Space Sci. 2002, 22, 65–71. [Google Scholar] [CrossRef]
  18. Reid, I.M. MF and HF radar techniques for investigating the dynamics and structure of the 50 to 110 km height region: A review. Prog. Earth Planet. Sci. 2015, 2, 33. [Google Scholar] [CrossRef]
  19. Hocking, W.K.; Kolomiyets, S.V. Radio meteor physics—A comparison between techniques from 1945 to the mid-1970’s. Radiotekhnika 2020, 201, 78–90. [Google Scholar] [CrossRef]
  20. May, P.T. Statistical errors in the determination of wind velocities by the spaced antenna technique. J. Atmos. Terr. Phys. 1988, 50, 21–32. [Google Scholar] [CrossRef]
  21. Doviak, R.J.; Zhang, G.; Cohn, S.A.; Brown, W.O. Comparison of spaced-antenna baseline wind estimators: Theoretical and simulated results. Radio Sci. 2004, 39, 1–21. [Google Scholar] [CrossRef]
  22. Zhang, G.; Doviak, R.J.; Vivekanandan, J.; Brown, W.O.; Cohn, S.A. Performance of correlation estimators for spaced-antenna wind measurement in the presence of noise. Radio Sci. 2004, 39, 1–16. [Google Scholar] [CrossRef]
  23. Zhao, L.; Chen, J.; Li, N.; Zhao, Z.; Wu, J. MF radar in Kunming and its preliminary observation results. Chin. J. Space Sci. 2011, 31, 27–33. [Google Scholar] [CrossRef]
  24. Holdsworth, D.A. Signal Analysis with Applications to Atmospheric Radars. Ph.D. Thesis, University of Adelaide, Adelaide, Australia, 1995. [Google Scholar]
  25. Briggs, B.H.; Phillips, G.J.; Shinn, D.H. The analysis of observations on spaced receivers of the fading of radio signals. Proc. Phys. Soc. Sect. B 1950, 63, 106. [Google Scholar] [CrossRef]
  26. Li, M.; Deng, Y.; Harding, B.J.; England, S. Climatology of dayside E-region zonal neutral wind shears from ICON-mighti observations. Space Weather 2024, 22, e2023SW003670. [Google Scholar] [CrossRef]
  27. Quan, L.; Cai, B.; Hu, X.; Xu, Q.; Li, L. Study of ionospheric D region changes during solar flares using MF radar measurements. Adv. Space Res. 2021, 67, 715–721. [Google Scholar] [CrossRef]
  28. Richards, M.A. Fundamentals of Radar Signal Processing; Mcgraw-Hill: New York, NY, USA, 2005. [Google Scholar]
Figure 1. Correlation index τ , τ 0.5 , τ , and τ .
Figure 1. Correlation index τ , τ 0.5 , τ , and τ .
Atmosphere 15 00899 g001
Figure 2. Autocorrelation function after adding additive white Gaussian noise with different SNR (left); cross-correlation function after adding additive white Gaussian noise with different SNR (right).
Figure 2. Autocorrelation function after adding additive white Gaussian noise with different SNR (left); cross-correlation function after adding additive white Gaussian noise with different SNR (right).
Atmosphere 15 00899 g002
Figure 3. Quadratic fitting denoising results of autocorrelation function under different SNR.
Figure 3. Quadratic fitting denoising results of autocorrelation function under different SNR.
Atmosphere 15 00899 g003
Figure 4. The change of mean square error of SNR with the increase of fitting points under different sampling frequencies.
Figure 4. The change of mean square error of SNR with the increase of fitting points under different sampling frequencies.
Atmosphere 15 00899 g004
Figure 5. The change in mean square error of SNR with the increase in fitting points under a sampling frequency of 10 Hz.
Figure 5. The change in mean square error of SNR with the increase in fitting points under a sampling frequency of 10 Hz.
Atmosphere 15 00899 g005
Figure 6. Under complex conditions, the mean square error with different sampling frequency and SNR of 1dB varies with the number of fitting points.
Figure 6. Under complex conditions, the mean square error with different sampling frequency and SNR of 1dB varies with the number of fitting points.
Atmosphere 15 00899 g006
Table 1. The definitions and calculation formulas of the correlation functions, the influence of noise on the correlation functions, and the corresponding noise reduction methods.
Table 1. The definitions and calculation formulas of the correlation functions, the influence of noise on the correlation functions, and the corresponding noise reduction methods.
Autocorrelation FunctionCross-Correlation Function
Definition without noise ρ α ( τ k ) = [ f α ( t ) f α ( t + τ k ) ] [ | f α ( t ) | 2 ] ρ α b ( τ k ) = [ f α ( t ) f b ( t + τ k ) ] [ | f α ( t ) | 2 ] [ | f b ( t ) | 2 ]
Calculation without noise ρ α ( τ k ) = n = 1 N r | k | f α ( t ) f α ( t + τ k ) ( N r | k | ) n = 1 N r | f α ( t ) | 2 ρ α b ( τ k ) = n = 1 N r | k | f α ( t ) f b ( t + τ k ) n = 1 N r | f α ( t ) | 2 n = 1 N r | f b ( t ) | 2
With noise ρ ( τ ) = 1 , τ = 0 [ s i ( t ) s i ( t + τ ) ] [ | s i ( t ) | 2 ] + [ | n i ( t ) | 2 ] , τ 0 ρ i j ( τ ) = [ s j ( t ) s j ( t ) ] + [ n c i ( t ) n c j ( t + τ ) ] ( [ | s i ( t ) | 2 ] + [ | n i ( t ) | 2 ] ) ( [ | s j ( t ) | 2 ] + [ | n j ( t ) | 2 ] ) , ( τ = 0 ) s i ( t ) s j ( t ) ( [ | s i ( t ) | 2 ] + [ | n i ( t ) | 2 ] ) ( [ | s j ( t ) | 2 ] + [ | n j ( t ) | 2 ] ) , ( τ 0 )
Noise reduction ρ ( τ ) = ρ ( τ ) F i = [ s ( t ) s i ( t + τ k ) ] [ | s ( t ) | 2 ] F i = [ | s i ( t ) | 2 ] + [ | n i ( t ) | 2 ] [ | s i ( t ) | 2 ] ρ i j ( τ ) = ρ i j ( τ ) F i F j = [ s i ( t ) s j ( t + τ k ) ] [ | s i ( t ) | 2 ] [ | s j ( t ) | 2 ] F i F j = [ | s i ( t ) | 2 ] + [ | n i ( t ) | 2 ] [ | s i ( t ) | 2 ] [ | s j ( t ) | 2 ] + [ | n j ( t ) | 2 ] [ | s j ( t ) | 2 ]
Table 2. The fitting time complexity of various fitting methods with the number of fitting points of 10.
Table 2. The fitting time complexity of various fitting methods with the number of fitting points of 10.
Fitting MethodsRuntime (s)
x 2 fitting0.000019
x 4 fitting0.000021
x 6 fitting0.000023
x 8 fitting0.000026
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, J.; Zhang, Y.; Wang, L.; Kang, G.; Li, N.; Wei, J. Polynomial Fitting-Based Noise Reduction for Correlation Functions in Medium-Frequency Radar. Atmosphere 2024, 15, 899. https://doi.org/10.3390/atmos15080899

AMA Style

Chen J, Zhang Y, Wang L, Kang G, Li N, Wei J. Polynomial Fitting-Based Noise Reduction for Correlation Functions in Medium-Frequency Radar. Atmosphere. 2024; 15(8):899. https://doi.org/10.3390/atmos15080899

Chicago/Turabian Style

Chen, Jinsong, Yang Zhang, Liming Wang, Guoqin Kang, Na Li, and Junfeng Wei. 2024. "Polynomial Fitting-Based Noise Reduction for Correlation Functions in Medium-Frequency Radar" Atmosphere 15, no. 8: 899. https://doi.org/10.3390/atmos15080899

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop