State Estimation (GPS) Part 2

Estimator Design by Pole Placement

This section discusses estimator design by pole placement (i.e., Ackerman’s method) for two reasons.Second, understanding of the concepts of this section will aid the understand of the results derived in Section 3.6.3.

Section 3.1.3 derived the controllable canonical form state space representation for the strictly proper transfer function in eqn. (3.6). For observer design by pole placement, it is also useful to define the observable canonical form. In discrete-time, the observable canonical form, state space representation is

tmp3AC154_thumb[2]

where

tmp3AC155_thumb[2]

andtmp3AC156_thumb[2]The observable canonical form will be convenient for selection of the observer gain matrix, assuming that we have method for transforming an arbitrary state space representation to observable canonical form.


The observer gain matrix design process has three steps. First, we find a similarity transformtmp3AC157_thumb[2]from the original state space representation

tmp3AC160_thumb[2]

to observable canonical form wheretmp3AC161_thumb[2] andtmp3AC162_thumb[2]Second, we select an observer gain vectortmp3AC163_thumb[2]using pole placement for the system in observable canonical form. Third, we use the equationtmp3AC164_thumb[2]to transform the observer gaintmp3AC165_thumb[2]back to the original state space representation. Throughout this section, we assume that the system is observable. If the system was not observable, then there would in general be no solution to the full state estimation problem.

The first step is to derive the similarity transformtmp3AC166_thumb[2]to transform the original state space representation of eqn. (3.79) to observable canonical form of eqn. (3.78). To achieve this, we begin by computing the observability matrixtmp3AC167_thumb[2]as defined in eqn. (3.73). Because the system is observable,tmp3AC168_thumb[2]is nonsingular; therefore,tmp3AC169_thumb[2]exists. Let the last column of tmp3AC170_thumb[2]be denoted as p so thattmp3AC181_thumb[2]

Using the fact thattmp3AC182_thumb[2] and focusing only on the last column yields

tmp3AC183_thumb[2]

Definingtmp3AC184_thumb[2]we havetmp3AC185_thumb[2]which shows that the output matrix transforms correctly. Next, consider the rows of

tmp3AC188_thumb[2]

Using the fact that

tmp3AC189_thumb[2]

we have that fortmp3AC191_thumb[2]

Definingtmp3AC192_thumb[2]we are ready to consider the product

tmp3AC193_thumb[2]

which has the desired form for the observable canonical form. Therefore, the similarity transform defined by U transforms the original observable system to observable canonical form.

In the observable canonical form, it is straightforward to design the observer gaintmp3AC194_thumb[2]because

tmp3AC196_thumb[2]

which has the characteristic equation

tmp3AC197_thumb[2]

Therefore, if a set of desired discrete-time pole locationstmp3AC198_thumb[2]are selected, then the desired characteristic equation can be computed. The vectortmp3AC199_thumb[2]to cause eqn. (3.80) to match the desired characteristic equation can be computed andtmp3AC200_thumb[2]

Example 3.16 Consider the system described in eqn. (3.4) with m =1, b = 2 and k = 0. Defining the state vector astmp3AC201_thumb[2]continuous-time state space model for the system is

tmp3AC206_thumb[2]

For a sample period of T = 0.1 seconds, the equivalent discrete-time state space model given by eqn. (3.65) is

tmp3AC207_thumb[2]

The observability matrix is

tmp3AC208_thumb[2]

so that

tmp3AC209_thumb[2]

The resulting similarity transform uses

tmp3AC210_thumb[2]

With the closed loop estimator poles (i.e. the eigenvalues oftmp3AC211_thumb[2]and oftmp3AC212_thumb[2]specified to be attmp3AC213_thumb[2]The observer gain matrices aretmp3AC218_thumb[2]

 

 

State estimator simulation for Example 3.16. Top - State estimation error versus time. The solid line isand the dashed line is Bottom - State and state estimates versus time. The solid lines are the states and the dashed lines are the estimates.

Figure 3.6: State estimator simulation for Example 3.16. Top – State estimation error versus time. The solid line istmp3AC220_thumb[2]and the dashed line istmp3AC221_thumb[2] Bottom – State and state estimates versus time. The solid lines are the states and the dashed lines are the estimates.

The convergence of the resulting state estimator is shown in Figure 3.6 for the initial conditionstmp3AC224_thumb[2]The top plot shows the estimation error with xi as a solid line and X2 as a dashed line. The bottom plot shows state (solid) and the estimated state (dashed). Both plots clearly show the convergence of the estimation error to zero. The plots are shown as continuous lines, but the state estimates are actually only defined at the sampling instants. A

There are clearly a variety of tradeoffs in the selection of the state estimation gain vector. If L is small then the convergence will be slow, but measurement noise will have a small effect on the state estimates. Alternatively, when L is large convergence of the error will be rapid, but measurement noise will have a relatively large effect on the state estimates. When it is possible to propagate estimates of the accuracy of the state estimates and of the measurement accuracy, then it is interesting to consider the time-varying gain vector that provides an optimal tradeoff between the accuracy of the state and measurements.

Next post:

Previous post: