Chaos control via Mathieu-Van der Pol system and Linear Optimal Control design with a non-ideal excitation and parametric uncertainties

<!-- metadata commented in wiki content
==Chaos Control via Mathieu-Van der Pol System and Linear Optimal Control Design with a Non-ideal Excitation and Parametric Uncertainties==

'''Joselma Pinheiro Gonçalves Vicente<sup>a</sup>, Fábio Roberto Chavarette<sup>b</sup>, Luiz Gustavo Pereira Roefero<sup>c</sup>'''
-->

==Abstract==

In this work, a non-ideal system based on the Mathieu-Van der Pol System is modeled considering its parameters as uncertain, which makes it closer to reality.

Numerical simulations are presented demonstrating that the system has a chaotic behavior with three positive Lyapunov exponents. Given such unstable and unpredictable behavior, the linear optimal control design is proposed to reduce the chaotic movement of this system to a fixed point. The simulation results show that the identification by Linear Optimal Control is very effective.

Keywords: Non-Ideal System, Uncertain Parameters, Optimal Linear Control

==1. Introduction==

Humanity is going through a very fast technological revolution, which makes well established companies and start-ups strive daily to invest in the search for introducing new, customizable, autonomous and more efficient products in the work market. For such purposes, mechanical engineering contributes with projects of advanced technological products, such as the manufacturing of electric and autonomous vehicles in the automotive sector. Therefore, with a simple software update, costs related to maintenance for changing oil, filters and other parts of cars powered by internal combustion engines, will be replaced.

The challenge of creating a new product requires developing projects that are straightforward concerning the objectives to be achieved and all stages to be completed, from their beginning, development, until their end. During the execution phase of a mechanical project, the probability of mistakes happening is high, leading to undesirable effects, that is, negative effects, or even chaotic ones. Considering this, managing the risks effectively may decrease the occurrence of a chaotic event that hinders the development of a project. Formally speaking, a chaotic event, or simply chaos, represents the irregularity of a dynamical system. Its essential characteristics are: sensitivity to initial conditions, which means that its evolution in time may be altered by small perturbations, and unpredictability, which implies that it is impossible to predict the posterior evolution of the system, even knowing its state during a long-time interval.

A priori, it is possible to state that a chaotic response can be understood based on a series of contraction-expansion-bending transformations, and among the tools applied to characterize chaos, the Lyapunov exponents stand out, which are one of the invariants of the system, and may be used even for evaluating other invariants [1].

Due to the possibility of a chaotic event happening during the development of projects, the alternative is to seek to mitigate this disorder state, that is, to manage the chaos that is inherent to the system, checking the weaknesses that unexpectedly arise, so as not to harm the production stage. Nowadays, Chaos Engineering contributes to carrying out experiments in the infrastructure and pointing out the weaknesses related to the system before they arise. This empirical verification process leads to invulnerable systems, and provides reliability for a complex system to operate under turbulent situations during the operational process of production. In 1990, Grebogi, Edward Ott and James Yorke demonstrated the possibility of controlling chaos [2]. Ever since then, controlling chaotic systems has been of great interest in mechanical engineering, and implies in designing state feedback control laws which stabilize the chaotic system around the instable equilibrium points. Many control techniques may be used to control a system, and each of them presents peculiarities and different results, especially in practice. To stabilize chaotic systems, in many works feedback linear control was used[10][11] [18]. In the present paper, we apply the Optimal Linear Control (OLC) methodology, proposed by Rafikov and Balthazar [18], which ensures the application of linear control to nonlinear systems. Among the main characteristics of Optimal Linear Control, we highlight the use of the state-space representation model, the obtaining of optimal control signals through the Algebraic Ricatti Equation (ARE), and the objective function, parameterized by the matrices Q and R to weigh the state and control vectors, respectively [18].

This work aims at proposing a new numerical model for the Mathieu-Van der Pol system coupled to a source with a rotative unbalance [4][5][6][7] and considering uncertain parameters; afterwards, a control technique is applied to stabilize the whole proposed system.

Nowadays, the applications in mechanical engineering require that the systems operate with optimal parameters under certain operation conditions, which demonstrates the need for developing realistic numerical models which consider uncertainties in the parameters and inputs of the system properly. In mechanical engineering, the uncertainties have been considered in several areas, such as structure dynamics [13], rotor dynamics [14], robot manipulators  [15], among others. In the process of modeling dynamical systems, uncertainties may be considered based on a probabilistic approach, through predefined intervals, or even the possibility theory [16].

==2. Mathematical Model==

The Mathieu-Van der Pol autonomous system is a model with four state variables and three positive Lyaponuv exponents, proposed by Shih-Yu et al.[3]. It is a hyper-chaotic system that was modeled via linear coupling of the non-linear and non-autonomous systems of Mathieu and Van der Pol [3], whose equations of motion are as follows::

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\dot{x}_1=x_2</math>
|-
| style="text-align: center;" | <math> \dot{x}_2=-(a+bx_3)x_1-(a+bx_3)x_1^3-cx_2+dx_3 </math>
|-
| style="text-align: center;" | <math> \dot{x}_3=x_4</math>
|-
| style="text-align: center;" | <math> \dot{x}_4=-ex_3+f(1-x_3^2)x_4+gx_1 </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (1)
|}

where <math display="inline">x_1</math>, <math display="inline">x_2</math>, <math display="inline">x_3</math> and <math display="inline">x_4</math> are the four states of the system and <math display="inline">a,b,c,d,e,f</math> and <math display="inline">g</math> are the parameters of the Mathieu-Van der Pol System [3][24].

The next section deals with non-ideal excitations, as well as presents a proposal for a new model called Non-Ideal Mathieu-Van der Pol System.

===2.1 Non-ideal System===

A system is said to be ideal when the excitation does not influence the response of the system, ie, the ideal excitation has an ideal power source [8][9][10][11][12]. On the other hand, when the excitation is influenced by the response of the system, the excitation is said to be non-ideal [19]. A non-ideal excitation can be modeled as shown in  Fig [[#img-1|1]] [4][11].

<div id='img-1'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Draft_Chavarette_796969493-motor_nonideal.png|441x441px|Non-ideal power source [11]   ]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figure 1:''' Non-ideal power source [11]   
|}

This representation on the based in Kononenko model and can be mathematically modeled according to equations [[#eq-2|2]] where <math display="inline">m_1</math> is a DC motor mass, <math display="inline">m_0</math> is unbalanced mass, <math display="inline">l</math> is the dumping, <math display="inline">k</math> is the rigidity, <math display="inline">r</math> is the eccentricity distance from the unbalanced mass to the torque source, <math display="inline">(I +m_0r^2)</math> is the moment of inertia of the unbalanced mass, and the state variables are <math display="inline">X</math> representation for the beam tip position and <math display="inline">\varphi </math> for the unbalanced mass angular position. The net torque is a function of the angular velocity <math display="inline">\dot{\varphi }</math> and described as for <math display="inline">S(\dot{\varphi }) = a-b \dot{\varphi }</math> [6][11][17], where a is the net torque applied by the DC motor and b is the resistive net torque constant.

<span id="eq-2"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>{} (m_1+m_0)\ddot{X}+l\dot{X}+kX=m_0r(\dot{\varphi }^2cos\varphi +\ddot{\varphi }sen \varphi ))</math>
|-
| style="text-align: center;" | <math> (I+m_0r^2)\ddot{\varphi }-m_0r\ddot{X}sen\varphi =s(\dot{\varphi }) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (2)
|}

The first characteristic of a non-ideal system is that it presents an equation that describes the interaction of the energy source with the ideal dynamic system, and the second characteristic is that the dynamics of a non-ideal system approaches the ideal case as the supplied power increases. For the present work the adimensional non-ideal mathematical model configuration [6][9] of the non-ideal power source is given by equations [[#eq-3|3]]

<span id="eq-3"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>{} \ddot{x}_1=-(a+bx_2)x_1-(a+bx_2)x_1^3-c\dot{x}_1+dx_2+\mu (\dot{z}^2 cosz+\ddot{z}senz)</math>
|-
| style="text-align: center;" | <math> \ddot{x}_2=ex_2+f(1-x_2^2)\dot{x}_2+gx_1</math>
|-
| style="text-align: center;" | <math> \ddot{z}=\varepsilon \ddot{x}_1senz+\alpha - \beta \dot{z} </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3)
|}

The angular position <math display="inline">z</math> is the non-ideal excitation response, <math display="inline">\alpha </math>, <math display="inline">\beta </math> are motor torque constants, <math display="inline">\epsilon </math> is the distance from the unbalanced mass to the rotation center of the DC motor and <math display="inline">\mu </math> is related to the moment of inertia of the system. All parameters are dimensionless positive constants. The term <math display="inline">\mu (\dot{z}^2 cosz+\ddot{z}senz)</math> is due to the interaction between the dynamical system and an energy source, for example, a DC motor with limited power supply. <math display="inline">\mu (\dot{z}^2 cosz+\ddot{z}senz)</math> is responsible for the Non-idealization of the system. The parameter <math display="inline">\alpha </math> is the applied torque constant and depends on the initial conditions, and <math display="inline">\beta </math> is the resistive net torque and has no influence of initial conditions and is considered as internal damping of the DC Motor.

Making <math display="inline">y_1=x_1</math>, <math display="inline">y_2=\dot{x_1}</math>,<math display="inline">y_3=z</math>,<math display="inline">y_4=\dot{z}</math>, <math display="inline">y_5=x_2</math> and <math display="inline">y_6=\dot{x_2}</math>, <math display="inline">y_5=z</math> and , has

<span id="eq-4"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>{} \dot{y}_1 = y_2,</math>
|-
| style="text-align: center;" | <math> \dot{y}_2 = - \frac{-(a+by_3)y_1-(a+by_3)y_1^3-cy_2+dy_3+\mu y_4^2cosy_3 +(\alpha{-\beta}y_4)\mu seny_3}{1-\mu \epsilon (seny_3)^2},</math>
|-
| style="text-align: center;" | <math> \dot{y}_3 = y_4,</math>
|-
| style="text-align: center;" | <math> \dot{y}_4 = -\frac{\epsilon (-(a+by_3)y_1-(a+by_3)y_1^3-cy_2+dy_3)seny_3+\mu \epsilon y_4^2 cosy_3 seny_3+\alpha -\beta y_4}{1-\mu \epsilon (sen y_3)^2},</math>
|-
| style="text-align: center;" | <math> \dot{y}_5 = y_6,</math>
|-
| style="text-align: center;" | <math> \dot{y}_6 =-ey_5+f(1-y_5^2)y_6+gy_1.</math>
|-
| style="text-align: center;" | 
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (4)
|}

Figure&nbsp;[[#img-2|2]] and Figure&nbsp;[[#img-3|3]] illustrate the behavior of the adopted dynamical model, using numerical values <math display="inline">a=91.7, b=5.023, c=0.01, d=91, e=87.001, f=18</math> and <math display="inline">g=9.5072</math> [24] for the Mathieu-Van der Pol System and <math display="inline">\alpha=0.8</math>, <math display="inline">\beta=1.5</math>, <math display="inline">\epsilon=0.3</math> and <math display="inline">\mu=0.2</math> for the non-ideal excitation [6].

<div id='img-2'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Draft_Chavarette_796969493-fhistorico.png|560x560px|Time History. (a)  displacement for the Mathieu system; (b)  velocity for the Mathieu system; (c)displacement for the non-ideal excitation; (d) velocity for the non-ideal excitation; (e) displacement for Van der Pol system, and (f) velocity for the Van der Pol system   ]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figure 2:''' Time History. (a)  displacement for the Mathieu system; (b)  velocity for the Mathieu system; (c)displacement for the non-ideal excitation; (d) velocity for the non-ideal excitation; (e) displacement for Van der Pol system, and (f) velocity for the Van der Pol system   
|}

<div id='img-3'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Draft_Chavarette_796969493-fdisplac.png|560x560px|Phase Portrait. (a) Mathieu system, (b) Non-ideal excitation, and (c) Van der Pol system   ]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figure 3:''' Phase Portrait. (a) Mathieu system, (b) Non-ideal excitation, and (c) Van der Pol system   
|}

Figure&nbsp;[[#img-4|4]] illustrates the six Lyapunov exponents, <math display="inline">\lambda _1=8.9985</math>; <math display="inline">\lambda _2=8.9539</math>; <math display="inline">\lambda _3=-6.4645</math>; <math display="inline">\lambda _4=-1.4999</math> and <math display="inline">\lambda _5=-3.203</math>, <math display="inline">\lambda _6=-3.2329</math>, demonstrating the presence of chaos with one positive Lyapunov exponent. The total time for computing the Lyapunov exponents was <math display="inline">\Delta \tau  =100.000</math> with a time-step of <math display="inline">\tau=1</math>.

<div id='img-4'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Draft_Chavarette_796969493-flyapunov.png|560x560px|Dynamics of Lyapunov expoents   ]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figure 4:''' Dynamics of Lyapunov expoents   
|}

Aiming to minimize vibrations and reduce the oscillatory motion caused in the system , in the following section we propose the application of Linear Optimal Control to reduce this chaotic behavior to a fixed point.

To consider the effect of the parameter uncertainties on the performance of the controller they are added to the state. The nominal values of the parameters are given below.

The real unknown parameters of system are supposed to be as follows:<math display="inline">\bar{a}=73.36+36.68r(t)</math>, <math display="inline">\bar{b}=4.024+2.012(t)</math>, <math display="inline">\bar{c}=0.008+0.002(t)</math>, <math display="inline">\bar{d}=72.8+36.4r(t)</math>, <math display="inline">\bar{e}=69.6008+34.8004r(t)</math>, <math display="inline">\bar{f}=14.4+7.2r(t)</math>, <math display="inline">\bar{g}=7.5928+3.8288r(t)</math>, <math display="inline">\bar{\alpha }=0.64+0.32r(t)</math>, <math display="inline">\bar{\mu }=0.16+0.08r(t)</math>, <math display="inline">\bar{\beta }=1.2+0.6r(t)</math> and <math display="inline">\bar{eta}=0.24+0.12r(t)</math>  where <math display="inline">r(t)</math> are normally distributed random functions [20]. Numerical simulation results are shown in the next section, where the robustness of the methods is assessed considering that the system parameters present uncertainties.

==3. Control Design==

The Optimal Linear Control design for the Non-ideal Mathieu-Van der Pol system reduces the oscillatory movement to a small stable orbit. Due to its simplicity in terms of configuration and implementation, the linear state feedback control is especially attractive [18].

It should be stressed that this approach is analytical and it can be used without dropping any non-linear term. Rewriting the equations of motion [[#eq-4|4]], in a state form:

<span id="eq-5"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>{} \dot{x} = A{x} + g({x}) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (5)
|}

Considering a vector function <math display="inline">\tilde{x}</math>, witch characterizes the desired trajectory, and the control vector <math display="inline">U</math>, consisting of two parts: <math display="inline"> \tilde{u} </math> for the feedforward and <math display="inline">u_f</math> for a linear feedback, in such way that,

<span id="eq-6"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>{} u_f = Bu </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (6)
|}

where <math display="inline">B</math> is a constant matrix, and driving the deviation of the trajectory of system  [[#eq-6|6]] to the desired one <math display="inline">y=x-\tilde{x}</math>.

The control is optimal, in order to transfer the non-linear system [[#eq-4|4]] from any initial to final state <math display="inline">y(t_f)=0</math>, minimizing the functional<math display="inline"> \tilde{J}=\frac{1}{2}\int ^\infty _{0} (y^T\tilde{Q}y+u^tRu)dt</math>, where the symmetric matrix <math display="inline">P(t)</math> is evaluated through the solution of the matrix Ricatti differential equation [19][20][21][22].

<span id="eq-7"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>{} PA+A^TP-PBR^{-1}B^TP+Q=0 </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (7)
|}

satisfying the final condition <math display="inline">P(t_f)=0</math>.

In addition, with the feedback control [[#eq-4|4]], there exists a neighborhood <math display="inline">\Gamma _0 \subset \Gamma </math>,<math display="inline">\Gamma \subset \Re ^n</math>, of the origin such that if <math display="inline">x_0 \in \Gamma _0</math>, the solution <math display="inline">x(t)=0</math>,<math display="inline">t\geq{0}</math>  of the controlled system [[#eq-7|7]] is locally asymptotically stable, and <math display="inline">J_{min}=x^T_0P(0)x_0</math>. Finally, if <math display="inline">\Gamma =\Re ^n</math> then the solution <math display="inline">y(t)=0</math>,<math display="inline">t>0</math> of the controlled system [[#eq-7|7]] is globally asymptotically stable.

Using the theorem [20] the dynamic error y can be minimized (<math display="inline">y \rightarrow 0</math>).

===3.1 Application of the Optimal Linear Control to the System===

From equation [[#eq-4|4]]:

<span id="eq-8"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>{} \dot{y}_1 = y_2+U,</math>
|-
| style="text-align: center;" | <math> \dot{y}_2 = - \frac{-(a+by_3)y_1-(a+by_3)y_1^3-cy_2+dy_3+\mu y_4^2cosy_3 +(\alpha{-\beta}y_4)\mu seny_3}{1-\mu \epsilon (seny_3)^2},</math>
|-
| style="text-align: center;" | <math> \dot{y}_3 = y_4,</math>
|-
| style="text-align: center;" | <math> \dot{y}_4 = -\frac{\epsilon (-(a+by_3)y_1-(a+by_3)y_1^3-cy_2+dy_3)seny_3+\mu \epsilon y_4^2 cosy_3 seny_3+\alpha -\beta y_4}{1-\mu \epsilon (sen y_3)^2},</math>
|-
| style="text-align: center;" | <math> \dot{y}_5 = y_6,</math>
|-
| style="text-align: center;" | <math> \dot{y}_6 =-ey_5+f(1-y_5^2)y_6+gy_1.</math>
|-
| style="text-align: center;" | 
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (8)
|}

where the function of control <math display="inline">U</math> is defined in equation [[#eq-5|5]].

Rewriting the dynamical system [[#eq-4|4]] in the space-state form:

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\dot{y}_1 = y_2 +U,</math>
|-
| style="text-align: center;" | <math> \dot{y}_2 = - \frac{-k_1y_1-k_2(y_1-y_3)+dy_6^2cos(y_5)+dsin(y_5)(a-by_6)}{m_1-dr(sin^2y_5)}</math>
|-
| style="text-align: center;" | <math> \dot{y}_3 = y_4,</math>
|-
| style="text-align: center;" | <math> \dot{y}_4 = -\frac{y_2}{m_2}(y_3-y_1),</math>
|-
| style="text-align: center;" | <math> \dot{y}_5 = y_6,</math>
|-
| style="text-align: center;" | <math> \dot{y}_6 = \frac{rsin(y_5)(-k_1y_1-k_2(y_1-y_3))+dy_6^2rcos(y_5)sin(y_5)+drsin(y_5)^2+(m_1-dr(sin^2(y_5)))(a-by_6)}{m_1-dr(sin^2(y_5))}.</math>
|-
| style="text-align: center;" | 
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (9)
|}

And then

<math display="inline"> B = \begin{pmatrix} 1 \\ 1 \\1 \\1 \\1 \\1 \end{pmatrix}</math>, <math display="inline">y = \begin{pmatrix} x_1-\tilde{x}_1 \\ x_2-\tilde{x}_2 \\ x_3-\tilde{x}_3 \\ x_4-\tilde{x}_4 \\ x_5-\tilde{x}_5 \\ x_6-\tilde{x}_6     \end{pmatrix}</math>, <math display="inline">\tilde{x}=\begin{pmatrix} 0 \\ 0 \\ 0 \\0 \\0\\0 \end{pmatrix}</math>,<math display="inline">Q=I_6</math>, <math display="inline">A=\begin{pmatrix} 0&1&0&0&0&0 \\ -91.77&-0.01&91.10&0&0&0\\9.50&0&0&1&0&0 \\ -0.09&0&0.09&-1.49&0&0 \\0&0&0&1&0&1 \\9.50&0&0&0&-87.00&17.99 \end{pmatrix}</math>

where the controllability matrix <math display="inline">R</math> of the system for the pair [<math display="inline">A,B</math>] is obtained by <math display="inline">R=[B|AB|A^2B|A^3B|A^4B|A^5B]</math>.

Thus, <math display="inline">R= \begin{pmatrix} 1\end{pmatrix}</math>. Then the Matrix

<math display="inline">P(t)</math> is given by <math display="inline">P= 10^4\begin{pmatrix}  0.216 & 0.031 &  0.406 & -0.030 & -0.647 &  0.040 \\ 0.031 & 0.004 &  0.058 & -0.004 & -0.093 &  0.005 \\ 0.406 & 0.058 &  0.767 & -0.058 & -1.218 &  0.071 \\ -0.030& -0.004& -0.058 &  0.004 &  0.092 & -0.005 \\ -0.647& -0.093& -1.218 &  0.092 &  1.944 & -0.123 \\ 0.040 & 0.005 &  0.076 & -0.005 & -0.123 &  0.008  \end{pmatrix}</math>

and (an optimal control)

<math display="inline">u= 10^2(1.6765 x_1 + 0.2493 x_2 + 3.2667 x_3 - 0.2353 x_4 - 4.6331 x_5 + 0.2116 x_6</math>).

The trajectories of the system with control can be seen in Figures [[#img-5|5]] and [[#img-6|6]]. According to the optimal control verification [23], the function   is numerically calculated by <math display="inline">L(t)=y^T\tilde{Q}y</math>,where <math display="inline">L(t)</math> is positive definite.

<div id='img-5'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Draft_Chavarette_796969493-fhistcontrol.png|560x560px|Time History - the blue line is related to the system with uncertain parameters; the red line, with fixed parameters; the black line, in turn, is related to the controlled system. (a)  displacement for the Mathieu system; (b)  velocity for the Mathieu system; (c)displacement for the non-ideal excitation; (d) velocity for the non-ideal excitation; (e) displacement for Van der Pol system, and (f) velocity for the  Van der Pol system   ]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figure 5:''' Time History - the blue line is related to the system with uncertain parameters; the red line, with fixed parameters; the black line, in turn, is related to the controlled system. (a)  displacement for the Mathieu system; (b)  velocity for the Mathieu system; (c)displacement for the non-ideal excitation; (d) velocity for the non-ideal excitation; (e) displacement for Van der Pol system, and (f) velocity for the  Van der Pol system   
|}

<div id='img-6'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:Draft_Chavarette_796969493-fdisplacontrol.png|560x560px|Phase Portrait - the blue line is related to the system with uncertain parameters; the red line, with fixed parameters; the black line, in turn, is related to the controlled system. (a) Mathieu system, (b) Non-ideal excitation, and (c) Van der Pol system   ]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" | '''Figure 6:''' Phase Portrait - the blue line is related to the system with uncertain parameters; the red line, with fixed parameters; the black line, in turn, is related to the controlled system. (a) Mathieu system, (b) Non-ideal excitation, and (c) Van der Pol system   
|}

==4. Conclusions==

In this work, the dynamics of a new model called the Non-Ideal Mathieu-Van Der Pol System with chaotic behavior is proposed and analyzed.

We proposed the use of a Linear Optimal Control strategy to reduce the chaotic motion of this system to a fixed point.

Figures [[#img-5|5]] and [[#img-6|6]]  illustrated the effectiveness of the control strategy applied to solve these problems.

==Acknowledgements==

The authors would like to thank Universidade Estadual Paulista ”Júlio de Mesquita Filho”, especially the Faculty of Engineering of Ilha Solteira, for the space provided to carry out the work, as well as the Laboratory of Complex Systems. We would also like to thank the Universidade Federal do Estado de Mato Grosso for the support given to the work.

==References==
[1] Nayfeh, A. H., Mook, D.T. Non-linear Oscillations, John Wiley e Sons, New York, 1979.

[2] Ott., E., Grebogi, C., Yorque, J. A. Controlling Chaos, Phys. Rev. Lett. 66:1196-1199, 1990.

[3] Shih-Yu, L., Sheng-Chieh, H., Cheng-Hsiung, Y., Zheng-Ming, G. Generating tri-chaos attractors with three positive Lyapunov exponents in new four order system via linear coupling, Nonlinear Dyn, 69:805-816, 2012.

[4] Balthazar, J. M., Mook, D. T., Weber, H.I., Brasil, R. M. L. R. F., Fenili, A., Belato, D., Felix, J.L.P. An Overview on Non-Ideal Vibrations. Meccanica. 38:613-621, 2003.

[5] Dantas, M. J. H., Balthazar, J. M. On local analysis of oscillations of a non-ideal and non-linear mechanical model. Meccanica. 39:313-330, 2004.

[6] Chavarette, F. R. On an optimal linear control of a chaotic non-ideal duffing system. Applied Mechanics and Materials. 138-139:50-55, 2012.

[7] Souza, S. L. T., Caldas, I. L., Viana, R. L., Balthazar, J. M. Control and chaos for vibro-impact and non-ideal oscillators. Journal of Theoretical and Applied Mechanics. 46-3:641-664, 2008.

[8] Belato, D., Weber, H. I., Balthazar, J. M., Mook, D.T. Chaotic vibrations of a nonideal electro-mechanical system. International Journal of Solids and Structures. 38:1699-1706, 2001.

[9] Piccirillo, V., Balthazar, J. M., Pontes JR, B. R., Felix, J. L. P. On a nonlinear and chaotic non-ideal vibrating system with shape memory alloy (SMA). Journal of Theoretical and Applied Mechanics. 46-3:597-620, 2008.

[10] Chavarette, F. R., Balthazar, J. M., Felix, J. L. P. Remarks on an optimal linear control design applied to a nonideal and an ideal structure coupled to an essentially nonlinear oscillator. Journal of Computational and Nonlinear Dynamics. 5:1-9, 2010.

[11] Chavarette, F. R. Control design applied to a non-ideal structural system with behavior chaotic. International Journal of Pure and Applied Mathematics. 86-3: 487-500, 2013.

[12] Tsuchida, M., Guilherme, K. L., Balthazar, J. M. On chaotic vibrations of a non-ideal system with two degrees of freedom: 1 : 2 resonance and sommerfeld effect. Journal of Sound and Vibration. 282: 1201-1207, 2005.

[13] Ghanem, R. G., Spanos, P. D. Stochastic Finite Elements: A Spectral Approach. Courier Dover Publications, 1991.

[14] Koroishi, E. H., Cavalini Jr, A. A., de Lima, A. M. G. , Stefen Jr, V. Stochastics modeling of exible rotors. Journal of the Brazilian Society of Mechanical Sciences and Engineering , 34:597-603, 2012.

[15] Lara-Molina, F. A., Rosario, J. M., Dumur, D., Wenger, P. Generalized predictive control of parallel robots. Proceeding of Robot Motion and Control, 422, Lecture Notes in Control and Information Sciences, 159-169, 2012.

[16] Moller, B., Beer, M. Fuzzy Randomness, Uncertainty in Civil Engineering and Computational Mechanics. Springer-Verlag, 2004.

[17] Tusset, A. M., Balthazar, J. M., Chavarette, F. R., Felix, J. L. P. On energy transfer phenomena, in a nonlinear ideal and nonideal essential vibrating systems, coupled to a (MR) magneto-rheological damper. Nonlinear Dyn. 69-4:1-22, 2012.

[18] Rafikov, M., Balthazar, J. M. On control and synchronization in chaotic and hyperchaotic system via linear control feedback. Communications in Nonlinear Science and Numerical Simulation. 13:1246-1255, 2008.

[19] Chavarette, F. R., Peruzzi, N. J., Balthazar, J. M., Barbanti, L., Damasceno, B. C. On an optimal linear control applied to a non-ideal load transportation system, modeled whit periodic coefficients. Applied Mechanics and Materials.  52-54:13-18, 2011.

[20] Chavarette, F. R. Optimal linear control to parametric uncertainties in a micro electro mechanical system. International Journal of Pure and Applied Mathematics. 83-4:539-548, 2013.

[21] Frahm, H. Adevice for damping vibrations of bodies, US Patent 989958, 1911.

[22] Viquié, R., Kerschen, G. Nonlinear vibration absorber coupled to a nonlinear primary system: A tuning methodology. Journal of Sound and Vibration, 326:780-793, 2009.

[23] Cveti-anin, L. Dynamics of the non-ideal mechanical systems: a review. Journal of the Serbian Society for Computational Mechanics, 4-2:75-86, 2010.

[24] Chavarette, F. R., Peruzzi, N. J., Lopes, M. L. M., Cossi, A. M. A Note of Optimal Linear Control Technical Applied in the Chaos Stabilization in the Mathieu-van der Pol Autonomous Oscillator. International Journal of Applied Mathematics, 25:861-869, 2012.