Optimal State Estimation

Chapter 9.0 - Optimal Smoothing

CHAPTER 9

Optimal smoothing

In a post mortem (after the fact) analysis, it is possible to wait for more observations to accumulate. In that case, the estimate can be improved by smoothing.

—Andrew Jazwinski [Jaz70, p. 143]


In previous chapters, we discussed how to obtain the optimal a priori and aposteriori state estimates. The a priori state estimate at time k, 09_00_Optimal_State_Estimation-1.jpg, is the state estimate at time k , based on all the measurements up to (but not including) time k. The a posteriori state estimate at time 09_00_Optimal_State_Estimation-1.jpg, is the state estimate at time k based on all the measurements up to and including time k:

09_00_Optimal_State_Estimation-3.jpg

There are often situations in which we want to obtain other types of state estimates. We will define09_00_Optimal_State_Estimation-1.jpgk, j as the estimate of xk given all measurements up to and including time j. With this notation, we see that

09_00_Optimal_State_Estimation-4.jpg

Now suppose, for example, that we have recorded measurements up to time index 54 and we want to obtain an estimate of the state at time index 33. Our theory in the previous chapters tells us how to obtain 09_00_Optimal_State_Estimation-5.jpg or 09_00_Optimal_State_Estimation-5.jpg, but those estimates only use the measurements up to and including times 32 and 33, respectively. If we have more measurements (e.g., measurements up to time 54) it stands to reason that we should be able to get an even better estimate of x33. This chapter discusses some ways of obtaining better estimates.

In another scenario, it may be that we are interested in obtaining an estimate of the state at a fixed time j. As measurements keep rolling in, we want to keep updating our estimate09_00_Optimal_State_Estimation-5.jpgj .In other words, we want to obtain 09_00_Optimal_State_Estimation-5.jpgj, j+1, 09_00_Optimal_State_Estimation-5.jpgj, j+2, … . This could be the case, for example, if a satellite takes a picture at time j. In order to more accurately process the photograph at time j we need an estimate of the satellite state (position and velocity) at time j. As the satellite continues to orbit, we may obtain additional range measurements of the satellite, so we can continue to update the estimate of xj and thus improve the quality of the processed photograph. This situation is called fixed-point smoothing because the time point for which we want to obtain a state estimate (time j in this example) is fixed, but the number of measurements that are available to improve that estimate continually changes. Fixed-point smoothing is depicted in Figure 9.1 and is discussed in Section 9.2.

09_00_Optimal_State_Estimation-8.jpg

Figure 9.1 Fixed-point smoothing. We desire an estimate of x4.Up until k = 4, the standard Kalman filter operates. At k = 4, we have 09_00_Optimal_State_Estimation-1.jpg= 09_00_Optimal_State_Estimation-1.jpg4, 4, which is the estimate of x4 based on measurements up to and including y3. As time progresses, we continue to refine our estimate of x4 based on an increasing number of measurements. At time k = N, we have 09_00_Optimal_State_Estimation-1.jpg4, N which is the estimate ofx4 based on measurements up to and including time N - 1.


Another type of smoothing is fixed-lag smoothing. In this situation, we want to obtain an estimate of the state at time (k - N) given measurements up to and including time k, where the time index k continually changes as we obtain new measurements, but the lag N is a constant. In other words, at each time point we have N future measurements available for our state estimate. We therefore want to obtain 09_00_Optimal_State_Estimation-1.jpgk-N, k for k=N, N+1, …, where N is a fixed positive integer. This could be the case, for example, if a satellite is continually taking photographs that are to be displayed or transmitted N time steps after the photograph is taken. In this case, since the photograph is processed N time steps after it is taken, we have N additional measurements after each photograph that are available to update the estimate of the satellite state and hence improve the quality of the photograph. Fixed-lag smoothing is depicted in Figure 9.2 and is discussed in Section 9.3.

The final type of smoothing is fixed-interval smoothing. In this situation, we have a fixed interval of measurements (y1, y2, …, yM) that are available, and we want to obtain the optimal state estimates at all the times in that interval. For each state estimate we want to use all of the measurements in the time interval. That is, we want to obtain 09_00_Optimal_State_Estimation-1.jpg0, M , 09_00_Optimal_State_Estimation-1.jpg1, M , …,09_00_Optimal_State_Estimation-1.jpgM, M .This is the case when we have recorded

09_00_Optimal_State_Estimation-10.jpg

Figure 9.2 Fixed-lag smoothing. We desire an estimate of the state at each time step based on measurements two time steps ahead. After processing y2, we form the estimate 09_00_Optimal_State_Estimation-1.jpg0, 2, which is the estimate of x0 based on measurements up to and including y2. Similarly, 09_00_Optimal_State_Estimation-1.jpg1, 3 is the estimate of x1 based on measurements up to and including y3.


some data that are available for post-processing. For example, if a manufacturing process has run over the weekend and we have recorded all of the data, and now we want to plot a time history of the best estimate of the process state, we can use all of the recorded data to estimate the states at each of the time points. Fixed-interval smoothing is depicted in Figure 9.3 and is discussed in Section 9.4.

09_00_Optimal_State_Estimation-11.jpg

Figure 9.3 Fixed-interval smoothing. We desire an estimate of the state at each time step based on all of the measurements in some interval. After processing all of the measurements from y1 to yM, we form the estimate 09_00_Optimal_State_Estimation-1.jpg0, M which is the estimate of x0 based on all the measurements. Similarly, 09_00_Optimal_State_Estimation-1.jpg1, M is the estimate of x1 based on all the measurements.


Our derivation of these optimal smoothers will be based on a form for the Kalman filter different than we have seen in previous chapters. Therefore, before we can discuss the optimal smoothers, we will first present an alternate Kalman filter form in Section 9.1.

UNLIMITED FREE
ACCESS
TO THE WORLD'S BEST IDEAS

SUBMIT
Already a GlobalSpec user? Log in.

This is embarrasing...

An error occurred while processing the form. Please try again in a few minutes.

Customize Your GlobalSpec Experience

Category: Engineering Consulting Services
Finish!
Privacy Policy

This is embarrasing...

An error occurred while processing the form. Please try again in a few minutes.