Research Article | Open Access

# Receding Horizon Unbiased FIR Filters and Their Application to Sea Target Tracking

**Academic Editor:**Ai-Guo Wu

#### Abstract

Finite impulse response (FIR) state estimation algorithms have been much discussed in literature lately. It is well known that they allow overcoming the Kalman filter divergence caused by modeling uncertainties. In this paper, new receding horizon unbiased FIR filters ignoring noise statistics for time-varying discrete state-space models are proposed. They have the following advantages. First, the proposed filters use only known means of state vector components at starting points of sliding windows. This allows us to take into account priory statistical information (on average) about specified movements of the system. Second, the iterative version of the filter has a Kalman-like form. Besides, its initialization does not include a training cycle in a batch form. Such filters may have a wide range of applications. In this paper, position and speed estimation of sea targets using angle measurements in azimuth and elevation is considered as an example.

#### 1. Introduction

Finite impulse response (FIR) state estimation algorithms can be considered as an alternative to the Kalman filter (KF). These algorithms allow us to overcome its divergence caused by model inaccuracy (incorrectly chosen models, temporary perturbations, and errors in the setting of the noise statistics) [1â€“3]. Divergence prevention is achieved by reducing the impact of the data or rejection the data out of the sliding window. In this context, the KF can be considered as an estimation algorithm with infinite impulse response (IIR).

Quite a number of methods have been offered to construct FIR algorithms. In [4], the FIR filter for a linear homogeneous system is proposed. The optimal unbiased FIR (OUFIR) and receding horizon algorithms are proposed in [2, 3, 5â€“7]. Initial conditions at starting points of sliding windows are assumed to be diffuse random variables or unknown arbitrary values. In [7], the receding horizon optimal unbiased FIR (RHOUFIR) filter which obtains information about the window initial conditions from finite observations is suggested. Another approach ensuring optimality and unbiasedness in a finite number of steps is described in [8, 9]. The RHOUFIR filter suggested by the authors uses known statistical information for parts of state vector components at starting points of sliding windows. Within the framework of covariance analysis this allows us to take into account priory statistical information about random biases, trends, and specified movements of the system. Training tasks of neural and neurofuzzy networks are another example of possible application of this filter [9].

In [10â€“14], the unbiased FIR (UFIR) filters ignoring the noise statistics of the process, measurements, and the statistics of state vector initial conditions are developed. Computer modeling shows that the UFIR filters can be more robust than the KF or the OUFIR filter if the size of the sliding window is properly fitted. There are several general approaches to select based on a heuristic choice, via the mean square errors, via real measurements, and via the use of bandlimited property of signals [13]. However, the proposed recursive filters have two disadvantages. First, a butch form of the algorithm is needed for the initialization of the recursive UFIR filter during the learning cycle which may be inconvenient or even problematic in some cases (e.g., in case of gaps in the observations, estimation parameters of nonlinear systems, and nonstationary processes). Second, these filters do not use priorly known information for initialization at starting points of sliding windows.

In this paper, the receding horizon UFIR (RHUFIR) filters using only known means of state vector components at starting points of sliding windows are suggested. They allow us to take into account (on average) priory statistical information about the system. Furthermore, for the iterative filter initialization no training cycle in a batch form is required. While the implementation of the RHOUFIR filter described in [8, 9] may be difficult due to the lack of sufficiently accurate statistical information about system, the algorithms under consideration may produce quite acceptable estimates. Theoretically this can be justified by a selection of the optimal or sub-optimal values of using the algorithms described in [13]. Another important feature is that the RHUFIR filter is computationally less demanding in comparison with the RHOUFIR. This is important, for example, in the tasks of neural and neurofuzzy networks training.

As a possible application of the proposed filters, the position and speed estimation of sea targets using angle measurements in azimuth and elevation from the observer-shipâ€™s video camera is considered. The general approach to solve this problem is to use Kalman filtering methods [15]. The difficulties arising due to the specificity of both the problem itself and the methods commonly used are well known. First, measurement models may be highly nonlinear near the state estimate [16]. Second, the object state may not be observable if the observer does not maneuver in a special way [17, 18]. Third, there is a large initial uncertainty about the position, speed, and acceleration of the observed object. All this can lead to a divergence of the KF and its various modifications [19â€“21].

The work is organized as follows. The RHUFIR state estimation problem is formulated in Section 2. Sections 3 and 4 concern the batch and iterative RHUFIR filters, respectively. Section 5 contains a comparative analysis of the RHUFIR and RHOUFIR filters. A diffuse approach to FIR estimation and its connection with the RHOUFIR filter is considered in Section 6. The position and speed estimation of sea objects is considered in Section 7. The conclusions are presented in Section 8.

#### 2. Problem Statement

Let us consider a linear discrete time-variant system of the formwhere is the state vector, is the measured output vector, and are random processes with zero means, i.e., , , , , , are known matrices of appropriate dimensions, .

Consider the discrete intervals , (the sliding windows), where is the horizon length. Let us assume that the following conditions are met.

A1: Means for arbitrary components of the state vector at the starting points of the sliding windows are known. A priori information on remaining components of , is absent and they are either unknown constants or random variables statistical characteristics of which are unknown. The value is used further if a priori information about all components of is absent.

A2: If then without loss of generality it is assumed that the state vector elements are arranged so that means , , are known, where , .

It is required to find the linear state estimate of the convolution-based batch formfrom measurements and its recursive representation for calculation, where the matrices and do not depend on a priori information about the unknown vector , , , and , if . The estimate is found under the conditions that it is unbiased

The first component in (2) allows to take into account (on average) a priori information about the system (1) at starting points of sliding windows , **. **This may be necessary if the pair is not observable but the estimate must be unbiased after processing the finite number of observations. Typical examples are models including random biases, trends and specified movements of the system (see the example in Section 7.3). The relation , , , is the simplest case describing a random bias with a known mean. Note that the inference of the UFIR filter in [11] for is based on the assumption that the pair is observable. Note also that this assumption is not required for the KF. The need for the first component in (3) may also arise in the tasks of neural and neurofuzzy networks training which have a separable structure [22]. It is assumed that the linearly incoming parameters are unknown and in addition to the training set a priori information only about the nonlinearly incoming parameters is given [9]. This information is obtained from the distribution of a generating sample, a training set, or some linguistic information. The RHUFIR filter is used to estimate the parameters by linearizing networks relations in the vicinity of the last estimate.

#### 3. Batch Receding Horizon UFIR Filter

The RHUFIR filter in a butch form is specified by the following theorem.

Theorem 1. *Let the condition be fulfilled. Then the unbiased estimate of the system state (1) in the batch form is determined by the expression (2), where, is arbitrary matrix, is the Moore Penrose inversion of , is th unit vector, and is the identity matrix of the size .*

*Proof. *First, we show that the estimation problem is equivalent to a dual control problem. Consider an auxiliary linear system of the formwhere is the state matric and is a control. Suppose that there is a control bringing the system (12) in a state satisfying the condition for any . Let us show that then there is an unbiased state estimate of the system (1) of the form (2), where and is determined by (6).

Using the identity(1) and (12) giveTaking into account this expression, the estimation error at the moment can be presented in the form where , and is the vector row with zero elements.

Averaging the left and right sides of this expression and using the condition get Now we find the control providing a solution to the dual problem. Iterating (12) giveswhere the matrix is determined by the system Using the boundary condition , we obtain with help of (19) the linear equations system We find a solution of this system in the form where is a constant unknown matrix and is arbitrary matric function satisfying the condition Substituting (22) into (21) and taking in account (23) gives provided , where Using the expressionswe obtainBut as and (10) is the general solution of the linear homogeneous system (23) then it implies (5) and (6).

*Comment 1. *The condition is observable for , and implies (3).

*Comment 2. *The filter takes a particularly simple form for and . We find from Theorem 1 thatwhere These relations are equal to the algorithm from [11] up to the notation used if to ignore the one-step prediction.

#### 4. Iterative Form of Batch Receding Horizon UFIR Filter

We will need the following assertion.

Lemma 2. *Let the condition (4) be true. Then the control bringing the system (12) in a state satisfying the condition can be presented in the form where*

*Proof. *Let us show that the control (29) indeed solves our problem. Consider the system (12) under the action of the controls (22) with . We show that . Substituting with in (12) gives From comparison of right parts of (31) and (32), it follows that the solutions of the systems really coincide if , i.e.,Iterating the system (32), we successively findThe expressions (34) and (35) follow from the identitiesLet us transform (35) using the orthogonal decomposition , where Let be linearly independent columns of the matrix . Using skeletal decomposition [23] yields , where , , and are a set of matrices that have rank . Since is the Gram matrix generated by linearly independent rows of the matrix and thenSubstituting these expressions in (35) and using the representation , we findIt follows from (37) that , where is some rectangular matrix. Substituting this expression in (39) givesThe RHUFIR filter is specified by the following theorem.

Theorem 3. *The state of the system is the unbiased estimate of the system state (1) for , , where *

*Proof. *Put in (2) where the matric functions are identified in Theorem 1. This implies the unbiasedness of the estimate at the moment and the relation From Lemma 2 it follows that In view of this where By (48) we findUsing the identitieswe transform this expression to the form (41) Putting in (48) gives .

*Comment 1. *The filter is equal to the algorithm from [11] up to the used notation for if to ignore the one-step prediction and to use Theorem 1 to initialize it:where Consider the system for the state transition matrix of the homogeneous part of the UFIR filter (41)

*Consequence 1. *The following representations are true:Sincethe statement follows from (34) and (35).

*Consequence 2. *Let us show using (56) that the estimate is unbiased indeed. Consider the system of equations for the expectation of the estimation error where is the arbitrary vector. We find from (56) for . Similarly, it is verified that for , , and , where is the arbitrary vector.

*Consequence 3. *Let the system (1) be time-invariant, , and . Then the representation for the state transition matrix has the formThe statement can be obtained from (56) with the help of elementary calculations.

#### 5. Comparison of RHUFIR and RHOUFIR Filters

Let us assume that the following conditions are held.

B1: and are the uncorrelated random processes with zero means and known covariance matrices , are known matrices of appropriate dimensions.

B2: Means and covariances are known for arbitrary components of the state vector at the starting points of the sliding windows. A priori information on remaining components of , is absent and they are either unknown constants or random variables the statistical characteristics of which are unknown. The value is used further if a priori information about all components of is absent.

B3: If then without loss of generality it is assumed that the state vector elements are arranged so that means and covariancesare known, where , .

B4: If are random variables then they are uncorrelated with , , and for .

It is required to find the unbiased linear state estimate of (1) minimizing the criterion . This estimate is called the optimal estimate of the system state (1).

Theorem 4 ([8, 9]). *(1) The optimal estimate of (the RHOUFIR filter) is determined by the following relations for , :**Prediction: **Correction:**(2) The estimation error covariance matrix for is given by the expression *

*Comment 1. *Formally, the RHUFIR filter follows from the RHOUFIR filter by specifying , , and . At the same time, the statement of the problem and the derivation of the UFIR filter do not rely on these assumptions.

*Comment 2. *The relation (71) is the decomposition of the error covariance matrix into two terms. The first one takes into account the known statistical information about initial conditions, process, and measurement noises covariance matrices on the estimate accuracy. The second one additionally reflects the impact of unknown initial conditions at starting points of sliding windows.

*Comment 3. *Similarly to the covariance matrix, the decomposition is true for the UOFIR filter gain. Indeed, it follows from (64) that , where the term coincides with the KF gain taking in account the known statistical information and the term reflects the additional effect connected with unbiasedness of the RHOUFIR filter.

#### 6. Receding Horizon FIR Filter with Diffuse Initialization

Consider an alternative approach to constructing the FIR estimator which is important in view of possible applications. Now, we assume that the condition B1 from the previous section and the following conditions are met.

C1: Means and covariances are known for arbitrary components of the state vector at the starting points of the sliding windows.

C2: If then are treated as random variables with zero mean and covariance matrix proportional to the large parameter , i.e., , , where is an arbitrary positive definite matrix.

C3: If then without loss of generality it is assumed that the state vector elements are arranged so that means and covariancesare known. A priori information on remaining components of , , is absent and they are treated as random variables with zero mean and covariance matrix proportional to the large parameter , i.e., where , is a large parameter.

C4: The random vector is not correlated with and for , .

C5: The random variables for , are uncorrelated with variables , .

It is required to find the limit relations for the KF as and to study their properties. In [8, 9], it was shown that as these relations are coincide with (63)â€“(70).

The conditions C1â€“C5 formalize the standard approach to the implementation of the KF for large initial uncertainty of the system state. At the same time, it is well known [9, 24] that the large values of can lead to a divergence of the KF. In [8, 9], the influence of the large values of on the KF divergence was studied. In particular, let be the error connected with calculations of . Then For the matrix becomes proportional to the large parameter and so divergence is possible even if the continuity condition of the matrix pseudoinversion finding is performed. Thus, the use of the limiting relations for the filter allows avoiding divergence of the FIR filters with the diffuse initialization.

Note that, for the special case , the relations (63)â€“(70) (the IIR filter) can be used for the state estimation in absence or incompleteness of the statistical information about the initial conditions of (1). Following [8, 9], we call them the KF with diffuse initial conditions (DKF).

#### 7. Application of Receding Horizon FIR Filters to Sea Target Tracking Problem

##### 7.1. Models of Objects and Measurements

To illustrate the capability of the proposed FIR algorithms, we use the following scenario. With the help of the camera installed on the ship, the bearing (), and elevation () angles to the target (Figure 1) are measured.

The expressions for and taking into account measurement errors have the following form:where and are the positions of the target, is altitude of the camera, and are the centered uncorrelated white noises with the variances and , respectively, and is the 4-quadrant arctangent.

It is assumed that the variability of is caused by a rolling with the dominant frequency and amplitude. More exactly, it is assumed that , where is approximately known and are unknown values to the observer and, is the sampling step. The function is interpreted as an uncontrolled disturbance in the elevation angle measurement channel.

The ship and the measurements model are described by a nonlinear discrete system of the form where is the state vector, is the measured output vector, and are random processes, is the control, and and are known vector functions.

The movement of the target to the observer is unknown and we rely on the common approach to describe its behavior consisting in the following [15, 24]. It is assumed that the target movement can be specified by means of a stochastic system of the form where is the state vector and is the random process. A priori information about the initial position of the target is absent and the right-hand parts of the system (79) and the state vector and the random process are chosen using the specifics of the problem under consideration.

Special cases of models (77), (78), and (79) are as follows [15, 24].

The ship moves at a constant fixed velocity in accordance with where the linear speed () and the heading angle () of the ship are known values and is the sampling period. The coordinates of the ship in the horizontal plane (,) are measured and the observation models have the formwhere and are the centered uncorrelated white noises with the known variances and , respectively, , , , and .

The motion of the target is described bywhere and are the target velocity projections on coordinates axes and , , and are the accelerations projections (the centered uncorrelated white noises with the variances , respectively).

##### 7.2. Models of Pseudo Measurements

Let us find expressions for the pseudo measurements model using the relations (75) and (76). The general idea [25] is to represent the nonlinear measurement model in the next pseudo linear form , where is the pseudo measurement vector, is the known function of the true measurements of, and is the pseudo measurement errors. The KF is used with , and , where is the predicted state value of x.

Rewrite (75) and (76) in the following form: We find from (84)Substitution of this expression into (85) givesand using (86) and (87), we getAs thenLinearizing the right-hand parts of (90) in a neighborhood of the points , , we findwhere and are correlated random processes defined by expressions It follows the relations for the pseudo measurements model whereLinearizing the right-hand sides of (93) and (94) in a neighborhood of the point , we also get