Optimal State Estimation

Chapter 9.5 - Summary

9.5 SUMMARY

In this chapter we derived the optimal smoothing filters. These filters, sometimes called retrodiction filters [BarOl], include the following variants.

  • j, k = E (xj y1, …, yk-1 ) (k j) is the output of the fixed-point smoother. In this filter we find the estimate of the state at the fixed time j when measurements continue to arrive at the filter at times greater than j. The time index j is fixed while k continues to increase as we obtain more measurements.

  • k-N, k = E (xk-N y1, …, yk ) for a fixed N is the output of the fixed-lag smoother. In this filter we find the estimate of the state at each time k while using measurements up to and including time (k + N). The time index k varies while N remains fixed.

  • k, N = E (xk y1, …, yN ) for a fixed N is the output of the fixed-interval smoother. In this filter we find the estimate of the state at each time k while using measurements up to and including time N. The time index k varies while the total number of measurements N is fixed. The two formulas we derived for this type of smoothing included the forward-backward smoother and the RTS smoother.

Just as steady-state filters can be used for standard filtering, we can also derive steady-state smoothers to save computational effort [Gel74]. An early survey of smoothing algorithms is given in [Med73].

PROBLEMS

Written exercises

9.1     Prove or disprove the following conjecture: The trace of the inverse of a matrix is equal to the inverse of the trace of the matrix.

9.2     Show that (A + B)-1 = B -1(AB -1+ I )-1.

9.3     Derive Equation (9.83).

9.4     Consider a scalar system with F = 1, H = 1, and R = 2Q.

  1. What is the steady-state value of the a priori estimation-error covariance ?

  2. Suppose that after the Kalman filter has reached steady state, the fixed-point smoother begins to operate. Find a closed-form solution to the covariance of the smoothed estimate Πk as a function of the time index k. What is the limiting value of Πk as k → ∞?

9.5      Repeat Problem 9.4 for the case R = 12Q. What is the percent improvement in the estimation-error covariance due to smoothing? Explain why the percent improvement due to smoothing for this case differs in the way that it does from the results of Problem 9.4.

9.6      Consider a scalar system with F = 1, H = 1, and R = 2Q. Suppose that the fixed-lag smoother for this system is in steady state so that = , Lk+1, i = Lk, i , = , = , for i = 1, …, N+1. Find closed-form expressions for, Lk, i , , and as functions of i. What is the limit as i → ∞ of Lk, i , , and?

9.7      Suppose you have a fixed-lag smoother as shown in Equation (9.43) that is in steady state. How do the eigenvalues of the fixed-lag smoother relate to the eigenvalues of the standard Kalman filter? What do you conclude about the stability of the fixed-lag smoother?

9.8      Solve Equation (9.10) for (yk - Hk) [assuming that ρ( Lk ) = r, where r is the number of measurements in the system]. Substitute the resulting expression for (yk - Hk) in the fixed-lag smoother equation for k+1-i, kto show that the smoothed state estimate can be driven by the state estimates without any input from the measurements [And79].

9.9 Suppose that f and b are unbiased estimates of x, and = Kff + Kbb . Show that ifis an unbiased estimate of x, then we must have Kf + Kb = I.

9.10     Consider a scalar system with F = 1, if = 1, and R = 2Q. Use the forward-backward smoother of Section 9.4.1 to find the steady-state value of the covariance of the smoothed state estimate.

9.11     Consider a scalar system with F = 1, H = 1, and R = 2Q. Use the RTS smoother of Section 9.4.2 to find the steady-state value of the covariance of the smoothed state estimate.

9.12     Consider a scalar system with F = 1, H = 1, and R = 2Q. Suppose that the forward filter has reached steady state. Use the RTS smoother of Section 9.4.2 to find the covariance of the smoothed state estimate for k =N, N - 1, N - 2, N - 3, and N - 4. At what point does the covariance of the smoothed state estimate get within 1% of its steady-state value?

9.13     Repeat Problem 9.12 for R = 12Q. How do you intuitively explain the quicker convergence of Pkto steady state?

9.14    Use the RTS smoother equations to show that constant states are not smoothable. That is, if F = I and Q = 0, then Pk = for all k.

Computer exercises

9.15 Consider the second-order system

09_05_Optimal_State_Estimation-16.jpg

where ω = 6 rad/s is the natural frequency of the system, and ζ = 0.16 is the damping ratio. The input w(t) is continuous-time white noise with a variance of 0.01. Measurements of the first state are taken every 0.5 s:

09_05_Optimal_State_Estimation-17.jpg

where v(tk) is discrete-time white noise with a variance of 10-4. The initial state, estimate, and covariance are

09_05_Optimal_State_Estimation-18.jpg
  1. Discretize the system equation.

  2.  Implement the discrete-time Kalman filter and the RTS smoother for 10 s (20 time steps). Plot the variance of the estimation error of the first state for the forward filter and for the RTS smoother on a single plot. Do the same for the second state. Why is the second state more smoothable than the first state?

9.16     Repeat Problem 9.15 with the continuous-time process noise w(t) having a variance of 1. How does this change the smoothability of the states?

9.17     Design a fixed-interval smoother for the system described in Problem 5.11 to estimate the state at each time on the basis of measurements at all 10 time steps.

  1. Plot the a posteriori covariance of the forward state estimate and the covariance of the smoothed state estimate as a function of time for both states.

  2. What are the percent improvements in the estimation-error variances due to smoothing for the two states at the initial time? Why is there so much more improvement for one state than for the other state?

  3. Simulate the system and smoother a hundred times or so, each simulation with a different noise history. On the basis of your simulations, derive a numerical estimate of the smoother estimation-error variances of the two states at the initial time. How do your numerical variances compare with the theoretical variances obtained in part (b)?

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.