1. Introduction
Frequency-modulated continuous-wave (FMCW) radar is a widely utilized radar system, esteemed for its high measurement accuracy and robust anti-interference capabilities. Its performance is minimally affected by adverse weather conditions such as rain and snow, rendering it highly suitable for applications in military, traffic management, automotive collision avoidance, and liquid level measurement [
1]. The FMCW radar generates a beat signal by mixing the transmitted signal with the reflected echo, thereby deriving information on the target’s distance and velocity [
2]. However, the beat signal acquired by FMCW radar is contaminated by noise, leading to significant frequency errors when using direct FFT. Therefore, the accurate estimation of the beat signal’s frequency is a crucial problem in FMCW radar signal processing.
The purpose of frequency estimation is to estimate the distance or velocity of a target object by using the phase difference or time delay of the beat signal, and its accuracy directly affects the accuracy of the distance and velocity measurements [
3]. Due to the non-integer relationship between the sampling frequency and the actual frequency, spectral leakage [
4] and the fence effect [
5] can result in frequency estimation biases. In practical engineering, FFT-based interpolation is typically employed to enhance the accuracy of frequency estimation. Rife et al. [
6] proposed the Rife algorithm, which employs double spectral line interpolation. This algorithm uses the ratio between the maximum spectral line and the adjacent second-largest spectral line to interpolate and estimate the signal, thereby improving the accuracy of frequency estimation. Nevertheless, when the actual frequency is close to the quantized frequency point, the frequency estimation accuracy decreases. Dou et al. [
7] proposed an improved Quinn and Rife Automatic Segmentation algorithm, which adaptively selects the optimal estimation method for the frequency interval, thereby improving the stability of the algorithm. Nevertheless, the complexity is not considered. Nian et al. [
8] proposed an anticipated Rife interpolation algorithm. This algorithm improved the original single-discriminant method of Rife, enhancing the estimation stability within any frequency range. Nevertheless, the algorithm’s noise suppression capability is relatively weak.
To reduce estimation errors caused by noise and spectral leakage, and to balance both anti-interference performance and complexity, a FMCW beat signal frequency estimation algorithm based on spectral refinement and phase interpolation (SP–Rife) is proposed. This algorithm achieves a more accurate frequency estimation by improving the signal resolution and utilizing phase interpolation. The main steps are as follows:
Spectral Refinement: The post-FFT spectral signal is refined to narrow the frequency search range and enhance frequency resolution, effectively separating the noise signal.
Frequency Deviation Factor is Defined: By defining the frequency deviation factor as the ratio of the maximum spectral amplitude to the sum of the maximum and second-largest spectral amplitudes, the relationship of spectral proportions is established. This approach enhances the algorithm’s capability to correct frequency deviations.
Phase Angle Interpolation: Leveraging the stability of phase angles, the algorithm converts spectral amplitudes into phase angles to determine the frequency deviation factor and interpolation direction, thereby increasing the estimation accuracy.
The rest of this paper is organized as follows: Firstly, a comprehensive review of existing frequency estimation algorithms is presented in
Section 2. Then, the SP–Rife algorithm is proposed in
Section 3. Furthermore, the simulation and experimental results are provided to validate the performance of the proposed algorithm in
Section 4. Finally, the research work of this paper is summarized in
Section 5.
2. Related Works
The frequency measurement algorithms with a high accuracy can be categorized as autocorrelation-function-based time-domain frequency estimation algorithms and Fourier-transform-based frequency-domain estimation algorithms. Basic time-domain frequency estimation algorithms mainly include the Fitz algorithm [
9] and the Kay algorithm [
10] and their improved algorithms for the signal-to-noise ratio threshold and constraint parameters. These methods estimate the frequency value from the phase relationship. Frequency-domain frequency estimation algorithms primarily rely on spectral peaks for signal frequency estimation. Examples include the Rife algorithm and the Quinn algorithm [
11], as well as improved algorithms designed to address errors that arise when estimated frequency values are close to quantized frequency points. Time-domain estimation algorithms have higher requirements on the signal, such as the signal length and signal-to-noise ratio threshold. In contrast, frequency estimation algorithms are more universally applicable [
12].
In response to the fence effect and spectral leakage effects caused by the Fourier transform, Wang et al. [
13] proposed a fast and accurate frequency estimation method for monotone signals based on DFT interpolation. Firstly, the frequency deviation direction was accurately identified by the three-sample DFT phase difference, and then a new unbiased estimator was developed for an arbitrary number of samples, which solved the problems of the accuracy of the existing methods in noisy conditions and the bias in a small number of samples. However, the method was still limited by the signal-to-noise ratio threshold in finding the peak DFT in the initial coarse estimation. Serbes et al. [
14] proposed the HAQSE algorithm, which employed a hybrid approach combining the half-shift and Q-shift DFT interpolation methods. The first iteration of the algorithm was modified using an unbiased A&M interpolator and then refined using the QSE algorithm to obtain the final frequency estimate. The algorithm required only two iterations, which reduced the computational cost, but the algorithm limited the computational length of the signal and performs poorly above a certain threshold.
Among frequency estimation algorithms, the Rife algorithm, which relied on the ratio of spectral amplitude values, was widely utilized in signal processing due to its relatively low computational complexity and minimal requirements for signal length. However, misjudging the direction of the second-largest spectral value leads to substantial estimation errors when the signal-to-noise ratio is relatively low. Additionally, while the Rife algorithm provided a high accuracy when the estimated frequency was near the center of the quantization frequency, its performance deteriorated significantly and errors increased when the estimated frequency was close to the quantization frequency.
In order to improve the accuracy of frequency interpolation algorithms in the frequency domain, Chen et al. [
15] proposed a non-iterative accurate frequency estimator using two DFT spectral lines. The method calculated the Fourier coefficients on the two spectral lines located midway between the two spectral lines adjacent to the maximum value and used the properties of sine function fitting to compute the amplitude ratios of the two intermediate spectral lines to the maximum spectral line, in order to estimate the position of the actual maximum magnitude Fourier coefficient. The complexity of the proposed method was lower than that of traditional methods. However, the method could only achieve CRLB performance under high signal-to-noise ratio conditions, and the estimation bias was larger when the frequency was close to the band edges.
The primary sine frequency estimation method was only suitable for signals with zero padding, resulting in a decrease in frequency estimation performance. Jiang et al. [
16] proposed a novel method for the frequency estimation of non-stationary discrete-time sinusoidal signals. They first estimated the signal parameters through a linear regression model, and then constructed a frequency estimation algorithm using the standard gradient descent approach. The method ensured the global exponential convergence of the frequency estimation and had a certain degree of robustness. However, the method relied on the persistently exciting (PE) assumption of the signal, and its performance could be affected when the noise level was high or the signal did not meet the PE condition. Zhang et al. [
17] introduced a frequency estimation method based on the amplitude ratio of two DFT samples from zero-filled signals. This method employed the analytical expression of the DFT signal after zero-filling, defining the amplitude ratio function between the two DFT samples nearest to the actual frequency. Utilizing the monotonicity of this ratio function, the residual frequency was estimated using either the Newton–Raphson method or a leveling search. While this method markedly improved the estimation accuracy, it also increased the complexity of signal processing and the overall algorithm.
To tackle the issue of incorrect interpolation direction in estimation algorithms that rely on spectral value ratios, Li et al. [
18] proposed a frequency estimation method that uses the peak of the sidelobe. It was not always possible to determine the peak of the first sidelobe from the DFT spectrum due to the fence effect. The method zero-filled the signal and then performed a Fourier transform to accurately locate the peaks of the main flap and the corresponding peaks of the first para flap. However, the method was significantly affected by noise. It was suitable for signals with high signal-to-noise ratios, and a sufficiently large number of samples were required to ensure that the method contributed to the identification of the peak. Borkowski et al. [
19] proposed frequency estimation in the interpolated discrete Fourier transforms with a generalized maximum sidelobe attenuation window, defining the Fourier analysis concerning the generalized MSD (GMSD) window. By double-zero padding, the frequency values were estimated using a three-point interpolation algorithm, which provided a greater flexibility for noisy signals but had specific requirements on the signal-to-noise ratio of the processed signals.
Cheng et al. [
20] introduced a phase angle interpolation-based improved Rife frequency estimation algorithm called PAI–Rife, which utilized frequency shifting techniques. First, they redefined the existing interpolation method. Then, frequencies outside the frequency shift threshold range were moved to the optimal estimation space. This proposed algorithm reduced the impact of the quantization frequency distribution on estimation accuracy and improved the frequency estimation precision. However, the algorithm’s noise resistance was limited, being effective only within the range of −10 dB to 10 dB, making it unsuitable for harsher environments.
Based on the frequency bias for the corrected frequency estimation, Shi et al. [
21] proposed the Golden–MacLeod (GM) algorithm. By improving the MacLeod algorithm and employing an iterative method to obtain the frequency bias, followed by the application of the golden ratio algorithm for frequency estimation, this method significantly reduced estimation errors and enhanced stability. However, the iterative process was influenced by the number of sample points, which led to an increased computational load.
An analysis of the discussed algorithms indicates that the performance of frequency estimation methods is predominantly limited by the trade-off between accuracy and efficiency, especially under low signal-to-noise ratio conditions. The key research objective is to reconcile the reduced computational complexity with the improved accuracy. This paper tackles this issue by combining phase spectrum refinement and phase angle interpolation with the Rife algorithm, aiming to optimize the overall performance of frequency estimation.
3. SP–Rife Algorithm
The SP–Rife algorithm is proposed to enhance the precision of frequency estimation and improve the algorithm’s noise resistance by combining spectral refinement and phase angle interpolation. The algorithm is primarily divided into three parts: spectral refinement, the definition of the frequency deviation factor, and phase angle interpolation estimation.
3.1. Spectral Refinement
Due to the limitation in frequency resolution, it is not possible to obtain the frequencies between adjacent discrete points, which may cause some useful frequencies to be missed. Increasing the number of sampling points will solve this problem, but it will result in higher hardware costs. The spectral refinement technique can improve the signal’s resolution, making it easier to identify and suppress noise signals. Therefore, to enhance the accuracy of frequency estimation at low signal-to-noise ratios, spectral refinement processing is performed on the signal before conducting interpolation estimation in the SP–Rife algorithm.
Let the index corresponding to the maximum value in the spectrum be denoted by
, and the index corresponding to the second-largest value be denoted by
. The corresponding spectra are represented as
and
where
is the signal length,
is the number of sampling points, and
is the sine signal.
Due to the fence effect, the frequency between points
and
has not been detected. Therefore, it is necessary to further refine the frequency resolution by a factor of
between the original adjacent spectral lines. In other words, the spectral amplitudes of the frequencies at the following frequency points need to be obtained:
After performing spectral refinement, it is only necessary to calculate the spectra corresponding to the mentioned frequency points, without computing the entire spectrum. For the frequency spectrum of each frequency point
between
and
, use Equation (4) to calculate:
where
is an integer from 0 to
L − 1 indicating the deviation of the refined frequency point with respect to
and
. In this way, a more refined spectral representation can be obtained between
and
, thus enabling a more accurate estimation of the frequency.
Subsequently, the search for the maximum and second-largest spectral values is performed again. In the SP–Rife algorithm, a 100-fold spectral refinement is conducted. Additionally, to find the true maximum spectral peak, the real part of the spectrum is used for determination.
3.2. Defining the Frequency Deviation Factor
The Rife algorithm is known for its high estimation accuracy when the estimated frequency is situated between two quantized frequency points. However, its accuracy diminishes as the estimated frequency approaches these quantized points. To address the limitations of the Rife algorithm and mitigate the effects of spectral leakage, the original sinusoidal signal is processed using the Rife–Vincent (I) window, and the frequency deviation factor is redefined.
The Rife–Vincent (I) window function is expressed as follows [
20]:
where
is the order of the window function,
is the sampling period, and
is the rectangular window.
Given the symmetry of the DFT for a real sequence, only the positive frequency component of the spectrum is considered in this context. The representation of the frequency spectrum for the windowed signal
is given by:
where
is the amplitude of the sinusoidal signal, and
is the index of the quantized frequency point corresponding to the true frequency of the signal.
The main frequency components of the signal can be localized by . However, relying on the maximum spectrum alone is not sufficient to accurately estimate the true frequency because quantization errors lead to deviations in the signal frequency. By considering both and , the energy distribution of the signal in the vicinity of the quantized frequency point can be obtained. This distribution usually reflects the true location of the signal frequency. If the true frequency of the signal is closer to the quantization point , then will be significantly larger than ; conversely, if the true frequency is closer to the quantization point , then will be relatively large. By calculating the ratio of these two spectra, it is possible to quantify this difference in energy distribution and thus estimate the degree of deviation of the signal frequency relative to the quantization point.
By utilizing Equation (6), the values of
and
can be obtained:
and
When
, the Rife–Vincent (I) window function does not exhibit any sidelobes, resulting in minimal impact on the signal [
6].
The true frequency of the signal always lies between and , with an interval of between these two indices. Additionally, it is known that the sampling frequency and the actual frequency do not have an integer relationship in reality, resulting in a deviation between and , known as the frequency deviation factor. It satisfies the relation . Based on the relationship between the spectra, the frequency deviation factor can be determined using the ratio of and .
Therefore, the redefined frequency deviation factor can be expressed as
The provides a way to quantify the frequency deviation. When is close to 1, it means that the energy of the signal is mainly concentrated near , which suggests that the true frequency of the signal may be closer to the quantization frequency point , whereas, when is close to 0, it means that the energy of the signal is more uniformly distributed between and , which may suggest that the true frequency of the signal is located in between these two quantization points. By using as a compensation for interpolation, the energy distribution of the signal near the quantization points is taken into account, allowing for the more accurate estimation between the quantization frequency points.
After substituting the refined spectrum into Equation (9), it is possible to perform interpolation correction among the refined frequency points.
3.3. Phase Angle Interpolation and Criterion
The phase angle varies within the range of 0 to , offering greater stability compared to the amplitude. Therefore, in the SP–Rife algorithm, frequency estimation is conducted through phase angle interpolation.
The spectrum of a sinusoidal signal can be further expressed in terms of phase and amplitude as follows:
From the above equations, it can be concluded that the magnitude of the phase angle corresponding to the spectrum is determined by
and
. Given
, when
corresponds to the quantized frequency point with the maximum spectrum value,
. It can be concluded that
and
are positive, while
is negative. As shown in
Figure 1, the denominator of
is always positive when
. When
, both the denominator and numerator of
are positive, so
is positive. When
, the numerator of
is positive and the denominator is negative, so
is negative.
Assuming
is in the first quadrant, it can be derived that
, and both
and
are in the third quadrant. From Equation (10), the following is obtained:
From the above analysis, it can be inferred that the phase angles corresponding to and are situated in the first quadrant, while the phase angle corresponding to is in the third quadrant. Given the substantial phase angle variation at , interpolation correction will be applied between and . The direction of frequency deviation correction is also determined based on the extent of the phase angle variation. Similarly, for other quadrant conditions of , this analysis allows for analogous conclusions.
Combining Equation (9), the frequency deviation factor
based on the phase angle can be obtained:
where
is the phase angle function, and
represents the correction direction. The determination of the correction direction is based on the following conditions:
If , the true frequency point is to the right of the maximum spectral position . If , the true frequency point is to the left of .
After frequency spectrum refinement and phase angle interpolation, the corrected frequency is calculated using Equation (16).
where
is the signal sampling frequency.
3.4. SP–Rife Algorithm
The SP–Rife algorithm, as presented in Algorithm 1, is proposed based on the comprehensive content from
Section 3.1,
Section 3.2,
Section 3.3. The computation process of the algorithm involves spectral refinement and phase angle interpolation estimation, aiming to optimize the issues of the spectral leakage, direction correction determination, and interpolation method.
Algorithm 1: SP–Rife Algorithm |
Input: | A sinusoidal signal |
Output: | Frequency estimation |
Initialize: | Windowing
|
Step 1: | Obtain the position of the maximum spectral
|
Step 2: | Refine the spectrum at by a factor of 100 and update |
Step 2: | Calculate and |
Step 3: | Calculate the phase angles , and at , and with phase angle function |
Step 4: | Calculate correction direction and frequency deviation If Else End of if |
Step 5: | Return |
The SP–Rife algorithm can be summarized as follows:
- 1.
The signal is windowed, and the position
corresponding to the maximum value in the spectrum is detected.
- 2.
The spectrum between
and
is refined 100 times, and the maximum and second-largest values in the refined spectrum are found.
- 3.
The phase function
is used to calculate the phase angle
corresponding to the maximum value in the spectrum and the phase angles
and
corresponding to the adjacent two frequencies.
- 4.
By comparing the values of
and
, the correction direction
and the frequency deviation factor
are determined. Assuming
, if
, interpolation is performed on the left side of the maximum spectrum, and
; if
, it is concluded that the actual frequency point is on the right side of the maximum spectrum, and .
- 5.
The final frequency estimation value is computed using Equation (16).
4. Simulation Analysis
The simulation and evaluation of the SP–Rife algorithm, along with the comparison to other algorithms (Rife [
6], ASIQ–Rife [
7], and PAI–Rife [
20]), are performed using the MATLAB platform. The evaluation metrics include computational complexity, and mean and standard deviation. To evaluate the accuracy of the estimation, additive Gaussian white noise with a zero mean is introduced into the signal, with the Cramer–Rao lower bound (CRLB) used as the benchmark. The simulation parameters are outlined in
Table 1.
4.1. Computational Complexity
In the Rife algorithm, a single N-point FFT operation requires
complex multiplications and
complex additions.
where
is the base-2 logarithmic operation.
During the algorithm execution, finding the maximum and second maximum values in the spectrum requires complex multiplications and complex additions.
Table 2 shows the computational complexities of different Rife algorithms. While the ASIQ–Rife algorithm avoids frequency deviation, its computational complexity is substantially higher due to the repeated use of the Rife and Quinn algorithms. The PAI–Rife algorithm recalculates the spectrum with a 40% probability following frequency shifting. On the other hand, the SP–Rife algorithm reduces the computational load by applying 100-fold spectrum interpolation to selected frequency points, limiting the number of computation points to 200. As a result, the SP–Rife algorithm presents a lower computational demand than the other two improved algorithms and exhibits only a 5.2% increase in complexity compared to the original Rife algorithm.
4.2. Estimated Performance at Different Frequencies
In this subsection, the performance of the SP–Rife algorithm is evaluated across various frequencies. The quantization frequency is set, is the quantization interval, and 11 quantization frequency points are taken between . The SNR is −4 dB.
The results comparing the SP–Rife algorithm with other Rife algorithms are illustrated in
Figure 2 and
Figure 3. The Rife algorithm exhibits significant estimation errors near quantized frequencies. Despite improvements in the ASIQ–Rife algorithm aimed at reducing estimation errors near quantized frequencies, its fundamental spectral interpolation method has not been altered. Consequently, the estimation accuracy remains consistent within the optimal range. The PAI–Rife algorithm changes the criteria for determining the correction direction and the interpolation method, resulting in a more accurate estimation of the frequency deviation factor, but there are still fluctuations in errors near quantized frequencies. In contrast, SP–Rife achieves the most stable frequency estimation performance by increasing the frequency resolution through spectral refinement.
Based on the statistical analysis of the standard deviation, the mean standard deviation of the Rife algorithm across 11 frequency points is approximately 0.068. For the ASIQ–Rife algorithm, the mean standard deviation is about 0.045, while the PAI–Rife algorithm has a mean standard deviation of approximately 0.018. In contrast, the SP–Rife algorithm exhibits a mean standard deviation of around 0.013. This indicates a significant improvement in frequency estimation accuracy with the SP–Rife algorithm, achieving an 80.9% reduction compared to the Rife algorithm and a 27.8% reduction compared to the PAI–Rife algorithm. The SP–Rife algorithm consistently delivers a superior frequency estimation accuracy across various frequency points relative to other algorithms.
4.3. Estimated Performance at SNRs
In addition to evaluating the algorithm’s estimation performance at various frequencies, the effects of noise on frequency estimation are also assessed. The quantized frequency point is set to .
Figure 4 and
Figure 5 present the comparison results, where (a) displays the data for SNR values ranging from −20 dB to 20 dB, and (b) offers an expanded view of (a) for SNR values between −10 dB and 10 dB. The ASIQ–Rife and PAI–Rife algorithms show an improved estimation accuracy compared to the original Rife algorithm under certain SNR conditions. Nevertheless, substantial errors are still evident at lower SNR levels. Conversely, the SP–Rife algorithm capitalizes on the stability of phase angles and spectral refinement techniques, effectively reducing the impact of noise on the frequency estimation. Consequently, the SP–Rife algorithm consistently achieves the most stable mean frequency estimation in low-SNR scenarios.
The calculations reveal that, within an SNR range of −10 dB to 10 dB, the Rife algorithm has an average standard deviation of approximately 0.063. Among the two enhanced Rife algorithms, the PAI–Rife algorithm performs better, with an average standard deviation of about 0.013. The SP–Rife algorithm delivers the best performance, with an average standard deviation of approximately 0.008. Compared to the Rife algorithm, the PAI–Rife algorithm shows a notable improvement of 87.3% in noise resilience, while the SP–Rife algorithm surpasses PAI–Rife with a 38.5% improvement in noise resilience. In low-SNR conditions, this enhancement is even more significant, with the standard deviation nearing the CRLB.
4.4. Correction Direction Misjudgment Rate
Frequency estimation errors can be attributed to two main categories: one is the inaccurate estimation of frequency deviation, leading to deviations in interpolation positions; the other is the misjudgment of the correction direction, resulting in increased frequency estimation errors. To address this, a comparative analysis of the misjudgment of the correction direction for each algorithm under different SNR conditions is shown in
Figure 6.
As previously analyzed, the Rife algorithm tends to exhibit a relatively high misjudgment rate for the correction direction when the frequency point is near the quantized frequency, leading to increased frequency errors after deviation correction. To address this issue, the ASIQ–Rife algorithm reassesses the correction direction following signal frequency shifting, which reduces the misjudgment rate. The PAI–Rife algorithm enhances the accuracy by shifting the signal to the optimal estimation range, thereby improving the correction direction judgment; however, significant errors persist under low-SNR conditions. In contrast, the SP–Rife algorithm effectively mitigates the effects of spectral leakage and noise by refining the frequency resolution through spectrum refinement. Additionally, it enhances the accuracy of correction direction judgment using phase angle interpolation. Consequently, the SP–Rife algorithm achieves a misjudgment rate for the correction direction that is significantly closer to zero and much lower compared to both the Rife and ASIQ–Rife algorithms, especially in high-noise environments.
4.5. Analysis of Actual Measurement Data
To validate the practical performance of the SP–Rife algorithm, actual data are collected using FMCW radar at five different distances. Subsequently, the SP–Rife algorithm is employed to estimate the frequency of the received echo signals containing distance information. The estimated frequency values are then used in Equation (16) to obtain the distance estimates.
where
is the sampling period,
is the speed of electromagnetic waves, and
is the signal bandwidth.
The experimental data acquisition environment is shown in
Figure 7. The metallic panel on the wall serves as the target for measurement. The radar-to-target distance is set using the launching platform, and then the equipment is fixed on the track. The FMCW radar transmits the signal, which reflects off the metallic panel in front and is received by the device’s receiving panel, generating beat signals. Before the measurement, the sampling frequency of the FMCW radar is set to
, the signal bandwidth
, and the sampling number
. The actual distances measured by the radar are 2378 mm, 2420 mm, 2561 mm, and 2653 mm.
The beat frequency of the FMCW radar’s beat signal is estimated using various frequency estimation algorithms, and the corresponding distances are calculated.
Table 3 displays the distance results obtained with different beat frequency estimation algorithms. For a more intuitive comparison, measurement errors from these algorithms are plotted in
Figure 8. The analysis reveals that the Rife algorithm provides the least accurate estimates across all datasets. The PAI–Rife algorithm reduces errors by optimizing the signal’s estimation range. Conversely, the SP–Rife algorithm achieves the most precise beat frequency estimates through spectrum refinement and phase interpolation, resulting in the smallest distance measurement errors.
5. Conclusions
Frequency leakage and the fence effect pose significant challenges in the estimation of FMCW signals. The commonly used Rife algorithm, while prevalent, is vulnerable to noise and has limitations concerning the frequency range. To address these issues, this paper proposes the SP–Rife algorithm, which leverages spectrum refinement and phase interpolation techniques. The algorithm enhances the frequency resolution through spectrum refinement, improving the signal’s robustness against noise. Additionally, by employing amplitude–phase conversion followed by interpolation estimation, the SP–Rife algorithm achieves a greater accuracy and precision in frequency estimation.
The simulation results indicate that the SP–Rife algorithm reduces the computational complexity by 30.1% compared to the Rife algorithm, demonstrating a significant advantage over the other two improved Rife algorithms. In comparison to the PAI–Rife algorithm, which shows a relatively strong estimation performance, the SP–Rife algorithm provides a 38.5% improvement in noise robustness and a 27.8% enhancement in frequency estimation accuracy. Overall, the SP–Rife algorithm achieves an effective balance between computational efficiency and accuracy, outperforming the Rife, ASIQ–Rife, and PAI–Rife algorithms in terms of overall estimation performance. Despite these advancements, the algorithm’s effectiveness is contingent upon the appropriate setting of the spectrum refinement factor, which is sensitive to varying environmental conditions. In future work, adaptive methods for adjusting the spectral refinement factor will be explored to enable automatic calibration based on the characteristics of signals and environmental noise, thereby enhancing the generality and robustness of the algorithm. The current work has solved the estimation accuracy, which provides a foundation for future parameter adaptation.