Mathematica 9 is now available

 Documentation /  Advanced Numerical Methods /  Function Index /  Pole Assignment and Stabilization by State Feedback /

NaturalFrequencyRegionQRDecomposition

5.4 Lyapunov Feedback Stabilization

It is possible to stabilize a system by state feedback using a certain Lyapunov equation solution. The following method, due to Armstrong (1975), stabilizes a continuous-time system with a damping factor .

1. Choose the damping factor Alpha such that the eigenvalues of have positive real parts.

2. Solve the Lyapunov equation for

.

3. Obtain the stabilizing feedback matrix .

For the discrete-time system, the method has the same format. The Lyapunov equation to be solved in this case is and the state feedback matrix is given by , where Beta is chosen such that . The feedback matrix is such that (see Datta (2003) for details). The method in the discrete-time case was developed by Armstrong and Rublein (1976).

The aforementioned algorithms are accessible via the function StateFeedbackGains using the method LyapunovShift. Additionally, Advanced Numerical Methods provides the method PartialLyapunovShift that reassigns only those eigenvalues that lie outside the specified region, which can be given only as DampingFactorRegion[] or SettlingTimeRegion[, ]. This typically results in feedback matrices with a smaller norm. However, the algorithm in this case requires computing the ordered Schur decomposition before performing the Lyapunov shift and is computationally more expensive.

Methods for Lyapunov feedback stabilization.

Make sure the application is loaded.

In[1]:=

Load the collection of test examples.

In[2]:=

This is the typical state-space system arising in the magnetic tape control problem.

In[3]:=

Out[3]=

In[4]:=

Out[4]=

The partial Lyapunov feedback stabilization method is applied to the same state-space system.

In[5]:=

Out[5]=

The system has been stabilized with the damping factor .

In[6]:=

Out[6]=

NaturalFrequencyRegionQRDecomposition



Any questions about topics on this page? Click here to get an individual response.Buy NowMore Information
THIS IS DOCUMENTATION FOR AN OBSOLETE PRODUCT.