Next Article in Journal
The Uniform Convergence Property of Sequence of Fractal Interpolation Functions in Complicated Networks
Next Article in Special Issue
CNN-Based Temporal Video Segmentation Using a Nonlinear Hyperbolic PDE-Based Multi-Scale Analysis
Previous Article in Journal
The Relation Dimension in the Identification and Classification of Lexically Restricted Word Co-Occurrences in Text Corpora
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Alternative Method to Estimate the Fourier Expansions and Its Rate of Change

by
Johnny Rodríguez-Maldonado
,
Cornelio Posadas-Castillo
and
Ernesto Zambrano-Serrano
*
Facultad de Ingeniería Mecánica y Eléctrica, Universidad Autónoma de Nuevo León, San Nicolás de los Garza 66455, NL, Mexico
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(20), 3832; https://doi.org/10.3390/math10203832
Submission received: 2 September 2022 / Revised: 7 October 2022 / Accepted: 14 October 2022 / Published: 17 October 2022

Abstract

:
This paper presents a methodology to obtain the Fourier coefficients (FCs) and the derivative Fourier coefficients (DFCs) from an input signal. Based on the Taylor series that approximates the input signal into a trigonometric signal model through the Kalman filter, consequently, the signal’s and successive derivatives’ coefficients are obtained with the state prediction and the state matrix inverse. Compared to discrete Fourier transform (DFT), the new class of filters provides noise reduction and sidelobe suppression advantages. Additionally, the proposed Taylor–Kalman–Fourier algorithm (TKFA) achieves a null-flat frequency response around the frequency operation. Moreover, with the proposed TKFA method, the decrement in the inter-harmonic amplitude is more significant than that obtained with the Kalman–Fourier algorithm (KFA), and the neighborhood of the null-flat frequency is expanded. Finally, the approximation of the input signal and its derivative can be performed with a sum of functions related to the estimated coefficients and their respective harmonics.
MSC:
65T40; 42A16; 41A58

1. Introduction

Signal analysis is implemented in many areas because it allows for obtaining characteristics or relevant information about systems. It has been employed in different areas, such as power systems, vibration analysis, speech recognition, energy processing, radar applications, time series modeling, biomedical engineering, and digital communications. For instance, [1] gave the generalized sampling expansions (GSE) in the realm of a one-dimensional quaternion Fourier transform. Ref. [2] employed the Fourier coefficients (FCs) to fit the GPS precise ephemeris. In [3,4], the Fourier series coefficients were used in a tracking control trajectory. In the both of the last cases, obtaining more information could help better control GPS and tracking control, not only regarding the trajectory, but also its speed.
The most commonly used tool to obtain information about systems is the Fourier transform. In the literature, it is possible to find some strategies to obtain better precision in the estimation of the FCs in different conditions. In [5], the discrete Fourier transform (DFT) is used to estimate the periodic signal’s FCs. On the other hand, the impact of the data windows in the estimation of FCs was analyzed by [6]. The recursive discrete Fourier transform (RDFT) [7,8,9] is another strategy that considers a structure of the DFT with a multi-rate sampling. One disadvantage of the last method is that it takes a sinusoidal signal as a model. Therefore, the estimates can be downgraded when the input is different from the sinusoidal signal, for example, a quasi-periodic signal. Another disadvantage is that obtained estimates provided a delay due to the implemented windows in the algorithms. To solve the first disadvantage, when the measurement signal does not satisfy the ideal conditions (when the input signal is not sinusoidal), this occurs when the measurement signal is quasi-periodic. In [10,11], the authors implemented a quasi-periodic signal model considering an optimal FIR filter, whereas in [12], the authors improved the approximation of a sufficiently smooth nonperiodic function defined on a compact interval by proposing alternative forms of Fourier series expansions. In [13], an algorithm was proposed that was derived from the least mean squares (LMS) that can estimate the FCs of a sinusoidal and quasi-periodic signal. Then, to solve the delay problem in the estimates, in [14,15], the FCs were obtained with a Kalman filter.
Systems can be modified through time by different circumstances. For example, the evolution in the systems can be due to wear or operation conditions. For this reason, is not enough to estimate the FCs. The derivative of the Fourier coefficients (DFCs) provides information on the evolution of systems over time. In this work, the FCs are estimated without delay because the estimations are achieved with a quasi-periodic dynamic signal model applied to the Kalman filter. Additionally, at the same time, the DFCs are estimated. With the Taylor polynomial approximation to the signal model and its derivatives, it can develop a state space model, where the first row in the state matrix is the approximation to the signal model, the second is the approximation to the first derivative, and so on. Thereby, with the state space model implemented in the Kalman filter, it is possible to obtain the FC and DFC estimates in synchronous time. Additionally, the harmonics and their derivatives can be estimated with an extension of the state matrix with its harmonics. Recently, some techniques have been proposed to estimate the FCs. Ref. [16] determined the coefficient in the representation of two polynomials as a linear combination of an arbitrary polynomial sequence to achieve the above derivative of the Fourier series. Ref. [17] considered unbiased Fourier coefficients for reducing systematics in applications of large cosmological data sets in grid size. In [18], the authors proposed an algorithm for the identification of significant frequencies and the estimation of the Fourier coefficients. However, they did not estimate the DFCs. In [19], the authors employed the Fourier cosine series expansion (COS) to value the guaranteed minimum death benefit products.
The aim of the proposed method is based on two stages. In the first stage, with the Taylor polynomial implemented in the Kalman filter, the signal and its derivatives are estimates. In the second stage, the coefficients and their derivatives are calculated by considering the estimations obtained from the Kalman filter and the state-proposed matrix. Here we observe that when a zero-order Taylor polynomial is employed as an analytic function in the Kalman filter, the results are equivalent to the traditional Fourier series. With the expansion of the Taylor polynomial with an order greater than zero, the estimation results of the Fourier coefficients are improved. Moreover, it can be confirmed in the magnitude response of the proposed method, where the magnitude in harmonic frequencies is flat.
The remainder of this paper is organized as follows. First, the state space observer approach is described in Section 2. The Kalman filter applied and the structure implemented to obtain the coefficients are illustrated in Section 3. The performance of the proposed method is shown in Section 4. A summary of our conclusions closes the paper in Section 5.

2. The Spectral Observer Approach

It is common to implement the Fourier series as a signal model [20,21] that is defined as follows
z k = a k 0 + = 1 N a k c o s 2 π k 2 N + 1 + b k s i n 2 π k 2 N + 1 ,
where the number of samples per period T is 2 N + 1 , and represent the harmonic components. However, here, we propose a quasi-periodic signal inspired by (1). In its first stage, it is approximate around t 0 by a M-th order Taylor polynomial as
ζ k = m = 0 M f 2 π k N 1 ( m ) m ! τ ( m ) ,
with t 0 = [ n 1 ] τ at time t = n τ , and  τ = t t 0 . The function f ( · ) in (2) is defined as
f 2 π k N 1 = c o s 2 π k N 1 + s i n 2 π k N 1 .
Hence, the real part of the function in (2) is denoted as
Δ r , k = m = 0 M a m c o s 2 π k N 1 ( m ) m ! τ ( m ) ,
and the imaginary part is
Δ i , k = m = 0 M b m s i n 2 π k N 1 ( m ) m ! τ ( m ) ,
where N 1 is the number of samples per cycle in the input signal, and M is the order of the polynomial approximation.
With (4) and its successive derivatives, the following matrix is obtained:
Λ r = α 1 , 1 α 1 , 2 α 1 , m α 2 , 1 α 2 , 2 α 2 , m α m , 1 α m , 2 α m , m c o s ( ω τ ) ,
where the coefficients α i , j with, i = 1 , , m and j = 1 , , m are given as follows:
α 1 , 1 = 1 , α 1 , 2 = τ , α 1 , m = τ m m ! ;
α 2 , 1 = j ω , α 2 , 2 = ( j ω τ + 1 ) , α 2 , m = j ω τ m m ! + τ m 1 ( m 1 ) ! ;
α m , 1 = ( j ω ) m ;
α m , 2 = ( j ω ) m τ + m ( j ω ) m 1 ;
α m , m = ( j ω ) m τ m m ! + m ( j ω ) m 1 τ m 1 ( K 1 ) ! + + 1 ;
with ω = 2 π N 1 and τ = 1 N 1 . On the other hand, the imaginary part (5) and its derivatives can be represented as
Λ i = j α 1 , 1 α 1 , 2 α 1 , m α 2 , 1 α 2 , 2 α 2 , m α m , 1 α m , 2 α m , m s i n ( ω τ ) .
Then, the signal model defined by (2) and its derivative is given as
ζ k ζ ¯ k = Λ r Λ i Λ i Λ r a b ,
where ζ ¯ k is the complex conjugate of ζ k . Therefore,
a b = Λ r Λ i Λ i Λ r 1 ζ k ζ ¯ k .
Note that a = [ a a 1 a 2 a m ] T , b = [ b b 1 b 2 b m ] T , and  ζ k , = [ ζ k ζ ˙ k ζ ¨ k ζ k ( m ) ] . As the coefficients a and b correspond to ζ = [ ζ k , ζ k ¯ ] T , consequently, the coefficients a 1 and b 1 correspond to ζ k ˙ ; therefore, its coefficients are called DFCs.
In the next section, the state transition matrix is extended to all sets of harmonics and is implemented to extract the coefficients.

3. The Kalman Filter Applied and Coefficient Estimates

This section illustrates the construction of the signal model and its implementation in the Kalman filter, where two methodologies are considered. The first method is when the signal model is approximate with a zero-order Taylor polynomial defined in (2). The second method is obtained when the signal model is approximate with a higher-than-zero-order.

3.1. Kalman Filter Algorithm (KFA)

With the matrix (6) and (14), we developed a signal model, as can be seen in (8), that is implemented in the Kalman filter. Ehen the order of the signal model is defined as m = 1 and l = 1 , the algorithm is called a Kalman filter algorithm (KFA). Thus, the state prediction and the measurement equation are defined as.
x k + 1 = Ψ x k + G w k ,
ζ k = H x k + 1 ,
where
Ψ = Λ r , Λ i , Λ i , Λ r ,
h = [ 1 0 , , 0 ] ( 2 ( K + 1 ) ) ,

3.2. Taylor–Kalman–Fourier Algorithm (TKFA)

The matrix (15) can be expanded with the array of Δ r and Δ i , where the functions c o s and s i n have the respective harmonics as
Λ r , = α 1 , 1 α 1 , 2 α 1 , m α 2 , 1 α 2 , 2 α 2 , m α m , 1 α m , 2 α m , m c o s ( ω τ ) .
and
Λ i , = j α 1 , 1 α 1 , 2 α 1 , m α 2 , 1 α 2 , 2 α 2 , m α m , 1 α m , 2 α m , m s i n ( ω τ ) .
with = 1 , , N .
On the other hand, when the transition matrix is expanded with > 1 , and the approximation to the model signal is greater than one ( m > 1 ), the filter is named a Taylor–Kalman–Fourier algorithm (TKFA), and the state transition matrix and the measurement equation are defined as
Ψ = B l o c k   d i a g Λ r , Λ i , Λ i , Λ r , ; = 1 , , N } ,
H = [ h h , , h ] N ( 2 ( K + 1 ) ) ,
Remark 1.
When the parameters m and ℓ of matrix Ψ in (15) are set as m = 1 and = 1 , then the algorithm is called a KFA.
With the discrete state space model (10a) and the truncated discrete signal model (10b), the standard discrete Kalman filter [22,23] is applied to obtain the state estimation x ^ k + 1 according to Algorithm 1.
Algorithm 1 Pseudocode of the Kalman filter algorithm.
1: procedure
2:       –Inputs
3:        s [ n ] Input current signal for  n = 1 up to ( 2 N + 1 ) ,
4:        x ^ [ n 1 ] Guess initial internal variables,
5:        ω ^ [ 0 ] Guess initial frequency,
6:        P [ n 1 ] Guess initial covariance,
7:        σ v 2 Variance of process noise,
8:        σ w 2 Variance of measurement noise,
9:        Ψ , H Signal model parameters,
10:       for  n = 1 to ( 2 N + 1 )  do
11:            –State prediction
12:             x ^ [ n ] Ψ x ^ [ n 1 ] ,
13:             P [ n ] Ψ P [ n 1 ] ( Ψ ) H + Γ Γ H σ v 2 ,
14:            –Measurement update
15:             K [ n ] P [ n ] H T ( H P [ n ] H T + σ w 2 ) 1 ,
16:             x ^ [ n ] x ^ [ n ] + K [ n ] ( s [ n ] H x ^ [ n ] ) ,
17:             P [ n ] ( I K [ n ] H ) P [ n ] ,
18:      end for
19: end procedure
However, with the last Kalman filter algorithm, it is not possible to obtain the coefficients a and b . Therefore, in the next subsection, the equation to obtain the coefficients and the approximation to the signal and its derivatives are illustrated.

3.3. Coefficients and Signal Approximation

With the state update implemented in the next equation, it is possible to obtain the coefficients and their derivatives
a b = Ψ 1 x ^ , x ^ C 1 × ( 2 ( K + 1 ) ) .
With the estimated coefficients a and b within a and b , respectively, the input signal and its derivatives are approximated as
ζ ^ ( 0 ) = = 1 N ζ ^ ( 0 ) [ 1 ] + a c o s ω k + b s i n ( ω k ) ; ζ ^ ( 1 ) = = 1 N ζ ^ ( 1 ) [ 1 ] +
a ˙ ( ω ) s i n ω k + b ˙ ( ω ) c o s ( ω k ) τ ; ζ ^ ( 2 ) = = 1 N ζ ^ ( 2 ) [ 1 ] +
a ¨ ( ω ) 2 c o s ω k b ¨ ( ω ) 2 s i n ( ω k ) τ 2 ;
ζ ^ ( N )
where the initial conditions ζ ^ ( 0 ) [ 0 ] , ζ ^ ( 1 ) [ 0 ] , , ζ ^ ( N ) [ 0 ] are defined as equal to zero. Algorithm 2 illustrates the implementation of Equation (18a)–(18c) to obtain the coefficients and its derivatives, as well as the approximation of the input signal and its respective derivative.
All set algorithm is illustrated in Figure 1, where Figure 1 (a) is the Kalman algorithm implemented as can be seen in Algorithm 1, where the input into the algorithm is the measurement signal, and it is possible to implement the model signals of Section 3.1 and Section 3.2. Note that the output is the states.
Algorithm 2 Pseudocode to coefficient estimates and signal approximation.
1: procedure
2:       –Inputs
3:        ζ ^ ( 0 ) 0, Initial value,
4:        ζ ^ ( 1 ) 0, Initial value,
5:       for  n = 1 to  do,
6:              [ c 1 , c 2 ] T Ψ 1 x ^ ,
7:              [ a , b ] T | [ c 1 , c 2 ] T | ,
8:              ζ ^ ( 0 ) , Compute (18a),
9:              ζ ^ ( 1 ) , Compute (18b),
10:             ζ ^ ( 2 ) , Compute (18c),
11:       end for
12: end procedure
The output of the Kalman algorithm is the input of a second stage, as seen in Figure 1 (b), which is divided into two subsections, Figure 1 (c) and (d). The coefficients are obtained by the block in Figure 1 (c) with the state matrix Ψ 1 and x ^ (as is illustrated in (17). Finally, with the coefficient estimates a and b implemented in Equation (18a)–(18c) by the block of Figure 1 (d), the approximations to the signal ζ ^ ( 0 ) and its derivatives ζ ^ ( 1 ) , ζ ^ ( 2 ) , , ζ ^ ( M ) , are obtained and are the output of the block diagram. Note that the last block in Figure 1 represents the procedure code that is shown in Algorithm 2, and thus the process is completed.

4. Performance of the Harmonic Filters

In this section, the performance of the proposed algorithm is shown. When the signal model is approximated with M = 0 in the polynomial form, defined in (2), the signal model is equivalent to the model proposed in [20] and here is named a Kalman filter algorithm (KFA). When the order signal approximation in the model is M > 0 and > 1 , the proposed algorithm is named a Taylor–Kalman–Fourier algorithm (TKFA). On the other hand, the filters are tested with a square wave signal as an input and are compared with DFT, and KFA. The results are illustrated in Figure 2, where the behavior between DFT, KFA and TKFA are very similar. Additionally, the zoom shows in detail the behavior of the estimates, and it can be seen that the precision of the methods is equivalent. To obtain more precise results in the comparison of the algorithms, the RMS Error defined by
R M S E = n = 1 N s [ n ] ζ ^ [ n ] N 2 ,
is obtained with respect to the measurement input signal s [ n ] and its estimate ζ [ n ] . The obtained results are shown in Table 1 for a different number of harmonics in the signal model by KFA and TKFA. The number of harmonics implemented in the model is represented by ω . Hence, ω f represents the fundamental model, and therefore, the transition matrix Ψ is developed only with = 1 . Note that the results obtained by the three algorithms are equivalent. However, with the DFT and KFA, the DFCs cannot be obtained. Thus, the results of the RMS error in Table 1 are obtained only using the coefficients a and b.
A second signal is tested and is illustrated in Figure 3. Additionally, its approximation is obtained with DFT, KFA, and TKFA. The zoom shows in detail the behavior of the estimates, and it can be seen that the precision of the methods is equivalent between DFT and KFA, but TKFA has a litter discrepancy in the estimation provided by the signal model, where the harmonics are included.
We propose to select this signal because its derivative is intuitive and easy to obtain. On the other hand, with the coefficients a 1 and b 1 obtained with (17) and implemented in (18b), we achieved approximation to the derivative of the signal, and the results are presented in Figure 4 where the path trajectory is the expected result because the derivative of a sawtooth signal is a square signal. The zoom shows in detail the behavior of the estimates. As the coefficients a 1 and b 1 are not possible to obtain with DFT and KFA, the approximation to the derivative ζ ^ ( 1 ) in Figure 4 only represents the estimates with TKFA.
With the proposal to make an analysis of the noise effect, Figure 5 shows the frequency responses of DFT, KFA and TKFA. The zoom shows in detail the behavior of the estimates. Using KFA and TKFA, the frequency response can be obtained using
G ( z ) = I + K H I Ψ z 1 1 K .
The results illustrated in Figure 5 were developed using the next parameters in the KFA and TKFA, Q = 10 3 , R = 1 and P r = 10 10 .
In Figure 6 the frequency response is illustrated with Q = 10 1 , P r = 10 10 and R = 1 . The zoom shows in detail the behavior of the estimates.
The comparison between Figure 5 and Figure 6 show how the inter-harmonic amplitude decreases in KFA and TKFA when the value in Q decreases. However, the decrement in the inter-harmonic amplitude is more significant in TKFA. A similar study of the noise was made in [20], where a complementary study was made of the spectral and short-time analysis. As TKFA achieves a null-flat frequency response around the frequency operation and its harmonics, the proposed method can achieve better results in applications where frequency and harmonics have a few distortions. For example, when the rotor bar has litter damage, the current has a frequency shift [24].

5. Conclusions

The aim of this work is the development of a method to estimate the coefficients of an input signal and the coefficients of the derivative input signal. To achieve the coefficient estimates, we develop a method that consists of two stages. First, a proposed signal model is implemented in the Kalman filter, and with this model, we are able to estimate the derivatives of the input signal. Second, with the state prediction and the inverse of the state matrix, the coefficients of the signal and the coefficients of the successive derivatives are obtained. On the other hand, in comparison to the DFT, the new class of filters provides advantages in the sense of noise reduction and the sidelobe suppression. Additionally, the proposed TKFA achieves a null-flat frequency response around the frequency operation. Moreover, with the proposed TKFA method, the decrement in the inter-harmonic amplitude is greater than that obtained with KFA, and the neighborhood of the null-flat frequency is expanded. The approximation to the input signal and its derivative can be achieved with a sum of functions that are related by the estimate coefficients and their respective harmonics. Thereby, the proposed method is related to the Fourier series. However, with the new filters, it is possible to estimate the derivative Fourier series with the derivative coefficient estimates.

Author Contributions

Conceptualization, J.R.-M., C.P.-C. and E.Z.-S.; methodology, J.R.-M., C.P.-C. and E.Z.-S.; software, J.R.-M., C.P.-C. and E.Z.-S.; validation, J.R.-M., C.P.-C. and E.Z.-S.; formal analysis, J.R.-M., C.P.-C. and E.Z.-S.; investigation, J.R.-M., C.P.-C. and E.Z.-S.; writing—original draft preparation, J.R.-M., C.P.-C. and E.Z.-S.; writing—review and editing, J.R.-M., C.P.-C. and E.Z.-S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by CONACyT-Mexico grant number 166654 & A1-S-31628.

Data Availability Statement

Not applicable.

Acknowledgments

All the authors acknowledge “Departamento de Electrónica y Automatización” from Facultad de Ingeniería Mecánica y Eléctrica at UANL and CONACyT-Mexico (research grant no. 166654 & A1-S-31628). E. Zambrano-Serrano and C. Posadas-Castillo acknowledge UANL (research grants 596-IT-2022 and 492-IT-2022).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Saima, S.; Li, B.; Adnan, S.M. New Sampling Expansion Related to Derivatives in Quaternion Fourier Transform Domain. Mathematics 2022, 10, 1217. [Google Scholar] [CrossRef]
  2. Ruibo, J.; Xiaoqiang, L. Using Fourier series to fit The GPS precise ephemeris. In Proceedings of the 2010 International Conference on Computer Application and System Modeling (ICCASM 2010), Taiyuan, China, 22–24 October 2010; Volume 11, pp. V11-96–V11-99. [Google Scholar] [CrossRef]
  3. Labiod, S.; Boubertakh, H.; Guerra, T.M. Fourier series-based adaptive tracking control for robot manipulators. In Proceedings of the 3rd International Conference on Systems and Control, Algiers, Algeria, 29–31 October 2013; pp. 968–972. [Google Scholar] [CrossRef]
  4. Gholipour, R.; Fateh, M.M. Adaptive task-space control of robot manipulators using the Fourier series expansion without task-space velocity measurements. Measurement 2018, 123, 285–292. [Google Scholar] [CrossRef]
  5. Van den Bos, A. Best linear unbiased estimation of the Fourier coefficients of periodic signals. IEEE Trans. Instrum. Meas. 1993, 42, 49–51. [Google Scholar] [CrossRef] [Green Version]
  6. Harris, F.J. On the use of windows for harmonic analysis with the discrete Fourier transform. Proc. IEEE 1978, 66, 51–83. [Google Scholar] [CrossRef]
  7. Hostetter, G. Recursive discrete Fourier transformation. IEEE Trans. Acoust. Speech Signal Process 1980, 28, 184–190. [Google Scholar] [CrossRef]
  8. Varkonyi-Koczy, A.R. A recursive fast Fourier transformation algorithm. IEEE Trans. Circuits Syst. II Analog Digit. Signal Process 1995, 42, 614–616. [Google Scholar] [CrossRef]
  9. Fan, C.P.; Su, G.A. Novel recursive discrete Fourier transform with compact architecture. In Proceedings of the 2004 IEEE Asia-Pacific Conference on Circuits and Systems, Tainan, Taiwan, 6–9 December 2004; Volume 2, pp. 1081–1084. [Google Scholar] [CrossRef]
  10. Park, S.H.; Kim, M.J.; Kwon, W.H.; Kwon, O.K. Short-time harmonic analysis via the state-space optimal FIR filter. In Proceedings of the SICE ’95, 34th SICE Annual Conference, Hokkaido, Japan, 26–28 July 1995; International Session Papers. pp. 1205–1210. [Google Scholar] [CrossRef]
  11. Park, S.H.; Kwon, W.H.; Kwon, O.K.; Kim, M.J. Short-time Fourier analysis via optimal harmonic FIR filters. IEEE Trans. Signal Process 1997, 45, 1535–1542. [Google Scholar] [CrossRef]
  12. Li, W. Alternative Fourier Series Expansions with Accelerated Convergence. Appl. Math. 2016, 7, 1824–1845. [Google Scholar] [CrossRef] [Green Version]
  13. Xiao, Y.; Tadokoro, Y.; Shida, K. Adaptive algorithm based on least mean p-power error criterion for Fourier analysis in additive noise. IEEE Trans. Signal Process 1999, 47, 1172–1181. [Google Scholar] [CrossRef]
  14. Nishi, K. Kalman filter analysis for quasi-periodic signals. In Proceedings of the 2001 IEEE International Conference on Acoustics, Speech, and Signal Processing, Proceedings (Cat. No.01CH37221), Salt Lake City, UT, USA, 7–11 May 2001; Volume 6, pp. 3937–3940. [Google Scholar] [CrossRef]
  15. Gruber, P.; Todtli, J. Estimation of quasiperiodic signal parameters by means of dynamic signal models. IEEE Trans. Signal Process 1994, 42, 552–562. [Google Scholar] [CrossRef]
  16. Kim, T.; Kim, D.; Jang, L.C.; Jang, G.W. Fourier Series for Functions Related to Chebyshev Polynomials of the First Kind and Lucas Polynomials. Mathematics 2018, 6, 276. [Google Scholar] [CrossRef] [Green Version]
  17. Sefusatti, E.; Crocce, M.; Scoccimarro, R.; Couchman, H.M.P. Accurate estimators of correlation functions in Fourier space. Mon. Not. R. Astron. Soc. 2016, 460, 3624–3636. [Google Scholar] [CrossRef] [Green Version]
  18. Gilbert, A.C.; Guha, S.; Indyk, P.; Muthukrishnan, S.; Strauss, M. Near-Optimal Sparse Fourier Representations via Sampling. In Proceedings of the Annual ACM Symposium on Theory of Computing, Montreal, QC, Canada, 19–21 May 2002; pp. 152–161. [Google Scholar] [CrossRef]
  19. Yu, W.; Yong, Y.; Guan, G.; Huang, Y.; Su, W.; Cui, C. Valuing Guaranteed Minimum Death Benefits by Cosine Series Expansion. Mathematics 2019, 7, 835. [Google Scholar] [CrossRef] [Green Version]
  20. Bitmead, R.; Tsoi, A.; Parker, P. A Kalman filtering approach to short-time Fourier analysis. IEEE Trans. Acoust. Speech Signal Process 1986, 34, 1493–1501. [Google Scholar] [CrossRef]
  21. Zhang, B.; Pong, M.H. Dynamic model and small signal analysis based on the extended describing function and Fourier series of a novel AM ZVS direct coupling DC/DC converter. In Proceedings of the Power Processing and Electronic Specialists Conference, Atlantic City, NJ, USA, 22–23 May 1997; Volume 1, pp. 447–452. [Google Scholar] [CrossRef]
  22. Manolakis, D.G.; Ingle, V.K.; Kogon, S.M. Statistical and Adaptive Signal Processing: Spectral Estimation, Signal Modeling, Adaptive Filtering and Array Processing; Artech House Publishers: Boston, MA, USA, 2005. [Google Scholar]
  23. Simon, D. Optimal State Estimation: Kalman, H Infinity, and Nonlinear Approaches; John Wiley & Sons: Hoboken, NJ, USA, 2006. [Google Scholar]
  24. Trujillo-Guajardo, L.A.; Rodriguez-Maldonado, J.; Moonem, M.A.; Platas-Garza, M.A. A Multiresolution Taylor–Kalman Approach for Broken Rotor Bar Detection in Cage Induction Motors. IEEE Trans. Instrum. Meas. 2018, 67, 1317–1328. [Google Scholar] [CrossRef]
Figure 1. Block diagram describing the proposed methodology to estimate the FCs, DFCs, signal approximation and its derivatives. In (a), the Kalman filter is implemented, (b) contains two subsections, (c) and (d). In (c) the coefficients with the state estimates provided by Kalman filter and with the state matrix Ψ 1 are obtained. Finally in block (d) Equation (18a)–(18d) is implemented to obtain the approximation signal.
Figure 1. Block diagram describing the proposed methodology to estimate the FCs, DFCs, signal approximation and its derivatives. In (a), the Kalman filter is implemented, (b) contains two subsections, (c) and (d). In (c) the coefficients with the state estimates provided by Kalman filter and with the state matrix Ψ 1 are obtained. Finally in block (d) Equation (18a)–(18d) is implemented to obtain the approximation signal.
Mathematics 10 03832 g001
Figure 2. Reconstruction signal with DFT, KFA, and TKFA algorithms. The model considers 15 harmonic components.
Figure 2. Reconstruction signal with DFT, KFA, and TKFA algorithms. The model considers 15 harmonic components.
Mathematics 10 03832 g002
Figure 3. Reconstruction signal with DFT, KFA and TKFA algorithms. In this case, the model of KFA and TKFA was constructed with the full set of harmonics components (64 harmonics).
Figure 3. Reconstruction signal with DFT, KFA and TKFA algorithms. In this case, the model of KFA and TKFA was constructed with the full set of harmonics components (64 harmonics).
Mathematics 10 03832 g003
Figure 4. Derivative approximation to the input signal with TKFA algorithm.
Figure 4. Derivative approximation to the input signal with TKFA algorithm.
Mathematics 10 03832 g004
Figure 5. Magnitude response of DFT, KFA and TKFA algorithms, with Q = 10 3 , P r = 10 10 and R = 1 .
Figure 5. Magnitude response of DFT, KFA and TKFA algorithms, with Q = 10 3 , P r = 10 10 and R = 1 .
Mathematics 10 03832 g005
Figure 6. Magnitude response of DFT, KFA and TKFA algorithms with Q = 10 1 , P r = 10 10 and R = 1 .
Figure 6. Magnitude response of DFT, KFA and TKFA algorithms with Q = 10 1 , P r = 10 10 and R = 1 .
Mathematics 10 03832 g006
Table 1. RMS Error (RMSE) obtained with DFT, KFA and TKFA with respect to the input signal and the approximation to the measurement signal ζ ^ ( 0 ) .
Table 1. RMS Error (RMSE) obtained with DFT, KFA and TKFA with respect to the input signal and the approximation to the measurement signal ζ ^ ( 0 ) .
RMSE
ω DFTKFATKFA ω DFTKFATKFA
Square signal ω f 0.430 0.430 0.430 ω 8 0.187 0.187 0.187
ω 1 0.430 0.430 0.430 ω 9 0.187 0.187 0.187
ω 2 0.308 0.308 0.308 ω 10 0.167 0.167 0.167
ω 3 0.308 0.308 0.308 ω 11 0.167 0.167 0.167
ω 4 0.249 0.249 0.249 ω 12 0.151 0.151 0.151
ω 5 0.249 0.249 0.249 ω 13 0.151 0.151 0.151
ω 6 0.213 0.213 0.213 ω 14 0.137 0.137 0.137
ω 7 0.213 0.213 0.213 ω 15 0.137 0.137 0.137
Sawtooth Signal ω f 0.069 0.069 0.069 ω 8 85 × 10 4 85 × 10 4 85 × 10 4
ω 1 0.069 0.069 0.069 ω 9 85 × 10 4 85 × 10 4 85 × 10 4
ω 2 0.028 0.028 0.028 ω 10 68 × 10 4 68 × 10 4 68 × 10 4
ω 3 0.028 0.028 0.028 ω 11 68 × 10 4 68 × 10 4 68 × 10 4
ω 4 0.016 0.016 0.016 ω 12 57 × 10 4 57 × 10 4 57 × 10 4
ω 5 0.016 0.016 0.016 ω 13 57 × 10 4 57 × 10 4 57 × 10 4
ω 6 0.011 0.011 0.011 ω 14 49 × 10 4 49 × 10 4 49 × 10 4
ω 7 0.011 0.011 0.011 ω 15 49 × 10 4 49 × 10 4 49 × 10 4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rodríguez-Maldonado, J.; Posadas-Castillo, C.; Zambrano-Serrano, E. Alternative Method to Estimate the Fourier Expansions and Its Rate of Change. Mathematics 2022, 10, 3832. https://doi.org/10.3390/math10203832

AMA Style

Rodríguez-Maldonado J, Posadas-Castillo C, Zambrano-Serrano E. Alternative Method to Estimate the Fourier Expansions and Its Rate of Change. Mathematics. 2022; 10(20):3832. https://doi.org/10.3390/math10203832

Chicago/Turabian Style

Rodríguez-Maldonado, Johnny, Cornelio Posadas-Castillo, and Ernesto Zambrano-Serrano. 2022. "Alternative Method to Estimate the Fourier Expansions and Its Rate of Change" Mathematics 10, no. 20: 3832. https://doi.org/10.3390/math10203832

APA Style

Rodríguez-Maldonado, J., Posadas-Castillo, C., & Zambrano-Serrano, E. (2022). Alternative Method to Estimate the Fourier Expansions and Its Rate of Change. Mathematics, 10(20), 3832. https://doi.org/10.3390/math10203832

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