Abstract
This paper presents a study of linear control systems based on exact feedback linearization and approximate feedback linearization. As exact feedback linearization is applied, a linear controller can perform the control objectives. The approximate feedback linearization is required when a nonlinear system presents a noninvolutive property. It uses a Taylor series expansion in order to compute a nonlinear transformation of coordinates to satisfy the involutivity conditions.
1. Introduction
The control theory of nonlinear system has been receiving increasing attention in recent years, both for its technical importance as well as for its impact in various fields of application. In strategic areas, such as aerospace, chemical and petrochemical industries, bioengineering, and robotics, a new practical application for this tool appears every day.
System nonlinearity is characterized when at least one component or subsystem is nonlinear [1]. The classic methods used in the study of nonlinear systems, particularly frequency analysis, are not usually applied to the nonlinear systems. It is necessary to use other methods to study the control of these systems.
When the nonlinear effects of systems become significant, nonlinear control techniques generally fail to produce the desired performance. In these cases, there are several nonlinear control techniques that have proven to supply better results.
Among the possible techniques used to deal with the nonlinear control problem, this paper studies the approximate feedback linearization technique for controlling nonlinear systems. The applications of these techniques are illustrated through some examples.
2. Analysis of Necessary and Sufficient Conditons for Exact Feedback Linearization
The objective of this section is to detail the necessary and sufficient conditions for input-state linearization. These conditions are limiting to the implementation of the exact feedback linearization technique. However, as a solution to this problem, the approximate linearization technique will be described.
Initially, we consider a single input for a nonlinear system, given by This system is known as linear in the control and linear in the input, among other denominations. The functions and are vectors with field in .
The system in the form (2.1) is input-state linearizable if there is a region in that meets the following conditions:(i) the vector field is linearly independents in ,(ii) the set of vectors for is involutive,
where indicates the Lie bracket and is defined as [2]: The first condition indicates the system controllability in which the vector fields are equivalent to the controllability matrix for linear systems [1].
The second involutivity condition enables to find out a new vector of linear state through the states feedback.
3. Exact Feedback Linearization
If the aforementioned conditions were satisfied, then it is possible to find a scalar function such that where defines a state transformation given by [1] A control input is then obtained through such that the closed-loop system has a new set of variables in linear state where and are in the canonical form of Brunovski however, there is no loss of generality; as every linear system may be written in the companion form by a state transformation [1].
3.1. Inverted Pendulum Linearization
Figure 1 shows an overview of a single pendulum. The pendulum may be described as a fixed end connected to a sphere through an inflexible rod and can be swung when stimulated by a force or has a nonzero initial condition.
The model of a pendulum is given by where , , , , , and are, respectively, gravitational acceleration, length of rod, mass of the sphere, friction coefficient, rotation angle, and angular speed.
3.1.1. Calculation of Controllability and Involutivity of the Simple Pendulum
As it is a second-order system, the vector field for verifying the controllability will be formed by, where, Hence, the vector field will be given by The vector fields are linearly independent; therefore, the simple pendulum is controllable.
In order to calculate the involutivity of the system, it is necessary to calculate the Lie bracket of and verify if it may be written as a linear combination of and The Lie bracket of is equal to zero, as the column vectors of the matrix are constant. Therefore, the system is involutive, which satisfy the conditions to apply input-state linearization.
3.1.2. Calculation of the New State Vector and the Input Variable Using Differential Geometry
The first component of the state vector should be obtained through the expression (3.1). As the order of the system is equal to 2, equation becomes resulting in Therefore, From (3.2), solving Therefore, In order to satisfy the last equation, has to be chosen.
We use (3.3) to compute the second component of the new state vector This gives the new state vector The new input variable is calculated from (3.4), where Therefore, This results in a linear system The system block diagram may be represented as shown in Figure 2.
The control structure, according to Figure 2, has two feedbacks. The first feedback is responsible for system linearization, eliminating the existing linearity. The second feedback is the controller project based on the state feedback.
3.2. The Need for Differential Geometry to Calculate the Exact Linearization
The linearization result shown previously might be obtained directly from When this procedure is applied, the linearization calculation would be significantly less. The problem occurs when the control law obtained by this means does not permit the exact linearization.
As an example of the aforementioned problem, consider a nonlinear system given by Applying (3.22), we will obtain: This result does not make equation (3.23) a linear system. It requires a variable transformation by inference or by differential geometry. An approach through differential geometry is more generic and can be implemented in software through simple programming.
3.2.1. The Use of Differential Geometry to Linearize the Previous Example
The first component of the state vector is obtained through As , we will have Therefore, If , obtained through the expression (3.2), then Substituting the result of (3.28) in the matrix , we obtain Consequently, This expression indicates that ; therefore, we should choose the new state in such way that it satisfies this condition.
Choosing , we obtain With the new state space given by ,
4. Approximate Feedback Linearization
There are situations in the design of controllers where the nonlinear terms are significant to the extent that failure to consider these terms in the system model could cause it to perform badly or to act restricted to a narrow operation band.
For systems where the nonlinear terms cannot be ignored, however, a physical model very close to reality can be made. The solution is to use exact feedback linearization. In this case, it is possible to linearize the system and apply linear control methods. Unfortunately exact linearization cannot be used for noninvolutive systems, as they do not satisfy the aforementioned conditions.
To solve this problem, the approximate feedback linearization, formulated by Krener [3], has been adopted using differential geometry techniques.
The approximate feedback linearization may be considered a near equivalent version of the exact feedback linearization. On the other hand, it offers the possibility to control noninvolutive systems, which was not possible using exact linearization. This new approach has recently been used in several applications.
As an example to validate this technique, the inverted pendulum model (Figure 3) with its characteristic nonlinear dynamic and an unstable open loop will be used. It is a classic problem in control theory and is widely used as a benchmark for testing of control algorithms (adaptive control, robust control, fuzzy, etc.).
4.1. Approximate Linearization Calculation
The approximate feedback linearization has been applied in cases where exact linearization is not feasible or has not yielded satisfactory results [4]. This process consists of finding approximate output functions that satisfy the involutivity condition up to determined system order.
The basic idea is [5], given where .
Considering a state transformation in the control variable, we obtain so that where , , , .
We apply a state transformation in the control variable such as where Differentiating (4.4), we obtain Substituting (4.1) in (4.7), Expanding in a Taylor series, Substituting (4.4) and (4.5) in (4.3), we obtain Equalizing the second order term of (4.9) and (4.10) so that the unknown constants cancel out the second-order terms , we obtain Thus, there is an approximate output , such as the set of vectors are involutive until certain order [6]. In other words, the gradient multiplied by the Lie bracket of and is equal to the superior terms . The following equation illustrates this: In order to find vector , it is necessary to expand the involutive distribution in a Taylor series on an equilibrium point . The expansion of the involutive distribution in can be represented through the following equation: where and
This expansion is carried out in order to obtain an involutivity until a certain order , where indicates the order in which the involutive distribution is expanded in a Taylor series to satisfy the involutivity condition. Therefore, the calculation of the approximate output function is done through the equation where, Once is known, the equation of the remaining state is given by And the input transformation The state equations are represented as It is important to emphasize that the choice for the representation of the linearized system was based on the Brunovski canonical form.
4.2. Inverted Pendulum Linearization
The proposed inverted pendulum consists of a rigid rod of negligible weight, attached to a mobile base called car and with a fixed mass at the other end. The inferior part of the pendulum rod is in direct contact with the car, which aims to keep the rod in the upright position according to the movement performance of the car control position. In the context of this paper, a pendulum that performs its movement in only one plane was conidered. In other words, it has only one degree of freedom.
We consider the following nonlinear function: where , , , and , and the movement equation may be described as [7] For the involutive order 1, there should be a , not trivial, that satisfies the equation below [8]: where resulting in integrating Once has been calculated, the remaining states and the control variable are obtained from [6]. The last column on the right side represents the liberalized system neglected terms.
The input variable is given by where resulting in system approximately linear given by where
5. Simulation Results
In the simulation, a comparison of a controller for an inverted pendulum using the Taylor linearization and the approximate feedback linearization techniques was carried out. The pole placement control technique was used, where the close loop poles in both linearizations were placed in , , , and . It can be seen from Figures 4 and 5 that the controller based on the approximate feedback linearization performed better for the initial positions , , and . From the simulation results, it can be seen that for the Taylor linearization, the controller managed to stabilize the inverted pendulum for , while with the approximate linearization, the controller managed the stabilization for , which represents an increase in the inverted pendulum operation region.
The simulation shows the effectiveness of the proposed controller. Another important factor was the response of the controller as regards its position. According to Figures 6 and 7, the displacement of the inverted pendulum was less for the controller based on the approximate feedback linearization.
Another simulation was performed increasing the weight of the fixed mass attached to the car which contains the weight five times.
6. Conclusions
This work describes feedback linearization techniques used as a tool for differential geometry calculus. The existing restrictions to the use of the exact input-state feedback linearization technique were analysed, and the impossibility of using exact linearization in involutive systems was observed. As a solution to this, we proposed approximate linearization.
To illustrate the problem of exact linearization, a simple pendulum was used, where the calculation of the input variable as well as the state transformations were performed using differential geometry as the main calculus tool.
In the approximate linearization technique, an inverted pendulum was used, as it is an example of great interest. The oscillating behaviour of the inverted pendulum rod reproduces the stabilization problem common in situations such as the trajectory of a projectile or the movement of a satellite. In addition, it is an unstable in open-loop, nonlinear, noninvolutive, and an underactuated system, which makes it attractive for application in advanced control techniques.
As a demonstration of the effectiveness of the method, a simulation was carried out to compare the Taylor linearization with the approximate feedback linearization, with the latter technique increasing the inverted pendulum operation region.
According to the results obtained with a five-kilogram weight, shown in Figure 8, on the top of the inverted pendulum flexible rod, it can be inferred that there is the possibility possess a good robustness.
Noninvolutivity, an important detail that prevents the application of the exact feedback linearization technique, was solved through coordinate transformation. This new state space, in turn, is involutive, enabling, therefore, exact linearization. In other words, as this coordinate transformation generates small errors, the linearization is no longer exact, making it an approximate linearization. However, it has been noted that despite these errors, the results indicate that the technique achieved a superior performance when compared to the Taylor linearization.