# State Fusion of Decentralized Optimal Unbiased FIR Filters.

1. IntroductionState estimation plays an important role in many applications such as control, moving target tracking, and timekeeping and clock synchronization [1]. The existing filtering methods can be classified by two types: the infinite impulse response (IIR) filter and the finite impulse response (FIR) filter [2]. Specifically, the former uses all the historical measurements, and a special case is the Kalman filter (KF) [3], while the latter utilizes limited memory over the most recent time interval [4]. It is due to the difference of structure that FIR-type filters exhibit some useful engineering features such as bounded input/bounded output (BIBO) stability, round off errors [5], and better robustness against temporary uncertainties [6].

With higher requirements of flexibility and accuracy, information fusion filtering theory for multisensor systems has been studied and widely applied [7, 8]. For example, [9,10] proposed the weighted fusion estimator according to the maximum likelihood and weight least square respectively with local filtering error cross covariance zero or not. In [11], an optimal decentralized fusion weighted by matrix was proposed to fuse local KF in the linear unbiased minimum variance (LUMV) sense. Furthermore, the situation of correlated multiplicative and additive noise was analyzed in [8], and the decentralized KF fusion problem for cross-correlated measurement noises was investigated in [12] which also discussed the cases with feedback from the fusion center to local sensors. Reference [13] focused on the decentralized fusion estimation problem for networked systems with random delays and packet losses, while [14] aimed at missing measurements and correlated noises. Reference [15] applied decentralized fusion to the maneuvering target tracking with multirate sampling and uncertainties in wireless sensor networks.

The fusion filtering methods nowadays take the KF as local filters mostly, so that they inherit the properties of KF, optimal but not robust under some uncertainties. Meanwhile there are few results considering decentralized FIR-type filter to deal with multisensor systems. Basically optimal FIR (OFIR) filter was obtained in [16] by minimum mean square error sense. Besides unbiased FIR (UFIR) filter was derived for real-time models and realized iteratively in [17], which can ignore the noise statistics but does not guarantee optimality. In between, [18] represents two forms of optimal unbiased FIR (OUFIR) filter, minimum variance unbiased FIR and embedded unbiasedness (EU) constraint on OFIR, and proved the identity of the two further. Contrary to the KF, the OUFIR filter is more robust and performs lower sensitivity to the initial condition [19]. In addition, there are a few new study and improved solutions on FIR filtering [2022] developed nowadays. Some practical applications were reported [23, 24] as well.

In this correspondence, we firstly develop OUFIR filter to an extension to handle the system with or without the control inputs universally. Then considering it as local filters, cross covariance between any two is determined by batch and iterative form under the LUMV sense. Finally an optimal decentralized state fusion filter is proposed, which shows better immunity against temporary model uncertainties than traditional KF. The rest of this paper is organized as follows. In Section 2, we describe the preliminaries of system model. Local OUFIR filters are developed in Section 3. Main results about fusion and cross covariance are given in Section 4. The simulation and conclusions are provided in Sections 5 and 6, respectively.

2. System Model and Preliminaries

Consider the following linear discrete time-invariant system with multiple sensors:

[x.sub.n] = A[x.sub.n-1] + E[u.sub.n] + B[w.sub.n], (1)

[y.sup.[i].sub.n] = [C.sup.[i]][x.sub.n] + [v.sup.[i].sub.n], (2)

where i = 1, 2, ..., r denotes the index of sensor, n is the discrete time index, [x.sub.n] [member of] [R.sup.k] and [mathematical expression not reproducible] denote the state and measurement vectors, respectively, [u.sub.n] [member of] [R.sup.p] is the known control signal, and the system matrices A, E, B, and [C.sup.[i]] are known with appropriate dimensions. In this paper, the process noise [w.sub.n] [member of] [R.sup.q] and measurement noise [mathematical expression not reproducible] are assumed to be correlated white noises with zero mean and

[mathematical expression not reproducible], (3)

[mathematical expression not reproducible] (4)

where [[delta].sub.nm] is the Kronecker delta function, i.e., if n = m, [[delta].sub.nm] = 1, otherwise, [[delta].sub.nm] = 0. For notational convenience, we assign [R.sub.ii] = [R.sub.i], [Q.sub.ii] = [Q.sub.i] when i = j. That is, the process noise [w.sub.n] is correlated with each sensor noise [v.sup.[i].sub.n] which is also correlated with each other.

Firstly the original state-space model (1) and (2) needs to be transformed to be a form with uncorrelated noise. In this way, we add the right-hand side of (1) with a zero term [y.sup.[i].sub.n] - [C.sup.[i]][x.sub.n] - [v.sup.[i].sub.n] as

[mathematical expression not reproducible], (5)

where [mathematical expression not reproducible], and [J.sub.i] is a coefficient matrix. To remove the relevance between [[bar.w].sub.[i].sub.n] and [v.sub.[i].sub.n], the expectation term E{[[bar.w].sub.[i].sub.n] [([v.sup.[i].sub.n]).sup.T]} should be zero by definition. That is,

[mathematical expression not reproducible], (6)

which further gives us [J.sub.i] = [BS.sub.i][([C.sup.[i]]B[S.sub.i] + [R.sub.i]).sup.-1]. It is not difficult to find that E{[[bar.w].sub.[i].sub.n]} = E{(I - [J.sub.i][C.sup.[i]])B[w.sub.n] - [J.sub.i][v.sup.[i].sub.n]} = 0, which means that the new-defined noise term [[bar.w].sup.[i].sub.n] is zero mean also. Accordingly, the cross state noise variance can be computed as

[mathematical expression not reproducible]. (7)

As can be seen, with an auxiliary matrix [J.sub.i], there is no correlation between [[bar.w].sup.[i].sub.m] and [v.sup.[i].sub.n], and we consider the state transition equation (5) and measurement equation (2) subsequently.

By defining an estimation horizon as [m, n], where m [??] N - N +1 is the initial time step and N is the horizon length, models (5) and (2) within the estimation horizon can be extended as [16] similarly,

[mathematical expression not reproducible], (8)

[mathematical expression not reproducible], (9)

where the initial state [x.sub.m] is assumed to be known, and the extended vectors are specified as

[mathematical expression not reproducible]. (10)

The involved extended matrices [mathematical expression not reproducible], and [mathematical expression not reproducible] are all N-dependent, which have the forms of

[mathematical expression not reproducible], (11)

[mathematical expression not reproducible], (12)

[C.sup.[i].sub.n,m] = [[bar.C].sup.[i].sub.n,m] [A.sup.[i].sub.n,m], (13)

[G.sup.[i].sub.n,m] = [[bar.C].sup.[i].sub.n,m] [E.sup.[i].sub.n,m], (14)

[F.sup.[i].sub.n,m] = [[bar.C].sup.[i].sub.n,m] [B.sup.[i].sub.n,m], (15)

where we assign

[mathematical expression not reproducible], (16)

and [E.sup.[i].sub.n,m] = [B.sup.[i].sub.n,m].

Once the extended state-space models (8) and (9) are available, the problem considered in this paper can be formulated as follows. Given the state-space model (1) and (2), we would like to design local OUFIR filter for each sensor.

3. Local Optimal Unbiased FIR Filter

A demonstration of the FIR filter structure is provided in Figure 1. From the figure, FIR estimators explicitly employ N most recent measurements unlike KF. It is in this way that some nice properties like better robustness are achieved. Based on [i.sup.th] for ith subsystem represented as (8) and (9), we design a local OUFIR filter in this section.

3.1. Batch Form. The linear FIR filter for each sensor can be expressed as the following batch form [25] generally:

[mathematical expression not reproducible] (17)

for gain matrices [H.sup.[i].sub.n] and [L.sup.[i].sub.n]. It is noted that the filter defined as (17) handles all the measurements and input values collected within the estimation horizon at one time.

Introduce unbiasedness constraint (or deadbeat constraint) [26]:

E {[[??].sup.[i].sub.n]} = E {[x.sub.n]} (18)

in the estimation (17) where [x.sub.n] can be specified as

[mathematical expression not reproducible], (19)

and [[bar.E].sup.[i].sub.n,m], [[bar.B].sup.[i].sub.n,m] are the first vector rows of [E.sup.[i].sub.n,m], [B.sup.[i].sub.n,m] respectively. By substituting (19) and (17) into (18), replacing the term [Y.sup.[i].sub.n,m] with (9), and providing the averaging, one arrives at two following constraints

[H.sup.[i].sub.n][C.sup.[i].sub.n,m] = [([[bar.A].sub.[i]]).sup.n-m], (20)

[L.sup.[i].sub.n] = [[bar.E].sup.[i].sub.n,m] - [H.sup.[i].sub.n] [G.sup.[i].sub.n,m]. (21)

Substituting (21) into (17) yields

[mathematical expression not reproducible]. (22)

Provided [[??].sup.[i].sub.n], the instantaneous filtering error can be defined as [e.sup.[i].sub.n] = [x.sub.n] - [[??].sup.[i].sub.n], then the OUFIR filter is derived by solving the following optimization problem

[mathematical expression not reproducible], (23)

where [H.sup.[i].sub.n] should satisfy the constrain (20).

By replacing [x.sub.n] with (19) and [[??].sup.[i].sub.n] with (22), and substituting the unbiasedness constrain (20), filtering error is obtained

[mathematical expression not reproducible]. (24)

Employing the trace operation and providing the averaging, the optimization problem (23) can be rewritten as

[mathematical expression not reproducible]. (25)

It is not difficult to get

[mathematical expression not reproducible], (26)

where [Q.sub.ij] is given in (7), and [R.sub.ij] is defined in (4). Referring to [19], we can obtain the following gain:

[mathematical expression not reproducible], (27)

where [Z.sub.w+v,n] = [Z.sub.w,n] + [R.sub.n,m] and

[Z.sub.w,n] = [F.sup.[i].sub.n,m] [Q.sup.[ii].sub.n,m] [([F.sup.[i].sub.n,m]).sup.T]. (28)

At this point, the local OUFIR filter with batch form is specified by (22) with filter gain [H.sup.[i].sub.n] provided as (26). As mentioned, this structure is not suitable for the subsequent fusion step, as N measurements are operated at one time and there is no error covariance available to quantify the estimation accuracy. To address these issues, an equivalent iterative computational formula is given below.

3.2. Equivalent Iterative Computation. Introducing [mathematical expression not reproducible], (22) can be divided into two parts shown as follows:

[[??].sup.[i].sub.1,n] = [[bar.E].sup.[i].sub.n,m] [U.sup.[i].sub.n,m], (29)

[[??].sup.[i].sub.2,n] = [H.sup.[i].sub.n] ([Y.sup.[i].sub.n,m] - [G.sup.[i].sub.n,m][U.sup.[i].sub.n,m]), (30)

and consider each components separately.

By former transformations, for (29), recursive algorithm is represented easily as

[mathematical expression not reproducible] (31)

Assigning [[??].sup.[i].sub.n,m] = [Y.sup.[i].sub.n,m] - [G.sup.[i].sub.n,m][U.sup.[i].sub.n,m], this part becomes [mathematical expression not reproducible]. For this estimate, we notice a iterative form proposed in [19], which is represented similarly

[mathematical expression not reproducible] (32)

where [[??].sup.[i].sub.n] and [[bar.H].sup.[i].sub.n] are gain matrices. Besides, [[??].sup.[i].sub.n] is the first vector rows of [[??].sup.[i].sub.n,m], substituting (31), which is given by

[mathematical expression not reproducible] (33)

[mathematical expression not reproducible] (34)

Then putting (33) into (32), iterative form is

[mathematical expression not reproducible]. (35)

Combining (31) and (35), we arrive at

[mathematical expression not reproducible]. (36)

which suggests easily that the prior estimate is [mathematical expression not reproducible]. Referring to [27], we finally come up with a general iterative OUFIR filtering algorithm similarly whose pseudocode is summarized as Algorithm 1 in which l is an auxiliary variable replaced with n. Unlike the OFIR or KF, it is because of unbiasedness condition that the filter does not require initial information. Instead that can be defined at s = m + k using a short batch form given in Table 1; k is the number of system states. In the Algorithm 1, variable l ranges from s to n. The true filtering result corresponds to l = n.

Once the local estimates [[??].sup.[i].sub.n] are corrected, i = 1, 2, ..., r, an optimal fusion strategy based on cross covariance is employed to get the overall estimates, which is introduced below.

4. Decentralized Fusion Filter

In this section, we will investigate the decentralized fusion filters and first deduce the corresponding filtering error cross covariance between any two local subsystems with FIR structure. Due to inefficient computation, then we find its fast iterative form. Based on these, a multisensor decentralized fusion OUFIR filtering algorithm is finally presented.

4.1. Cross Covariance Matrices and Iterative Computation. For the [i.sup.th] sensor subsystem, filtering error [e.sup.[i].sub.n] is written as (24) where filtering gain [H.sup.[i].sub.n] is specified as (27). So the cross covariance matrix between the [i.sup.th] and the [j.sup.th] subsystem can be showed as

[mathematical expression not reproducible] (37)

where [Q.sup.[ij].sub.n,m], [R.sup.[ij].sub.n,m] are defined in (26).

Algorithm 1: General iterative OUFIR filtering algorithm. 1: for n = N - 1, N, ... do 2: m = n-N+ 1, s = m + k Table 1: Computation of the initial state 3: for l = s : n do 4: [mathematical expression not reproducible] 5: [mathematical expression not reproducible] 6: [mathematical expression not reproducible] 7: [mathematical expression not reproducible] 8: [mathematical expression not reproducible] 9: [mathematical expression not reproducible] 10: [mathematical expression not reproducible] 11: [mathematical expression not reproducible] 12: end for 13: [mathematical expression not reproducible] 14: end for

The batch cross covariance given by (37) is exact but takes a large amount of computation. To find a recursive formula, one may combine [[??].sup.[i].sub.l] given in Algorithm 1 and (5) with respect to an auxiliary variable l, so [e.sup.[i].sub.l] can be rewritten as

[mathematical expression not reproducible]. (38)

Substitute given by (5) and (2) into (38), then transform it into

[mathematical expression not reproducible]. (39)

So [P.sup.[ij].sub.l] can now be computed iteratively over the recursive form of

[mathematical expression not reproducible] (40)

where l ranges from s to n, and the output is taken when l = n similarly. The initial value [p.sup.[ij].sub.s-1] is given in a short batch form by (37) as

[mathematical expression not reproducible]. (41)

Clearly the pseudocode is shown in Algorithm 2.

4.2. Decentralized OUFIR Filtering Algorithm. Based on the local filters [[??].sup.[i].sub.n] in Algorithm 1 and the filtering error cross covariance matrices between any two local filters in Algorithm 2, we can get the decentralized fusion filter weighted by matrices in the LUMV sense [7]. Introduce the synthetically unbiased estimator

[mathematical expression not reproducible], (42)

where a; are arbitrary matrices. From the unbiasedness assumption, we have [mathematical expression not reproducible]. Taking the expectation in (42) yields

[[alpha].sub.1] + [[alpha].sub.2] + ... + [[alpha].sub.r] = [bar.I] (43)

where [bar.I] = [[[I.sub.k] [I.sub.k] ... [I.sub.k]].sup.T.sub.krxk] and [I.sub.k] denotes k x k identity matrix. Then fusion estimation error is [e.sub.n] = [x.sub.n] - [[??].sup.o.sub.n] = [[summation].sup.r.sub.i=1] [[alpha].sub.i][e.sup.[i].sub.n]. Let [LAMBDA] = [[[[alpha].sub.1] [[alpha].sub.2] ... [[alpha].sub.r]].sub.krxk]; the error variance matrix of the fusion estimator is

P = E{[e.sub.n][e.sup.T.sub.n]} = [LAMBDA][summation][[LAMBDA].sup.T] (44)

where [summation] [??] [[[P.sup.[ij].sub.n]].sub.i,j=1,2,...,r] is an kr x kr symmetric positive definite matrix where [mathematical expression not reproducible] is from (40).

Algorithm 2: Iterative form for [P.sup.[ij].sub.n]. 1: for n = N - 1, N, ... do 2: m = n - N + 1, s = m + k 3: for i = 1 : r do 4: for j = 1 : r do 5: [mathematical expression not reproducible] 6: for l = s : n do 7: [mathematical expression not reproducible] 8: end for 9 [P.sub.[ij].sub.n] = [P.sub.[ij].sub.l] 10: end for 11: end for 12: end for Algorithm 3: Optimal decentralized fusion OUFIR filter. 1: [J.sub.i] = B[S.sub.i] [([C.sup.[i]]B[S.sub.i] + [R.sub.i]).sup.-1] 2: [[bar.A].sub.[i]] = (I - [J.sub.i][C.sub.[i]])A 3: [[bar.u].sub.[i].sub.n] = (I - [J.sub.i][C.sub.[i]])E[u.sub.n] + [J.sub.i][y.sup.[i].sub.n] 4: [[bar.w].sup.[i].sub.n] = (I - [J.sub.i][C.sub.[i]])B[w.sub.n] - [J.sub.i][v.sup.[i].sub.n] 5: Algorithm 1: General Iterative OUFIR Filtering Algorithm 6: Algorithm 2: Iterative Form for [P.sup.[ij].sub.n] 7: [summation] = [[P.sup.[ij].sub.n]], i = 1, 2 ... r 8: [bar.[LAMBDA]] = [([[bar.I].sup.T][[summation].sup.-1] [bar.I]).sup.-1] [[bar.I].sup.T][[summation].sup.-1] 9: [mathematical expression not reproducible]

To find the optimal weights [LAMBDA] under the constrain (43), applying the Lagrange multiplier approach, the target function is defined as

J = tr ([LAMBDA][summation][[LAMBDA].sup.T]) + 2[GAMMA]([LAMBDA][bar.I] - [I.sub.k]) (45)

where [GAMMA] = [([[lambda].sub.ij]).sub.kxk]. Set [partial derivative]J/[partial derivative][LAMBDA] = 0 and note that [[summation].sup.T] = [summation]; we have

2[LAMBDA][summation] + 2[GAMMA][[bar.I].sup.T] = 0 (46)

Combining (43) and (46) yields the matrix equation as

[mathematical expression not reproducible]. (47)

Noted [summation] is a symmetric positive definite matrix and using the formula of the inverse matrix [28] we have

[mathematical expression not reproducible] (48)

Pseudocode is summarized by Algorithm 3, and the corresponding structure is shown in Figure 2. As can be seen, each sensor subsystem estimates the state independently at time n. Next fusion layer determines optimal matrix weights. As a basic requirement, the proposed decentralized fusion filter will have better accuracy than any local one. Different from fusion KF, the proposed method does not relay on the initial information about [[??].sub.0], [P.sub.0], and [P.sup.[ij].sub.0], and it performs better robustness.

5. Simulation

In this section, we consider two basic models of moving target tracking system and 1 degree of freedom (1-DOF) torsion system. Fusion OUFIR filter is applied to above models under various conditions such as temporary modeling uncertainties and unknown noise statistics. Compared to fusion Kalman filter (conveniently FOUF and FKF denote as fusion OUFIR filter and fusion KF filter respectively next), it demonstrates better performance.

5.1. Moving Target Tracking System. With three sensors, a simplified model can be specified by (1) and (2) with E = [[1 1].sup.T], B = [[0.1 1].sup.T], and

[mathematical expression not reproducible] (49)

where i = 1, 2, 3, [tau] = 0.1, control signal [u.sub.n] = 0.1, and [C.sub.1] = [C.sub.2] = [C.sub.3] = [1 0]. And the measurement noise [v.sup.[i].sub.n] is represented as [v.sup.[i].sub.n] = [b.sub.i][w.sub.n] + [[xi].sup.[i].sub.n], which are correlated with the white process noise [[sigma].sup.2.sub.w] = 0.1. The correlation coefficients are [b.sub.1] = 0.5, [b.sub.2] = 0.8, and [b.sub.3] = 0.4. [[xi].sup.[i].sub.n] are white Gaussian noise with zero mean and variances [mathematical expression not reproducible], and [R.sub.[xi][3]] = 10, respectively, and are independent of [w.sub.n]. Easily we can get [R.sub.ij] and [S.sub.i] as

[mathematical expression not reproducible] (50)

[S.sub.i] = E {[w.sub.n] [([v.sup.[i].sub.n]).sup.T]} = Q[b.sub.i.sup.T] (51)

5.1.1. Estimation Accuracy with Correct Noise Variances. In this section, we assume that all the modeling parameters are known accurately and show that FOUF outperforms any local filter. The process was simulated over 500 points, and the estimation horizon length used is N = 20. The average root mean square error (RMSE) based on 100 Monte-Carlo (MC) runs is shown in Figure 3. The result proves that the proposed approach has better accuracy exactly.

5.1.2. Sensitivity to Incorrect Noise Variance. The knowledge of noise is typically not known exactly in practice, which does not guarantee optimality. To figure out the effect of noise variances, we give an idea about the FKF and FOUF immunity to noise statistics' errors in the worst case. Introducing error coefficient [beta], designed covariances are substituted here with [Q.sup.*] = [beta]Q.

[square root of tr([J.sub.n])] of FKF and FOUF over 100 MC runs for [beta] = 0.01, 0.1, 1, 10, 100 are shown in Figure 4 and also listed in Table 2. As can be seen, errors in the FKF grow rapidly with [beta] and even become unacceptable, while the proposed FOUF has much higher immunity against noise error when [beta] [less than or equal to] 1. Otherwise, performance of FOUF exhibits closely to FKFs.

5.1.3. Immunity to Errors in the Noise Statistics. To investigate the effect of temporary inaccurate statistics, the state noise variance is set to vary as follows: [[sigma].sup.2.sub.w] = 0.5 when 100 [less than or equal to] n [less than or equal to] 200 and [[sigma].sup.2.sub.w] = 0.1 otherwise. And invariable value [[sigma].sup.2.sub.w] = 0.1 is used in all the algorithms. The average RMSEs based on 100 MC runs are provides by Figure 5. One can see that FOUF has much higher immunity against errors of noise statistic than FKF. Besides it is noted that the first state is more obvious about robustness performance.

5.1.4. Robustness against Model Uncertainties. To identify the robustness of each tracker against temporary model uncertainties, we set

[mathematical expression not reproducible] (52)

where d = 0.1 when 100 [less than or equal to] n [less than or equal to] 300 and d = 0 otherwise. The process is simulated over 500 subsequent points with correct noise statistics. Comparing the FKF as a benchmark, Figures 6(a) and 6(b) illustrate the typical response under the model mismatch introduced. From figures, one infers that the FOUF's peak value is smaller and recovers to normal faster than the FKF's, particularly with the first state.

5.2. 1-DOF Torsion System. We consider 1-DOF torsion system [29] as the object whose schematic and parameters are provided by Figure 7 and Table 3. Referring from [29] of the equations of motion, the continuous-time model can be described by

[[partial derivative]x/[partial derivative]t] = Ax + Bu, (53)

where [mathematical expression not reproducible] and

[mathematical expression not reproducible]. (54)

Discretize this process with sampling time [T.sub.s] = 0.01s, and the discrete state-space model is achieved as

[x.sub.n+1] = G[x.sub.n] + F[u.sub.n] + [w.sub.n], (55)

where [w.sub.n] is introduced to describe the stochastic properties of states, [u.sub.n] is input signal, [mathematical expression not reproducible], and [mathematical expression not reproducible]. The angular positions of two loads are measured using three groups of sensors, and the observation equation is given by [y.sup.[i].sub.n] = [C.sup.[i]][x.sub.n] + [v.sup.[i].sub.n], where i = 1, 2, 3,

[mathematical expression not reproducible]. (56)

The correlation between state and measurement noises is described as [v.sup.[i].sub.n] = [b.sup.[i]][w.sub.n] + [[xi].sup.[i].sub.n], where the coefficients are

[mathematical expression not reproducible], (57)

The measurement noises variances are assigned as [mathematical expression not reproducible].

5.2.1. Imprecisely Defined Noise Statistics. With inaccurateness, the state noise variance is set to vary as follows: [mathematical expression not reproducible] when 150 [less than or equal to] n [less than or equal to] 300 and [mathematical expression not reproducible] otherwise. Meanwhile invariable value [mathematical expression not reproducible], is used in the filter. The average RMSEs based on 100 MC runs are provided by Figure 8. It suggests again that FOUF performs closely to the FKF under the correct noise, but it has better robustness against noise errors than FKF.

5.2.2. Temporary Model Uncertainty. Similarly, we assume that [B.sub.2] changes unpredictably as [B.sub.2] = 0.5 when 150 [less than or equal to] n [less than or equal to] 300 and [B.sub.2] = 0.015 otherwise, which makes A vary simultaneously. Average RMSEs produced by FKF and FOUF over 100 Monte-Carlo simulations are shown Figure 9.

Clearly the result obtained in this case coincides with that in the previous model. At this point, one can conclude that proposed FOUF is less influenced by noise and more robust against temporary model mismatch.

6. Conclusion

Based on optimal fusion weighted by matrices in the LUMV sense, this paper presents a decentralized fusion OUFIR filter for discrete-invariant control system. Fundamentally the FOUF has higher accuracy than each local filter. The simulation results show that the proposed algorithm maintains the advantages of OUFIR filter and suggests lowlier sensitive to noise statistics and more robust against modeling uncertainties.

https://doi.org/10.1155/2018/1505137

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This paper is supported by the National Natural Science Foundation of China (NSFC 61773183).

References

[1] D. Simon, Optimal State Estimation: Kalman, H<x>, and Nonlinear Approaches, Wiley-Interscience, 2006.

[2] V. R. Algazi, M. Suk, and C.-S. Rim, "Design of almost minimax FIR filters in one and two dimensions by WLS techniques," IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 33, no. 6, pp. 590-596,1986.

[3] R. E. Kalman, "A new approach to linear filtering and prediction problems," Journal of Basic Engineering, vol. 82, no. 1, pp. 34-45, 1960.

[4] S. Zhao, Y. S. Shmaliy, B. Huang, and F. Liu, "Minimum variance unbiased FIR filter for discrete time-variant systems," Automatica, vol. 53, pp. 355-361, 2015.

[5] 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, no. 3, pp. 545-551, 2002.

[6] A. H. Jazwinski, Stochastic Process and Filtering Theory, Academic Press, New York, NY, USA, 1970.

[7] S.-L. Sun and Z.-L. Deng, "Multi-sensor optimal information fusion Kalman filter," Automatica, vol. 40, no. 6, pp. 1017-1023 (2005), 2004.

[8] T. Tian, S. Sun, and N. Li, "Multi-sensor information fusion estimators for stochastic uncertain systems with correlated noises," Information Fusion, vol. 27, pp. 126-137, 2016.

[9] N. A. Carlson, "Federated square root filter for decentralized parallel processes," IEEE Transactions on Aerospace and Electronic Systems, vol. 26, no. 3, pp. 517-525, 1990.

[10] K. C. Chang, R. K. Saha, and Y. Bar-Shalom, "On optimal track-to-track fusion," IEEE Transactions on Aerospace and Electronic Systems, vol. 33, no. 4, pp. 1271-1276, 1997.

[11] S.-L. Sun, "Multi-sensor optimal information fusion Kalman filters with applications," Aerospace Science and Technology, vol. 8, no. 1, pp. 57-62, 2004.

[12] E. Song, Y. Zhu, J. Zhou, and Z. You, "Optimal Kalman filtering fusion with cross-correlated sensor noises," Automatica, vol. 43, no. 8, pp. 1450-1456, 2007.

[13] N. Li, S. Sun, and J. Ma, "Multi-sensor distributed fusion filtering for networked systems with different delay and loss rates," Digital Signal Processing, vol. 34, pp. 29-38, 2014.

[14] R. Caballero-Aguila, I. Garcia-Garrido, and J. Linares-Perez, "Optimal fusion filtering in multisensor stochastic systems with missing measurements and correlated noises," Mathematical Problems in Engineering, vol. 2013, Article ID 418678, 14 pages, 2013.

[15] X. Yang, W.-A. Zhang, L. Yu, and K. Xing, "Multi-rate distributed fusion estimation for sensor network-based target tracking," IEEE Sensors Journal, vol. 16, no. 5, pp. 1233-1242, 2016.

[16] Y. S. Shmaliy, "Linear optimal FIR estimation of discrete time-invariant state-space models," IEEE Transactions on Signal Processing, vol. 58, no. 6, pp. 3086-3096, 2010.

[17] Y. S. Shmaliy, "An iterative Kalman-like algorithm ignoring noise and initial conditions," IEEE Transactions on Signal Processing, vol. 59, no. 6, pp. 2465-2473, 2011.

[18] S. Zhao, Y. S. Shmaliy, F. Liu, O. Ibarra-Manzano, and S. H. Khan, "Effect of embedded unbiasedness on discrete-time optimal FIR filtering estimates," EURASIP Journal on Advances in Signal Processing, vol. 2015, no. 1, Article ID 83, 2015.

[19] S. Zhao, Y. S. Shmaliy, and F. Liu, "Fast Kalman-like optimal unbiased FIR filtering with applications," IEEE Transactions on Signal Processing, vol. 64, no. 9, pp. 2284-2297, 2016.

[20] C. K. Ahn, S. Zhao, Y. S. Shmaliy, and H. Li, "On the 12-1ot and Hot performances of continuous-time deadbeat H2 FIR filter," IEEE Transactions on Circuits and Systems II: Express Briefs, vol. PP, no. 99, pp. 1-1, 2017.

[21] C. K. Ahn, Y. S. Shmaliy, P. Shi, and Y. Zhao, "Receding-Horizon 12-1ot FIR filter with embedded deadbeat property," IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 64, no. 2, pp. 211-215, 2017.

[22] C. K. Ahn and Y. S. Shmaliy, "New receding horizon fir estimator for blind smart sensing of velocity via position measurements," IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 65, no. 1, pp. 135-139, 2018.

[23] Y. Xu, Y. S. Shmaliy, Y. Li, and X. Chen, "UWB-based indoor human localization with time-delayed data using EFIR filtering," IEEE Access, vol. 5, pp. 16676-16683, 2017.

[24] R. Rasoulzadeh and A. M. Shahri, "Accuracy improvement of a multi-MEMS inertial measurement unit by using an iterative UFIR filter," in Proceedings of the 25th European Navigation Conference, ENC 2017, pp. 279-286, che, May 2017.

[25] S. H. Han, W H. Kwon, and P. S. Kim, "Quasi-deadbeat minimax filters for deterministic state-space models," Institute of Electrical and Electronics Engineers Transactions on Automatic Control, vol. 47, no. 11, pp. 1904-1908, 2002.

[26] W. H. Kwon and S. H. Han, Receding Horizon Control: Model Predictive Control for State Models, Springer Science & Business Media, 2006.

[27] Y. S. Shmaliy, S. H. Khan, S. Zhao, and O. Ibarra-Manzano, "General unbiased FIR filter with applications to GPS-based steering of oscillator frequency," IEEE Transactions on Control Systems Technology, vol. 25, no. 3, pp. 1141-1148, 2017.

[28] N. Xu, Stochastic Signal Estimation and System Control, Beijing Industry University Press, Beijing, China, 2001.

[29] G. Fuchs, V. Satek, V. Vopenka, J. Kunovsky, and M. Kozek, "Application of the Modern Taylor Series Method to a multitorsion chain," Simulation Modelling Practice and Theory, vol. 33, pp. 89-101, 2013.

Xuefeng Fan and Fei Liu [ID]

Key Laboratory of Advanced Process Control for Light Industry, Jiangnan University, Wuxi 214122, China

Correspondence should be addressed to Fei Liu; fliu@jiangnan.edu.cn

Received 2 November 2017; Accepted 1 February 2018; Published 3 June 2018

Academic Editor: Rajesh Khanna

Caption: Figure 1: Operation time diagrams of the FIR structures.

Caption: Figure 2: The structure of decentralized OUFIR filter.

Caption: Figure 3: RMSEs produced by local and fusion filters: (a) the first state and (b) the second state.

Caption: Figure 4: Average [square root of tr([J.sub.n])] against system noises uncertainties [beta] = 0.01, 0.1, 1, 10, 100 between FKF and FOUF.

Caption: Figure 5: RMSEs produced by FKF and FUF with imprecisely defined noise statistics: (a) the first state and (b) the second state.

Caption: Figure 6: Estimation errors by FKF and FOUF filter against temporary model uncertainty in a gap of 100 [less than or equal to] n [less than or equal to] 150: (a) the first state and (b) the second state.

Caption: Figure 7: Schematic of 1-DOF Torsion system.

Caption: Figure 8: RMSEs produced by FKF and FOUF with imprecisely defined noise statistics: (a) the first state and (b) the second state.

Caption: Figure 9: RMSEs produced by FKF and FOUF with temporary model uncertainty: (a) the first state and (b) the second state.

Table 1: Computation of the initial state. 1: m = n- N+1, s = m + k 2: [M.sub.s-2] = [[bar.B].sup.[i].sub.s-2,m], [Q.sub.s-2,m] [([F.sup.[i].sub.s-2,m]).sup.T] 3: [M.sub.s-1] = [[bar.B].sup.[i].sub.s-1,m], [Q.sub.s-1,m] [([F.sup.[i].sub.s-1,m]).sup.T] 4: [U.sub.s-1] = [[bar.B].sup.[i].sub.s-1,m], [Q.sub.s-1,m] [([[bar.B].sup.[i].sub.s-1,m]).sup.T] 5: [P.sub.s-1] = [U.sub.s-1,m] - [[bar.A].sup.[i]][M.sub.s-2] [Z.sup.-1.sub.w+v,s-2][([M.sub.s-2]).sup.T] [([[bar.A].sup.[i]]).sup.T] 6: [N.sub.s-1] = [[[(C.sup.[i].sub.s-1,m]).sup.T] [Z.sup.-1.sub.w+v,s-1][C.sup.[i].sub.s-1,m]].sup.-1] 7: [G.sub.s-1] = [A.sup.[i].sub.s-1,m] - [[bar.A].sup.[i]][M.sub.s-2][Z.sup.-1.sub.w+v,s-2] [C.sup.[i].sub.s-1,m] 8: [mathematical expression not reproducible] Table 2: Performance of FOUF and FKF with respect to different [beta]. [beta] 0.01 0.1 1 10 100 FKF 4.42 1.86 1.24 1.28 2.53 FOUF 1.58 1.47 1.28 1.29 2.58 Table 3: Physical parameters of 1-DOF torsion system, where EMI is equal to equivalent moment of inertia, and EVD is equivalent viscous damping. Constant Value Unit [K.sub.s] 1.0 N.m/rad [J.sub.1] 2.2e-4 Kg.[m.sup.2] [J.sub.2] 5.45e-4 Kg.[m.sup.2] [B.sub.1] 0.015 N.m.s/[rad.sup.2] [B.sub.2] 0.015 N.m.s/rad [[theta].sub.1] -- rad [[theta].sub.1] -- rad Constant Description [K.sub.s] Stiffness of flexible coupling [J.sub.1] EMI at the rigid load shaft [J.sub.2] EMI at the torsion load shaft [B.sub.1] EVD of rigid load [B.sub.2] EVD of torsion load [[theta].sub.1] angular positions of rigid load [[theta].sub.1] angular positions of torsion load

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Fan, Xuefeng; Liu, Fei |

Publication: | Journal of Electrical and Computer Engineering |

Date: | Jan 1, 2018 |

Words: | 5745 |

Previous Article: | A Mixed-Signal Programmable Time-Division Power-On-Reset and Volume Control Circuit for High-Resolution Hearing-Aid SoC Application. |

Next Article: | Uniform Local Binary Pattern for Fingerprint Liveness Detection in the Gaussian Pyramid. |

Topics: |