Feedback linearization - davidar/scholarpedia GitHub Wiki

Feedback linearization is the process of determining a feedback law and a change of coordinates that transform a nonlinear system into a linear and controllable one.

Table of Contents Problem formulation Notations Relative degree and normal forms Conditions for feedback linearization Example References See also

Consider a single-input nonlinear system (without output)

<math>\label{sys_wo_y}

\dot x = f(x) + g(x)u </math>

with internal state <math>x \in \mathbb{R}^n</math> and input <math>u \in \mathbb{R}\ .</math> Given a point <math>x^\circ \in \mathbb{R}^n\ ,</math> the problem of feedback linearization consists in finding a neighborhood <math>U</math> of <math>x^\circ\ ,</math> a coordinate transformation <math>z=\phi(x)</math> defined on <math>U\ ,</math> and a feedback law

<math>

u=\alpha(x)+\beta(x)v </math> where <math>x \in U</math> and <math>v \in \mathbb{R}</math> is an external input, such that the resulting closed loop system

<math>

\dot x = f(x) + g(x)\alpha(x)+g(x)\beta(x)v </math> expressed in the <math>z</math> coordinates is linear and controllable.

Let <math>\lambda</math> be a real-valued function defined on a subset <math>U</math> of <math>\mathbb{R}^n\ ;</math> the differential of <math>\lambda\ ,</math> denoted <math>d\lambda(x)\ ,</math> is the row vector

<math>

d\lambda(x) := \left( \frac{\partial \lambda}{\partial x_1}(x) \ \frac{\partial \lambda}{\partial x_2}(x) \ \cdots \ \frac{\partial \lambda}{\partial x_n}(x) \right):= {\partial \lambda \over \partial x}(x)\,.</math> Let <math>f</math> be an <math>n</math>-vector-valued function defined on <math>U\ .</math> The function <math>L_f\lambda</math> is the real-valued function defined as

<math>L_f\lambda(x):=\sum_{i=1}^n{\partial\lambda\over\partial

x_i}f_i(x) = {\partial\lambda\over\partial x}f(x)\,.</math> This function is sometimes called derivative of <math>\lambda</math> along <math>f\ .</math> If <math>\lambda</math> is being differentiated <math>k</math> times along <math>f\ ,</math> the notation <math>L_f^k \lambda</math> is used; in addition <math>L_f^0 \lambda=\lambda\ .</math>

Let <math>g</math> be another <math>n</math>-vector-valued function defined on <math>U\ .</math> The Lie product (or bracket) of <math>f</math> and <math>g</math> denoted by <math>[f,g]</math> is the <math>n</math>-vector-valued function defined by

<math>

[f,g](x):=\frac{\partial g}{\partial x}f(x)-\frac{\partial f}{\partial x}g(x)\;. </math> If <math>g</math> is bracketed with <math>f</math> several times the notation <math>ad^k_fg</math> is used, which is defined as <math>ad^k_fg=[f,ad^{k-1}_fg]\ ,</math> with <math>ad^0_f g=g\ .</math>

Let <math>f_i \ i=1,\ldots,d</math> be <math>n</math>-vector-valued functions defined on <math>U\ .</math> The assignment, with each <math>x \in U\ ,</math> of the subspace

<math>

\Delta(x)=\mathrm{span} \{f_1(x), \ldots, f_d(x) \} </math> defines a distribution <math>\Delta\ .</math> It is said that <math>f \in \Delta</math> if <math>f(x) \in \Delta(x)\ \forall x \in U\ .</math> The distribution <math>\Delta</math> is involutive if <math> \tau_1~\in~\Delta</math> and <math>\tau_2 \in~\Delta \Rightarrow [\tau_1,] \in \Delta\;. </math> It can be proved that <math>\Delta</math> is involutive if and only if <math> [f_i,] \in \Delta</math> for all <math>i,j=1,\ldots,d \ .</math>

The point of departure of the whole analysis is the notion of relative degree of a system, which is formally described in the following way. The single-input single-output nonlinear system

<math>\label{411}

\begin{array}{rcl} \dot x&=& f(x)+g(x)u\\ y&=& h(x) \end{array} </math>

with internal state <math>x \in \mathbb{R}^n\ ,</math> input <math>u \in \mathbb{R}\ ,</math> and output <math>y \in \mathbb{R}</math> is said to have relative degree  <math>r</math> at a point <math>x^\circ</math> if

(i) <math>L_gL_f^kh(x)=0</math> for all <math>x</math> in a neighborhood of <math>x^\circ</math> and all <math>k =0, \ldots, r-2</math>

(ii) <math>L_gL_f^{r-1}h(x^\circ )\ne 0\;.</math>

It is possible to show that the functions <math>h(x)\ ,</math> <math>L_fh(x),\ldots,L_f^{r- 1}h(x)</math> can be used in order to define, at least partially, a local coordinates transformation around <math>x^\circ\ .</math> This fact is based on the following property.

Lemma 1. If system \eqref{411} has relative degree <math>r</math> at <math>x^\circ\ ,</math> then the row vectors

<math>dh(x^\circ ),dL_fh(x^\circ

),\ldots,dL_f^{r-1}h(x^\circ )</math> are linearly independent.

Lemma 1 shows that necessarily <math>r \le n</math> and that the <math>r</math> functions <math>h(x), L_fh(x), \ldots, L_f^{r-1}h(x)</math> qualify as a partial set of new coordinate functions around the point <math>x^\circ\ .</math> If <math>r</math> is strictly less than <math>n\ ,</math> the set of new coordinates can be completed as described in the following.

Lemma 2. Suppose that system \eqref{411} has relative degree <math>r</math> at <math>x^\circ\ .</math> Then <math>r\le n\ .</math> If <math>r</math> is strictly less than <math>n\ ,</math> it is always possible to find <math>n-r</math> more functions <math>\psi_{1}(x),\ldots,\psi_{r-r}(x)</math> such that the mapping

<math>\Phi(x)=\left(\begin{array}{c} \psi_1(x)\\ \ldots\\ \psi_{n-r}(x)\\ h(x)\\ L_fh(x)\\ \ldots\\ L_f^{r-1}h(x) \end{array} \right)</math>

qualifies as a local coordinates transformation in a neighborhood of <math>x^\circ\ .</math> Moreover, it is always possible to choose <math>\psi_{1}(x), \ldots, \psi_{n-r}(x)</math> in such a way that <math>L_g \psi_i(x)=0</math> for all <math>i = 1, \ldots, n-r</math> and all <math>x</math> around <math>x^\circ.</math>

Setting

<math>

z = \left( \begin{array}{c} \psi_{1}(x)\\ \psi_{2}(x) \\ \cdots \\ \psi_{n-r}(x) \end{array} \right), \quad \xi = \left( \begin{array}{c} h(x)\\ L_fh(x) \\ \cdots \\ L_f^{r-1}h(x) \end{array} \right)\;, </math> the description of the system in the new coordinates is as follows

<math>\label{417a}

\begin{array}{rcl} \dot z &=& f_0(z,\xi) \\ \dot \xi_1&=& \xi_2\\ \dot \xi_2&=& \xi_3\\ &&\cdots\\ \dot \xi_{r-1}&=&\xi_r\\ \dot \xi_r&=& q(z,\xi)+b(z,\xi)u\\ y &=& \xi_1 \end{array} </math>

where

<math>

\begin{array}{rcl} f_0(z,\xi) &=& \left( \begin{array}{c} L_f\psi_1(\Phi^{-1}(z,\xi))\\ \cdots\\ L_f\psi_{n-r}(\Phi^{-1}(z,\xi))\end{array} \right) \\ q(z,\xi)&=&L_f^rh(\Phi^{-1}(z,\xi))\\ b(z,\xi)&=&L_gL_f^{r-1}h(\Phi^{-1}(z,\xi))\;. \end{array} </math> Note that <math>b(z^\circ,\xi^\circ) \ne 0</math> where <math>(z^\circ,\xi^\circ)=\Phi(x^\circ )\ .</math> Thus, the function <math>b(z,\xi)</math> is nonzero for all <math>(z,\xi)</math> in a neighborhood of <math>(z^\circ,\xi^\circ)\ .</math>

Equations \eqref{417a} are said to be in normal form.

In this section conditions and constructive procedures are given for finding a solution to the feedback linearization problem.

Consider the nonlinear system with output \eqref{411} and suppose that at some point <math>x^\circ</math> the system has relative degree equal to the dimension of the state, i.e. <math>r=n\ .</math> In this case, the change of coordinates that puts the system in normal form is given by

<math>\label{lin_coord}

\Phi(x)= \left( \begin{array}{c} h(x)\\ L_f h(x)\\ \cdots\\ L_f^{n-1}h(x) \end{array} \right) \;. </math>

Performing the change of coordinates <math>\xi=\Phi(x)</math> and neglecting the output, the system is described by

<math>

\begin{array}{rcl} \dot \xi_1&=& \xi_2\\ \dot \xi_2&=& \xi_3\\ &&\cdots\\ \dot \xi_{n-1}&=&\xi_n\\ \dot \xi_n&=& q(\xi)+b(\xi)u \end{array} </math> where

<math>

\begin{array}{rcl} q(\xi)&=&L_f^rh(\Phi^{-1}(\xi))\\ b(\xi)&=&L_gL_f^{r-1}h(\Phi^{-1}(\xi))\;. \end{array} </math> Recall that in a neighborhood of <math>\xi^{\circ}=\Phi(x^{\circ})</math> the function <math>b(\xi)</math> is nonzero.

Apply now the following feedback control law

<math>

u=\frac{1}{b(\xi)}(-q(\xi)+v) </math> which is well defined in a neighborhood of <math>\xi^{\circ}\ .</math> The resulting closed loop system is given by

<math>

\begin{array}{rcl} \dot \xi_1&=& \xi_2\\ \dot \xi_2&=& \xi_3\\ &&\cdots\\ \dot \xi_{n-1}&=&\xi_n\\ \dot \xi_n&=& v \end{array} </math> which is linear and controllable. Thus, we conclude that any nonlinear system with relative degree <math>n</math> at some point <math>x^{\circ}</math> can be transformed, using a coordinate transformation and a feedback, into a linear and controllable system.

In general, the two transformations used in order to obtain a linear and controllable system can be interchanged: one can first apply a feedback and then change the coordinates without modifying the result. Note that the feedback just used expressed in the <math>x</math> coordinates is given by

<math>

u=\frac{1}{b(\Phi(x))}(-q(\Phi(x))+v)=\frac{1}{L_gL_f^{n-1}h(x)}(-L_f^nh(x)+v)\;. </math> The above feedback is called linearizing feedback and the coordinates defined by \eqref{lin_coord} are called linearizing coordinates.

Of course, the basic feature that made feedback linearization possible was that system \eqref{411} had relative degree <math>n</math> (at <math>x^{\circ}</math>). Thus, the system without output \eqref{sys_wo_y} can be feedback linearized if it is possible to find an "output" function <math>h(x)</math> such that the corresponding system with output \eqref{411} has relative degree <math>n</math> (at <math>x^{\circ}</math>). The above condition turns out to be also necessary for solving the feedback linearization problem as stated in the following.

Lemma 3. The feedback linearization problem is solvable if and only of there exists a neighborhood <math>U</math> of <math>x^{\circ}</math> and a real-valued function <math>h(x)</math> defined on <math>U\ ,</math> such that system \eqref{411} has relative degree <math>n</math> at <math>x^{\circ}\ .</math>

The question now arises of when, given a system of the form \eqref{sys_wo_y} (namely a pair of <math>n</math>-vector-valued functions <math>f(x)</math> and <math>g(x)</math>) there exists an "output function" <math>h(x)</math> such that the resulting system \eqref{411} has relative degree <math>n</math> at <math>x^{\circ}\ .</math> This is answered by the following result.

Theorem 1. Let <math>f(x)\ ,</math> <math>g(x)\ ,</math> and <math>x^{\circ}</math> be given. The feedback linearization problem is solvable (i.e. there exists an "output" function <math>h(x)</math> for which system \eqref{411} has relative degree <math>n</math> at <math>x^{\circ}</math>) if and only if the following conditions are satisfied

(i) the matrix <math>\left( g(x^{\circ}) \ ad_f g(x^{\circ}) \cdots \ ad_f^{n-2} g(x^{\circ}) \ ad_f^{n-1} g(x^{\circ}) \right)</math> is nonsingular

(ii) the distribution <math>\mathrm{span} \{ g,\ ad_f g, \cdots, ad_f^{n-2}g \}</math> is involutive near <math>x^{\circ}\ .</math>

It can be shown that the "output" function <math>h(x)</math> for which the system has relative degree <math>n</math> at <math>x^{\circ}</math> is a solution of the following system of first order partial differential equations

<math>

L_gh(x)=L_{ad_{f}g}h(x)= \ldots = L_{ad^{n-2}_{f}g}h(x)=0\;. </math> If the assumptions of Theorem 1 are satisfied, then it is guaranteed that there exist solutions to the above system of equations. Once <math>h(x)</math> has been determined, the linearizing feedback and the linearizing coordinates can be obtained as discussed before.

Consider the system

<math>

\dot x = \left(\begin{array}{c} x_3(1+x_2) \\x_1 \\x_2(1+x_1) \end{array}\right) +\left(\begin{array}{c}0 \\1+x_2 \\-x_3 \end{array}\right) u\;. </math>

In order to check whether or not this system is feedback linearizable around <math>x=0\ ,</math> we compute the functions <math>ad_f g(x)</math> and <math>ad^2_f g(x)</math> and test the conditions of Theorem 1.

Appropriate calculations show that

<math>

ad_f g(x)=\left(\begin{array}{c}0 \\x_1 \\-(1+x_1)(1+2 x_2) \end{array}\right) </math> and that

<math>

ad^2_f g(x)=\left(\begin{array}{c}(1+x_2)(1+2 x_2)(1+x_1)-x_3 x_1 \\ x_3(1+x_2) \\-x_3(1+x_2)(1+2 x_2) \end{array}\right) </math> The matrix

<math>

\left( g(0) \ ad_f g(0) \ ad_f^{2} g(0) \right)= \left(\begin{array}{ccc}0 & 0 & 1 \\1 & 0 & 0 \\0 & -1 & 0\end{array}\right) </math> has rank 3 and therefore the condition (i) is satisfied. It is also easily checked that the product <math>[g,ad_f](x)</math> has the form

<math>

[g,ad_f](x)=\left(\begin{array}{c}0 \\ \star \\ \star \end{array}\right) </math> and therefore also the condition (ii) is satisfied, because the matrix

<math>

\left( g(x) \ ad_f g(x) \ [g,ad_f](x) \right) </math> has rank 2 for all <math>x</math> near <math>x=0\ .</math>

In the present case, it is easily seen that a function <math>h(x)</math> that solves the system of equations

<math>

L_g h(x)=L_{ad_{f}g}h(x)= 0 </math> is given by

<math>

h(x)=x_1\;. </math>

From our previous discussion, we know that considering this as "output" will yield a system having relative degree 3 at the point <math>x=0\ .</math> Then, locally around <math>x=0\ ,</math> the system will be transformed into a linear and controllable one by means of the state feedback

<math>

u=\frac{-L_f^3h(x)+v}{L_g L_f^2 h(x)}= \frac{-x_3^2(1+x_2)-x_2 x_3(1+x_2^2)-x_1(1+x_1)(1+2x_2)-x_1x_2(1+x_1)+v} {(1+x_1)(1+x_2)(1+2x_2)-x_1x_3} </math> and the coordinate tranformation

<math>

\begin{array}{rcl} \xi_1 &=& h(x) = x_1\\ \xi_2 &=& L_f h(x)=x_3(1+x_2)\\ \xi_3 &=& L_f^2 h(x)=x_3 x_1+(1+x_1)(1+x_2)x_2\;. \end{array} </math>

Control theory

Category:Control Theory Category:Multiple_Curators