Computationally Efficient Robust Chinese Remainder Theorem for Frequency Estimation from Undersampled Waveforms
We present an efficient algorithm based on the robust Chinese remainder theorem (CRT) to perform single frequency determination from multiple undersampled waveforms. The optimal estimate of common remainder in robust CRT, which plays an important role in the final frequency estimation, is first discussed. To avoid the exhausted searching in the optimal estimation, we then provide an improved algorithm with the same performance but less computation. Besides, the sufficient and necessary condition of the robust estimation was proposed. Numerical examples are also provided to verify the effectiveness of the proposed algorithm and related conclusions.
Frequency estimation from noisy sinusoid is a hot research topic in engineering, which is applied in almost all engineering fields, such as communications, radar, astronomy, medicine, measurements, and geophysical exploration [1–6]. At present, many of the strategies are based on the Nyquist sampling theorem, for instance, discrete Fourier transform (DFT) algorithm [7–9], the least-error squares algorithm [10–13], and equally spaced samples algorithm . However, one of the main drawbacks of these methods is that the sampling frequency is usually proportional to the original frequency, which requires much higher speed analog-to-digital convertor. Recently, frequency estimation from undersampled waveforms has been studied in [15–17]. The basic idea of these methods is to fuse multiple remnant information provided by multiple sub-Nyquist sampled sequences. Then the original frequency can be reconstructed by using the Chinese remainder theorem (CRT) [18–20]. The biggest merit of these methods is that they greatly decrease the rigorous requirement for hardware; that is, a very high frequency can be determined by several lower rate sample devices. In addition, the method is also used in synthetic aperture radar [21, 22] and wireless sensor network .
The problem of frequency determination in undersampled condition is summarized as follows.
Assume that a single frequency signal is modelled as a complex sinusoid form: where is an unknown frequency we are interested in the signal and it may be high, is the signal amplitude, and . In practical applications, noises are unavoidable and thus lead to the received signal having the following expression: where denotes noise.
Let the sampling rates be , . Then we can get the corresponding sampled signal as follows: If the sampling duration is seconds, then
The -point DFT of each of the aforementioned sampled signals may provide an estimate of the original frequency. Becauseare much smaller than in practical applications, the -point DFT only provides a folded remainder of modulo . These remainders may be erroneous when the signal is contaminated by the additive noise. So, the above frequency estimation problem is equivalent to determining the original number from the erroneous remainders , .
This problem has been studied in  where a fast robust CRT algorithm is proposed. However, the number of searches is high and no closed-form solution is provided. It has been validated that this method has a poorer performance than the improved closed-form robust CRT (ICRT), which is proposed in . The latter method improved the performance because of the remainders’ proper choosing. However, the estimation of common remainder requires searching, which leads to a high computation complexity. Motivated by this, we propose a computational efficient algorithm which has the same performance as the ICRT method but with much less computation.
The remaining of this paper is organized as follows. In Section 2, we present the CRT and analyse its error sensitivity. In Section 3, we propose the robust CRT frequency estimation algorithm. Besides, we give a sufficient and necessary condition of the robust estimation based on the proposed algorithm. In Section 4, we apply the proposed algorithm for signal in additive noise with low sample rates. We end with concluding remarks in Section 5.
2. The Error Sensitivity of CRT
In this section, we first briefly describe the extended CRT where all the moduli have a common factor larger than 1. We then briefly describe the error sensitivity of CRT.
For convenience, we denote by the remainder of modulo .
Let be the moduli, and let , be the remainders of the positive integer ; that is, where . It is not hard to see that the original number can be uniquely reconstructed from its remainders if and only if is less than the least common multiple of the remainders. In the case when all the moduli are coprime, the CRT has a simple formula. In this paper, we consider the case when any paired moduli have the same greatest common remainder (gcd) and the remaining integers factorized by are coprime integers ; that is, It follows from (6) that One can see from (7) that all of the remainder moduli have the same value, which is called the common remainder, denoted by .
Let , and let . Then we obtain from (5) that According to the CRT, if and only if , can be uniquely reconstructed as where , is the modular multiplicative inverse of modulo . Therefore, can be uniquely reconstructed by
Now, we consider the case when the remainders have errors.
Let the th erroneous remainder be where denotes the error. In order to reconstruct the original number from its erroneous remainders , , we need to determine and according to (10). Suppose that the estimate of the common remainder is . Then can be estimated by where and stands for the rounding integer.
According to (12), if the error of the common remainder is constrained by then . Consequently, can be accurately reconstructed. It follows from (10) that the estimate of the original number is Obviously, the estimate error is .
Now, we consider condition (13), which leads to the correct estimate of . It can be rewritten as
Note that the fact that does not satisfy (15) will lead to the error estimate of by (12). As is usually a large integer, this may cause a large error of . For example, let , , , and . Then we have , , and . Suppose that the contaminated data are , , and . Then can be reconstructed as 2661. However, in the case when all errors are the same, the error of reconstruction is reduced significantly. Suppose that the estimate of is for . Then we have , which means that the estimate error of is . Therefore, the estimate of the common remainder is significant to the estimation.
3. Robust CRT Frequency Estimation Algorithm
In this section, we give the optimal estimate of the common remainder firstly; then we propose the efficient frequency estimation algorithm. Before introducing our method, we introduce a kind of circular distance function.
Definition 1. For real numbers and , the circular distance of to for cycle is defined as where and .
Putting then we can conclude that may be different from each other due to errors . The optimal estimate is the point which minimizes the summation of the deviation between and variable , where . Based on the circular distance defined above, we can obtain the summation of deviation as follows: Then, the optimal estimate can be determined by Next, we give the optimal estimate in the following theorem.
Theorem 2. The optimal estimate in (19) can be determined by where .
Proof. According to the definition of the circular distance, we have
Let . One can easily see that is continuous in field of real number. Now, we consider the differentiability of at which satisfies
For , where is a small positive number, we have . It follows that the left-sided derivative of at is .
For , we obtain . It follows that the right-sided derivative of at is . Thus, is not differentiable at .
Let . Then we can obtain where is the number of remainders satisfying (22). Thus, does not have a local minimum value at . Note that is differentiable at all , so the minimum point must satisfy , that is, in (19), such that After equating the result, we obtain Since and , it follows that belongs to Hence, the optimal estimate in (19) can be simplified as Combining (26) and (27), we can draw the conclusion.
Theorem 2 tells us that the estimate of common remainder belongs to the set . It is not difficult to find that the cardinality of the set equals the number of the remainders. So, the number of searches is reduced sharply when compared with other searching methods. After optimal estimate is determined, the original frequency can be reconstructed. For convenience, we conclude the proposed efficient algorithm in Algorithm 1.
Based on the proposed algorithm, we have a sufficient and necessary condition of the robust estimate as shown in Algorithm 1.
Theorem 3. Let . Assume that and . is a real number in the range . Then holds if and only if
Proof. We proof the sufficiency.
From (**), we have Since and , it follows that We make this into three cases as follows.
Case A. .
According to (30), we obtain . It follows from (12) that for . From (9), we have Here, we use the conclusion that . Consequently,
Case B. .
According to (30), we obtain . It follows from (12) that . By using CRT, we have . Consequently, .
Case C. .
According to (30), we obtain . It follows from (12) that . From (9), we have . Consequently, . This completes the proof of the sufficiency.
We then proof the necessity.
Combination of (10) and (****) yields This leads to Since , it follows that Note that both of and are integers; we have three possible cases: , , and .
Case A. .
By (34), we have . That is, .
Case B. .
By (34), we have .
Case C. .
By (34), we have . That is, .
Therefore, we have Combing (**) and (37), we can draw the conclusion.
4. Simulation Results
In the simulations, the coprime numbers from to are 5, 7, 11, 13, 19, 23, and 29, respectively. Two possibilities of are considered: and . The frequencies which are to be estimated is taken integers randomly and uniformly distributed in the range for each . The additive noise in (4) is additive white Gauss noise. The sampling frequencies are, which are much smaller than the unknown frequencies. The observation of the time duration is 1 s, and the number of the trials is 10000 for each signal-to-noise ratio (SNR).
We compare the proposed method with the ICRT algorithm which is the most effective one of the existing methods. To assess the performance, we take the trail-fail rate (TFR) as the measurement. In each trial, if the estimate of the frequency is then the trial passed; otherwise, the trail failed.
Figure 1 gives the curve of the TFR on different SNR. It shows that the probabilities of robust estimation are almost equal when the two methods have the same . That is to say, the two methods have nearly the same performance. The figure also shows that the performance of the estimation is better when than that when , which is in accordance with the theory. From (13) we can conclude that the larger can lead to the larger error tolerance.
Figure 2 presents the curve of the total running time over the number of remainders. The two methods are carried out under the same operational conditions. One can see that our proposed algorithm has nearly equal time no matter is, while the ICRT algorithm has essentially linear running time. This is because the number of searches for the common remainder is for the proposed algorithm. However, the number of searches is for the ICRT method, where is the searching step. In this simulation, the searching step is , and .
Figure 3 gives the robust estimation probability of the common remainder for searching methods when is 10, 1, and 0.1, respectively. It shows that the estimate of the common remainder can achieve an approximation of the optimal estimate whenis small enough. Clearly, decreasing the searching step can improve the performance at the cost of increasing the computation amount.
In this paper, we have proposed an efficient algorithm based on the robust CRT to perform single frequency estimation from multiple undersampled waveforms. We transformed the problem into the robust CRT firstly. Then we presented the optimal estimate of the common remainder and the original frequency. Based on the proposed algorithm, we presented the sufficient and necessary condition of the robust estimation. We finally applied the efficient algorithm to estimate the original frequency from multiple undersampled waveforms. Compared with the improved robust CRT algorithm, it has nearly the same performance but less computation. After using the undersampled waveforms, the sampling rate reduced sharply.
The authors would like to thank the referee for his (her) valuable comments and suggestions which improved this work to a great extent. This work was supported by the National Natural Science Foundation of China (Grants nos. 61172092, 61302069) and the Research Fund for the Doctoral Programs of Higher Education of China (no. 20130201110014).
B. G. Quinn and E. J. Hannan, The Estimation and Tracking of Frequency, Cambridge University Press, New York, NY, USA, 2001.
E. Lavopa, P. Zanchetta, M. Sumner, and F. Cupertino, “Real-time estimation of fundamental frequency and harmonics for active shunt power filters in aircraft electrical systems,” IEEE Transactions on Industrial Electronics, vol. 56, no. 8, pp. 2875–2884, 2009.View at: Publisher Site | Google Scholar
T. Grandke, “Interpolation algorithms for discrete fourier transforms of weighted signals,” IEEE Transactions on Instrumentation and Measurement, vol. 32, no. 2, pp. 350–355, 1983.View at: Google Scholar
B. Krishna, H. C. Krishna, and K. Y. Lin, Computational Number Theory and Digital Signal Processing: Fast Algorithms and Error Control Techniques, CRC, Boca Raton, Fla, USA, 1994.
C. Ding, D. Pei, and A. Salomaa, Chinese Remainder Theorem: Applications in Computing, Coding, Cryptography, World Scientific, Singapore, 1999.
K. H. Rosen, Elementary Number Theory and Its Applications, Addison-Wesley, Mass, 5th edition, 2010.
Z. H. Yuan, Y. k. Deng, F. Li et al., “Multichannel InSAR DEM reconstruction through improved closed-form robust Chinese remainder theorem,” IEEE Geoscience and Remote Sensing Letters, 2013.View at: Google Scholar