Fast and accurate estimation and adjustment of a local
|
|
- Belinda Hawkins
- 6 years ago
- Views:
Transcription
1 86 ieee transactions on ultrasonics ferroelectrics and frequency control vol. 53 no. 5 may 006 An Unbiased FIR Filter for TIE Model of a Local Clock in Applications to GPS-Based Timekeeping Yuriy S. Shmaliy Senior Member IEEE Abstract An unbiased finite impulse response (FIR) filter is proposed to estimate the time-interval error (TIE) K-degree polynomial model of a local clock in global positioning system (GPS)-based timekeeping in the presence of noise that is not obligatory Gaussian. Generic coefficients for the unbiased FIRs are derived. The low-degree FIRs and noise power gains are given. An estimation algorithm is proposed and examined for the TIE measurements of a crystal clock in the presence of the uniformly distributed sawtooth noise induced by the multichannel GPS timing receiver. Based upon this algorithm we show that the unbiased FIR estimates are consistent with the reference (rubidium) measurements and fit them better than the standard Kalman filter. I. Introduction Fast and accurate estimation and adjustment of a local clock performance making possible for a variety of modern digital systems to operate in common time with minimum slips is of importance for the global positioning system (GPS)-based timekeeping. Here the time interval error (TIE) between the GPS time and the local clock time is measured in the presence of noise induced by the GPS timing receiver. The TIE then is estimated and the correction signal adjusts the clock for the GPS time. The standard deviation of the noise using commercially available receivers is about 30 ns can reach 10 0 ns [1] and may be improved by removal of systematic errors to no less than 3 5 ns [1] []. Having such a large noise the measured data is usually not appropriate for clock correction and a TIE tracking filter with a large time constant is applied. To obtain filtering in an optimum way a TIE model of a local clock must be known for the filter memory. Such a model was proposed in [3] as the second degree Taylor polynomial and is now basic in timekeeping being practically proven. In the discrete time it may be written as: x 1 (n) =x 1 (0) + x (0)τn + x 3(0) τ n + w 1 (n τ) (1) Manuscript received April 6 005; accepted May This material was discussed December at the Precise Time and Time Interval (PTTI) Meeting Washington DC. The work was supported by the CONACyT Project J38818-A. Y. S. Shmaliy is with the Guanajuato University Salamanca Gto. Mexico ( shmaliy@salamanca.ugto.mx). He also is with the Kharkiv National University of Radio Electronics Kharkiv Ukraine. Fig. 1. Typical TIE 1 PPS sawtooth noise induced by a GPS timing sensor SynPaQ III. where n = ; τ = t n t n 1 is a time step multiple to the 1 s; t n is a discrete time; x 1 (0) is an initial time error; x (0) is an initial fractional frequency offset of a local clock from the reference frequency; x 3 (0) is an initial linear fractional frequency drift rate; and w 1 (n τ) is a random component caused by the oscillator noise and environment. In GPS-based measurements (1) is observed via the mixture: ξ 1 (n) =x 1 (n)+v 1 (n) () in which v 1 (n) is a noisy component induced at the receiver (noise of a measurement set is usually small). In modern receivers such as the Motorola M1+ (see [4]) and SynPaQ III GPS Sensor (Synergy Systems LLC San Diego CA) a random variable v 1 (n) is uniformly distributed owing to the sawtooth noise caused by a principle of the 1 PPS (one pulse per second) signal formation. Fig. 1 shows a typical structure of v 1 (n) in a short time; that is the modulo v max Brownian TIE associated with a phase of a receiver local oscillator where v max is the noise bound (in SynPaQ III it is about 50 ns. If the sawtooth correction is used the noise is reduced by the factor of about 5 and becomes near Gaussian). In long-term measurements v 1 (n) exhibits nonstationary excursions due to uncertainties in the GPS time [5] with different satellite sets in a view. To estimate optimally the states of different clocks (atomic and crystals) several filtering algorithms have been examined during a couple of decades. For the state space equivalent of (1) the problem was formulated by Allan and Barnes in [6] to apply Kalman filtering. The solutions then were given in [7] [8]. Thereafter the Kalman algorithm was used by many authors [9] [13] and its applications in time scales were recently outlined in [14]. The problem with the Kalman filter arises when noise is not white and the model demonstrates temporary uncertainties; thus the estimate may become biased and noisy. To overcome advanced Kalman algorithms were proposed in [15] [16] being however not yet adopted for GPS applications /$0.00 c 006 IEEE
2 shmaliy: proposed finite impulse response filter 863 An alternative approach is known as finite impulse response (FIR) filtering allowing for noise of arbitrary distribution. In contrast to the infinite impulse response (IIR) structures (including Kalman filters) FIR structures have inherent properties such as a bounded input/bounded output (BIBO) stability and robustness against temporary model uncertainties and round-off errors [17]. They may be used independently or combined with Kalman filters [1] [18]. A general optimal FIR filtering algorithm with embedded unbiasedness for state space models was recently proposed in [19]. Especially for GPS-based timekeeping and a linear TIE polynomial model an unbiased FIR filter was designed and studied in [0]. In this paper we propose a new unbiased FIR filter and algorithm intended for real time estimating the TIE K- degree polynomial model in the presence of a GPS noise of arbitrary distribution (with or without the sawtooth correction). The rest of the paper is organized as follows. In Section II we formulate the problem. In Section III a design of an unbiased FIR filter is given (the necessary coefficients are postponed in Appendix I). The low-degree FIRs and noise-power gains are derived in Section IV. An unbiased FIR-filtering algorithm for a single multichannel GPS timing receiver is described in Section V and applied for a local crystal clock without the sawtooth correction. The FIR estimates are compared here to the reference measurement (rubidium) and to those obtained using the 3-state Kalman filter (Appendix II). Concluding remarks are drawn in Section VI. II. Problem Formulation Most commonly the TIE polynomial model projects aheadonahorizonofn points from the start point n =0 with the K-degree Taylor polynomial: x 1 (n) = K p=0 x p+1 τ p n p p! + w 1 (n τ) = x 1 + x τn + x 3 τ n + x 4 6 τ 3 n 3 + w 1 (n τ) (3) where x l+1 x l+1 (0) l [0K] are initial states of the clock and w 1 (n τ) is a noise with known properties. By extending the time derivatives of the TIE model to the Taylor series (3) and () become respectively: λ(n) =A(n)λ(0) + w(n τ) (4) ξ(n) =Cλ(n)+v(n) (5) where λ(n) =[x 1 (n)x (n)...x K+1 (n)] T is a vector [(K + 1) 1] of the clock states and a time-varying clock matrix [(K +1) (K +1)]is: 1 τn τ n /... (τn) K /K! 0 1 τn... (τn) K 1 /(K 1)! A(n) = (τn) K /(K )! (6) The most common situation that may be assumed in timekeeping is when all or several clock states are observable by (5) via M (independent or dependent) measurements in the presence of correlated or uncorrelated noises. Then the observation vector is ξ(n) = [ξ 1 (n)ξ (n)...ξ M (n)] T and: c C = 0 c (7) c MM... 0 is a measurement matrix [M (K + 1)] with typically unit components c uu u [1M]. The clock noise vector is w(n τ) = [w 1 (n τ)w (n τ)...w K+1 (n τ)] T with the components caused by the oscillator noises. In view of Fig. 1 the noise vector v(n) =[v 1 (n)v (n)...v M (n)] T contains correlated or uncorrelated components that are not obligatory Gaussian. It is important that the GPS noise v(n) dominates on a horizon N; that is typically wu(n τ) N vl (n) N. Therefore w(n τ) maybeneglected in the averaging FIR procedure (it cannot be discarded in the Kalman filter). The problem now formulates as follows. Given a noisefree TIE model (4) w(n τ) = 0 that is observable via (5) in the presence of a mean zero noise v(n) =0distributed arbitrary with a known covariance v(n)v T (n). We want to derive a real-time unbiased FIR estimator of the clock states λ(n) supposing in this paper that M = K +1andC is an identity matrix. III. An Unbiased FIR Filter Using N points of the nearest past the FIR estimate ˆλ(n) =[ˆx 1 (n)ˆx (n)...ˆx K+1 (n)] T at n-th point is given by neglecting w(n τ): ˆλ(n) = N 1 H(i)ξ(n i) = q(n)γ = d(n)γ + r(n)γ (8a) where an auxiliary unit matrix (N 1) is Γ =[ ] T and the matrix [(K +1) (K + 1)] of unknown FIRs is: h K (i) H(i) = 0 h K 1 (i) (9) h 0 (i) in which the l-th FIR has inherent properties: h l (i) = { h l (i) 0 i N 1 0 otherwise and N 1 h l(i) = 1. The es-
3 864 ieee transactions on ultrasonics ferroelectrics and frequency control vol. 53 no. 5 may 006 timates vector and its deterministic and random constituents all of [(K +1) N] are respectively: q(n) =[H(0)ξ(n) H(1)ξ(n 1)...H(N 1)ξ(n N +1)] d(n) =[H(0)Cλ(n) H(1)Cλ(n 1)...H(N 1)Cλ(n N +1)] r(n) =[H(0)v(n) H(1)v(n 1)...H(N 1)v(n N +1)]. The mean square error (MSE) of ˆλ(n) is written as: J(n) = [λ(n) ˆλ(n)] T [λ(n) ˆλ(n)] =[λ(n) d(n)γ] T [λ(n) d(n)γ] + [r(n)γ] T [r(n)γ] producing the estimate bias and variance respectively: (10) (11) (1) (13) bias[ˆλ(n)] = λ(n) d(n)γ (14) var[ˆλ(n)] = [r(n)γ] T [r(n)γ]. (15) A. Generic Coefficients for the FIR of an Unbiased Filter The necessary and sufficient condition for the unbiased estimate follows straightforwardly from (14); that is: λ(n) =d(n)γ (16) providing the rule to derive the FIRs for the clock states: x 1 (n) WK T x (n)... = λ 1(n) WK 1 T λ (n)... (17) x K+1 (n) W0 T λ K+1 (n) where W l =[h l (0)h l (1)...h l (N 1)] T (18) x K+1 l (n) λ K+1 l (n) = x K+1 l (n 1)... x K+1 l (n N +1). (19) For the clock (K +1 l)-th state (17) thus yields a relation: x K+1 l (n) =W T l λ K+1 l(n). (0) It is seen that (17) is valid for arbitrary n. Setting n = 0 expressing the components in (19) with the l- degree polynomial by (3) and involving the property N 1 h l (i) = 1 we go from (17) to the unbiasedness (or deadbeat) constraint: FW l = G (1) in which G =[ ] T is of [(l +1) 1] and N 1 F = (N 1) () 0 1 l... (N 1) l We notice that (1) is inherent for any other unbiased estimators e.g. for the minimum variance unbiased (MVU) and best linear unbiased estimator (BLUE). It follows from Kalman filtering that an optimum estimate is achieved if the model degree is equal to a number of the filter states minus one. This means that the components in (18) also may be substituted by the l-degree polynomial such as: h l (i) = l a jl i j (3) j=0 where a jl are still unknown coefficients. Embedded (3) the constraint (1) becomes: DB = G (4) where the FIR coefficients matrix is B =[a 0l a 1l... a ll ] T and an auxiliary quadratic matrix [(l +1) (l +1)] is given by: d 0 d 1 d... d l d 1 d d 3... d l+1 D = d d 3 d 4... d l (5) d l d l+1 d l+... d l where the generic component d m = N 1 im m [0 l] is determined by the Bernoulli polynomials (Appendix I). An analytic solution of (4) yields the generic coefficients for (3): a jl =( 1) j M (j+1)1 (6) D in which D and M (j+1)1 are the determinant and minor of (5) respectively. Determined a jl and h l (n) the unbiased FIR estimate of the clock (K +1 l)-th state becomes 1 : ˆx K+1 l (n) = N 1 h l (i)ξ K+1 l (n i) (7a) = Wl T Ξ K+1 l (n) (7b) = G T (FF T ) 1 FΞ K+1 l (n) (7c) 1 The second reviewer mentioned that (3) can be written as W l = F T B; therefore (4) follows from (1) with D = FF T.Wenotice that (4) then may be solved by B =(FF T ) 1 G and (7b) rewritten as (7c).
4 shmaliy: proposed finite impulse response filter 865 where: ξ K+1 l (n) Ξ K+1 l (n) = ξ K+1 l (n 1).... (8) ξ K+1 l (n N +1) B. Estimate Noise The estimate noise variance (15) now may be rewritten as: var[ˆλ(n)] = [r(n)γ] T [r(n)γ] K = WK p T R p+1(n)w K p p=0 (9) where the autocorrelation matrix R l (n) of(n N) has a generic component R l (i j) = v l (i)v l (j) i j [n n N + 1]. Accordingly the estimate variance associated with (7b) calculates: σ K+1 l (n) =WT l R K+1 l(n)w l. (30) It is important that by large N the sawtooth noise becomes delta-correlated. This degenerates R l (n) tothe diagonal form with the components R l (i i) =σvl (i) and tends (30) toward: σ K+1 l (n) =σ v(k+1 l) (n)wt l W l (31) where σv(k+1 l) (n) is a variance of the noise perturbing the (K +1 l)-th clock state at nth point. IV. Applications to the Clock TIE Polynomial Model In applications K is identified for the filter memory on a horizon [0 N 1]) by the clock precision. Typically it is assumed that (3) fits atomic clocks by K 1 and crystal clocks by K. However K =3mayberequiredfor low-precision crystal clocks. Below we derive the unbiased FIRs for all these cases. A. Low-Order FIRs for Unbiased Filters Setting l =0 1 3 and using the coefficient d m (Appendix I) we first calculate (6). Then (3) leads to the unique FIRs namely: h 1 (i) = h 0 (i) = 1 N (3) (N 1) 6i (33) N(N +1) The sawtooth noise produced by SynPaQ III becomes practically delta-correlated if N = 1800 with τ =1sorN = 180 with τ =10s that corresponds to 0.5 hour of averaging. Fig.. Unique FIRs of the unbiased filter: h 0 (i) by (3); h 1 (i) by (33); h (i) by (34); and h 3 (i) by (35). An example is given for N = 60. h (i) = 3(3N 3N +) 18(N 1)i +30i N(N +1)(N +) (34) h 3 (i) = 8(N 3 3N +7N 3) 0(6N 6N +5)i N(N +1)(N +)(N +3) + 10(N 1)i 140i 3 N(N +1)(N +)(N +3). (35) Fig. sketches (3) (35) for N = 60. A uniform FIR (3) corresponds to simple averaging and is optimal in a sense of a minimum-produced noise. This FIR is practically proven to be reasonable in GPS-based common view measurements [5]. A linear FIR (33) was derived in [0] by using linear regression to compensate a bias of simple averaging. Its kernel starts with a maximum h (0) = (N 1) N(N+1) = 4 N > 0andgoestoaminimum N 1 h (N 1) = (N ) N(N+1) = N 1 N < 0 having zero at n 0 = N 1 h 3. Its special peculiarity is r = lim (0) N h = (N 1) that allows one to synthesize a FIR by saving r = for arbitrary N. It is surprising that the FIR synthesized in such a way is equal to that derived in [1] for the 1-step linear prediction on a horizon [1 N]. B. Noise-Power Gains The noise-power gain corresponding to the l-degree FIR is specified by (31) to be g l (N) =σ K+1 l /σ v(k+1 l) = W T l W l. Its values associated with (3) (35) are given below respectively: g 0 (N) = 1 N (36) (N 1) g 1 (N) = N(N +1) (37) g (N) = 3(3N 3N +) N(N +1)(N +) (38) g 3 (N) = 8(N 3 3N +7N 3) N(N +1)(N +)(N +3). (39)
5 866 ieee transactions on ultrasonics ferroelectrics and frequency control vol. 53 no. 5 may 006 Fig. 4. Structure of the (K+1)-state unbiased FIR filtering algorithm for the K-degree TIE polynomial model observable with a single GPS timing receiver. Fig. 3. Noise gains of the unbiased FIR filters: l = 0 by simple averaging (36); l = 1 by (37); k = by (38); and k = 3 by (39). Dashed lines are the upper bounds calculated by (40). Fig. 3 illustrates (36) (39) manifesting that unbiasedness is achieved at increase of noise. Indeed the curves for l>0 trace above the lower bound 1/ N associated with simple averaging (l = 0) that produces minimum noise (among all filters). It also follows that by large N the noise gain is performed by g l (N) = (l +1)/ N and thus traces below the upper bound: { (l +1)/ N N (l +1) gl (N) 1 N < (l +1). (40) V. An Unbiased FIR Filtering Algorithm for a Single Multichannel GPS-Timing Receiver We now consider an important practical case in which the measurement ξ 1 (n) ofatiex 1 (n) is obtained with a single multichannel GPS-timing receiver. Here observations of the higher-order states may be formed by increments of the lower-order estimates. An unbiased FIR filtering algorithm then is written as: ˆλ(n) =q(n)γ (41) ξ k (n) = 1 τ [ˆx k 1(n) ˆx k 1 (n 1)] k > 1 (4) where q(n) is given by (10) and the observation components for (10) k>1 are formed by (4). The algorithm is illustrated in Fig. 4. The clock first state estimate ˆx 1 (n) is obtained with h K (i) atahorizon of N K points. The observation ξ (n) for the second state x (n) then is formed using (4) by increments of ˆx 1 (n). Accordingly ˆx (n) is achieved with h K 1 (i)atahorizonof N K 1 points. Inherently the first accurate value of ˆx (n) appears at (N K +N K 1 )th point starting from n =0.The last state estimate ˆx K+1 (n) is calculated with h 0 (i) ata horizon of N 0 points using ξ K+1 (n) thatisformedinthe same manner as ξ (n). The first correct value of ˆx K+1 (n) appears at (N K + N K N 0 )th point. For the quadratic TIE model (K =crystalclocks) the 3-state unbiased FIR batch algorithm becomes by (7a) and (4): ˆx (n) = 1 τ ˆx 1 (n) = 1 1 j=0 1 h (i)ξ 1 (n i) (43) h 1 (j)[ˆx 1 (n j) ˆx 1 (n j 1)] (44) ˆx 3 (n) = [ˆx (n r) ˆx (n r 1)] τn 0 r=0 (45) where h (i) andh 1 (i) are given by (34) for N = N and (33) for N = N 1 respectively. For the linear TIE model (K = 1 atomic clocks) (41) and (4) simplify to the - state form of: ˆx 1 (n) = 1 1 h 1 (i)ξ 1 (n i) (46) ˆx (n) = [ˆx 1 (n j) ˆx 1 (n j 1)] τn 0 j=0 (47) where h 1 (i) is given by (33) with N = N 1. Each state also may be calculated using the matrix forms (7b) or (7c). Below as an example of application we use this algorithm to estimate the TIE model of an oven crystal clock embedded to the Stanford Frequency Counter SR60 (Stamford Research Systems Inc. Sunnyvale CA). The measurement is done with SynPaQ III and SR60 for τ = 1 s (GPS measurement). Simultaneously to get a reference trend the TIE of the same crystal clock is measured by SR65 (Stamford Research Systems Inc.) for the rubidium clock (Rb-measurement). Initial time and frequency shifts between two measurements then are eliminated statistically and a transition to τ = 10 s is provided by the data thinning in time. At the early stage the TIE model was identified to be quadratic K =andn l are
6 shmaliy: proposed finite impulse response filter 867 TABLE I Average Error (error) and Allan Deviation (σ) ofthe Estimate (Est) for 9.7 Hours and τ =10s: F is FIR and K is Kalman. Errors are Given for Rb-Measurements. x ns y 10 1 D /s Est error σ x(10) error σ y(10) error σ D (10) F K K-F determined for each estimate in the minimum MSE sense 3. We also compare the unbiased FIR estimates to those obtained with the 3-state Kalman filter (Appendix II). A. Several Hours Measurements In this experiment a short-term measurement of the TIE has been done during several hours [Fig. 5(a)]. The algorithm then was run. The horizons were identified for τ = 10 s to be N 1 = 155 or 0.43 hours N = 950 or.64 hours and N 3 = 860 or.39 hours for the Rbmeasurements. Thereafter we set the values of q s in the Kalman filter (Appendix II) to obtain the minimum MSEs for the FIR estimates. Fig. 5 and Table I illustrate these studies showing that the unbiased FIR estimates ˆx 1 (n) ˆx (n) and ˆx 3 (n) and the relevant Kalman estimates ˆx(n) ŷ(n) and ẑ(n) respectively are consistent with however some differences. It follows from Table I that the FIR filter works accurately. Fig. 5(a) shows that ˆx 1 (n) andˆx(n) trackthemean value of the GPS measurement and that their offsets from the Rb-measurement are coursed mostly by the GPS time uncertainty. In this experiment a maximum estimate error of about 60 ns was indicated between the 8th and 9th hours when a time shift in the 1 PPS signal has occurred. It follows [Fig. 5(b)] that ˆx (n) andŷ(n) fitwellthe weighted by 1/τ increments of the Rb-measurement. Even so there are two special ranges (dashed). In range I the frequency shift of about has occurred in the span between the 7th and 8th hours and no appreciable error is indicated in a range of large time shifts (between the 8th and 9th hours) in Fig. 5(a). We associate it with the frequency shift in SR65. In range II the Kalman filter demonstrates a brightly pronounced instability caused likelybythetemporarymodel uncertainty but the FIR estimate still is consistent. We watch for a bit shifted trends of ˆx 3 (n) andẑ(n) in Fig. 5(c) that may be explained by some inconsistency between the q s and N l. It also is seen that ẑ(n) traces much upper ˆx 3 (n) after about 8.7 hours. We associate it with the Kalman filter instability like the case of a range II in Fig. 5(a). 3 To identify K ˆx 1 (n) is compared to the reference (rubidium) measurement in the MSE sense by changing N for K [0 3]. The minimum MSE identifies K and determines N for the TIE estimate. Given K the other N l are determined in the same manner. (a) (b) (c) Fig. 5. Short-time measurement and estimation of the crystal clock TIE model with the 3-state unbiased FIR algorithm and the 3-state Kalman filter: (a) TIE (b) fractional frequency offset and (c) linear fractional frequency drift rate. The experiment was repeated for τ =1s.Theresults are presented in Table II to mention that on the whole the picture (Fig. 5) remains the same. The only principle point to notice is that the Allan deviations of all estimates are reduced by a large number of the points. The FIR and Kalman estimates behave here closer to each othereven though the former is still more accurate with its lower error and much lower Allan variance. B. Long-Term Measurements The same crystal clock was later examined during about.5 days using only the unbiased FIR filter. The measure-
7 868 ieee transactions on ultrasonics ferroelectrics and frequency control vol. 53 no. 5 may 006 TABLE II Average Error (error) and Allan Deviation (σ) ofthe Estimate (Est) for 9.7 Hours and τ =1s: F is FIR and K is Kalman. Errors are Given for Rb-Measurements. x ns y 10 1 D /s Est error σ x(1) error σ y(1) error σ D (1) F K K-F ments inherently show oscillations caused by day s variations in temperature [Fig. 6(a)] and like the previous case all FIR estimates fit well the Rb-measurement. Using ˆx (n) the temperature drift [Fig. 6(b)] was estimated to be about (14to4 C) and ˆx 3 (n) calculatesthe aging rate by ˆx 3 (n) = /day [Fig. 6(c)]. (a) VI. Conclusions In this paper we presented an unbiased FIR filter for the TIE K-degree polynomial model of a local clock. In contrast to the standard Kalman filter the proposed solution does not require the noises to be white and does not involve any knowledge about noises in the algorithm. Instead the FIR filter needs a length N l that is determined for a given clock using a reference source in the minimum MSE sense. The filter produces a noise with a variance that reduces as a reciprocal of N l. We notice that timekeeping operates with large horizons N 1 and thus one should not expect appreciable discrepancy between the optimum andunbiasedestimates. The trade-off between the 3-state unbiased FIR algorithm and the 3-state Kalman algorithm has shown their consistency. However as it was demonstrated experimentally the FIR filter produces a smaller error and a lower Allan variance for the sawtooth noise. Moreover it may be advanced further to be optimal in a sense of a minimum MSE that is currently under investigation. Appendix A Coefficients of Matrix (5) The coefficients for (5) are calculated by: d m = N 1 i m = 1 m +1 [B m+1(n) B m+1 ] where B n (x) is the Bernoulli polynomial and B n = B n (0) is the Bernoulli coefficient. For low orders B n (x) maybe found in the reference books. For high orders the following recurrent relation is valid: B n (x) =n B n 1 (x)dx + B n. (b) (c) Fig. 6. Long-term measurement and estimation of the crystal clock TIE with the 3-state unbiased FIR algorithm: (a) TIE (b) fractional frequency offset and (c) linear fractional frequency drift rate. Several low order coefficients d m are given below: d 0 = N N(N 1) d 1 = d = N(N 1)(N 1) 6 d 3 = N (N 1) 4 d 4 = N(N 1)(N 1)(3N 3N 1) 30 d 5 = N (N 1) (N N 1) 1 d 6 = N(N 1)(N 1)(3N 4 6N 3 +3N +1). 4 For large horizon N 1 the coefficients d m may be calculated by d m N 1 = N m+1 m+1.
8 shmaliy: proposed finite impulse response filter 869 Appendix B Three-State Kalman Filtering Algorithm In the state space the TIE model (1) is given by: x(n) y(n) = 1 ττ / x(n 1) 01 τ y(n 1) + w 1(n τ) w (n τ) z(n) 00 1 z(n 1) w 3 (n τ) x(n) =Ax(n 1) + w(n τ) and (5) becomes assuming a single receiver ξ(n) = [ 100 ] x(n) y(n) + v(n) z(n) ξ(n) =Cx(n)+v(n). The noises w(n τ) andv(n) are mean zero and jointly uncorrelated. The sawtooth noise v(n) has a uniform distribution p(v) =1/v max and correlated increments. Its white Gaussian approximation has a variance V = σv = 1 vmax v max v max v dv = vmax/3. The autocorrelation matrix of the white Gaussian noise w(n) is given by [19]: q 1 + qτ 3 + q3τ 4 0 Ψ = τ q τ + q3τ 3 8 q + q3τ 3 q 3τ 6 q τ + q3τ 3 8 q 3τ 6 q 3τ q 3τ q 3 in which the diffusion coefficients q s namely q 1 q and q 3 specify the white FM noise (WHFM) white random walk FM noise (WRFM) and white random run FM noise (RRFM) respectively in the τ-domain power law. The linear Kalman filtering algorithm reads as follows. Enter the q s R n 1 andˆx n 1 then calculate recursively: R n = AR n 1 A T + Ψ K n = R n C T (C R n C T + V ) 1 ˆx n = Aˆx n 1 + K n (ξ n CAˆx n 1 ) R n =(I K n C) R n. Acknowledgment The author would like to thank Dr. Raymond Filler of the U.S. Army Research Development and Engineering Command (RDECOM) CERDEC; Dr. Charles Greenhall of the Jet Propulsion Laboratory (JPL) California Institute of Technology; Dr. Judah Levine of the National Institute of Standards and Technology (NIST); and two anonymous reviewers for valuable comments and remarks. References [1] W. Lewandowski G. Petit and C. Thomas Precision and accuracy of GPS time transfer IEEE Trans. Instrum. Meas. vol. 4 no. pp Apr [] F. Meyer Common-view and melting-pot GPS time transfer with the UT+ in Proc. 3nd Annu. Precise Time and Time Interval Meeting 000 pp [3] D. W. Allan J. E. Gray and H. E. Machlan The National Bureau of Standards Atomic Time Scale: Generation Stability Accuracy and Accessibility. NBS Monograph 140 Time and Frequency: Theory and Fundamentals National Institute of Standards and Technology 1974 pp [4] R. M. Hambly and T. A. Clark Critical evaluation of the Motorola M1+ GPS timing receiver vs. the master clock at the United States Naval Observatory Washington DC in Proc. 34th Annu. Precise Time and Time Interval Meeting 00 pp [5] J. Levine Time transfer using multi-channel GPS receivers IEEE Trans. Ultrason. Ferroelect. Freq. Contr. vol. 46 no. pp Mar [6] D. W. Allan and J. A. Barnes Optimal time and frequency transfer using GPS signals in Proc. 36th Annu. Freq. Contr. Symp. 198 pp [7] P. V. Tryon and R. H. Jones Estimation of parameters in models for cesium beam atomic clocks J. Res. National Bureau of Standards pp [8] J. W. Chaffee Relating the Allan variance to the diffusion coefficients of a linear stochastic differential equation model for precision oscillators IEEE Trans. Ultrason. Ferroelect. Freq. Contr. vol. 34 no. 6 pp Nov [9] S. R. Stein and R. L. Filler Kalman filter analysis for real time applications of clocks and oscillators in Proc. 4th Annu. Freq. Contr. Symp pp [10] J. A. Barnes R. H. Jones P. V. Tryon and D. W. Allan Stochastic models for atomic clocks in Proc. 14th Annu. Precise Time and Time Interval Meeting 198 pp [11] L. Breakiron Timescale algorithms combining cesium clocks and hydrogen masers in Proc. 3rd Annu. Precise Time and Time Interval Meeting 1991 pp [1] C. Greenhall Kalman plus weights: A time scale algorithm in Proc. 33rd Annu. Precise Time and Time Interval Meeting 001 pp [13] K. Senior P. Koppang and J. Ray Developing an IGS time scale IEEE Trans. Ultrason. Ferroelect. Freq. Contr. vol. 50 no. 6 pp Jun [14] L. Galleani and P. Tavella On the use of the Kalman filter in timescales Metrologia vol. 40 pp [15] C. J. Masreliez Approximate non-gaussian filtering with linear state and observation relations IEEE Trans. Automat. Contr. vol. 0 pp [16] H. Wu and G. Chen Suboptimal Kalman filtering for linear systems with Gaussian-sum type of noise Math. Comput. Model. vol. 9 no. 5 pp [17] O. K. Kwon W. H. Kwon and K. S. Lee FIR filters and recursive forms for discrete-time state-space models Automatica vol. 5 pp [18] W. H. Kwon P. S. Kim and P. Park A receding horizon Kalman FIR filter for discrete time-invariant systems IEEE Trans. Automat. Contr. vol. 44 no. 9 pp [19] W. H. Kwon P. S. Kim and S. H. Han A receding horizon unbiased FIR filter for discrete-time state space models Automatica vol. 38 pp [0] Y. S. Shmaliy A simple optimally unbiased MA filter for timekeeping IEEE Trans. Ultrason. Ferroelect. Freq. Contr. vol. 49 no. 6 pp Jun. 00. [1] P. Heinonen and Y. Neuvo FIR-median hybrid filters with predictive FIR structures IEEE Trans. Acoust. Speech Signal Processing vol. 36 no. 6 pp Jun Yuriy S. Shmaliy (M 96 SM 00) was born January He received the B.S. M.S. and Ph.D. degrees in and 198 respectively from Aviation Institute of Kharkiv Kharkiv Ukraine all in electrical engineering. In 199 he received the Doctor of Technical Sc. degree from the Railroad Academy of Kharkiv Kharkiv Ukraine. In March 1985 he joined the Kharkiv MilitaryUniversityKharkivUkraine.Heserved as professor beginning in In 1999 he joined the Kharkiv National University of
9 870 ieee transactions on ultrasonics ferroelectrics and frequency control vol. 53 no. 5 may 006 Radio Electronics Kharkiv Ukraine. Since November 1999 he has been with the Guanajuato University of Mexico Salamanca Gto. Mexico as a professor. Dr. Shmaliy has 178 papers and 80 patents. He was awarded a title Honorary Radio Engineer of the USSR in 1991; was listed in Marquis Who s Who in the World in 1998; and was listed in Outstanding People of the 0th Century Cambridge England in He is a member of several professional societies and organizing and program committees of International Symposia. His current interests include the statistical theory of precision resonators and oscillators optimal estimation and stochastic signal processing for frequency and time.
STUDIES OF THE UNBIASED FIR FILTER FOR THE TIME ERROR MODEL IN APPLICATIONS TO GPS-BASED TIMEKEEPING 1
STUDIES OF THE UNBIASED FIR FILTER FOR THE TIME ERROR MODEL IN APPLICATIONS TO GPS-BASED TIMEKEEPING 1 Yu. Shmaliy and O. Ibarra-Manzano Guanajuato University, FIMEE, Mexico E-mail: shmaliy@salamanca.ugto.mx
More informationImplementation of Digital Unbiased FIR Filters with Polynomial Impulse Responses
Circuits Syst Signal Process (2012) 31:611 626 DOI 10.1007/s00034-011-9330-9 Implementation of Digital Unbiased FIR Filters with Polynomial Impulse Responses Paula Castro-Tinttori Oscar Ibarra-Manzano
More informationCity, University of London Institutional Repository
City Research Online City, University of London Institutional Repository Citation: Zhao, S., Shmaliy, Y. S., Khan, S. & Liu, F. (2015. Improving state estimates over finite data using optimal FIR filtering
More informationTHE LONG-TERM STABILITY OF THE U.S. NAVAL OBSERVATORY S MASERS
THE LONG-TERM STABILITY OF THE U.S. NAVAL OBSERVATORY S MASERS Demetrios Matsakis, Paul Koppang Time Service Department U.S. Naval Observatory Washington, DC, USA and R. Michael Garvey Symmetricom, Inc.
More informationA KALMAN FILTER FOR ATOMIC CLOCKS AND TIMESCALES
33rdAnnual Precise Time and Time Interval (PTTI) Meeting A KALMAN FILTER FOR ATOMIC CLOCKS AND TIMESCALES Lee A. Breakiron U.S. Naval Observatory Washington, DC 2392542,USA Abstract In a test of whether
More informationSTUDIES OF NPL S CLOCK ENSEMBLE ALGORITHM
43 rd Annual Precise Time and Time Interval (PTTI) Systems and Applications Meeting STUDIES OF NPL S CLOCK ENSEMBLE ALGORITHM Setnam L. Shemar, John A. Davis, and Peter B. Whibberley National Physical
More informationA KALMAN FILTER CLOCK ALGORITHM FOR USE IN THE PRESENCE OF FLICKER FREQUENCY MODULATION NOISE
A KALMAN FILTER CLOCK ALGORITHM FOR USE IN THE PRESENCE OF FLICKER FREQUENCY MODULATION NOISE J. A. Davis, C. A. Greenhall *, and P. W. Stacey National Physical Laboratory Queens Road, Teddington, Middlesex,
More informationSTUDIES OF AN OPTIMALLY UNBIASED MA FILTER INTENDED FOR GPS-BASED TIMEKEEPING
33rdAnnual Precise Time and Time Interval (PTTI) Meeting STUDIES OF AN OPTIMALLY UNBIASED MA FILTER INTENDED FOR GPS-BASED TIMEKEEPING. Y Shmaliy, 0. Ibarra-Manzano, R. Rojas-Laguna, and R. Vazguez-Bautista
More informationLee A. Breakiron U.S. Naval Observatory Washington, DC ABSTRACT
A COMPARATIVE STUDY OF CLOCK RATE AND DRIFT ESTIMATION Lee A. Breakiron U.S. Naval Observatory Washington, DC 20392 ABSTRACT Five different methods of drift determination and four different methods of
More informationFIR Filters for Stationary State Space Signal Models
Proceedings of the 17th World Congress The International Federation of Automatic Control FIR Filters for Stationary State Space Signal Models Jung Hun Park Wook Hyun Kwon School of Electrical Engineering
More informationA KALMAN FILTER FOR ATOMIC CLOCKS AND TIMESCALES
33rdAnnual Precise Time and Time Interval (PTTI) Meeting A KALMAN FILTER FOR ATOMIC CLOCKS AND TIMESCALES Lee A. Breakiron U.S. Naval Observatory Washington, DC 2392-542,USA Abstract In a test of whether
More informationDesign of FIR Smoother Using Covariance Information for Estimating Signal at Start Time in Linear Continuous Systems
Systems Science and Applied Mathematics Vol. 1 No. 3 2016 pp. 29-37 http://www.aiscience.org/journal/ssam Design of FIR Smoother Using Covariance Information for Estimating Signal at Start Time in Linear
More informationEE482: Digital Signal Processing Applications
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu EE482: Digital Signal Processing Applications Spring 2014 TTh 14:30-15:45 CBC C222 Lecture 11 Adaptive Filtering 14/03/04 http://www.ee.unlv.edu/~b1morris/ee482/
More informationTHE LONG-TERM STABILITY OF THE U.S. NAVAL OBSERVATORY S MASERS
THE LONG-TERM STABILITY OF THE U.S. NAVAL OBSERVATORY S MASERS Demetrios Matsakis, Paul Koppang Time Service Department U.S. Naval Observatory Washington, DC, USA and R. Michael Garvey Symmetricom, Inc.
More informationThe Statistics of GPS
The Statistics of GPS Demetrios Matsakis U.S. Naval Observatory, 3450 Massachusetts Avenue NW, Washington, DC., USA, 20392 ABSTRACT The Global Positioning System (GPS) is an extremely effective satellite-based
More informationDiscrete Simulation of Power Law Noise
Discrete Simulation of Power Law Noise Neil Ashby 1,2 1 University of Colorado, Boulder, CO 80309-0390 USA 2 National Institute of Standards and Technology, Boulder, CO 80305 USA ashby@boulder.nist.gov
More informationThe Development of a New Kalman-Filter Time Scale at NIST
The Development of a ew Kalman-Filter Time Scale at IST Jian Yao, Thomas Parker, and Judah Levine Time and Frequency Division, ational Institute of Standards and Technology ABSTRACT We report on a preliminary
More informationCity, University of London Institutional Repository
City Research Online City, University of London Institutional Repository Citation: Granados-Cruz, M., Shmaliy, Y. S., Khan, S., Ahn, C. K. & Zhao, S. (15). New results in nonlinear state estimation using
More informationA NEW REALIZATION OF TERRESTRIAL TIME
A NEW REALIZATION OF TERRESTRIAL TIME G. Petit BIPM, Pavillon de Breteuil, 92312 Sèvres Cedex, France E-mail: gpetit@bipm.org Abstract Terrestrial Time TT is a time coordinate in a geocentric reference
More informationTHE FUTURE MODEL OF TA (TL)
THE FUTURE MODEL OF TA (TL) Shinn-Yan Lin, Hsin-Ming Peng, Weng-Hung Tseng, Hung-Tien Lin, and Chia-Shu Liao National Standard Time and Frequency Lab., TL, Chunghwa Telecom Co. Ltd. No. 12 Lane 551, Min-Tsu
More informationEE482: Digital Signal Processing Applications
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu EE482: Digital Signal Processing Applications Spring 2014 TTh 14:30-15:45 CBC C222 Lecture 11 Adaptive Filtering 14/03/04 http://www.ee.unlv.edu/~b1morris/ee482/
More informationINVESTIGATING THE BIASES IN ALLAN AND HADAMARD VARIANCES AS MEASURES OF M TH ORDER RANDOM STABILITY
43 rd Annual Precise Time and Time Interval (PTTI) Systems and Applications Meeting INVESTIGATING THE BIASES IN ALLAN AND HADAMARD VARIANCES AS MEASURES OF M TH ORDER RANDOM STABILITY Victor R. Reinhardt
More informationSYMMETRICOM TIME-SCALE SYSTEM
SYMMETRICOM TIME-SCALE SYSTEM Timothy Erickson, Venkatesan Ramakrishnan, Samuel R. Stein 3 Symmetricom, Inc., Boulder, CO, USA, tierickson@symmetricom.com Symmetricom, Inc., Santa Rosa, CA, USA, vramakrishnan@symmetricom.com
More informationModeling and Estimation of Stationary and Non-stationary Noises of Rubidium Atomic Clock
ISSN : 48-96, Vol. 4, Issue 7( Version ), July 014, pp.44-49 RESEARCH ARTICLE OPEN ACCESS Modeling and Estimation of Stationary and Non-stationary Noises of Rubidium Atomic Clock 1 Deepak Mishra, Alak
More informationSparse Least Mean Square Algorithm for Estimation of Truncated Volterra Kernels
Sparse Least Mean Square Algorithm for Estimation of Truncated Volterra Kernels Bijit Kumar Das 1, Mrityunjoy Chakraborty 2 Department of Electronics and Electrical Communication Engineering Indian Institute
More informationAcomplex-valued harmonic with a time-varying phase is a
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 46, NO. 9, SEPTEMBER 1998 2315 Instantaneous Frequency Estimation Using the Wigner Distribution with Varying and Data-Driven Window Length Vladimir Katkovnik,
More informationCONFIDENCE ON THE THREE-POINT ESTIMATOR OF FREQUENCY DRIFT*
CONFIDENCE ON THE THREE-POINT ESTIMATOR OF FREQUENCY DRIFT* Marc A. Weiss and Christine Hackman Time and Frequency Division National Institute of Standards and Technology 325 Broadway Boulder, CO 80303
More information5 Kalman filters. 5.1 Scalar Kalman filter. Unit delay Signal model. System model
5 Kalman filters 5.1 Scalar Kalman filter 5.1.1 Signal model System model {Y (n)} is an unobservable sequence which is described by the following state or system equation: Y (n) = h(n)y (n 1) + Z(n), n
More informationAn Ensemble of Atomic Fountains
An Ensemble of Atomic Fountains Steven Peil, Scott Crane, James Hanssen, Thomas B. Swanson and Christopher R. Ekstrom Clock Development Division United States Naval Observatory Washington, DC 39 Abstract
More informationGENERALIZATION OF CAMPBELL S THEOREM TO NONSTATIONARY NOISE. Leon Cohen
GENERALIZATION OF CAMPBELL S THEOREM TO NONSTATIONARY NOISE Leon Cohen City University of New York, Hunter College, 695 Park Ave, New York NY 10056, USA ABSTRACT Campbell s theorem is a fundamental result
More informationAn Analysis of Errors in RFID SAW-Tag Systems with Pulse Position Coding
An Analysis of Errors in RFID SAW-Tag Systems with Pulse Position Coding YURIY S. SHMALIY, GUSTAVO CERDA-VILLAFAÑA, OSCAR IBARRA-MANZANO Guanajuato University, Department of Electronics Salamanca, 36855
More informationLessons in Estimation Theory for Signal Processing, Communications, and Control
Lessons in Estimation Theory for Signal Processing, Communications, and Control Jerry M. Mendel Department of Electrical Engineering University of Southern California Los Angeles, California PRENTICE HALL
More informationTiming Applications and User Equipment R. Michael Garvey. Time and Frequency Services with Galileo Workshop. 5-6 December 2005
Timing Applications and User Equipment R. Michael Garvey Time and Frequency Services with Galileo Workshop 5-6 December 2005 Outline Atomic Clock Technologies Atomic Clock Vendors Rubidium Gas Cell Cesium
More informationGaussian, Markov and stationary processes
Gaussian, Markov and stationary processes Gonzalo Mateos Dept. of ECE and Goergen Institute for Data Science University of Rochester gmateosb@ece.rochester.edu http://www.ece.rochester.edu/~gmateosb/ November
More informationParameter Estimation in a Moving Horizon Perspective
Parameter Estimation in a Moving Horizon Perspective State and Parameter Estimation in Dynamical Systems Reglerteknik, ISY, Linköpings Universitet State and Parameter Estimation in Dynamical Systems OUTLINE
More informationLecture 10 Linear Quadratic Stochastic Control with Partial State Observation
EE363 Winter 2008-09 Lecture 10 Linear Quadratic Stochastic Control with Partial State Observation partially observed linear-quadratic stochastic control problem estimation-control separation principle
More informationAPPLICATION OF THE GSF-1 ALGORITHM TO THE NEAR-OPTIMAL TIMESCALE PREDICTION OF THE HYDROGEN MASER
APPLICATION OF THE GSF-1 ALGORITHM TO THE NEAR-OPTIMAL TIMESCALE PREDICTION OF THE HYDROGEN MASER Laurent-Guy Bernier METAS - Swiss Federal Office of Metrology and Accreditation Lindenweg 50, CH-3003 Bern-Wabern,
More informationA statistic for describing pulsar and clock stabilities
Astron. Astrophys. 326, 924 928 (1997) ASTRONOMY AND ASTROPHYSICS A statistic for describing pulsar and clock stabilities Demetrios N. Matsakis 1, J. H. Taylor 2, and T. Marshall Eubanks 1 1 U.S. Naval
More informationHere represents the impulse (or delta) function. is an diagonal matrix of intensities, and is an diagonal matrix of intensities.
19 KALMAN FILTER 19.1 Introduction In the previous section, we derived the linear quadratic regulator as an optimal solution for the fullstate feedback control problem. The inherent assumption was that
More informationUNCERTAINTY OF STABILITY VARIANCES BASED ON FINITE DIFFERENCES
UNCERTAINTY OF STABILITY VARIANCES BASED ON FINITE DIFFERENCES C. A. Greenhall W. J. Riley Jet Propulsion Laboratory Symmetricom, Inc. California Institute of Technology Abstract This work gives an algorithm
More informationA NEW REALIZATION OF TERRESTRIAL TIME
CCTF/04-17 A NEW REALIZATION OF TERRESTRIAL TIME G. Petit BIPM, Pavillon de Breteuil, 92312 Sèvres Cedex France- e-mail: gpetit@bipm.org ABSTRACT Terrestrial Time TT is a time coordinate in a geocentric
More informationVARIANCES BASED ON DATA WITH DEAD TIME BETWEEN THE MEASUREMENTS *
VARIANCES BASED ON DATA WITH DEAD TIME BETWEEN THE MEASUREMENTS * James A. Barnes Austron Incorporated Boulder, Colorado 80301 David W. Allan Time and Frequency Division National Bureau of Standards Boulder,
More informationResearch Article Time Synchronization and Performance of BeiDou Satellite Clocks in Orbit
Navigation and Observation Volume 213, Article ID 37145, 5 pages http://dx.doi.org/1.1155/213/37145 Research Article Time Synchronization and Performance of BeiDou Satellite Clocks in Orbit Han Chunhao,
More informationClock Modeling & Algorithms for Timescale Formation
Clock Modeling & Algorithms for Timescale Formation K. Senior U.S. Naval Research Laboratory (NRL) 26 July 2012 IGS Workshop 2012 University of Warmia and Mazury Olsztyn, Poland Outline: Motivation / history
More informationDEGREES OF FREEDOM AND THREE-CORNERED HATS
33rd Annual Precise Time and Time Interval (PTTI) Meeting DEGREES OF FREEDOM AND THREE-CORNERED HATS Christopher R. Ekstrom and Paul A. Koppang U.S. Naval Observatory, Washington, D.C. Abstract The three-cornered
More informationNOISE ROBUST RELATIVE TRANSFER FUNCTION ESTIMATION. M. Schwab, P. Noll, and T. Sikora. Technical University Berlin, Germany Communication System Group
NOISE ROBUST RELATIVE TRANSFER FUNCTION ESTIMATION M. Schwab, P. Noll, and T. Sikora Technical University Berlin, Germany Communication System Group Einsteinufer 17, 1557 Berlin (Germany) {schwab noll
More informationKalman Filters with Uncompensated Biases
Kalman Filters with Uncompensated Biases Renato Zanetti he Charles Stark Draper Laboratory, Houston, exas, 77058 Robert H. Bishop Marquette University, Milwaukee, WI 53201 I. INRODUCION An underlying assumption
More informationPROTOTYPE OF THE DLR OPERATIONAL COMPOSITE CLOCK: METHODS AND TEST CASES
PROTOTYPE OF THE DLR OPERATIONAL COMPOSITE CLOCK: METHODS AND TEST CASES Matthias Suess and Jens Hammesfahr E-mail: Matthias.suess@dlr.de German Aerospace Centre (DLR) Oberpfaffenhofen, Germany Abstract
More informationStochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions
International Journal of Control Vol. 00, No. 00, January 2007, 1 10 Stochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions I-JENG WANG and JAMES C.
More informationChapter 2 Wiener Filtering
Chapter 2 Wiener Filtering Abstract Before moving to the actual adaptive filtering problem, we need to solve the optimum linear filtering problem (particularly, in the mean-square-error sense). We start
More informationTitle without the persistently exciting c. works must be obtained from the IEE
Title Exact convergence analysis of adapt without the persistently exciting c Author(s) Sakai, H; Yang, JM; Oka, T Citation IEEE TRANSACTIONS ON SIGNAL 55(5): 2077-2083 PROCESS Issue Date 2007-05 URL http://hdl.handle.net/2433/50544
More informationQuadratic Extended Filtering in Nonlinear Systems with Uncertain Observations
Applied Mathematical Sciences, Vol. 8, 2014, no. 4, 157-172 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.12988/ams.2014.311636 Quadratic Extended Filtering in Nonlinear Systems with Uncertain Observations
More informationOptimal Mean-Square Noise Benefits in Quantizer-Array Linear Estimation Ashok Patel and Bart Kosko
IEEE SIGNAL PROCESSING LETTERS, VOL. 17, NO. 12, DECEMBER 2010 1005 Optimal Mean-Square Noise Benefits in Quantizer-Array Linear Estimation Ashok Patel and Bart Kosko Abstract A new theorem shows that
More informationAdaptiveFilters. GJRE-F Classification : FOR Code:
Global Journal of Researches in Engineering: F Electrical and Electronics Engineering Volume 14 Issue 7 Version 1.0 Type: Double Blind Peer Reviewed International Research Journal Publisher: Global Journals
More informationECE531 Lecture 11: Dynamic Parameter Estimation: Kalman-Bucy Filter
ECE531 Lecture 11: Dynamic Parameter Estimation: Kalman-Bucy Filter D. Richard Brown III Worcester Polytechnic Institute 09-Apr-2009 Worcester Polytechnic Institute D. Richard Brown III 09-Apr-2009 1 /
More informationA NEW SYSTEM FOR THE GENERATION OF UTC(CH)
A NEW SYSTEM FOR THE GENERATION OF UTC(CH) L.G. Bernier and G. Schaller METAS Swiss Federal Office of Metrology Lindenweg 50, Bern-Wabern, CH-3003, Switzerland laurent-guy.bernier@metas.ch Abstract A new
More informationAD-A REPORT DOCUMENTATION PAGE OM8 No UNCLASSIFIED
UNCLASSIFIED SECURITY CLASSIFICATION OF THIS PAGE Form Approved REPORT DOCUMENTATION PAGE OM8 No. 0704-0188 1b RESTRICTIVE MARKINGS AD-A212 005 3. DISTRIBUTION /AVAILABILITY OF REPORT Approved for Public
More informationThe Discrete Kalman Filtering of a Class of Dynamic Multiscale Systems
668 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL 49, NO 10, OCTOBER 2002 The Discrete Kalman Filtering of a Class of Dynamic Multiscale Systems Lei Zhang, Quan
More informationA Subspace Approach to Estimation of. Measurements 1. Carlos E. Davila. Electrical Engineering Department, Southern Methodist University
EDICS category SP 1 A Subspace Approach to Estimation of Autoregressive Parameters From Noisy Measurements 1 Carlos E Davila Electrical Engineering Department, Southern Methodist University Dallas, Texas
More informationRESPONSE OF A GEOSYNCHRONOUS SPACECRAFT S CRYSTAL OSCILLATOR TO SOLAR FLARES: RESULTS OF A SPACE EXPERIMENT
RESPONSE OF A GEOSYNCHRONOUS SPACECRAFT S CRYSTAL OSCILLATOR TO SOLAR FLARES: RESULTS OF A SPACE EXPERIMENT J. Camparo, A. Presser, S. Lalumondiere, and S. Moss The Aerospace Corporation PO Box 92957,
More informationON November 5, 2000 (GPS week 1087), a new set of
PREPRINT FOR PROCEEDINGS OF THE IEEE/EIA INTERNATIONAL FREQUENCY CONTROL SYMPOSIUM, JUNE, Developing an IGS Time Scale Ken Senior, Paul Koppang, Demetrios Matsakis, and Jim Ray Abstract Currently, the
More informationQuantization and Compensation in Sampled Interleaved Multi-Channel Systems
Quantization and Compensation in Sampled Interleaved Multi-Channel Systems The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation
More informationAdaptive Filter Theory
0 Adaptive Filter heory Sung Ho Cho Hanyang University Seoul, Korea (Office) +8--0-0390 (Mobile) +8-10-541-5178 dragon@hanyang.ac.kr able of Contents 1 Wiener Filters Gradient Search by Steepest Descent
More informationTHE RANGE COVERED BY A CLOCK ERROR IN THE CASE OF WHITE FM
! 90th Annual Precise Time and Time Interval (PTTI) Meeting THE RANGE COVERED BY A CLOCK ERROR IN THE CASE OF WHITE FM Patrizia Tavella* and Daniela MeoO JCIstituto Elettrotecnico Nazionale G. Ferraris
More informationProbability Space. J. McNames Portland State University ECE 538/638 Stochastic Signals Ver
Stochastic Signals Overview Definitions Second order statistics Stationarity and ergodicity Random signal variability Power spectral density Linear systems with stationary inputs Random signal memory Correlation
More informationPrediction of ESTSP Competition Time Series by Unscented Kalman Filter and RTS Smoother
Prediction of ESTSP Competition Time Series by Unscented Kalman Filter and RTS Smoother Simo Särkkä, Aki Vehtari and Jouko Lampinen Helsinki University of Technology Department of Electrical and Communications
More informationAN ALGORITHM FOR THE ITALIAN ATOMIC TIME SCALE
AN ALGORITHM FOR THE ITALIAN ATOMIC TIME SCALE F. Cordara, G. Vizio, P. Tavella, V. Pettiti Istituto Elettrotecnico Nazionale Galileo Ferraris Corso Massimo d9azeglio 42 10125 Torino - Italy Abstract During
More informationUsing the Kalman Filter to Estimate the State of a Maneuvering Aircraft
1 Using the Kalman Filter to Estimate the State of a Maneuvering Aircraft K. Meier and A. Desai Abstract Using sensors that only measure the bearing angle and range of an aircraft, a Kalman filter is implemented
More informationMOMENT functions are used in several computer vision
IEEE TRANSACTIONS ON IMAGE PROCESSING, VOL. 13, NO. 8, AUGUST 2004 1055 Some Computational Aspects of Discrete Orthonormal Moments R. Mukundan, Senior Member, IEEE Abstract Discrete orthogonal moments
More information2262 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 47, NO. 8, AUGUST A General Class of Nonlinear Normalized Adaptive Filtering Algorithms
2262 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 47, NO. 8, AUGUST 1999 A General Class of Nonlinear Normalized Adaptive Filtering Algorithms Sudhakar Kalluri, Member, IEEE, and Gonzalo R. Arce, Senior
More informationAdaptive Filtering Part II
Adaptive Filtering Part II In previous Lecture we saw that: Setting the gradient of cost function equal to zero, we obtain the optimum values of filter coefficients: (Wiener-Hopf equation) Adaptive Filtering,
More informationAN ADAPTIVE ALGORITHM TO EVALUATE CLOCK PERFORMANCE IN REAL TIME*
AN ADAPTIVE ALGORITHM TO EVALUATE CLOCK PERFORMANCE IN REAL TIME* Dr. James A. Barnes Austron Boulder, Co. Abstract Kalman filters and ARIMA models provide optimum control and evaluation techniques (in
More informationUpdate on the In-orbit Performances of GIOVE Clocks
Update on the In-orbit Performances of GIOVE Clocks Pierre Waller, Francisco Gonzalez, Stefano Binda, ESA/ESTEC Ilaria Sesia, Patrizia Tavella, INRiM Irene Hidalgo, Guillermo Tobias, GMV Abstract The Galileo
More informationUncertainties of drift coefficients and extrapolation errors: Application to clock error prediction
Vernotte et al - E 1319 - Final Version - Metrologia - February 1, 21 1 Uncertainties of drift coefficients and extrapolation errors: Application to clock error prediction F Vernotte 1, J Delporte 2, M
More informationLTI Systems, Additive Noise, and Order Estimation
LTI Systems, Additive oise, and Order Estimation Soosan Beheshti, Munther A. Dahleh Laboratory for Information and Decision Systems Department of Electrical Engineering and Computer Science Massachusetts
More informationChapter 2 Fundamentals of Adaptive Filter Theory
Chapter 2 Fundamentals of Adaptive Filter Theory In this chapter we will treat some fundamentals of the adaptive filtering theory highlighting the system identification problem We will introduce a signal
More informationUNIT 1. SIGNALS AND SYSTEM
Page no: 1 UNIT 1. SIGNALS AND SYSTEM INTRODUCTION A SIGNAL is defined as any physical quantity that changes with time, distance, speed, position, pressure, temperature or some other quantity. A SIGNAL
More information5 Generation and Dissemination of Time and Frequency Standard
5 Generation and Dissemination of Time and Frequency Standard 5-1 Algorithm of Ensemble Atomic Time HANADO Yuko, IMAE Michito, AIDA Masanori, HOSOKAWA Mizuhiko, ITO Hiroyuki, NAKAGAWA Fumimaru, and SHIMIZU
More informationKALMAN PLUS WEIGHTS: A TIME SCALE ALGORITHM*
33rdAnnual Precise Time and Time nterval (PlTl) Meeting KALMAN PLUS WEGHTS: A TME SCALE ALGORTHM* Charles A. Greenhall Jet Propulsion Laboratory, California nstitute of Technology 48 Oak Grove Drive, MS
More informationCOMP Signals and Systems. Dr Chris Bleakley. UCD School of Computer Science and Informatics.
COMP 40420 2. Signals and Systems Dr Chris Bleakley UCD School of Computer Science and Informatics. Scoil na Ríomheolaíochta agus an Faisnéisíochta UCD. Introduction 1. Signals 2. Systems 3. System response
More informationCONTROL SYSTEMS, ROBOTICS, AND AUTOMATION - Vol. V - Prediction Error Methods - Torsten Söderström
PREDICTIO ERROR METHODS Torsten Söderström Department of Systems and Control, Information Technology, Uppsala University, Uppsala, Sweden Keywords: prediction error method, optimal prediction, identifiability,
More informationAsymptotic Analysis of the Generalized Coherence Estimate
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 49, NO. 1, JANUARY 2001 45 Asymptotic Analysis of the Generalized Coherence Estimate Axel Clausen, Member, IEEE, and Douglas Cochran, Senior Member, IEEE Abstract
More informationDominant Pole Localization of FxLMS Adaptation Process in Active Noise Control
APSIPA ASC 20 Xi an Dominant Pole Localization of FxLMS Adaptation Process in Active Noise Control Iman Tabatabaei Ardekani, Waleed H. Abdulla The University of Auckland, Private Bag 9209, Auckland, New
More informationEfficient Use Of Sparse Adaptive Filters
Efficient Use Of Sparse Adaptive Filters Andy W.H. Khong and Patrick A. Naylor Department of Electrical and Electronic Engineering, Imperial College ondon Email: {andy.khong, p.naylor}@imperial.ac.uk Abstract
More informationTime Series Prediction by Kalman Smoother with Cross-Validated Noise Density
Time Series Prediction by Kalman Smoother with Cross-Validated Noise Density Simo Särkkä E-mail: simo.sarkka@hut.fi Aki Vehtari E-mail: aki.vehtari@hut.fi Jouko Lampinen E-mail: jouko.lampinen@hut.fi Abstract
More informationParametric Signal Modeling and Linear Prediction Theory 1. Discrete-time Stochastic Processes
Parametric Signal Modeling and Linear Prediction Theory 1. Discrete-time Stochastic Processes Electrical & Computer Engineering North Carolina State University Acknowledgment: ECE792-41 slides were adapted
More informationAdaptive MMSE Equalizer with Optimum Tap-length and Decision Delay
Adaptive MMSE Equalizer with Optimum Tap-length and Decision Delay Yu Gong, Xia Hong and Khalid F. Abu-Salim School of Systems Engineering The University of Reading, Reading RG6 6AY, UK E-mail: {y.gong,x.hong,k.f.abusalem}@reading.ac.uk
More informationEFFECTS OF ILL-CONDITIONED DATA ON LEAST SQUARES ADAPTIVE FILTERS. Gary A. Ybarra and S.T. Alexander
EFFECTS OF ILL-CONDITIONED DATA ON LEAST SQUARES ADAPTIVE FILTERS Gary A. Ybarra and S.T. Alexander Center for Communications and Signal Processing Electrical and Computer Engineering Department North
More informationA Framework for Optimizing Nonlinear Collusion Attacks on Fingerprinting Systems
A Framework for Optimizing onlinear Collusion Attacks on Fingerprinting Systems egar iyavash Department of Electrical and Computer Engineering University of Illinois at Urbana-Champaign Coordinated Science
More informationA New Approach for Computation of Timing Jitter in Phase Locked Loops
A New Approach for Computation of Timing Jitter in Phase ocked oops M M. Gourary (1), S. G. Rusakov (1), S.. Ulyanov (1), M.M. Zharov (1),.. Gullapalli (2), and B. J. Mulvaney (2) (1) IPPM, Russian Academy
More informationCooperative Communication with Feedback via Stochastic Approximation
Cooperative Communication with Feedback via Stochastic Approximation Utsaw Kumar J Nicholas Laneman and Vijay Gupta Department of Electrical Engineering University of Notre Dame Email: {ukumar jnl vgupta}@ndedu
More informationALLAN VARIANCE ESTIMATED BY PHASE NOISE MEASUREMENTS
ALLAN VARIANCE ESTIMATED BY PHASE NOISE MEASUREMENTS P. C. Chang, H. M. Peng, and S. Y. Lin National Standard Time & Frequency Lab., TL, Taiwan 1, Lane 551, Min-Tsu Road, Sec. 5, Yang-Mei, Taoyuan, Taiwan
More information798 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL. 44, NO. 10, OCTOBER 1997
798 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL 44, NO 10, OCTOBER 1997 Stochastic Analysis of the Modulator Differential Pulse Code Modulator Rajesh Sharma,
More informationOn Modern and Historical Short-Term Frequency Stability Metrics for Frequency Sources
On Modern and Historical Short-Term Frequency Stability Metrics for Frequency Sources Michael S. McCorquodale, Ph.D. Founder and CTO, Mobius Microsystems, Inc. EFTF-IFCS, Besançon, France Session BL-D:
More informationUCSD ECE250 Handout #20 Prof. Young-Han Kim Monday, February 26, Solutions to Exercise Set #7
UCSD ECE50 Handout #0 Prof. Young-Han Kim Monday, February 6, 07 Solutions to Exercise Set #7. Minimum waiting time. Let X,X,... be i.i.d. exponentially distributed random variables with parameter λ, i.e.,
More informationThe Kalman Filter. Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience. Sarah Dance
The Kalman Filter Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience Sarah Dance School of Mathematical and Physical Sciences, University of Reading s.l.dance@reading.ac.uk July
More informationIntelligent Embedded Systems Uncertainty, Information and Learning Mechanisms (Part 1)
Advanced Research Intelligent Embedded Systems Uncertainty, Information and Learning Mechanisms (Part 1) Intelligence for Embedded Systems Ph. D. and Master Course Manuel Roveri Politecnico di Milano,
More informationTM-Radiation From an Obliquely Flanged Parallel-Plate Waveguide
1534 IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 50, NO. 11, NOVEMBER 2002 TM-Radiation From an Obliquely Flanged Parallel-Plate Waveguide Jae Yong Kwon, Member, IEEE, Jae Wook Lee, Associate Member,
More informationMODERN video coding standards, such as H.263, H.264,
146 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS FOR VIDEO TECHNOLOGY, VOL. 16, NO. 1, JANUARY 2006 Analysis of Multihypothesis Motion Compensated Prediction (MHMCP) for Robust Visual Communication Wei-Ying
More informationOptimal Polynomial Control for Discrete-Time Systems
1 Optimal Polynomial Control for Discrete-Time Systems Prof Guy Beale Electrical and Computer Engineering Department George Mason University Fairfax, Virginia Correspondence concerning this paper should
More information