<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.1 20151215//EN" "http://jats.nlm.nih.gov/publishing/1.1/JATS-journalpublishing1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="1.1">
<front>
<journal-meta>
<journal-id journal-id-type="pmc">CMC</journal-id>
<journal-id journal-id-type="nlm-ta">CMC</journal-id>
<journal-id journal-id-type="publisher-id">CMC</journal-id>
<journal-title-group>
<journal-title>Computers, Materials &#x0026; Continua</journal-title>
</journal-title-group>
<issn pub-type="epub">1546-2226</issn>
<issn pub-type="ppub">1546-2218</issn>
<publisher>
<publisher-name>Tech Science Press</publisher-name>
<publisher-loc>USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">22993</article-id>
<article-id pub-id-type="doi">10.32604/cmc.2022.022993</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Article</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Intelligent Dynamic Inversion Controller Design for Ball and Beam System</article-title>
<alt-title alt-title-type="left-running-head">Intelligent Dynamic Inversion Controller Design for Ball and Beam System</alt-title>
<alt-title alt-title-type="right-running-head">Intelligent Dynamic Inversion Controller Design for Ball and Beam System</alt-title>
</title-group>
<contrib-group content-type="authors">
<contrib id="author-1" contrib-type="author" corresp="yes">
<name name-style="western"><surname>Mehedi</surname><given-names>Ibrahim M.</given-names></name><xref ref-type="aff" rid="aff-1">1</xref>
<xref ref-type="aff" rid="aff-2">2</xref><email>imehedi@kau.edu.sa</email>
</contrib>
<contrib id="author-2" contrib-type="author">
<name name-style="western"><surname>Aljohani</surname><given-names>Abdulah Jeza</given-names></name><xref ref-type="aff" rid="aff-1">1</xref>
<xref ref-type="aff" rid="aff-2">2</xref></contrib>
<contrib id="author-3" contrib-type="author">
<name name-style="western"><surname>Alam</surname><given-names>Md Mottahir</given-names></name><xref ref-type="aff" rid="aff-1">1</xref></contrib>
<contrib id="author-4" contrib-type="author">
<name name-style="western"><surname>Mahmoud</surname><given-names>Mohamed</given-names></name><xref ref-type="aff" rid="aff-3">3</xref></contrib>
<contrib id="author-5" contrib-type="author">
<name name-style="western"><surname>Abdulaal</surname><given-names>Mohammed J.</given-names></name><xref ref-type="aff" rid="aff-1">1</xref></contrib>
<contrib id="author-6" contrib-type="author">
<name name-style="western"><surname>Bilal</surname><given-names>Muhammad</given-names></name><xref ref-type="aff" rid="aff-1">1</xref>
<xref ref-type="aff" rid="aff-2">2</xref></contrib>
<contrib id="author-7" contrib-type="author">
<name name-style="western"><surname>Alasmary</surname><given-names>Waleed</given-names></name><xref ref-type="aff" rid="aff-4">4</xref></contrib>
<aff id="aff-1"><label>1</label><institution>Department of Electrical and Computer Engineering (ECE), King Abdulaziz University</institution>, <addr-line>Jeddah, 21589</addr-line>, <country>Saudi Arabia</country></aff>
<aff id="aff-2"><label>2</label><institution>Center of Excellence in Intelligent Engineering Systems (CEIES), King Abdulaziz University</institution>, <addr-line>Jeddah, 21589</addr-line>, <country>Saudi Arabia</country></aff>
<aff id="aff-3"><label>3</label><institution>Electrical and Computer Engineering, Electrical and Computer Engineering</institution>, <addr-line>TN, 38505</addr-line>, <country>United States</country></aff>
<aff id="aff-4"><label>4</label><institution>Computer and Information Systems, Umm Al-Qura University</institution>, <addr-line>Makkah</addr-line>, <country>Saudi Arabia</country></aff>
</contrib-group>
<author-notes>
<corresp id="cor1"><label>&#x002A;</label>Corresponding Author: Ibrahim M. Mehedi. Email: <email>imehedi@kau.edu.sa</email></corresp>
</author-notes>
<pub-date pub-type="epub" date-type="pub" iso-8601-date="2022-03-26"><day>26</day>
<month>03</month>
<year>2022</year></pub-date>
<volume>72</volume>
<issue>2</issue>
<fpage>2341</fpage>
<lpage>2355</lpage>
<history>
<date date-type="received"><day>25</day><month>8</month><year>2021</year></date>
<date date-type="accepted"><day>20</day><month>12</month><year>2021</year></date>
</history>
<permissions>
<copyright-statement>&#x00A9; 2022 Mehedi et al.</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Mehedi et al.</copyright-holder>
<license xlink:href="https://creativecommons.org/licenses/by/4.0/">
<license-p>This work is licensed under a <ext-link ext-link-type="uri" xlink:type="simple" xlink:href="https://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution 4.0 International License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p>
</license>
</permissions>
<self-uri content-type="pdf" xlink:href="TSP_CMC_22993.pdf"></self-uri>
<abstract>
<p>The Ball and beam system (BBS) is an attractive laboratory experimental tool because of its inherent nonlinear and open-loop unstable properties. Designing an effective ball and beam system controller is a real challenge for researchers and engineers. In this paper, the control design technique is investigated by using Intelligent Dynamic Inversion (IDI) method for this nonlinear and unstable system. The proposed control law is an enhanced version of conventional Dynamic Inversion control incorporating an intelligent control element in it. The Moore-Penrose Generalized Inverse (MPGI) is used to invert the prescribed constraint dynamics to realize the baseline control law. A sliding mode-based intelligent control element is further augmented with the baseline control to enhance the robustness against uncertainties, nonlinearities, and external disturbances. The semi-global asymptotic stability of IDI control is guaranteed in the sense of Lyapunov. Numerical simulations and laboratory experiments are carried out on this ball and beam physical system to analyze the effectiveness of the controller. In addition to that, comparative analysis of RGDI control with classical Linear Quadratic Regulator and Fractional Order Controller are also presented on the experimental test bench.</p>
</abstract>
<kwd-group kwd-group-type="author">
<kwd>Ball and beam system (BBS)</kwd>
<kwd>unstable system</kwd>
<kwd>internal model</kwd>
<kwd>PD controllers</kwd>
<kwd>IDI control</kwd>
<kwd>robust control</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="s1"><label>1</label><title>Introduction</title>
<p>The ball and beam system (BBS) is a standard laboratory experimental setup and serves as an essential benchmark to investigate and validate the performance of various control strategies [<xref ref-type="bibr" rid="ref-1">1</xref>]. The tilt angle of the beam is the most significant action through which the ball position is controlled systematically by varying the tilt angle. It is an under-actuated open-loop unstable system having two degrees of freedom. The fundamental idea of BBS can be applied to a wide range of applications such as stabilizing the aircraft in its horizontal position while landing, tanker on-road carrying liquid, aircraft yaw roll control, goods handling problem by industrial robots, etc.</p>
<p>The control of an under-actuated system is an active research field having fewer control inputs than the state variables that need to be controlled. The experimental setup of BBS deals with the two loops control strategy. The outer loop is responsible for ball position control, whereas the inner loop is engaged for servo angle control. Due to its double loop control strategies, it is not easy to control this system using linear controllers [<xref ref-type="bibr" rid="ref-2">2</xref>]. Therefore, several control techniques are proposed in the literature to investigate the tracking and stabilization problem of BBS [<xref ref-type="bibr" rid="ref-3">3</xref>&#x2013;<xref ref-type="bibr" rid="ref-5">5</xref>].</p>
<p>A feedback linearization method is applied for BBS, which produces singularity and does not produce good results [<xref ref-type="bibr" rid="ref-3">3</xref>]. A Fuzzy-based controller is investigated in [<xref ref-type="bibr" rid="ref-6">6</xref>], showing the tracking performance numerically without implementing it on a practical testbed. The Fuzzy and Neural Networks based algorithms are presented in [<xref ref-type="bibr" rid="ref-7">7</xref>,<xref ref-type="bibr" rid="ref-8">8</xref>]. The robust control technique based on the sliding mode approach is also investigated for the ball balancing problem, see [<xref ref-type="bibr" rid="ref-9">9</xref>,<xref ref-type="bibr" rid="ref-10">10</xref>]. A Backstepping controller is designed in [<xref ref-type="bibr" rid="ref-11">11</xref>] to simulate the tracking problem of BBS. An anti-windup compensator (AWC) control scheme is also demonstrated in [<xref ref-type="bibr" rid="ref-12">12</xref>]. Similarly, the optimal performance of Model Predictive Control (MPC) is investigated to evaluate the tracking characteristics for ball and beam equipment [<xref ref-type="bibr" rid="ref-5">5</xref>].</p>
<p>Recent work was proposed designing a controller comprising two embedded loops based on linear fractional-order calculus [<xref ref-type="bibr" rid="ref-2">2</xref>]. The outer-loop and inner-loop control the ball&#x0027;s location on the beam and angle of servo gear, respectively. The prime focus of this method of control is to demonstrate a fractional-order control scheme to enhance efficiency, which is reported in several journals [<xref ref-type="bibr" rid="ref-13">13</xref>&#x2013;<xref ref-type="bibr" rid="ref-17">17</xref>]. However, the performance is unsatisfactory because an integer order traditional controller is related to the non-integer order integrator, which is used in the outer loop.</p>
<p>In addition to that, Nonlinear Dynamic Inversion (NDI) was also implemented for this ball and beam underactuated system, see [<xref ref-type="bibr" rid="ref-18">18</xref>]. However, this control method requires the exact model of the plant dynamics.</p>
<p>Therefore, NDI control may possess certain shortcomings and limitations while developing the algorithms. This includes the simplified approximations to introduce the inverse model of plant dynamics, square dimensionality condition, useful nonlinearity cancellations, etc.</p>
<p>Recently a control law focused on the inversion principle is proposed known as Generalized Dynamic Inversion (GDI), which addresses the limitations of the classical NDI approach [<xref ref-type="bibr" rid="ref-19">19</xref>]. The control objectives of GDI are defined in the form of differential state constraints which are calculated along the solution paths of the system dynamics. The control law is obtained by inverting the suggested constraints using Moore-Penrose Generalized Inverse (MPGI) [<xref ref-type="bibr" rid="ref-20">20</xref>]. Further, the GDI algorithm is made intelligent by augmenting the switching term based on the concept of classical SMC principle to ensure robustness against parametric uncertainties and un-modeled dynamics while guaranteeing semi-global asymptotic stability conditions [<xref ref-type="bibr" rid="ref-21">21</xref>&#x2013;<xref ref-type="bibr" rid="ref-25">25</xref>].</p>
<p>A two-loop control architecture is employed in this paper to solve the tracking and stabilization problem for two degrees of freedom under actuated BBS. A proportional derivative (PD) control is implemented in the outer loop to generate rotary servo angle commands based on the ball position on the beam. An intelligent version of GDI, Intelligent Dynamic Inversion (IDI), is used to process the inner loop&#x0027;s angular command. Identical (baseline) control and switching (continuous) control elements make up this IDI control system. The baseline control enforces the constraint dynamics based on the deviation function of the rotary servo angle for stable attitude tracking, while the switching control ensures robustness against nonlinearities and external perturbations. As a result, the proposed control guaranteed a semi-global attitude tracking system and bounded tracking errors. Quanser&#x0027;s BBS is used to simulate and perform laboratory experiments on the controller performance and compare the controller performance to existing control strategies.</p>
<p>The remaining sections of the paper are described as follows. The modeling of BBS is presented in Section 2. The two-loop control structure is discussed in Section 3. The classical PD control design for the outer position loop is presented in Section 4. The design implementation of IDI control for tracking the servo angle command is given in Section 5. The detailed stability analysis that will guarantee semi-global asymptotic stability was illustrated in Section 6. Controller validation through computer simulations and experimental findings are reported in Sections 7 and 8. Finally, the concluding remarks are given in Section 9.</p>
</sec>
<sec id="s2"><label>2</label><title>Description of BBS</title>
<p>Quanser developed BBS, a benchmark laboratory test bench [<xref ref-type="bibr" rid="ref-1">1</xref>], to validate control algorithms that are defined relative degrees with nonlinear attributes. The two significant ball and beam setup components are the ball and beam unit and rotary servo unit. One end of the finite length beam is fixed by a shaft, whereas the other end is connected to the rotary servo unit, as shown in the schematic diagram in <xref ref-type="fig" rid="fig-1">Fig. 1</xref>. Therefore, this unstable architecture is well known as an open-loop system in control technology. The core objective is to orient and stabilize the ball placed on the beam according to the commanded reference input profile. In <xref ref-type="fig" rid="fig-1">Fig. 1</xref>, <italic>x</italic> is the linear displacement ball, <italic>m<sub>b</sub> </italic>is the mass, <italic>r<sub>b</sub> </italic>is the beam&#x0027;s radius, and <italic>J<sub>b</sub> </italic>is the ball&#x0027;s momentum, respectively. The relative angle of the beam concerning the horizontal plane is denoted by <italic>&#x03B1;,</italic> whereas the gravitational acceleration is mentioned by <italic>g</italic>. The expression of the linear acceleration of the ball is given as <xref ref-type="disp-formula" rid="eqn-1">Eq. (1)</xref>.
<disp-formula id="eqn-1"><label>(1)</label><mml:math id="mml-eqn-1" display="block"><mml:mrow><mml:mover><mml:mi>x</mml:mi><mml:mo>&#x00A8;</mml:mo></mml:mover></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mtext>t</mml:mtext></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mrow><mml:mi>b</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mspace width="thickmathspace" /><mml:mi>g</mml:mi><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mi>&#x03B1;</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:msubsup><mml:mi>r</mml:mi><mml:mrow><mml:mi>b</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup></mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mrow><mml:mi>b</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:msubsup><mml:mi>r</mml:mi><mml:mrow><mml:mi>b</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup><mml:mo>+</mml:mo><mml:mrow><mml:msub><mml:mi>J</mml:mi><mml:mrow><mml:mi>b</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:mfrac></mml:math></disp-formula></p>
<fig id="fig-1"><label>Figure 1</label><caption><title>Ball and beam system</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_22993-fig-1.png"/></fig>
<p>Similarly, the relation between the beam angle <italic>&#x03B1;</italic> and servo angle <italic>&#x03B8;<sub>l</sub></italic>(<italic>t</italic>) is expressed as
<disp-formula id="eqn-2"><label>(2)</label><mml:math id="mml-eqn-2" display="block"><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mi>&#x03B1;</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mrow><mml:msub><mml:mi>r</mml:mi><mml:mrow><mml:mi>a</mml:mi><mml:mi>r</mml:mi><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mrow><mml:mi>b</mml:mi><mml:mi>e</mml:mi><mml:mi>a</mml:mi><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:mfrac></mml:math></disp-formula>where <italic>L<sub>beam</sub> </italic>is the length of the beam and <italic>r<sub>arm</sub> </italic>is the gap between the rotary servo output gear and the coupled joint. To relate the linear displacement <italic>x</italic> with the rotary servo angle <italic>&#x03B8;</italic>, place the expression of <italic>sin&#x03B1;</italic>(<italic>t</italic>) in the linear acceleration equation given by <xref ref-type="disp-formula" rid="eqn-2">Eq. (2)</xref>, it implies
<disp-formula id="eqn-3"><label>(3)</label><mml:math id="mml-eqn-3" display="block"><mml:mrow><mml:mover><mml:mi>x</mml:mi><mml:mo>&#x00A8;</mml:mo></mml:mover></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mtext>t</mml:mtext></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mrow><mml:mi>b</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mi>g</mml:mi><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mrow><mml:mi>l</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mspace width="thickmathspace" /><mml:mrow><mml:msub><mml:mi>r</mml:mi><mml:mrow><mml:mi>a</mml:mi><mml:mi>r</mml:mi><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mspace width="thickmathspace" /><mml:msubsup><mml:mi>r</mml:mi><mml:mrow><mml:mi>b</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup></mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mrow><mml:mi>b</mml:mi><mml:mi>e</mml:mi><mml:mi>a</mml:mi><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mspace width="thickmathspace" /><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mi>b</mml:mi></mml:msub></mml:mrow><mml:msubsup><mml:mi>r</mml:mi><mml:mrow><mml:mi>b</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup><mml:mo>+</mml:mo><mml:mrow><mml:msub><mml:mi>J</mml:mi><mml:mrow><mml:mi>b</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mfrac></mml:math></disp-formula></p>
<p>The second major component of the ball and beam setup is the rotary servo base unit which comprises a DC motor with a gearbox, as shown in <xref ref-type="fig" rid="fig-2">Fig. 2</xref>. The subscript <italic>m</italic> represents the motor attributes in the given schematic, whereas <italic>l</italic> characterizes the load shaft variables. In the rotary servo unit, the resistance is denoted by <italic>R</italic>, <italic>L</italic> is the inductance, <italic>&#x03C9;</italic> is the angular speed, <italic>&#x03C4;</italic> is the applied torque, <italic>J</italic> is the moment of inertia, and <italic>B</italic> stands for the viscous friction. The variable <italic>e<sub>b</sub> </italic>represents the back-electromotive voltage, where <italic>k<sub>m</sub> </italic>is the back-emf constant. To realize the dynamical model of the rotary servo unit, the mechanical and electrical systems are combined to form the relationship between the angular speed of the motor load shaft <italic>&#x03C9;<sub>l</sub> </italic>and the applied controlled voltage <italic>V<sub>m</sub></italic>. The dynamical system is represented by the following equations
<disp-formula id="eqn-4"><label>(4)</label><mml:math id="mml-eqn-4" display="block"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>&#x03B8;</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mi>l</mml:mi></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:math></disp-formula>
<disp-formula id="eqn-5"><label>(5)</label><mml:math id="mml-eqn-5" display="block"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>&#x03C9;</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mi>l</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>J</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>q</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:msub><mml:mi>B</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>q</mml:mi><mml:mi>v</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:msub><mml:mi>A</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></disp-formula>where the equivalent moment of inertia <italic>J<sub>eq</sub> </italic>is expressed as
<disp-formula id="eqn-6"><label>(6)</label><mml:math id="mml-eqn-6" display="block"><mml:mrow><mml:msub><mml:mi>J</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>q</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:msub><mml:mi>&#x03B7;</mml:mi><mml:mi>g</mml:mi></mml:msub></mml:mrow><mml:msubsup><mml:mi>k</mml:mi><mml:mi>g</mml:mi><mml:mn>2</mml:mn></mml:msubsup><mml:mrow><mml:msub><mml:mi>J</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:msub><mml:mi>J</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:math></disp-formula></p>
<fig id="fig-2"><label>Figure 2</label><caption><title>Rotary servo unit</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_22993-fig-2.png"/></fig>
<p>and the equivalent damping term <italic>B<sub>eqv</sub> </italic>is given by
<disp-formula id="eqn-7"><label>(7)</label><mml:math id="mml-eqn-7" display="block"><mml:mrow><mml:msub><mml:mi>B</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>q</mml:mi><mml:mi>v</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B7;</mml:mi><mml:mi>g</mml:mi></mml:msub></mml:mrow><mml:msubsup><mml:mi>k</mml:mi><mml:mi>g</mml:mi><mml:mn>2</mml:mn></mml:msubsup><mml:mrow><mml:msub><mml:mi>&#x03B7;</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mi>t</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:msub><mml:mi>B</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>q</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mfrac></mml:math></disp-formula></p>
<p>In the expressions of the equivalent moment of inertia and damping, <italic>&#x03B7;<sub>g</sub> </italic>and <italic>&#x03B7;<sub>m</sub> </italic>represent the efficiency of the gearbox and motor, respectively. The motor torque constant is denoted by <italic>k<sub>t,</sub> </italic>and <italic>k<sub>g</sub> </italic>is the gear ratio. The actuator gain <italic>A<sub>m</sub> </italic>is computed as
<disp-formula id="eqn-8"><label>(8)</label><mml:math id="mml-eqn-8" display="block"><mml:mrow><mml:msub><mml:mi>A</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B7;</mml:mi><mml:mi>g</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mi>g</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B7;</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mi>t</mml:mi></mml:msub></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mfrac></mml:math></disp-formula></p>
<p>The numerical value of the BBS parameters is listed in <xref ref-type="table" rid="table-1">Tab. 1</xref>, see [<xref ref-type="bibr" rid="ref-1">1</xref>,<xref ref-type="bibr" rid="ref-26">26</xref>].</p>
<table-wrap id="table-1"><label>Table 1</label><caption><title>BBS specifications</title></caption>
<table frame="hsides">
<colgroup>
<col align="left"/>
<col align="left"/>
<col align="left"/>
</colgroup>
<thead>
<tr>
<th align="left">Parameters</th>
<th align="left">Value</th>
<th align="left">Unit</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left"><inline-formula id="ieqn-57"><mml:math id="mml-ieqn-57"><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mrow><mml:mi>b</mml:mi><mml:mi>e</mml:mi><mml:mi>a</mml:mi><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>, Length of beam</td>
<td align="left"><inline-formula id="ieqn-58"><mml:math id="mml-ieqn-58"><mml:mn>0.4255</mml:mn></mml:math></inline-formula></td>
<td align="left"><inline-formula id="ieqn-59"><mml:math id="mml-ieqn-59"><mml:mi>m</mml:mi></mml:math></inline-formula></td>
</tr>
<tr>
<td align="left"><inline-formula id="ieqn-60"><mml:math id="mml-ieqn-60"><mml:mrow><mml:msub><mml:mi>r</mml:mi><mml:mrow><mml:mi>a</mml:mi><mml:mi>r</mml:mi><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>, Distance between shaft and coupled joint</td>
<td align="left"><inline-formula id="ieqn-61"><mml:math id="mml-ieqn-61"><mml:mn>0.0254</mml:mn></mml:math></inline-formula></td>
<td align="left"><inline-formula id="ieqn-62"><mml:math id="mml-ieqn-62"><mml:mi>m</mml:mi></mml:math></inline-formula></td>
</tr>
<tr>
<td align="left"><inline-formula id="ieqn-63"><mml:math id="mml-ieqn-63"><mml:mrow><mml:msub><mml:mi>r</mml:mi><mml:mi>b</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, Radius of ball</td>
<td align="left"><inline-formula id="ieqn-64"><mml:math id="mml-ieqn-64"><mml:mn>7.68</mml:mn><mml:mrow><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula></td>
<td align="left"><inline-formula id="ieqn-65"><mml:math id="mml-ieqn-65"><mml:mi>N</mml:mi><mml:mi>m</mml:mi><mml:mrow><mml:mo>/</mml:mo></mml:mrow><mml:mi>A</mml:mi></mml:math></inline-formula></td>
</tr>
<tr>
<td align="left"><inline-formula id="ieqn-66"><mml:math id="mml-ieqn-66"><mml:mrow><mml:msub><mml:mi>J</mml:mi><mml:mi>b</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, Ball inertia</td>
<td align="left"><inline-formula id="ieqn-67"><mml:math id="mml-ieqn-67"><mml:mn>4.129</mml:mn><mml:mo>&#x22C5;</mml:mo><mml:mrow><mml:msup><mml:mn>10</mml:mn><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>6</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula></td>
<td align="left"><inline-formula id="ieqn-68"><mml:math id="mml-ieqn-68"><mml:mi>k</mml:mi><mml:mi>g</mml:mi><mml:mo>.</mml:mo><mml:mrow><mml:msup><mml:mi>m</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula></td>
</tr>
<tr>
<td align="left"><inline-formula id="ieqn-69"><mml:math id="mml-ieqn-69"><mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mi>b</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, Mass of ball</td>
<td align="left"><inline-formula id="ieqn-70"><mml:math id="mml-ieqn-70"><mml:mn>0.064</mml:mn></mml:math></inline-formula></td>
<td align="left"><inline-formula id="ieqn-71"><mml:math id="mml-ieqn-71"><mml:mi>k</mml:mi><mml:mi>g</mml:mi></mml:math></inline-formula></td>
</tr>
<tr>
<td align="left"><inline-formula id="ieqn-72"><mml:math id="mml-ieqn-72"><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, Motor armature resistance</td>
<td align="left"><inline-formula id="ieqn-73"><mml:math id="mml-ieqn-73"><mml:mn>2.6</mml:mn></mml:math></inline-formula></td>
<td align="left"><inline-formula id="ieqn-74"><mml:math id="mml-ieqn-74"><mml:mrow><mml:mi>&#x03A9;</mml:mi></mml:mrow></mml:math></inline-formula></td>
</tr>
<tr>
<td align="left"><inline-formula id="ieqn-75"><mml:math id="mml-ieqn-75"><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, Motor armature inductance</td>
<td align="left"><inline-formula id="ieqn-76"><mml:math id="mml-ieqn-76"><mml:mn>0.18</mml:mn></mml:math></inline-formula></td>
<td align="left"><inline-formula id="ieqn-77"><mml:math id="mml-ieqn-77"><mml:mi>m</mml:mi><mml:mi>H</mml:mi></mml:math></inline-formula></td>
</tr>
<tr>
<td align="left"><inline-formula id="ieqn-78"><mml:math id="mml-ieqn-78"><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mi>t</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, Motor current-torque constant</td>
<td align="left"><inline-formula id="ieqn-79"><mml:math id="mml-ieqn-79"><mml:mn>7.68</mml:mn><mml:mo>&#x22C5;</mml:mo><mml:mrow><mml:msup><mml:mn>10</mml:mn><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula></td>
<td align="left"><inline-formula id="ieqn-80"><mml:math id="mml-ieqn-80"><mml:mi>N</mml:mi><mml:mi>m</mml:mi><mml:mrow><mml:mo>/</mml:mo></mml:mrow><mml:mi>A</mml:mi></mml:math></inline-formula></td>
</tr>
<tr>
<td align="left"><inline-formula id="ieqn-81"><mml:math id="mml-ieqn-81"><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, Motor back-emf constant</td>
<td align="left"><inline-formula id="ieqn-82"><mml:math id="mml-ieqn-82"><mml:mn>7.68</mml:mn><mml:mo>&#x22C5;</mml:mo><mml:mrow><mml:msup><mml:mn>10</mml:mn><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula></td>
<td align="left"><inline-formula id="ieqn-83"><mml:math id="mml-ieqn-83"><mml:mi>V</mml:mi><mml:mrow><mml:mo>/</mml:mo></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>r</mml:mi><mml:mi>a</mml:mi><mml:mi>d</mml:mi><mml:mrow><mml:mo>/</mml:mo></mml:mrow><mml:mi>s</mml:mi><mml:mi>e</mml:mi><mml:mi>c</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula></td>
</tr>
<tr>
<td align="left"><inline-formula id="ieqn-84"><mml:math id="mml-ieqn-84"><mml:mrow><mml:msub><mml:mi>J</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, Rotor moment of inertia</td>
<td align="left"><inline-formula id="ieqn-85"><mml:math id="mml-ieqn-85"><mml:mn>3.9</mml:mn><mml:mo>&#x22C5;</mml:mo><mml:mrow><mml:msup><mml:mn>10</mml:mn><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>7</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula></td>
<td align="left"><inline-formula id="ieqn-86"><mml:math id="mml-ieqn-86"><mml:mi>k</mml:mi><mml:mi>g</mml:mi><mml:mrow><mml:msup><mml:mi>m</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula></td>
</tr>
<tr>
<td align="left"><inline-formula id="ieqn-87"><mml:math id="mml-ieqn-87"><mml:mrow><mml:msub><mml:mi>J</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>q</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>, Gear moment of inertia</td>
<td align="left"><inline-formula id="ieqn-88"><mml:math id="mml-ieqn-88"><mml:mn>9.76</mml:mn><mml:mo>&#x22C5;</mml:mo><mml:mrow><mml:msup><mml:mn>10</mml:mn><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>5</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula></td>
<td align="left"><inline-formula id="ieqn-89"><mml:math id="mml-ieqn-89"><mml:mi>k</mml:mi><mml:mi>g</mml:mi><mml:mrow><mml:msup><mml:mi>m</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula></td>
</tr>
<tr>
<td align="left"><inline-formula id="ieqn-90"><mml:math id="mml-ieqn-90"><mml:mrow><mml:msub><mml:mi>B</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>q</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>, Gear viscous damping</td>
<td align="left"><inline-formula id="ieqn-91"><mml:math id="mml-ieqn-91"><mml:mn>0.015</mml:mn></mml:math></inline-formula></td>
<td align="left"><inline-formula id="ieqn-92"><mml:math id="mml-ieqn-92"><mml:mi>N</mml:mi><mml:mi>m</mml:mi><mml:mrow><mml:mo>/</mml:mo></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>r</mml:mi><mml:mi>a</mml:mi><mml:mi>d</mml:mi><mml:mrow><mml:mo>/</mml:mo></mml:mrow><mml:mi>s</mml:mi><mml:mi>e</mml:mi><mml:mi>c</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula></td>
</tr>
<tr>
<td align="left"><inline-formula id="ieqn-93"><mml:math id="mml-ieqn-93"><mml:mrow><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>, Maximum input current</td>
<td align="left"><inline-formula id="ieqn-94"><mml:math id="mml-ieqn-94"><mml:mn>1</mml:mn></mml:math></inline-formula></td>
<td align="left"><inline-formula id="ieqn-95"><mml:math id="mml-ieqn-95"><mml:mi>A</mml:mi></mml:math></inline-formula></td>
</tr>
<tr>
<td align="left"><inline-formula id="ieqn-96"><mml:math id="mml-ieqn-96"><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>, Maximum motor speed</td>
<td align="left"><inline-formula id="ieqn-97"><mml:math id="mml-ieqn-97"><mml:mn>628.3</mml:mn></mml:math></inline-formula></td>
<td align="left"><inline-formula id="ieqn-98"><mml:math id="mml-ieqn-98"><mml:mi>r</mml:mi><mml:mi>a</mml:mi><mml:mi>d</mml:mi><mml:mrow><mml:mo>/</mml:mo></mml:mrow><mml:mi>s</mml:mi><mml:mi>e</mml:mi><mml:mi>c</mml:mi></mml:math></inline-formula></td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3"><label>3</label><title>Cascade Control Architecture</title>
<p>A cascade control architecture based on the time separation principle is proposed to control the ball&#x0027;s position on the beam, as shown in <xref ref-type="fig" rid="fig-3">Fig. 3</xref>. Based on the measured ball position <italic>x</italic>, the outer positional loop computes the desired angle <italic>&#x03B8;<sub>ld</sub> </italic>of the servo load shaft by engaging the classical PD control. The inner fast attitude loop incorporates IDI control to track the desired angle <italic>&#x03B8;<sub>ld</sub> </italic>by generating the corresponding motor controlled voltage <italic>V</italic> such that the desired ball position on the beam is achieved in minimal time.</p>
<fig id="fig-3"><label>Figure 3</label><caption><title>Cascade control architecture</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_22993-fig-3.png"/></fig>
</sec>
<sec id="s4"><label>4</label><title>Outer Loop PD Control Design</title>
<p>Given the desired and actual ball&#x0027;s position, i.e., <italic>x<sub>d</sub> </italic>and <italic>x,</italic> the classical PD controller is engaged in the outer loop to generate the desired servo angle command <italic>&#x03B8;<sub>ld</sub> </italic>for the inner loop to minimize the positional error of the ball placed on the beam. The asymptotically stable virtual (desired) dynamics to compute <italic>&#x03B8;<sub>ld</sub> </italic>based on the positional error of the ball in the horizontal plane is expressed as
<disp-formula id="eqn-9"><label>(9)</label><mml:math id="mml-eqn-9" display="block"><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mi>d</mml:mi></mml:msub></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mrow><mml:mover><mml:mi>x</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mover><mml:mi>x</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mi>d</mml:mi></mml:msub></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mo>+</mml:mo><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mi>p</mml:mi></mml:msub></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>x</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mi>d</mml:mi></mml:msub></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>d</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></disp-formula>where <italic>k<sub>p</sub> </italic>and <italic>k<sub>d</sub> </italic>are controller gains chosen to calculate the suitable value of <italic>&#x03B8;<sub>ld</sub> </italic>such that <italic>x</italic> approached <italic>x<sub>d</sub> </italic>in minimum time. It is noteworthy to mention here that the position feedback of the ball is acquired from an analog sensor with its inherent noise. By taking its derivative, it will generate an amplified high-frequency signal transmitted to the inner loop via feedback which causes the grinding noise in the motor. The derivative is replaced by a first-order filter governed following dynamics to compensate for this effect.
<disp-formula id="eqn-10"><label>(10)</label><mml:math id="mml-eqn-10" display="block"><mml:mi>H</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>s</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>f</mml:mi></mml:msub></mml:mrow><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mi>s</mml:mi><mml:mo>+</mml:mo><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>f</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mfrac></mml:math></disp-formula>where <italic>&#x03C9;<sub>f</sub> </italic>represents the cutoff frequency, whose value is set to be 6.28 rad/s.</p>
</sec>
<sec id="s5"><label>5</label><title>Attitude Control Using IDI</title>
<p>This section discussed the design and implementation of IDI control for attitude control of the rotary servo unit of BBS. The proposed non-square inversion-based intelligent controller proves to be an efficient and effective control strategy to deal with complex nonlinear systems. The control law comprises two major components, i.e., an equivalent or baseline control and a switching control component. The general form of IDI control law is expressed as follows
<disp-formula id="eqn-11"><label>(11)</label><mml:math id="mml-eqn-11" display="block"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>I</mml:mi><mml:mi>D</mml:mi><mml:mi>I</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>q</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>r</mml:mi><mml:mi>b</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></disp-formula></p>
<p>The first component <italic>u<sub>eq</sub> </italic>represents the baseline controller, which enforces the prescribed constraint dynamics based on the attitude deviation function for stable attitude tracking. On the other hand, <italic>u<sub>rbt</sub> </italic>represents the discontinuous control element whose design principle is based on the classical SMC approach to address robustness against uncertain dynamics, parametric variations, and external disturbances.</p>
<sec id="s5_1"><label>5.1</label><title>Baseline Control Design</title>
<p>The equations of motion of rotary servo unit angular rate dynamics given by <xref ref-type="disp-formula" rid="eqn-5">Eq. (5)</xref> is expressed as
<disp-formula id="eqn-12"><label>(12)</label><mml:math id="mml-eqn-12" display="block"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>&#x03C9;</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mi>l</mml:mi></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mi>F</mml:mi><mml:mo>+</mml:mo><mml:mi>G</mml:mi><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow></mml:math></disp-formula>where <inline-formula id="ieqn-1"><mml:math id="mml-ieqn-1"><mml:mi>F</mml:mi><mml:mo>=</mml:mo><mml:msubsup><mml:mi>J</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>q</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msub><mml:mi>B</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>q</mml:mi><mml:mi>v</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula id="ieqn-2"><mml:math id="mml-ieqn-2"><mml:mi>G</mml:mi><mml:mo>=</mml:mo><mml:msubsup><mml:mi>J</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>q</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msub><mml:mi>A</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, and <inline-formula id="ieqn-3"><mml:math id="mml-ieqn-3"><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> represent the controlled voltage. The design is initiated by determining two constraint differential equations based on the attitude deviation function, which implies
<disp-formula id="eqn-13"><label>(13)</label><mml:math id="mml-eqn-13" display="block"><mml:msub><mml:mrow><mml:mover><mml:mi>&#x03B6;</mml:mi><mml:mo>&#x00A8;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mrow><mml:mi>l</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>&#x03B6;</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mrow><mml:mi>l</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B6;</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mrow><mml:mi>l</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:math></disp-formula>
<disp-formula id="eqn-14"><label>(14)</label><mml:math id="mml-eqn-14" display="block"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>&#x03B6;</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B6;</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:math></disp-formula>where <inline-formula id="ieqn-4"><mml:math id="mml-ieqn-4"><mml:mrow><mml:msub><mml:mi>&#x03B6;</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:msubsup><mml:mi>e</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow><mml:mn>2</mml:mn></mml:msubsup><mml:mo>=</mml:mo><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>d</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:msup><mml:mo stretchy="false">)</mml:mo><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> and <inline-formula id="ieqn-5"><mml:math id="mml-ieqn-5"><mml:mrow><mml:msub><mml:mi>&#x03B6;</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:msubsup><mml:mi>e</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow><mml:mn>2</mml:mn></mml:msubsup><mml:mo>=</mml:mo><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>d</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:msup><mml:mo stretchy="false">)</mml:mo><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> are the deviation functions and <inline-formula id="ieqn-6"><mml:math id="mml-ieqn-6"><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the angular speed. The order of the constraint equation is equal to the relative degree of the respective controlled variable. The coefficients <inline-formula id="ieqn-7"><mml:math id="mml-ieqn-7"><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula id="ieqn-8"><mml:math id="mml-ieqn-8"><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>, and <inline-formula id="ieqn-9"><mml:math id="mml-ieqn-9"><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> are selected such that the characteristic polynomials are given by <xref ref-type="disp-formula" rid="eqn-13">Eqs. (13)</xref>, and <xref ref-type="disp-formula" rid="eqn-14">(14)</xref> is Hurwitz, such that <inline-formula id="ieqn-10"><mml:math id="mml-ieqn-10"><mml:mrow><mml:msub><mml:mi>&#x03B6;</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">&#x2192;</mml:mo><mml:mn>0</mml:mn></mml:math></inline-formula> and <inline-formula id="ieqn-11"><mml:math id="mml-ieqn-11"><mml:mrow><mml:msub><mml:mi>&#x03B6;</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">&#x2192;</mml:mo><mml:mn>0</mml:mn></mml:math></inline-formula> as <inline-formula id="ieqn-12"><mml:math id="mml-ieqn-12"><mml:mi>t</mml:mi><mml:mo stretchy="false">&#x2192;</mml:mo><mml:mi mathvariant="normal">&#x221E;</mml:mi></mml:math></inline-formula> [<xref ref-type="bibr" rid="ref-19">19</xref>]. Computing and substituting the time derivatives in the constraint dynamics given by <xref ref-type="disp-formula" rid="eqn-13">Eqs. (13)</xref> and <xref ref-type="disp-formula" rid="eqn-14">(14)</xref> yields the following under-determined equivalent algebraic system.
<disp-formula id="eqn-15"><label>(15)</label><mml:math id="mml-eqn-15" display="block"><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:mrow><mml:mrow><mml:msub><mml:mi>A</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mrow><mml:msub><mml:mi>A</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:mrow><mml:mrow><mml:msub><mml:mi>B</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mrow><mml:msub><mml:mi>B</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo stretchy="false">&#x21D2;</mml:mo><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mrow><mml:mi mathvariant="script">B</mml:mi></mml:mrow></mml:mrow></mml:math></disp-formula>where the elements of the controls coefficient column vector function <inline-formula id="ieqn-13"><mml:math id="mml-ieqn-13"><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mo>:</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="bold">R</mml:mi></mml:mrow></mml:mrow><mml:mi>n</mml:mi></mml:msup></mml:mrow><mml:mo>&#x00D7;</mml:mo><mml:mo stretchy="false">[</mml:mo><mml:mrow><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow><mml:mo>,</mml:mo><mml:mi mathvariant="normal">&#x221E;</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">&#x2192;</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="bold">R</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mo>&#x00D7;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> is given by
<disp-formula id="ueqn-1">
<mml:math id="mml-ueqn-1" display="block"><mml:mrow><mml:msub><mml:mi>A</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>2</mml:mn><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mi>G</mml:mi></mml:math></disp-formula>
<disp-formula id="ueqn-2">
<mml:math id="mml-ueqn-2" display="block"><mml:mrow><mml:msub><mml:mi>A</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>2</mml:mn><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mi>G</mml:mi></mml:math></disp-formula>and the elements of the controls load vector function <inline-formula id="ieqn-14"><mml:math id="mml-ieqn-14"><mml:mrow><mml:mrow><mml:mi mathvariant="script">B</mml:mi></mml:mrow></mml:mrow><mml:mo>:</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="bold">R</mml:mi></mml:mrow></mml:mrow><mml:mi>n</mml:mi></mml:msup></mml:mrow><mml:mo>&#x00D7;</mml:mo><mml:mo stretchy="false">[</mml:mo><mml:mrow><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow><mml:mo>,</mml:mo><mml:mi mathvariant="normal">&#x221E;</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">&#x2192;</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="bold">R</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mo>&#x00D7;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> are given as
<disp-formula id="ueqn-3">
<mml:math id="mml-ueqn-3" display="block"><mml:mrow><mml:msub><mml:mi>B</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>2</mml:mn><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>&#x03C9;</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>l</mml:mi><mml:mi>d</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>2</mml:mn><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:msubsup><mml:mrow><mml:mover><mml:mi>e</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow><mml:mn>2</mml:mn></mml:msubsup><mml:mo>&#x2212;</mml:mo><mml:mn>2</mml:mn><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mi>F</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>2</mml:mn><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>e</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:msubsup><mml:mi>e</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow><mml:mn>2</mml:mn></mml:msubsup></mml:math></disp-formula>
<disp-formula id="ueqn-4">
<mml:math id="mml-ueqn-4" display="block"><mml:mrow><mml:msub><mml:mi>B</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>2</mml:mn><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>&#x03C9;</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>l</mml:mi><mml:mi>d</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>2</mml:mn><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mi>F</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:msubsup><mml:mi>e</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow><mml:mn>2</mml:mn></mml:msubsup></mml:math></disp-formula></p>
<p>The baseline or equivalent controlled voltage <inline-formula id="ieqn-15"><mml:math id="mml-ieqn-15"><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>q</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is computed by inverting the algebraic system given by <xref ref-type="disp-formula" rid="eqn-15">Eq. (15)</xref> using MPGI yields the following expression
<disp-formula id="eqn-16"><label>(16)</label><mml:math id="mml-eqn-16" display="block"><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>q</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mo>+</mml:mo></mml:msup></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="script">B</mml:mi></mml:mrow></mml:mrow><mml:mo>+</mml:mo><mml:mi>P</mml:mi><mml:mi>&#x03BB;</mml:mi></mml:math></disp-formula>where <inline-formula id="ieqn-16"><mml:math id="mml-ieqn-16"><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mo>+</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula> represents the MPGI of <inline-formula id="ieqn-17"><mml:math id="mml-ieqn-17"><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow></mml:math></inline-formula> defined as
<disp-formula id="eqn-17"><label>(17)</label><mml:math id="mml-eqn-17" display="block"><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mo>+</mml:mo></mml:msup></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo fence="false" stretchy="false">{</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mi>T</mml:mi></mml:msup></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:msup><mml:mo fence="false" stretchy="false">}</mml:mo><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mi>T</mml:mi></mml:msup></mml:mrow></mml:math></disp-formula></p>
<p><bold>Remark 1:</bold> The matrix <inline-formula id="ieqn-18"><mml:math id="mml-ieqn-18"><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mo>&#x00D7;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> appeared to be a tall matrix, which causes the value of <italic>P</italic> to become zero by using the property of pseudo-inverse <inline-formula id="ieqn-19"><mml:math id="mml-ieqn-19"><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mo>+</mml:mo></mml:msup></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:math></inline-formula>. Hence the under-determined nature of constraint dynamics does not qualify with the benefit of the null control vector due to the ineffectiveness of the null projection matrix. Hence, the null control is supposed to be incompetent and not considered in the present control design methodology.</p>
</sec>
<sec id="s5_2"><label>5.2</label><title>Switching Control Design</title>
<p>A suitable discontinuous control element is determined to solve the reachability problem, which ensures robustness against parametric uncertainties, system nonlinearities, and external disturbances. The resultant GDI-SMC or IDI control law is constructed to be of the following form
<disp-formula id="ueqn-5">
<mml:math id="mml-ueqn-5" display="block"><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>q</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mo>+</mml:mo></mml:msup></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="script">K</mml:mi></mml:mrow></mml:mrow><mml:mfrac><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mo stretchy="false">&#x2225;</mml:mo><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow><mml:mo stretchy="false">&#x2225;</mml:mo></mml:mrow></mml:mfrac></mml:math></disp-formula>
<disp-formula id="eqn-18"><label>(18)</label><mml:math id="mml-eqn-18" display="block"><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mo>+</mml:mo></mml:msup></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="script">B</mml:mi></mml:mrow></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mo>+</mml:mo></mml:msup></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="script">K</mml:mi></mml:mrow></mml:mrow><mml:mfrac><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mo stretchy="false">&#x2225;</mml:mo><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow><mml:mo stretchy="false">&#x2225;</mml:mo></mml:mrow></mml:mfrac></mml:math></disp-formula>where <inline-formula id="ieqn-20"><mml:math id="mml-ieqn-20"><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo stretchy="false">[</mml:mo><mml:mrow><mml:msub><mml:mi>s</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mrow><mml:msub><mml:mi>s</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:msup><mml:mo stretchy="false">]</mml:mo><mml:mi>T</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula> is an appropriate switching surface vector, whose vector elements are defined as
<disp-formula id="eqn-19"><label>(19)</label><mml:math id="mml-eqn-19" display="block"><mml:mrow><mml:msub><mml:mi>s</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>&#x03BE;</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03BE;</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:msup><mml:mrow><mml:mo largeop="false">&#x222B;</mml:mo></mml:mrow><mml:mrow></mml:mrow></mml:msup><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:msub><mml:mi>&#x03BE;</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:math></disp-formula>
<disp-formula id="eqn-20"><label>(20)</label><mml:math id="mml-eqn-20" display="block"><mml:mrow><mml:msub><mml:mi>s</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:msub><mml:mi>&#x03BE;</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow><mml:msup><mml:mrow><mml:mo largeop="false">&#x222B;</mml:mo></mml:mrow><mml:mrow></mml:mrow></mml:msup><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:msub><mml:mi>&#x03BE;</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:math></disp-formula></p>
<p>Moreover, <inline-formula id="ieqn-21"><mml:math id="mml-ieqn-21"><mml:mrow><mml:mrow><mml:mi mathvariant="script">K</mml:mi></mml:mrow></mml:mrow></mml:math></inline-formula> is <inline-formula id="ieqn-22"><mml:math id="mml-ieqn-22"><mml:mn>2</mml:mn><mml:mo>&#x00D7;</mml:mo><mml:mn>2</mml:mn></mml:math></inline-formula> diagonal gain matrix with positive valued diagonal elements <inline-formula id="ieqn-23"><mml:math id="mml-ieqn-23"><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula id="ieqn-24"><mml:math id="mml-ieqn-24"><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>, ensuring that the sliding mode dynamics are maintained and attained. Therefore, the controlled voltage <inline-formula id="ieqn-25"><mml:math id="mml-ieqn-25"><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> given by the <xref ref-type="disp-formula" rid="eqn-18">Eq. (18)</xref> is reasonably competent to attract the sliding surfaces given by <xref ref-type="disp-formula" rid="eqn-19">Eqs. (19)</xref> and <xref ref-type="disp-formula" rid="eqn-20">(20)</xref> and enforce the system trajectories to slide along the sliding surfaces.</p>
</sec>
</sec>
<sec id="s6"><label>6</label><title>Stability Analysis</title>
<p>The time derivatives of the sliding surfaces are evaluated as
<disp-formula id="eqn-21"><label>(21)</label><mml:math id="mml-eqn-21" display="block"><mml:msub><mml:mrow><mml:mover><mml:mi>s</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>&#x03BE;</mml:mi><mml:mo>&#x00A8;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mrow><mml:mi>l</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>&#x03BE;</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03BE;</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></disp-formula>
<disp-formula id="eqn-22"><label>(22)</label><mml:math id="mml-eqn-22" display="block"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>s</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>&#x03BE;</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03BE;</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></disp-formula></p>
<p>It has been observed that the asymptotic convergence of the sliding surfaces <inline-formula id="ieqn-26"><mml:math id="mml-ieqn-26"><mml:mrow><mml:msub><mml:mi>s</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula id="ieqn-27"><mml:math id="mml-ieqn-27"><mml:mrow><mml:msub><mml:mi>s</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> to zero implies the asymptotic realization of the constraint dynamics given by the <xref ref-type="disp-formula" rid="eqn-13">Eqs. (13)</xref> and <xref ref-type="disp-formula" rid="eqn-14">(14)</xref>, and its equivalent algebraic form given by <xref ref-type="disp-formula" rid="eqn-15">Eq. (15)</xref>. Hence, the sliding vector <inline-formula id="ieqn-28"><mml:math id="mml-ieqn-28"><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow></mml:math></inline-formula> is defined as
<disp-formula id="eqn-23"><label>(23)</label><mml:math id="mml-eqn-23" display="block"><mml:mrow><mml:mrow><mml:mover><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mo mathvariant="bold">&#x02D9;</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:mrow><mml:mi mathvariant="script">B</mml:mi></mml:mrow></mml:mrow></mml:math></disp-formula></p>
<p>To prove the asymptotic convergence of the error dynamics, place the control law <inline-formula id="ieqn-29"><mml:math id="mml-ieqn-29"><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mi>m</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> given by <xref ref-type="disp-formula" rid="eqn-18">Eq. (18)</xref> in sliding mode dynamics given by <xref ref-type="disp-formula" rid="eqn-23">Eq. (23)</xref>, resulting in
<disp-formula id="ueqn-6">
<mml:math id="mml-ueqn-6" display="block"><mml:mrow><mml:mrow><mml:mover><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mo mathvariant="bold">&#x02D9;</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mo>+</mml:mo></mml:msup></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="script">B</mml:mi></mml:mrow></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mo>+</mml:mo></mml:msup></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="script">K</mml:mi></mml:mrow></mml:mrow><mml:mfrac><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mo stretchy="false">&#x2225;</mml:mo><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow><mml:mo stretchy="false">&#x2225;</mml:mo></mml:mrow></mml:mfrac></mml:mrow><mml:mo>}</mml:mo></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:mrow><mml:mi mathvariant="script">B</mml:mi></mml:mrow></mml:mrow></mml:math></disp-formula>
<disp-formula id="eqn-24"><label>(24)</label><mml:math id="mml-eqn-24" display="block"><mml:mo>=</mml:mo><mml:mo fence="false" stretchy="false">{</mml:mo><mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mo>+</mml:mo></mml:msup></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>2</mml:mn><mml:mo>&#x00D7;</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mrow><mml:mo fence="false" stretchy="false">}</mml:mo><mml:mrow><mml:mrow><mml:mi mathvariant="script">B</mml:mi></mml:mrow></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mo>+</mml:mo></mml:msup></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="script">K</mml:mi></mml:mrow></mml:mrow><mml:mfrac><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mo stretchy="false">&#x2225;</mml:mo><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow><mml:mo stretchy="false">&#x2225;</mml:mo></mml:mrow></mml:mfrac></mml:math></disp-formula></p>
<p>Consider the following positive definite Lyapunov candidate function
<disp-formula id="eqn-25"><label>(25)</label><mml:math id="mml-eqn-25" display="block"><mml:mi>V</mml:mi><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mn>2</mml:mn></mml:mfrac><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow><mml:mi>T</mml:mi></mml:msup></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow></mml:math></disp-formula></p>
<p>The time derivative of <italic>V</italic> is computed as
<disp-formula id="ueqn-7">
<mml:math id="mml-ueqn-7" display="block"><mml:mrow><mml:mover><mml:mi>V</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow><mml:mi>T</mml:mi></mml:msup></mml:mrow><mml:mrow><mml:mrow><mml:mover><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mo mathvariant="bold">&#x02D9;</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow><mml:mi>T</mml:mi></mml:msup></mml:mrow><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mo>+</mml:mo></mml:msup></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>2</mml:mn><mml:mo>&#x00D7;</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mrow><mml:mrow><mml:mi mathvariant="script">B</mml:mi></mml:mrow></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mo>+</mml:mo></mml:msup></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="script">K</mml:mi></mml:mrow></mml:mrow><mml:mfrac><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mo stretchy="false">&#x2225;</mml:mo><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow><mml:mo stretchy="false">&#x2225;</mml:mo></mml:mrow></mml:mfrac></mml:mrow><mml:mo>}</mml:mo></mml:mrow></mml:math></disp-formula>
<disp-formula id="eqn-26"><label>(26)</label><mml:math id="mml-eqn-26" display="block"><mml:mo>=</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow><mml:mi>T</mml:mi></mml:msup></mml:mrow><mml:mo fence="false" stretchy="false">{</mml:mo><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mo>+</mml:mo></mml:msup></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>2</mml:mn><mml:mo>&#x00D7;</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mo fence="false" stretchy="false">}</mml:mo><mml:mrow><mml:mrow><mml:mi mathvariant="script">B</mml:mi></mml:mrow></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow><mml:mi>T</mml:mi></mml:msup></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mo>+</mml:mo></mml:msup></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="script">K</mml:mi></mml:mrow></mml:mrow><mml:mfrac><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mo stretchy="false">&#x2225;</mml:mo><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow><mml:mo stretchy="false">&#x2225;</mml:mo></mml:mrow></mml:mfrac></mml:math></disp-formula></p>
<p>The 2-(induced) norms of <inline-formula id="ieqn-30"><mml:math id="mml-ieqn-30"><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mo>+</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula> and <inline-formula id="ieqn-31"><mml:math id="mml-ieqn-31"><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mi mathvariant="script">A</mml:mi></mml:mrow></mml:mrow><mml:mo>+</mml:mo></mml:msup></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>2</mml:mn><mml:mo>&#x00D7;</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> are unity. It follows that <inline-formula id="ieqn-32"><mml:math id="mml-ieqn-32"><mml:mrow><mml:mover><mml:mi>V</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mo>&#x2264;</mml:mo><mml:mn>0</mml:mn></mml:math></inline-formula> is guaranteed if <inline-formula id="ieqn-33"><mml:math id="mml-ieqn-33"><mml:mrow><mml:msub><mml:mi>&#x03C3;</mml:mi><mml:mrow><mml:mrow><mml:mtext>min</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mrow><mml:mi mathvariant="script">K</mml:mi></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mo>&#x003E;&#x2225;</mml:mo><mml:mrow><mml:mrow><mml:mi mathvariant="script">B</mml:mi></mml:mrow></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="bold">x</mml:mi></mml:mrow></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mn>0</mml:mn><mml:mo stretchy="false">)</mml:mo><mml:mo>,</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">&#x2225;</mml:mo></mml:math></inline-formula>, where <inline-formula id="ieqn-34"><mml:math id="mml-ieqn-34"><mml:mrow><mml:msub><mml:mi>&#x03C3;</mml:mi><mml:mrow><mml:mrow><mml:mtext>min</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is the minimum singular value function. Moreover, because <inline-formula id="ieqn-35"><mml:math id="mml-ieqn-35"><mml:mrow><mml:mover><mml:mi>V</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:math></inline-formula> occurs only at <inline-formula id="ieqn-36"><mml:math id="mml-ieqn-36"><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:msub><mml:mn>0</mml:mn><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>, it follows Lassale&#x0027;s principle [<xref ref-type="bibr" rid="ref-27">27</xref>] that <inline-formula id="ieqn-37"><mml:math id="mml-ieqn-37"><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:msub><mml:mn>0</mml:mn><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> is attractive. Moreover, the discontinuous term in the right-hand side of <xref ref-type="disp-formula" rid="eqn-26">Eq. (26)</xref> implies that this condition on <inline-formula id="ieqn-38"><mml:math id="mml-ieqn-38"><mml:mrow><mml:mrow><mml:mi mathvariant="script">K</mml:mi></mml:mrow></mml:mrow></mml:math></inline-formula> guarantees finite-time convergence of <inline-formula id="ieqn-39"><mml:math id="mml-ieqn-39"><mml:mrow><mml:mrow><mml:mi mathvariant="bold">s</mml:mi></mml:mrow></mml:mrow></mml:math></inline-formula> and <inline-formula id="ieqn-40"><mml:math id="mml-ieqn-40"><mml:mrow><mml:mrow><mml:mover><mml:mi mathvariant="bold-italic">s</mml:mi><mml:mo mathvariant="bold">&#x02D9;</mml:mo></mml:mover></mml:mrow></mml:mrow></mml:math></inline-formula> to the zero vectors, which implies asymptotic convergence of <inline-formula id="ieqn-41"><mml:math id="mml-ieqn-41"><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula id="ieqn-42"><mml:math id="mml-ieqn-42"><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> to zero. Finally, because it is assumed that the matrix gain <inline-formula id="ieqn-43"><mml:math id="mml-ieqn-43"><mml:mrow><mml:mrow><mml:mi mathvariant="script">K</mml:mi></mml:mrow></mml:mrow></mml:math></inline-formula> can be increased arbitrarily such that the stability condition is satisfied for any initial state condition, it follows that the design guarantees semi-global asymptotic stability of <inline-formula id="ieqn-44"><mml:math id="mml-ieqn-44"><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>,</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mover><mml:mi>e</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>,</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mover><mml:mi>e</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula>.</p>
</sec>
<sec id="s7"><label>7</label><title>Numerical Simulation</title>
<p>This section investigates the performance of IDI control law through computer simulations on the dynamical simulator of BBS developed by considering the model of ball and beam unit given by <xref ref-type="disp-formula" rid="eqn-3">Eq. (3)</xref> and the model of rotary servo unit given by <xref ref-type="disp-formula" rid="eqn-4">Eqs. (4)</xref> and <xref ref-type="disp-formula" rid="eqn-5">(5)</xref>. The simulation environment is created in Simulink/Matlab, in which the closed-loop performance is evaluated by commanding the square wave profile having an amplitude of <inline-formula id="ieqn-45"><mml:math id="mml-ieqn-45"><mml:mo>&#x00B1;</mml:mo><mml:mn>20</mml:mn><mml:mi>m</mml:mi></mml:math></inline-formula> with <inline-formula id="ieqn-46"><mml:math id="mml-ieqn-46"><mml:mn>0.05</mml:mn><mml:mi>H</mml:mi><mml:mi>z</mml:mi></mml:math></inline-formula> frequency, assuming nominal system parameters. The initial condition of the system state vector is set as <inline-formula id="ieqn-47"><mml:math id="mml-ieqn-47"><mml:mrow><mml:mo>[</mml:mo><mml:mi>x</mml:mi><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mrow><mml:mover><mml:mi>x</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mover><mml:mi>&#x03B8;</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mi>l</mml:mi></mml:msub></mml:mrow><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mover><mml:mi>&#x03C9;</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mi>l</mml:mi></mml:msub></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mn>0</mml:mn><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mn>0</mml:mn><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mn>0</mml:mn><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mn>0</mml:mn><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mn>0</mml:mn><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mn>0</mml:mn><mml:mo>]</mml:mo></mml:mrow></mml:math></inline-formula> The time histories of the position tracking curve of the ball is shown in <xref ref-type="fig" rid="fig-4">Fig. 4</xref>, which is quite appealing. The inner loop rotary servo angle tracking performance using IDI control is depicted in <xref ref-type="fig" rid="fig-5">Fig. 5</xref>, which demonstrates faster convergence towards the angular reference command. Finally, the servo motor controlled voltage well within the saturation limit is displayed in <xref ref-type="fig" rid="fig-6">Fig. 6</xref>.</p>
<fig id="fig-4"><label>Figure 4</label><caption><title>Simulations: ball position</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_22993-fig-4.png"/></fig>
<fig id="fig-5"><label>Figure 5</label><caption><title>Simulations: beam angle</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_22993-fig-5.png"/></fig>
<fig id="fig-6"><label>Figure 6</label><caption><title>Simulations: controlled voltage</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_22993-fig-6.png"/></fig>
</sec>
<sec id="s8"><label>8</label><title>Experimental Results</title>
<p>To investigate the real-time performance of balancing the ball on a beam, the IDI control law developed in Section 5 is implemented on Quanser&#x0027;s ball and beam real hardware, as shown in <xref ref-type="fig" rid="fig-7">Fig. 7</xref>, which proves to be a benchmark testbed for testing and validating control techniques.</p>
<fig id="fig-7"><label>Figure 7</label><caption><title>Experimental setup</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_22993-fig-7.png"/></fig>
<p>The experimental setup comprised the power module VoltPAQ-X1 for supplying voltage to the rotary servo motor module. The servo motor is placed at the right corner of the beam to move the beam in an upward and downward direction. The ball placed on the beam rolls on the track accordingly by changing the slope of the beam. The ball position and the angular position feedback commands are sent to the computer through Quanser&#x0027;s Q8-USB data acquisition card. The communication between Simulink/MATLAB and the experimental hardware is accomplished through the QUARC software. The sampling time is chosen to be <inline-formula id="ieqn-48"><mml:math id="mml-ieqn-48"><mml:mn>2</mml:mn><mml:mi>m</mml:mi><mml:mi>s</mml:mi></mml:math></inline-formula> for all the experiments.</p>
<sec id="s8_1"><label>8.1</label><title>Step Response Analysis</title>
<p>In this analysis, the core objective is to present the comparative analysis of IDI control with classical Linear Quadratic Regulator (LQR) and Fractional Order Control (FOC) strategies. The ball is commanded to follow the step input profile in this setup. The initial position of system states are <inline-formula id="ieqn-49"><mml:math id="mml-ieqn-49"><mml:mi>x</mml:mi><mml:mo>=</mml:mo><mml:mo>&#x2212;</mml:mo><mml:mn>18</mml:mn><mml:mi>c</mml:mi><mml:mi>m</mml:mi></mml:math></inline-formula>, <inline-formula id="ieqn-50"><mml:math id="mml-ieqn-50"><mml:mrow><mml:mover><mml:mi>x</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mi>c</mml:mi><mml:mi>m</mml:mi></mml:math></inline-formula>, <inline-formula id="ieqn-51"><mml:math id="mml-ieqn-51"><mml:mrow><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mi>l</mml:mi></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mi>d</mml:mi><mml:mi>e</mml:mi><mml:mi>g</mml:mi></mml:math></inline-formula>, and <inline-formula id="ieqn-52"><mml:math id="mml-ieqn-52"><mml:mrow><mml:msub><mml:mrow><mml:mover><mml:mi>&#x03B8;</mml:mi><mml:mo>&#x02D9;</mml:mo></mml:mover></mml:mrow><mml:mi>l</mml:mi></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mi>d</mml:mi><mml:mi>e</mml:mi><mml:mi>g</mml:mi></mml:math></inline-formula>, whereas the desired step command is set to be <inline-formula id="ieqn-53"><mml:math id="mml-ieqn-53"><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mi>d</mml:mi></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>10</mml:mn><mml:mi>c</mml:mi><mml:mi>m</mml:mi></mml:math></inline-formula>. The three control methodologies were implemented on a practical BBS test bench, and the tracking curves of ball position, beam angle, and the controlled voltages are shown in <xref ref-type="fig" rid="fig-8 fig-9 fig-10">Figs. 8&#x2013;10</xref>. It is evident from <xref ref-type="fig" rid="fig-8">Fig. 8</xref> that the response curves of positional tracking with IDI control are faster in terms of rising time and settling time.</p>
<fig id="fig-8"><label>Figure 8</label><caption><title>Experimental results: ball position</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_22993-fig-8.png"/></fig>
<fig id="fig-9"><label>Figure 9</label><caption><title>Experimental results: attitude tracking</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_22993-fig-9.png"/></fig>
<fig id="fig-10"><label>Figure 10</label><caption><title>Experimental results: controlled voltage</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_22993-fig-10.png"/></fig>
<p>Furthermore, the FOC suffers from high overshoot, whereas many steady-state errors were observed in the LQR response. On the other hand, the IDI controller demonstrates very stable time-domain performance near reference input with less overshoot and steady-state error. The time evolution of desired and actual angular profiles and the time histories of the controlled voltages generated in response to the step input command are shown in <xref ref-type="fig" rid="fig-9">Figs. 9</xref> and <xref ref-type="fig" rid="fig-10">10</xref>, respectively.</p>
</sec>
<sec id="s8_2"><label>8.2</label><title>Disturbance Rejection Capability</title>
<p>This section highlighted the disturbance rejection capability of the IDI controller when the ball is subjected to some manual disturbances. The ball position on the beam in a steady state is stimulated by tapping the ball at different time instants. The ball position plot under the effect of disturbances is shown in <xref ref-type="fig" rid="fig-11">Fig. 11</xref>. It is evident that when disturbances are applied, the ball position is perturbed but immediately returns to its stable position by control reaction and re-stabilized. The jerks appeared in the angular profile due to manual perturbations, which are depicted in <xref ref-type="fig" rid="fig-12">Fig. 12</xref>. The controlled voltage under the supplied capacity is displayed in <xref ref-type="fig" rid="fig-13">Fig. 13</xref>, whose magnitude fluctuates in response to the applied disturbances. The experimental results demonstrate better disturbance rejection capability and intelligent attributes of IDI control.</p>
<fig id="fig-11"><label>Figure 11</label><caption><title>Experimental results: ball position</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_22993-fig-11.png"/></fig>
<fig id="fig-12"><label>Figure 12</label><caption><title>Experimental results: attitude tracking</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_22993-fig-12.png"/></fig>
<fig id="fig-13"><label>Figure 13</label><caption><title>Experimental results: controlled voltage</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_22993-fig-13.png"/></fig>
</sec>
<sec id="s8_3"><label>8.3</label><title>Square-wave Tracking</title>
<p>In this plot, the controllers&#x2019; capability is further analyzed by commanding a reference square wave profile with an amplitude of <inline-formula id="ieqn-54"><mml:math id="mml-ieqn-54"><mml:mo>&#x00B1;</mml:mo><mml:mn>20</mml:mn><mml:mi>m</mml:mi></mml:math></inline-formula> with <inline-formula id="ieqn-55"><mml:math id="mml-ieqn-55"><mml:mn>0.05</mml:mn><mml:mi>H</mml:mi><mml:mi>z</mml:mi></mml:math></inline-formula> frequency and the variation of 20% in the system dynamical parameters. The time evolution of the critical performance indices is illustrated in <xref ref-type="fig" rid="fig-14 fig-15 fig-16">Figs. 14&#x2013;16</xref>. The smooth tracking response of the ball&#x0027;s position placed on the beam in response to the square wave input command is shown in <xref ref-type="fig" rid="fig-14">Fig. 14</xref>. The rotary servo angular position tracking curve required to steer the ball to follow the square waveform command and the control effort in the form of generated controlled voltage are shown in <xref ref-type="fig" rid="fig-15">Figs. 15</xref> and <xref ref-type="fig" rid="fig-16">16</xref>, respectively.</p>
<fig id="fig-14"><label>Figure 14</label><caption><title>Experimental results: ball position</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_22993-fig-14.png"/></fig>
<fig id="fig-15"><label>Figure 15</label><caption><title>Experimental results: attitude tracking</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_22993-fig-15.png"/></fig>
<fig id="fig-16"><label>Figure 16</label><caption><title>Experimental results: controlled voltage</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_22993-fig-16.png"/></fig>
<p>While examining the different behaviors of the controller response in the simulation and experimental results, it is claimed that the proposed controller proves itself a reasonable control approach for stabilizing the under-actuated electro-mechanical systems.</p>
</sec>
</sec>
<sec id="s9"><label>9</label><title>Conclusion</title>
<p>This paper investigates the performance of IDI control for a standard nonlinear ball and beam unstable system. A two-loops structured control is implemented in which PD control is responsible for generating the desired servo angle command based on the positional error of the ball placed on the beam. Furthermore, IDI control is engaged in the inner loop for precise attitude tracking. The results are verified through computer simulations and experimental studies in a MATLAB environment. Furthermore, a comparative analysis is presented, which exhibits the superior performance of IDI control over LQR and FOC in terms of faster error convergence with lesser steady-state error. Hence the proposed methodology is quite appealing and feasible for practical electro-mechanical systems.</p>
</sec>
</body>
<back>
<ack>
<p>The authors extend their appreciation to the Deputyship for Research &#x0026; Innovation, Ministry of Education in Saudi Arabia for funding this research work through the Project Number (IFPRC-023-135-2020) and King Abdulaziz University, DSR, Jeddah, Saudi Arabia.</p>
</ack>
<fn-group>
<fn fn-type="other"><p><bold>Funding Statement:</bold> This research work was funded by Deputyship for Research &#x0026; Innovation, Ministry of Education in Saudi Arabia under Grant No. (IFPRC-023-135-2020).</p></fn>
<fn fn-type="conflict"><p><bold>Conflicts of Interest:</bold> The authors declare that they have no conflicts of interest to report regarding the present study.</p></fn>
</fn-group>
<ref-list content-type="authoryear">
<title>References</title>
<ref id="ref-1"><label>[1]</label><mixed-citation publication-type="book"><person-group person-group-type="author"><string-name><given-names>J.</given-names> <surname>Apkarian</surname></string-name>, <string-name><given-names>P.</given-names> <surname>Karam</surname></string-name>, <string-name><given-names>M.</given-names> <surname>Levis</surname></string-name> and <string-name><given-names>H.</given-names> <surname>Gurocak</surname></string-name></person-group>, <source>Quanser Innovative Edutech Laboratory Guide-Ball and Beam Experiment</source>, <publisher-loc>Canada</publisher-loc>: <publisher-name>Quanser Inc.</publisher-name>, <year>2012</year>.</mixed-citation></ref>
<ref id="ref-2"><label>[2]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>I. M.</given-names> <surname>Mehedi</surname></string-name>, <string-name><given-names>U. M.</given-names> <surname>Al-Saggaf</surname></string-name>, <string-name><given-names>R.</given-names> <surname>Mansouri</surname></string-name> and <string-name><given-names>M.</given-names> <surname>Bettayeb</surname></string-name></person-group>, &#x201C;<article-title>Two degrees of freedom fractional controller design: Application to the ball and beam system</article-title>,&#x201D; <source>Measurement</source>, vol. <volume>135</volume>, pp. <fpage>13</fpage>&#x2013;<lpage>22</lpage>, <year>2019</year>.</mixed-citation></ref>
<ref id="ref-3"><label>[3]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>J.</given-names> <surname>Hauser</surname></string-name>, <string-name><given-names>S.</given-names> <surname>Sastry</surname></string-name> and <string-name><given-names>P.</given-names> <surname>Kokotovic</surname></string-name></person-group>, &#x201C;<article-title>Nonlinear control via approximate input-output linearization: The ball and beam example</article-title>,&#x201D; <source>IEEE Transactions on Automatic Control</source>, vol. <volume>37</volume>, no. <issue>3</issue>, pp. <fpage>392</fpage>&#x2013;<lpage>398</lpage>, <year>2002</year>.</mixed-citation></ref>
<ref id="ref-4"><label>[4]</label><mixed-citation publication-type="conf-proc"><person-group person-group-type="author"><string-name><given-names>D. N.</given-names> <surname>Kouya</surname></string-name> and <string-name><given-names>F. A.</given-names> <surname>Okou</surname></string-name></person-group>, &#x201C;<article-title>A new adaptive state feedback controller for the ball and beam system</article-title>,&#x201D; in <conf-name>Canadian Conf. on Electrical and Computer Engineering (CCECE)</conf-name>, <conf-loc>Niagara Falls, Canada</conf-loc>, pp. <fpage>247</fpage>&#x2013;<lpage>252</lpage>, <year>2011</year>.</mixed-citation></ref>
<ref id="ref-5"><label>[5]</label><mixed-citation publication-type="conf-proc"><person-group person-group-type="author"><string-name><given-names>D.</given-names> <surname>Martinez</surname></string-name> and <string-name><given-names>F.</given-names> <surname>Ruiz</surname></string-name></person-group>, &#x201C;<article-title>Nonlinear model predictive control for a ball and beam</article-title>,&#x201D; in <conf-name>4th Colombian Workshop on Circuits and Systems (CWCAS)</conf-name>, <conf-loc>Barranquilla, Colombia</conf-loc>, pp. <fpage>1</fpage>&#x2013;<lpage>5</lpage>, <year>2012</year>.</mixed-citation></ref>
<ref id="ref-6"><label>[6]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>T.</given-names> <surname>Wu</surname></string-name> and <string-name><given-names>Y. T.</given-names> <surname>Juang</surname></string-name></person-group>, &#x201C;<article-title>Adaptive fuzzy sliding-mode controller of uncertain nonlinear systems</article-title>,&#x201D; <source>ISA Transactions</source>, vol. <volume>47</volume>, no. <issue>3</issue>, pp. <fpage>279</fpage>&#x2013;<lpage>285</lpage>, <year>2008</year>.</mixed-citation></ref>
<ref id="ref-7"><label>[7]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>A.</given-names> <surname>Nejadfard</surname></string-name>, <string-name><given-names>M. J.</given-names> <surname>Yazdanpanah</surname></string-name> and <string-name><given-names>I.</given-names> <surname>Hassanzadeh</surname></string-name></person-group>, &#x201C;<article-title>Riction compensation of double inverted pendulum on a cart using locally linear neuro-fuzzy model</article-title>,&#x201D; <source>Neural Computing and Applications</source>, vol. <volume>22</volume>, no. <issue>2</issue>, pp. <fpage>337</fpage>&#x2013;<lpage>347</lpage>, <year>2013</year>.</mixed-citation></ref>
<ref id="ref-8"><label>[8]</label><mixed-citation publication-type="conf-proc"><person-group person-group-type="author"><string-name><given-names>W.</given-names> <surname>Luhao</surname></string-name> and <string-name><given-names>S.</given-names> <surname>Zhanshi</surname></string-name></person-group>, &#x201C;<article-title>Lqr-fuzzy control for double inverted pendulum</article-title>,&#x201D; in <conf-name>Int. Conf. on Digital Manufacturing &#x0026; Automation</conf-name>, <conf-loc>Changcha, China</conf-loc>, pp. <fpage>900</fpage>&#x2013;<lpage>903</lpage>, <year>2010</year>.</mixed-citation></ref>
<ref id="ref-9"><label>[9]</label><mixed-citation publication-type="conf-proc"><person-group person-group-type="author"><string-name><given-names>P.</given-names> <surname>Feng</surname></string-name>, <string-name><given-names>L.</given-names> <surname>Lu</surname></string-name> and <string-name><given-names>X.</given-names> <surname>Dingyu</surname></string-name></person-group>, &#x201C;<article-title>Double inverted pendulum sliding mode variation structure control based on fractional order exponential approach law</article-title>,&#x201D; in <conf-name>Proc. CCDC</conf-name>, <conf-loc>Guiyang, China</conf-loc>, pp. <fpage>5077</fpage>&#x2013;<lpage>5080</lpage>, <year>2013</year>.</mixed-citation></ref>
<ref id="ref-10"><label>[10]</label><mixed-citation publication-type="conf-proc"><person-group person-group-type="author"><string-name><given-names>R. M.</given-names> <surname>Nagarale</surname></string-name> and <string-name><given-names>B. M.</given-names> <surname>Patre</surname></string-name></person-group>, &#x201C;<article-title>Decoupled neural fuzzy sliding mode control of nonlinear systems</article-title>,&#x201D; in <conf-name>Proc. FUZZ-IEEE</conf-name>, <conf-loc>Hyderabad, India</conf-loc>, pp. <fpage>1</fpage>&#x2013;<lpage>7</lpage>, <year>2013</year>.</mixed-citation></ref>
<ref id="ref-11"><label>[11]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>Z.</given-names> <surname>Song</surname></string-name> and <string-name><given-names>K.</given-names> <surname>Sun</surname></string-name></person-group>, &#x201C;<article-title>Adaptive backstepping sliding mode control with fuzzy monitoring strategy for a kind of mechanical system</article-title>,&#x201D; <source>ISA Transactions</source>, vol. <volume>53</volume>, no. <issue>1</issue>, pp. <fpage>125</fpage>&#x2013;<lpage>133</lpage>, <year>2014</year>.</mixed-citation></ref>
<ref id="ref-12"><label>[12]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>M.</given-names> <surname>Mehdi</surname></string-name>, <string-name><given-names>M.</given-names> <surname>Rehan</surname></string-name> and <string-name><given-names>F. M.</given-names> <surname>Malik</surname></string-name></person-group>, &#x201C;<article-title>A novel anti-windup framework for cascade control systems: An application to underactuated mechanical systems</article-title>,&#x201D; <source>ISA Transactions</source>, vol. <volume>3</volume>, no. <issue>1</issue>, pp. <fpage>802</fpage>&#x2013;<lpage>815</lpage>, <year>2014</year>.</mixed-citation></ref>
<ref id="ref-13"><label>[13]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>I. M.</given-names> <surname>Mehedi</surname></string-name>, <string-name><given-names>U. M.</given-names> <surname>Al-Saggaf</surname></string-name>, <string-name><given-names>R.</given-names> <surname>Mansouri</surname></string-name> and <string-name><given-names>M.</given-names> <surname>Bettayeb</surname></string-name></person-group>, &#x201C;<article-title>Stabilization of a double inverted rotary pendulum through fractional order integral control scheme</article-title>,&#x201D; <source>International Journal of Advanced Robotic Systems</source>, vol. <volume>16</volume>, no. <issue>4</issue>, pp. <fpage>1</fpage>&#x2013;<lpage>9</lpage>, <year>2019</year>.</mixed-citation></ref>
<ref id="ref-14"><label>[14]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>I. M.</given-names> <surname>Mehedi</surname></string-name></person-group>, &#x201C;<article-title>State feedback based fractional order control scheme for linear servo cart system</article-title>,&#x201D; <source>Journal of Vibroengineering</source>, vol. <volume>20</volume>, no. <issue>1</issue>, pp. <fpage>782</fpage>&#x2013;<lpage>792</lpage>, <year>2018</year>.</mixed-citation></ref>
<ref id="ref-15"><label>[15]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>U. M.</given-names> <surname>Al-Saggaf</surname></string-name>, <string-name><given-names>I. M.</given-names> <surname>Mehedi</surname></string-name>, <string-name><given-names>R.</given-names> <surname>Mansouri</surname></string-name> and <string-name><given-names>M.</given-names> <surname>Bettayeb</surname></string-name></person-group>, &#x201C;<article-title>Rotary flexible joint control by fractional order controllers</article-title>,&#x201D; <source>International Journal of Control Automation and Systems</source>, vol. <volume>15</volume>, no. <issue>6</issue>, pp. <fpage>2561</fpage>&#x2013;<lpage>2569</lpage>, <year>2017</year>.</mixed-citation></ref>
<ref id="ref-16"><label>[16]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>U. M.</given-names> <surname>Al-Saggaf</surname></string-name>, <string-name><given-names>I. M.</given-names> <surname>Mehedi</surname></string-name>, <string-name><given-names>R.</given-names> <surname>Mansouri</surname></string-name> and <string-name><given-names>M.</given-names> <surname>Bettayeb</surname></string-name></person-group>, &#x201C;<article-title>State feedback with fractional integral control design based on the bode&#x0027;s ideal transfer function</article-title>,&#x201D; <source>International Journal of Systems Science</source>, vol. <volume>47</volume>, no. <issue>1</issue>, pp. <fpage>149</fpage>&#x2013;<lpage>161</lpage>, <year>2016</year>.</mixed-citation></ref>
<ref id="ref-17"><label>[17]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>M.</given-names> <surname>Bettayeb</surname></string-name>, <string-name><given-names>C.</given-names> <surname>Boussalem</surname></string-name> and <string-name><given-names>R.</given-names> <surname>Mansouri</surname></string-name></person-group>, &#x201C;<article-title>Stabilization of an inverted pendulum-cart system by fractional pi-state feedback</article-title>,&#x201D; <source>ISA Transactions</source>, vol. <volume>53</volume>, no. <issue>2</issue>, pp. <fpage>508</fpage>&#x2013;<lpage>516</lpage>, <year>2014</year>.</mixed-citation></ref>
<ref id="ref-18"><label>[18]</label><mixed-citation publication-type="conf-proc"><person-group person-group-type="author"><string-name><given-names>A.</given-names> <surname>Young</surname></string-name>, <string-name><given-names>C.</given-names> <surname>Cao</surname></string-name> and <string-name><given-names>N.</given-names> <surname>Hovakimyan</surname></string-name></person-group>, &#x201C;<article-title>Control of a nonaffine double-pendulum system via dynamic inversion and time-scale separation</article-title>,&#x201D; in <conf-name>Proc. ACC</conf-name>, <conf-loc>Minneapolis, MN, USA</conf-loc>, pp. <fpage>6</fpage>&#x2013;<lpage>12</lpage>, <year>2006</year>.</mixed-citation></ref>
<ref id="ref-19"><label>[19]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>A. I.</given-names> <surname>Iskanderani</surname></string-name> and <string-name><given-names>I. M.</given-names> <surname>Mehedi</surname></string-name></person-group>, &#x201C;<article-title>Experimental application of robust and converse dynamic control for rotary flexible joint manipulator system</article-title>,&#x201D; <source>Mathematical Problems in Engineering</source>, vol. <volume>2021</volume>, Article ID 8917134, pp. <fpage>1</fpage>&#x2013;<lpage>9</lpage>, <year>2021</year>.</mixed-citation></ref>
<ref id="ref-20"><label>[20]</label><mixed-citation publication-type="book"><person-group person-group-type="author"><string-name><given-names>A.</given-names> <surname>Ben-Israel</surname></string-name> and <string-name><given-names>T. N.</given-names> <surname>Greville</surname></string-name></person-group>, <source>Generalized Inverses: Theory and Applications</source>, <publisher-loc>Piscataway</publisher-loc>, <publisher-loc>NJ</publisher-loc>, <publisher-loc>USA</publisher-loc>, <publisher-name>Springer</publisher-name>, <year>2003</year>.</mixed-citation></ref>
<ref id="ref-21"><label>[21]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>I. M.</given-names> <surname>Mehedi</surname></string-name>, <string-name><given-names>H. S. M.</given-names> <surname>Shah</surname></string-name>, <string-name><given-names>U. M.</given-names> <surname>Al-Saggaf</surname></string-name>, <string-name><given-names>R.</given-names> <surname>Mansouri</surname></string-name> and <string-name><given-names>M.</given-names> <surname>Bettayeb</surname></string-name></person-group>, &#x201C;<article-title>Adaptive fuzzy sliding mode control of a pressure-controlled artificial ventilator</article-title>,&#x201D; <source>Journal of Healthcare Engineering</source>, vol. <volume>2021</volume>, Article ID 1926711, pp. <fpage>1</fpage>&#x2013;<lpage>10</lpage>, <year>2021</year>.</mixed-citation></ref>
<ref id="ref-22"><label>[22]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>I. M.</given-names> <surname>Mehedi</surname></string-name>, <string-name><given-names>H. S. M.</given-names> <surname>Shah</surname></string-name>, <string-name><given-names>U. M.</given-names> <surname>Al-Saggaf</surname></string-name>, <string-name><given-names>R.</given-names> <surname>Mansouri</surname></string-name> and <string-name><given-names>M.</given-names> <surname>Bettayeb</surname></string-name></person-group>, &#x201C;<article-title>Fuzzy PID control for respiratory systems</article-title>,&#x201D; <source>Journal of Healthcare Engineering</source>, vol. <volume>2021</volume>, Article ID 7118711, pp. <fpage>1</fpage>&#x2013;<lpage>6</lpage>, <year>2021</year>.</mixed-citation></ref>
<ref id="ref-23"><label>[23]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>I. M.</given-names> <surname>Mehedi</surname></string-name>, <string-name><given-names>U.</given-names> <surname>Ansari</surname></string-name> and <string-name><given-names>U. M.</given-names> <surname>Al Saggaf UM</surname></string-name></person-group>, &#x201C;<article-title>Three degree of freedom rotary double inverted pendulum stabilization by using robust generalized dynamic inversion control: Design and experiments</article-title>,&#x201D; <source>Journal of Vibration and Control</source>, vol. <volume>26</volume>, no. <issue>23&#x2013;24</issue>, pp. <fpage>2174</fpage>&#x2013;<lpage>2184</lpage>, <year>2020</year>.</mixed-citation></ref>
<ref id="ref-24"><label>[24]</label><mixed-citation publication-type="conf-proc"><person-group person-group-type="author"><string-name><given-names>I. M.</given-names> <surname>Mehedi</surname></string-name>, <string-name><given-names>U. M.</given-names> <surname>Al-Saggaf</surname></string-name> and <string-name><given-names>S. S.</given-names> <surname>Ali</surname></string-name></person-group>, &#x201C;<article-title>Design and simulation of fractional order control scheme for rotary flexible link</article-title>,&#x201D; in <conf-name>Proc. ICIAS</conf-name>, <conf-loc>Kuala Lumpur, Malaysia</conf-loc>, pp. <fpage>1</fpage>&#x2013;<lpage>5</lpage>, <year>2018</year>.</mixed-citation></ref>
<ref id="ref-25"><label>[25]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>I. M.</given-names> <surname>Mehedi</surname></string-name></person-group>, &#x201C;<article-title>Full state-feedback solution for a flywheel based satellite energy and attitude control scheme</article-title>,&#x201D; <source>Journal of Vibroengineering</source>, vol. <volume>19</volume>, no. <issue>5</issue>, pp. <fpage>3522</fpage>&#x2013;<lpage>3532</lpage>, <year>2017</year>.</mixed-citation></ref>
<ref id="ref-26"><label>[26]</label><mixed-citation publication-type="book"><person-group person-group-type="author"><string-name><given-names>J.</given-names> <surname>Apkarian</surname></string-name>, <string-name><given-names>M.</given-names> <surname>Levis</surname></string-name> and <string-name><given-names>H.</given-names> <surname>Gurocak</surname></string-name></person-group>, <source>Quanser Innovative Edutech-Laboratory Guide-SRV02 Laboratory User Manual</source>, <publisher-loc>Canada</publisher-loc>: <publisher-name>Quanser Inc</publisher-name>, <year>2009</year>.</mixed-citation></ref>
<ref id="ref-27"><label>[27]</label><mixed-citation publication-type="book"><person-group person-group-type="author"><string-name><given-names>H.</given-names> <surname>Khalil</surname></string-name></person-group>, <source>Nonlinear Systems</source>, <publisher-loc>NJ</publisher-loc>, <publisher-loc>USA</publisher-loc>, <publisher-loc>Pearson</publisher-loc>: <publisher-name>Prentice-Hall</publisher-name>, <year>1996</year>.</mixed-citation></ref>
</ref-list>
</back>
</article>