LMIs in Control/pages/stabilization of second order systems

From testwiki
Jump to navigation Jump to search

LMIs in Control/pages/stabilization of second order systems

Stabilization is a vastly important concept in controls, and is no less important for second order systems. A second order system can be conceptualized most simply by the model of a mass-spring-damper. Velocity and position are of course chosen as the states for this system, and the state space model can be written as it is below. The goal of stabilization in this context is to design a control law that is made up of two controller gain matrices KpRr*mp, KdRr*md. These allow the construction of a stabilized closed loop controller.

The System

Here, we want to stabilize a second order system of the following form:

Mx¨+Dx˙+Kx=Bu,

where xRn and uRr are the state vector and the control vector, respectively, and M (called the "mass matrix"), D (called the "structural damping matrix"), K (called the "stiffness matrix"), and B are the system coefficient matrices of appropriate dimensions.

To make the system follow standard convention, we reformulate the system as:

A2x¨+A1x˙+A0x=Buyd=Cdx˙yp=Cpx,

where: xRn and uRr are the state vector and the control vector, respectively; ydRmd and ypRmp are the derivative output vector and the proportional output vector, respectively; and A2,A1,A0,B,Cd, and Cp are the system coefficient matrices of appropriate dimension. Note that A2 must be >0, and A0,A2 must be Sn.

To further define: x isRn and is the state vector, A0 is Rn*n and is the state matrix on x , A1 is Rn*n and is the state matrix on x˙ , A2 is Rn*n and is the state matrix on x¨, B is Rn*r and is the input matrix, u is Rr and is the input, Cd and Cp are Rm*n and are the output matrices, yd is Rm and is the output from Cd, and yp is Rm and is the output from Cp.


Rn

The Data

The matrices A2,A1,A0,B,Cd,Cp.

The Optimization Problem

For the system described, we choose the following control law

u=Kpyp+Kdyd=KpCpx^+KdCdx,

with KpRr*mp, KdRr*md, we obtain the closed-loop system as follows:

A2x¨+(A1BKpCp)x˙+(A0BKdCd)x=0.

We are tasked to design a state feedback control law such that the above system is hurwitz stable.

First, in order to solve this problem, we need to introduce a Lemma. This Lemma comes from Appendix A.6 in "LMI's in Control systems" by Guang-Ren Duan and Hai-Hua Yu. This Lemma states the following:

A2>0,A1+A1T>0,A0>0

The LMI: Stabilization of Second Order Systems

There is a solution if there exists matrices KpRr*mp and KdRr*md that satisfy the following LMIs:

A0BKdCd>0,

and

(A1BKpCp)+(A1BKpCp)T>0.

Conclusion:

Finally, having solved the LMI the optimization will produce two matrices, Kp and Kd that can be substituted into the system as

u=KpCpx˙+KdCdx

to obtain a stabilized second order system.

Implementation

This implementation requires Yalmip and Sedumi. https://github.com/rezajamesahmed/LMImatlabcode/blob/master/stab2ndorder.m

Robust Stabilization of Second-Order Systems

This LMI comes from

  • [1] - "LMIs in Control Systems: Analysis, Design and Applications" by Guang-Ren Duan and Hai-Hua Yu

Other resources:

References

Duan, G. (2013). LMIs in control systems: analysis, design and applications. Boca Raton: CRC Press, Taylor & Francis Group.

Return to Main Page:

Template:BookCat