<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Energy Res.</journal-id>
<journal-title>Frontiers in Energy Research</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Energy Res.</abbrev-journal-title>
<issn pub-type="epub">2296-598X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">887561</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2022.887561</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Energy Research</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Chaotic Oscillation Control Model of Power System Under Electromechanical Power Disturbance</article-title>
<alt-title alt-title-type="left-running-head">Hu</alt-title>
<alt-title alt-title-type="right-running-head">Chaos Oscillation of Power System</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Hu</surname>
<given-names>Qiankun</given-names>
</name>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1701801/overview"/>
</contrib>
</contrib-group>
<aff>
<institution>School of Intelligent Manufacturing</institution>, <institution>Huanghuai University</institution>, <addr-line>Zhumadian</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1155401/overview">Fadi AL-Turjman</ext-link>, Near East University, Cyprus</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1706264/overview">Auwalu Mubarak</ext-link>, Near East University, Cyprus</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1707583/overview">Min Cai</ext-link>, Xuzhou University of Technology, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1708119/overview">Zhongqiang Zhang</ext-link>, Xuzhou University of Technology, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Qiankun Hu, <email>huqiankun@huanghuai.edu.cn</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Smart Grids, a section of the journal Frontiers in Energy Research</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>14</day>
<month>04</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>10</volume>
<elocation-id>887561</elocation-id>
<history>
<date date-type="received">
<day>01</day>
<month>03</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>22</day>
<month>03</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2022 Hu.</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Hu</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>The influence of electromechanical power on the power system is controlled in order to stabilize the power system. The author establishes a fourth-order power system model with a power disturbance term based on the dissipative property; the possibility of the existence of a system chaotic attractor is analyzed using the Lyapunov exponent spectrum, bifurcation diagram, phase diagram, spectral entropy, etc., and the influence of the power disturbance term on the motion state of the system is studied. It can be seen that under the influence of the disturbance frequency, the system will exhibit sufficient dynamic behavior. The parameters of the power disturbance term are more sensitive to the influence of the system power angle, and when the disturbance amplitude reaches a certain value, the power angle will increase sharply, and eventually the system will become unstable. The experimental results show that when controller parameters <inline-formula id="inf1">
<mml:math id="m1">
<mml:mrow>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>60</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> are selected, then<inline-formula id="inf2">
<mml:math id="m2">
<mml:mi>&#x3c4;</mml:mi>
</mml:math>
</inline-formula> &#x3d; 0.01, <inline-formula id="inf3">
<mml:math id="m3">
<mml:mi>&#x3b7;</mml:mi>
</mml:math>
</inline-formula> &#x3d; 0.3, and <inline-formula id="inf4">
<mml:math id="m4">
<mml:mi>&#x3b5;</mml:mi>
</mml:math>
</inline-formula> &#x3d; 0.001. Moreover, when two groups of different control objectives <italic>r</italic> &#x3d; 1.2 &#x2b; 0.1 sin(<italic>t</italic>) and <italic>r</italic> &#x3d; 1.2 are chosen, it can be clearly seen that the power angle <inline-formula id="inf5">
<mml:math id="m5">
<mml:mi>&#x3b4;</mml:mi>
</mml:math>
</inline-formula> in the system is very unstable before the controller is connected and the fluctuations are large and irregular, whereas after the controller is added for 150&#xa0;s, the power angle <inline-formula id="inf6">
<mml:math id="m6">
<mml:mi>&#x3b4;</mml:mi>
</mml:math>
</inline-formula> tends to become stable, the fluctuation range becomes small and regular, and the system has almost no chattering. In order to reflect the superiority of the author&#x2019;s control method, under the same parameter conditions, the symbolic function <inline-formula id="inf7">
<mml:math id="m7">
<mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>sgn</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is used as a switching function of the controller and the tracking control numerical simulation is carried out for the same control objectives <italic>r</italic> &#x3d; 1.2 &#x2b; 0.1 sin(<italic>t</italic>) and <italic>r</italic> &#x3d; 1.2. The controller that adopts the relay characteristic function <inline-formula id="inf8">
<mml:math id="m8">
<mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mi>s</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mi>&#x3b5;</mml:mi>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula> as a switching function has a better control effect and the system is smooth and stable without chattering.</p>
</abstract>
<kwd-group>
<kwd>electromechanical power</kwd>
<kwd>power system</kwd>
<kwd>chaotic oscillation</kwd>
<kwd>control</kwd>
<kwd>model</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Chaos is the seemingly random motion that occurs in deterministic systems. It is the main direction of nonlinear research and generally exists in all macroscopic and microscopic systems of the universe (<xref ref-type="bibr" rid="B20">Wang et al., 2017</xref>). At present, chaos theory has penetrated mathematics, physics, chemistry, electronics, information science, biology, geology, meteorology, and cosmology, as well as economics and human brain science, and almost all natural and humanistic fields such as music, art, and sports. Therefore, the chaos theory has a high theoretical research value and practical application value (<xref ref-type="bibr" rid="B4">Bi et al., 2016</xref>). The power system is a dynamic system with strong coupling, high nonlinearity, and multiple parameters. Its dynamic behavior exhibits many complex nonlinear electromechanical oscillation phenomena. When the system is running normally, it exhibits periodic oscillation (<xref ref-type="bibr" rid="B7">Hou et al., 2021</xref>). With the rapid development of large-scale power systems characterized by large units and ultra-high voltage grids, there is a potential threat to the safe operation of the system. Various emergencies and uncertain factors cause continuous and irregular oscillations of the system operating parameters, and chaos often occurs in the actual operation of the power system. In severe cases, the power system may become unstable or even collapse, thereby causing large-scale power outages. The control strategies are aimed at precise integer-order mathematical models, and, in the actual power system, are influenced by external factors and the nature of the system itself. There are uncertainties in system parameters and external disturbances; in particular, new energy and distributed generation are increasingly being introduced into the grid, which further increases the uncertainty of the power system (<xref ref-type="bibr" rid="B22">Wang et al., 2019</xref>). These uncertainties easily cause chaotic oscillation of the system and seriously affect the stability of the power system (<xref ref-type="bibr" rid="B25">Zhao and Kamwa, 2020</xref>). On the other hand, integer-order systems are approximate idealizations of fractional-order systems. The fractional-order power systems are mainly based on chaotic synchronization. In particular, there is less research on chaos control under uncertain factors, there are not many researches based on chaos control. <xref ref-type="fig" rid="F1">Figure 1</xref> shows the physical model of power oscillation. Therefore, for integer-order and fractional-order power system models, considering the influence of internal parameters and uncertain factors of the system, it is necessary to study the chaos control strategy. This will provide the scientific basis and reference data for early warning and processing of the power system&#x27;s chaotic oscillation. Relative to the relationship between the state variables of the system and input variables, for a linear system that satisfies the superposition principle, the nonlinear system is a kind of irregular motion and is ubiquitous in nature, with the most important motion behavior being the chaotic motion (<xref ref-type="bibr" rid="B14">Preece and Milanovic, 2016</xref>). The definition of chaos has not yet reached a consensus in the industry and its expression is complex, with disorder and irregularity; it is often manifested as a nonlinear system under certain conditions and exhibits unpredictable phenomena. It is a manifestation of the fusion of variability and immutability, the presence or absence of sequences or rules. Chaotic motion is neither an invisible phenomenon nor a phenomenon that can be seen by the naked eye. It can be observed only with the help of advanced instruments. It is ubiquitous in the motion state of all things in the universe, and this phenomenon can be observed every day in life (<xref ref-type="bibr" rid="B2">Banerjee et al., 2018</xref>; <xref ref-type="bibr" rid="B10">Luo et al., 2020</xref>). For example, from the lit blue smoke that suddenly tumbles irregularly after rising and finally dissipates, to the calm stock market that becomes suddenly chaotic because of a stock change, and also the sputtering of water droplets, whose state trajectory is a chaotic phenomenon. Since the 1960s, there have been more and more studies on chaos. It has developed into a very huge discipline system, radiating to physics, economics, finance, meteorology, biology, sociology and other disciplines. At present, the research of various chaos theories will affect the development of modern discipline system, radiation to physics, economics and finance, meteorology, biology, sociology, and other disciplines, and now, every kind of chaos theory research will affect the further development of this modern disciplinary system.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Physical model of power oscillation.</p>
</caption>
<graphic xlink:href="fenrg-10-887561-g001.tif"/>
</fig>
<p>The author proposes a study on the chaotic oscillation control model of the power system under electromechanical power disturbance and establishes a fourth-order power system model with a power disturbance term based on the dissipative property; the possibility of the existence of a system chaotic attractor is analyzed using the Lyapunov exponent spectrum, bifurcation diagram, phase diagram, spectral entropy, etc., and the influence of the power disturbance term on the motion state of the system is studied. It can be seen that under the influence of the disturbance frequency, the system will exhibit sufficient dynamic behavior. The parameters of the power disturbance term are more sensitive to the influence of the system power angle, and when the disturbance amplitude reaches a certain value, the power angle will increase sharply and eventually the system will become unstable.</p>
</sec>
<sec id="s2">
<title>2 Literature Review</title>
<p>In response to this research problem, <xref ref-type="bibr" rid="B9">Liu et al. (2017)</xref> proposed the scientific concept of chaos, which is an important achievement in nonlinear science in recent years. <xref ref-type="bibr" rid="B23">Zhang et al. (2015)</xref> proposed a method to calculate the Lyapunov index of a system from a time series, moving chaos research from the theoretical stage to the practical application stage. <xref ref-type="bibr" rid="B6">Errouissi et al. (2017)</xref> proposed the OGY control method. Since then, researchers have continued to discover new chaotic systems and hyperchaotic systems, and chaos control has become a new hot spot in chaos research. <xref ref-type="bibr" rid="B3">Besselmann et al. (2016)</xref> established a simple power system model and proved that the system will exhibit chaotic oscillation. <xref ref-type="bibr" rid="B11">Mishra et al. (2021)</xref> studied in detail the case of considering and ignoring the damping winding, a bifurcation phenomenon in a three-node power system. Since then, the research on chaos control of power system has sprung up. As the world's largest energy producer and consumer, China's goal is the safe and stable operation of power system, which further promotes the in-depth study of power system control theory. Based on a simple interconnected power system model, <xref ref-type="bibr" rid="B17">Shen et al. (2020)</xref> analyzed the chaos generation mechanism of a three-parameter system. <xref ref-type="bibr" rid="B24">Zhang et al. (2020)</xref> focused on the analysis of its chaotic phenomenon for the classical generator rocking equation. <xref ref-type="bibr" rid="B12">Mousakazemi (2019)</xref> studied the chaotic phenomenon of a power system under different instability modes, and, for the first time, two kinds of instability phenomena were found in the ruptured state of the chaotic limit cycle: system angular instability as well as simultaneous voltage and angular instability. <xref ref-type="bibr" rid="B15">Qin et al. (2018)</xref> reviewed several bifurcation phenomena in power systems and introduced the application of the chaos theory in short-term load forecasting of power systems. <xref ref-type="bibr" rid="B21">Wang et al. (2020)</xref> adopted a nonlinear feedback control method, an active feedback control method, and control research on fractional-order power system chaos. <xref ref-type="bibr" rid="B13">Overbye and Klump (2015)</xref> used the fuzzy sliding mode control method to suppress the chaos of the power system and reduce chattering to a certain extent. The authors proposed a study on the chaotic oscillation control model of a power system under electromechanical power disturbance, established a fourth-order power system model with a power disturbance term based on the dissipative property, analyzed the possibility of the existence of a system chaotic attractor using the Lyapunov exponent spectrum, bifurcation diagram, phase diagram, spectral entropy, etc., and studied the influence of the power disturbance term on the motion state of the system, where it can be seen that under the influence of the disturbance frequency, the system will exhibit sufficient dynamic behavior. The parameters of the power disturbance term are more sensitive to the influence of the system power angle, and when the disturbance amplitude reaches a certain value, the power angle will increase sharply and eventually the system will become unstable. The experimental results show that when the controller parameters <inline-formula id="inf9">
<mml:math id="m9">
<mml:mrow>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>60</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> are selected, then <inline-formula id="inf10">
<mml:math id="m10">
<mml:mi>&#x3c4;</mml:mi>
</mml:math>
</inline-formula> &#x3d; 0.01, <inline-formula id="inf11">
<mml:math id="m11">
<mml:mi>&#x3b7;</mml:mi>
</mml:math>
</inline-formula> &#x3d; 0.3, and <inline-formula id="inf12">
<mml:math id="m12">
<mml:mi>&#x3b5;</mml:mi>
</mml:math>
</inline-formula> &#x3d; 0.001. Moreover, when two groups of different control objectives r &#x3d; 1.2 &#x2b; 0.1 sin(<italic>t</italic>) and <italic>r</italic> &#x3d; 1.2 are selected, it can be clearly seen that in a system before the controller is connected, the power angle <inline-formula id="inf13">
<mml:math id="m13">
<mml:mi>&#x3b4;</mml:mi>
</mml:math>
</inline-formula> is very unstable, the fluctuation range is large, and there is no regularity. However, after the controller is added for 150&#xa0;s, the power angle <inline-formula id="inf14">
<mml:math id="m14">
<mml:mi>&#x3b4;</mml:mi>
</mml:math>
</inline-formula> tends to become stable, the fluctuation range becomes small and regular, and the system has almost no chattering. The system is smooth and stable without chattering.</p>
</sec>
<sec id="s3">
<title>3 Methods</title>
<sec id="s3-1">
<title>3.1 System Modeling</title>
<p>The output voltage <inline-formula id="inf15">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> of the excitation controller is constrained by the controller, according to the magnitude of the input voltage <inline-formula id="inf16">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> that varies within the finite interval <inline-formula id="inf17">
<mml:math id="m17">
<mml:mrow>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, and the system equation can be described as<disp-formula id="e1">
<mml:math id="m18">
<mml:mrow>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>&#x3b4;</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mi>&#x3c9;</mml:mi>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>&#x3c9;</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>d</mml:mi>
<mml:mi>&#x3c9;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mi>E</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:msup>
<mml:mi>x</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msup>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>E</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:msup>
<mml:mi>X</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>X</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:msup>
<mml:mi>E</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:msup>
<mml:mi>x</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mi>d</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:msup>
<mml:mi>X</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>X</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:msup>
<mml:mi>T</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>E</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>In the above equation, <inline-formula id="inf18">
<mml:math id="m19">
<mml:mi>&#x3b4;</mml:mi>
</mml:math>
</inline-formula> is the power angle of the generator, <inline-formula id="inf19">
<mml:math id="m20">
<mml:mi>&#x3c9;</mml:mi>
</mml:math>
</inline-formula> is the angular frequency of the generator, <inline-formula id="inf20">
<mml:math id="m21">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the fundamental frequency of the synchronous motor, <italic>M</italic> is the inertia of the generator rotor, <italic>d</italic> is the damping factor, <inline-formula id="inf21">
<mml:math id="m22">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the generator of the input power, <inline-formula id="inf22">
<mml:math id="m23">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>G</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is delivered by the electromagnetic power of the generator, <italic>x</italic> is the reactance of the transmission line, <inline-formula id="inf23">
<mml:math id="m24">
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the reactance of the generator, <inline-formula id="inf24">
<mml:math id="m25">
<mml:mrow>
<mml:msub>
<mml:msup>
<mml:mi>x</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mi>d</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the transient reactance of the generator, <inline-formula id="inf25">
<mml:math id="m26">
<mml:mrow>
<mml:msub>
<mml:msup>
<mml:mi>T</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the time constant of the generator stator winding, <inline-formula id="inf26">
<mml:math id="m27">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the excitation voltage, <inline-formula id="inf27">
<mml:math id="m28">
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the infinite voltage, and <italic>V</italic> is the terminal voltage of the generator and can be expressed as<disp-formula id="e2">
<mml:math id="m29">
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:msup>
<mml:mi>x</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mi>d</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:msqrt>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:msup>
<mml:mi>x</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>x</mml:mi>
<mml:msup>
<mml:mi>E</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:msup>
<mml:mi>E</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:msqrt>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>In the power system, the excitation link of the terminal amplitude limiting method is generally used to protect the equipment of the system, and the output voltage <inline-formula id="inf28">
<mml:math id="m30">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> of the excitation controller is expressed as<disp-formula id="e3">
<mml:math id="m31">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3e;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>When the input voltage <inline-formula id="inf29">
<mml:math id="m32">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is in the interval <inline-formula id="inf30">
<mml:math id="m33">
<mml:mrow>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, the output voltage <inline-formula id="inf31">
<mml:math id="m34">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> (<xref ref-type="bibr" rid="B16">Ramanathan and Vittal, 2015</xref>), and when <inline-formula id="inf32">
<mml:math id="m35">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> exceeds the output upper limit (or lower limit), the output voltage <inline-formula id="inf33">
<mml:math id="m36">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> remains at the upper limit <inline-formula id="inf34">
<mml:math id="m37">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> (or lower limit <inline-formula id="inf35">
<mml:math id="m38">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>) and no longer changes with <inline-formula id="inf36">
<mml:math id="m39">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, until <inline-formula id="inf37">
<mml:math id="m40">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> returns to the range of the interval <inline-formula id="inf38">
<mml:math id="m41">
<mml:mrow>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> (<xref ref-type="bibr" rid="B8">Liu et al., 2016</xref>).</p>
<p>The parameter values are shown in <xref ref-type="table" rid="T1">Table 1</xref>. The dissipation characteristics of Equation (1) can be calculated as shown in Equation (4):<disp-formula id="e4">
<mml:math id="m42">
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mo>&#x2207;</mml:mo>
<mml:mi>V</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>&#x3b4;</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>&#x3c9;</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>&#x3c9;</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>E</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msup>
<mml:mi>E</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>E</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:msup>
<mml:mi>x</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:msup>
<mml:mi>T</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mn>10</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mn>7</mml:mn>
<mml:mn>6</mml:mn>
</mml:mfrac>
<mml:mo>.</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
<label>(4)</label>
</disp-formula>Corresponding to time <italic>t, V</italic> is always shrinking during the movement, as shown in Equation (5). This shows the possibility of the existence of chaotic attractors in the system.<disp-formula id="e5">
<mml:math id="m43">
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>V</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mn>10</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mn>7</mml:mn>
<mml:mn>6</mml:mn>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>At the same time, by substituting the system parameter values in <xref ref-type="table" rid="T1">Table 1</xref>, the algebraic equation shown in <xref ref-type="disp-formula" rid="e6">Equation 6</xref> can be obtained (<xref ref-type="bibr" rid="B19">Stankovic et al., 2015</xref>).<disp-formula id="e6">
<mml:math id="m44">
<mml:mrow>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>120</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
<mml:mi>&#x3c9;</mml:mi>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>d</mml:mi>
<mml:mi>&#x3c9;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:msup>
<mml:mi>E</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mrow>
<mml:mn>0.4</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.5</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>10</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.5</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>0.4</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.5</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:msup>
<mml:mi>E</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>0.4</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>0.4</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.5</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>10</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1.05</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>System parameter values.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Parameter name</th>
<th align="center">Numerical value</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Synchronous motor fundamental frequency (<inline-formula id="inf39">
<mml:math id="m45">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>)</td>
<td align="center">60</td>
</tr>
<tr>
<td align="left">Transmission line reactance (X)</td>
<td align="center">0.4</td>
</tr>
<tr>
<td align="left">Generator dynamic reactance (<inline-formula id="inf40">
<mml:math id="m46">
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>)</td>
<td align="center">1</td>
</tr>
<tr>
<td align="left">Excitation limiter voltage reference (<inline-formula id="inf41">
<mml:math id="m47">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>)</td>
<td align="center">2</td>
</tr>
<tr>
<td align="left">Equivalent moment of inertia (<italic>h</italic>)</td>
<td align="center">4.9</td>
</tr>
<tr>
<td align="left">Generator transient reactance (<inline-formula id="inf42">
<mml:math id="m48">
<mml:mrow>
<mml:msub>
<mml:msup>
<mml:mi>x</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mi>d</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>)</td>
<td align="center">0.4</td>
</tr>
<tr>
<td align="left">Bus voltage reference (<inline-formula id="inf43">
<mml:math id="m49">
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>)</td>
<td align="center">1.05</td>
</tr>
<tr>
<td align="left">Excitation limiter voltage minimum (<inline-formula id="inf44">
<mml:math id="m50">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>)</td>
<td align="center">0</td>
</tr>
<tr>
<td align="left">Infinite bus voltage (<inline-formula id="inf45">
<mml:math id="m51">
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>)</td>
<td align="center">1</td>
</tr>
<tr>
<td align="left">Generator stator time constant (<inline-formula id="inf46">
<mml:math id="m52">
<mml:mrow>
<mml:msub>
<mml:msup>
<mml:mi>T</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>)</td>
<td align="center">10</td>
</tr>
<tr>
<td align="left">Excitation time constant (<inline-formula id="inf47">
<mml:math id="m53">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>)</td>
<td align="center">1</td>
</tr>
<tr>
<td align="left">Excitation limiter voltage max. (<inline-formula id="inf48">
<mml:math id="m54">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>)</td>
<td align="center">5</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Among these, the output voltage <inline-formula id="inf49">
<mml:math id="m55">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is expressed as<disp-formula id="e7">
<mml:math id="m56">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mn>5</mml:mn>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3e;</mml:mo>
<mml:mn>5</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mn>5</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<p>The voltage <italic>V</italic> is expressed as<disp-formula id="e8">
<mml:math id="m57">
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mn>0.5</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.4</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x22c5;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>0.4</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.5</mml:mn>
<mml:msup>
<mml:mi>E</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>0.5</mml:mn>
<mml:msup>
<mml:mi>E</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:msqrt>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>By setting a set of system parameters <italic>d</italic> &#x3d; 0.5, <inline-formula id="inf50">
<mml:math id="m58">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1.3</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, and<inline-formula id="inf51">
<mml:math id="m59">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>150</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, and substituting them into Equation (6) for iterative calculation, a set of equilibrium points can be obtained: <inline-formula id="inf52">
<mml:math id="m60">
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>&#x3b4;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>&#x3c9;</mml:mi>
<mml:mo>,</mml:mo>
<mml:msup>
<mml:mi>E</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1.0409</mml:mn>
<mml:mo>,</mml:mo>
<mml:mn>0,1.3559</mml:mn>
<mml:mo>,</mml:mo>
<mml:mn>1.9229</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>.<disp-formula id="e9">
<mml:math id="m61">
<mml:mrow>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>&#x3b4;</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mi>&#x3c9;</mml:mi>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>&#x3c9;</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>d</mml:mi>
<mml:mi>&#x3c9;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mi>E</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:msup>
<mml:mi>x</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>cos</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>cos</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msup>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>E</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:msup>
<mml:mi>x</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:msup>
<mml:mi>E</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:msup>
<mml:mi>x</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mi>d</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:msup>
<mml:mi>x</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:msup>
<mml:mi>T</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>E</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mi>A</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>To set the simulated environment parameters, system parameters <italic>D</italic> &#x3d; 2, <inline-formula id="inf53">
<mml:math id="m62">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 1.2, <inline-formula id="inf54">
<mml:math id="m63">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 190 and the disturbance term parameters <inline-formula id="inf55">
<mml:math id="m64">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 0.2, <inline-formula id="inf56">
<mml:math id="m65">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 0.2, <inline-formula id="inf57">
<mml:math id="m66">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 0, <inline-formula id="inf58">
<mml:math id="m67">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 0 are selected; at this time, the system behaves in a state of random chaotic oscillation. In order to suppress the chaotic oscillation phenomenon, a controller to observe the control effect after the system runs for 150&#xa0;s is added. Here, the controller parameters <inline-formula id="inf59">
<mml:math id="m68">
<mml:mrow>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>60</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf60">
<mml:math id="m69">
<mml:mi>&#x3c4;</mml:mi>
</mml:math>
</inline-formula> &#x3d; 0.01, <inline-formula id="inf61">
<mml:math id="m70">
<mml:mi>&#x3b7;</mml:mi>
</mml:math>
</inline-formula> &#x3d; 0.3, and <inline-formula id="inf62">
<mml:math id="m71">
<mml:mi>&#x3b5;</mml:mi>
</mml:math>
</inline-formula> &#x3d; 0.001 are selected. Then, two groups of different control objectives r &#x3d; 1.2 &#x2b; 0.1 sin(<italic>t</italic>) and <italic>r</italic> &#x3d; 1.2 are chosen to control the effect.</p>
</sec>
<sec id="s3-2">
<title>3.2 Influence of the Power Disturbance Term on the System</title>
<p>We will substitute <inline-formula id="inf63">
<mml:math id="m72">
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>&#x3b4;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>&#x3c9;</mml:mi>
<mml:mo>,</mml:mo>
<mml:msup>
<mml:mi>E</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1.0409</mml:mn>
<mml:mo>,</mml:mo>
<mml:mn>0,1.3559</mml:mn>
<mml:mo>,</mml:mo>
<mml:mn>1.9229</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> as the initial value of the system into <xref ref-type="disp-formula" rid="e9">Equation 9</xref> for iterative operation (<xref ref-type="bibr" rid="B18">Silva et al., 2017</xref>). Using the bifurcation diagram, Lyapunov exponent spectrum, and phase diagram, the effects of these four parameters on the evolution process of the state attractor of the system are shown respectively (<xref ref-type="bibr" rid="B18">Silva et al., 2017</xref>).</p>
<sec id="s3-2-1">
<title>3.2.1 Disturbance Amplitude <inline-formula id="inf64">
<mml:math id="m73">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</title>
<p>The presence of the electromagnetic disturbance term may cause the system to appear ultra-high voltage, which affects the stable operations of the system, therefore it is necessary to study the electromagnetic disturbance term. When <inline-formula id="inf65">
<mml:math id="m74">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mi>H</mml:mi>
<mml:mi>z</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, and <inline-formula id="inf66">
<mml:math id="m75">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0.2</mml:mn>
<mml:mi>H</mml:mi>
<mml:mi>z</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> are selected and when <inline-formula id="inf67">
<mml:math id="m76">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>0,0.2255</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, the <inline-formula id="inf68">
<mml:math id="m77">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> bifurcation diagram and Lyapunov exponent spectrum of the system are obtained (<xref ref-type="bibr" rid="B1">Alam et al., 2020</xref>).</p>
<p>When <inline-formula id="inf69">
<mml:math id="m78">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>0,0.1256</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, the two largest Lyapunov exponents (LEs) of the system within the parameter range are kept close to 0&#xa0;at the same time, that is, the distribution of the LE value of the system in this range is (0,0,-,-,-) by the LE value and the bifurcation diagram. The specific value of LE at the boundary point <inline-formula id="inf70">
<mml:math id="m79">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0.1256</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> of the quasi-periodic state and the chaotic state of the system is (0,0,&#x2212;0.0944,&#x2212;0.2502,&#x2212;1.1088). When <inline-formula id="inf71">
<mml:math id="m80">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> continues to increase beyond 0.1256, the maximum LE of the system will also increase rapidly and becomes greater than 0, and the system motion realizes the transition process from quasi-periodic to chaotic. In the range of <inline-formula id="inf72">
<mml:math id="m81">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>0.1256</mml:mn>
<mml:mo>,</mml:mo>
<mml:mn>0.2255</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, the maximum LE of the system is greater than 0; combined with the bifurcation diagram, it can be seen that the system is running in a chaotic motion state (<xref ref-type="bibr" rid="B5">Dong et al., 2018</xref>). <xref ref-type="disp-formula" rid="e9">Equation 9</xref> cannot be in a stable state for a long time and causes voltage collapse under excessive disturbance. <xref ref-type="table" rid="T2">Table 2</xref> gives the LE and operating states of the system for some <inline-formula id="inf73">
<mml:math id="m82">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> values.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>LE and system status at different <inline-formula id="inf74">
<mml:math id="m83">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> values.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">
<inline-formula id="inf75">
<mml:math id="m84">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
<th align="center">LE</th>
<th align="center">System operating status</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">0.1237</td>
<td align="center">(0,0,&#x2212;0.0646,&#x2212;0.2723,&#x2212;1.0364)</td>
<td align="left">Quasi-period</td>
</tr>
<tr>
<td align="left">0.2250</td>
<td align="center">(0.07620,&#x2212;0.0217,&#x2212;0.2882,&#x2212;0.570)</td>
<td align="left">Chaos</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3-2-2">
<title>3.2.2 Electromagnetic Disturbance Frequency <inline-formula id="inf76">
<mml:math id="m85">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</title>
<p>Frequency, another important parameter of power disturbance, will also directly affect the motion state of the power system. Without considering the load disturbance and choosing in the range of <inline-formula id="inf77">
<mml:math id="m86">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x2208;(0.7780,0.8205) <inline-formula id="inf78">
<mml:math id="m87">
<mml:mo>&#x222a;</mml:mo>
</mml:math>
</inline-formula> (0.8740,0.9315), the maximum LE is greater than 0, indicating that the system is moving in a chaotic state. When <inline-formula id="inf79">
<mml:math id="m88">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x2208;(0,0.7780) <inline-formula id="inf80">
<mml:math id="m89">
<mml:mo>&#x222a;</mml:mo>
</mml:math>
</inline-formula> (0.9315,2), both the largest and second-largest LEs of the system tend to 0, and it can be seen from the bifurcation diagram that the system is in a quasi-periodic motion state at this time. The change in the motion state of the system from periodic to chaotic corresponds to the above analysis. The corresponding system LEs and operating states of some <inline-formula id="inf81">
<mml:math id="m90">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> values are shown in <xref ref-type="table" rid="T3">Table 3</xref>.</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>LE and system states at different <inline-formula id="inf82">
<mml:math id="m91">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> values.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">
<inline-formula id="inf83">
<mml:math id="m92">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
<th align="center">LE</th>
<th align="center">System operating status</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">0.4200</td>
<td align="center">(0,&#x2212;0.0469,&#x2212;0.0491,&#x2212;0.2336,&#x2212;1.0371)</td>
<td align="left">Cycle</td>
</tr>
<tr>
<td align="left">0.8000</td>
<td align="center">(0.0838,0,&#x2212;0.0329,&#x2212;0.2612.&#x2212;1.1564)</td>
<td align="left">Chaos</td>
</tr>
<tr>
<td align="left">0.8500</td>
<td align="center">(0,&#x2212;0.0517,&#x2212;0.1310,&#x2212;0.1389,&#x2212;1.0451)</td>
<td align="left">Cycle</td>
</tr>
<tr>
<td align="left">0.9000</td>
<td align="center">(0.0773,0,&#x2212;0.0196&#x2212;0.2138,&#x2212;1.2107)</td>
<td align="left">Chaos</td>
</tr>
<tr>
<td align="left">1.2680</td>
<td align="center">(0,&#x2212;0.0266,&#x2212;0.1972,&#x2212;0.2050,&#x2212;0.9380)</td>
<td align="left">Cycle</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3-2-3">
<title>3.2.3 Load Disturbance Amplitude <inline-formula id="inf84">
<mml:math id="m93">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</title>
<p>The presence of load disturbances may cause problems such as harmonics and voltage fluctuations that impair power quality and threaten the stability of the power system; therefore, it is necessary to study the load disturbance term. Without considering the electromagnetic disturbance, <inline-formula id="inf85">
<mml:math id="m94">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 0, <inline-formula id="inf86">
<mml:math id="m95">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 0, and <inline-formula id="inf87">
<mml:math id="m96">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 0.2 are selected. The bifurcation graph and Lyapunov exponent spectrum of the system are obtained with respect to <inline-formula id="inf88">
<mml:math id="m97">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> when <inline-formula id="inf89">
<mml:math id="m98">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x2208;(0,0.2555). In the range of <inline-formula id="inf90">
<mml:math id="m99">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x2208;(0, 0.1370), both the maximum LE and the second-largest LE of the system are close to 0 and the system is now in a quasi-periodic motion state. When <inline-formula id="inf91">
<mml:math id="m100">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x2208; (0.1370, 0.2055), the maximum LE of the system remains 0, the next largest LE is less than or equal to 0, and the system runs in a periodic state and the number of cycles is large; however, when <inline-formula id="inf92">
<mml:math id="m101">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is changed in the above, no trace of periodic motion is found. When <inline-formula id="inf93">
<mml:math id="m102">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x2208;(0.2055,0.2555), the maximum LE of the system within the parameter range is obviously positive; at this time, the motion state of the system is in a typical chaotic state. And when the <inline-formula id="inf94">
<mml:math id="m103">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> value exceeds 0.2555, the system will collapse due to the phenomenon of power angle divergence due to excessive disturbance. <xref ref-type="table" rid="T4">Table 4</xref> shows the LE and operating states of the system for some <inline-formula id="inf95">
<mml:math id="m104">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> values.</p>
<table-wrap id="T4" position="float">
<label>TABLE 4</label>
<caption>
<p>LE and system status at different <inline-formula id="inf96">
<mml:math id="m105">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> values.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">
<inline-formula id="inf97">
<mml:math id="m106">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
<th align="center">LE</th>
<th align="center">State of motion</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">0.0903</td>
<td align="center">(0,0,&#x2212;0.1003,&#x2212;0.2182,&#x2212;1.0739)</td>
<td align="left">Quasi-period</td>
</tr>
<tr>
<td align="left">0.1703</td>
<td align="center">(0,&#x2212;0.0628,&#x2212;0.1287,&#x2212;0.2227,&#x2212;0.9660)</td>
<td align="left">Cycle</td>
</tr>
<tr>
<td align="left">0.2180</td>
<td align="center">(02167,0,&#x2212;0.1242,&#x2212;0.2461,&#x2212;1.0505)</td>
<td align="left">Chaos</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3-2-4">
<title>3.2.4 Load Disturbance Frequency <inline-formula id="inf98">
<mml:math id="m107">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</title>
<p>Similar to the above analysis method, without considering electromagnetic disturbance, <inline-formula id="inf99">
<mml:math id="m108">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 0.02, <inline-formula id="inf100">
<mml:math id="m109">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; O, and <inline-formula id="inf101">
<mml:math id="m110">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 0 are chosen, and the bifurcation graph and Lyapunov exponent spectrum of the system are obtained with respect to <inline-formula id="inf102">
<mml:math id="m111">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> when <inline-formula id="inf103">
<mml:math id="m112">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x2208;(0,2).</p>
<p>When compared to the situation of changing A, changing B brings about similar changes to the system and there are only a few minor differences. When <inline-formula id="inf104">
<mml:math id="m113">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is in the range of (0.7615, 0.8130) <inline-formula id="inf105">
<mml:math id="m114">
<mml:mo>&#x222a;</mml:mo>
</mml:math>
</inline-formula> (0.8855, 0.9575), the maximum LE of the system is greater than 0, which shows that the motion state of the system is chaotic. When <inline-formula id="inf106">
<mml:math id="m115">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x2208;(0.7135,0.7615) <inline-formula id="inf107">
<mml:math id="m116">
<mml:mo>&#x222a;</mml:mo>
</mml:math>
</inline-formula> (0.8130,0.8855) <inline-formula id="inf108">
<mml:math id="m117">
<mml:mo>&#x222a;</mml:mo>
</mml:math>
</inline-formula> (1.2565,1.2775) <inline-formula id="inf109">
<mml:math id="m118">
<mml:mo>&#x222a;</mml:mo>
</mml:math>
</inline-formula> (1.6820,1.7040), the maximum LE of the system is 0 and the rest of the LEs are obviously less than 0; when <inline-formula id="inf110">
<mml:math id="m119">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x2208;(0,0.7135) <inline-formula id="inf111">
<mml:math id="m120">
<mml:mo>&#x222a;</mml:mo>
</mml:math>
</inline-formula> (0.9575,1.2565) <inline-formula id="inf112">
<mml:math id="m121">
<mml:mo>&#x222a;</mml:mo>
</mml:math>
</inline-formula> (1.2775,1.6820) <inline-formula id="inf113">
<mml:math id="m122">
<mml:mo>&#x222a;</mml:mo>
</mml:math>
</inline-formula> (1.7040,2), the two largest LE values in the system approach 0 at the same time, which means that the system operates in a quasi-periodic state within this parameter range. From the phase diagrams of the system at <inline-formula id="inf114">
<mml:math id="m123">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 0.75 and <inline-formula id="inf115">
<mml:math id="m124">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 0.90, it can be clearly seen that the motion state of the system changes from periodic to chaotic, corresponding to the above analysis. The corresponding system LEs and operating states of some <inline-formula id="inf116">
<mml:math id="m125">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> values are shown in <xref ref-type="table" rid="T5">Table 5</xref>.</p>
<table-wrap id="T5" position="float">
<label>TABLE 5</label>
<caption>
<p>LE and system status at different <inline-formula id="inf117">
<mml:math id="m126">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> values.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">
<inline-formula id="inf118">
<mml:math id="m127">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
<th align="center">LE</th>
<th align="center">Operating status</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">0.5000</td>
<td align="center">(0,0,&#x2212;0.0686,&#x2212;0.4534,&#x2212;0.8472)</td>
<td align="left">Quasi-period</td>
</tr>
<tr>
<td align="left">0.7500</td>
<td align="center">(0,&#x2212;0.0620,&#x2212;0.0724,&#x2212;0.1066,&#x2212;1.1257)</td>
<td align="left">Cycle</td>
</tr>
<tr>
<td align="left">0.8000</td>
<td align="center">(0.0637,0,&#x2212;0.0274,&#x2212;0.2505,&#x2212;1.1525)</td>
<td align="left">Chaos</td>
</tr>
<tr>
<td align="left">0.8500</td>
<td align="center">(0,&#x2212;0.1122,&#x2212;0.1144,&#x2212;0.1686,&#x2212;0.9716)</td>
<td align="left">Cycle</td>
</tr>
<tr>
<td align="left">0.9000</td>
<td align="center">(0.1004,0,&#x2212;0.0281,&#x2212;0.2503,&#x2212;1.1887)</td>
<td align="left">Chaos</td>
</tr>
<tr>
<td align="left">1.2700</td>
<td align="center">(0,&#x2212;0.1305,&#x2212;0.1636,&#x2212;0.1679,&#x2212;0.9047)</td>
<td align="left">Cycle</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</sec>
</sec>
<sec id="s4">
<title>4 Results and Analysis</title>
<p>The system parameters D &#x3d; 2, <inline-formula id="inf119">
<mml:math id="m128">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 1.2, and <inline-formula id="inf120">
<mml:math id="m129">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 190 and the disturbance term parameters <inline-formula id="inf121">
<mml:math id="m130">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 0.2, <inline-formula id="inf122">
<mml:math id="m131">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 0.2, and <inline-formula id="inf123">
<mml:math id="m132">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 0, <inline-formula id="inf124">
<mml:math id="m133">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 0 are selected; at this time, the system behaves in a state of random chaotic oscillation. In order to suppress the chaotic oscillation phenomenon, a controller to observe the control effect after the system runs for 150&#xa0;s is added. Here, the controller parameter <inline-formula id="inf125">
<mml:math id="m134">
<mml:mrow>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>60</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, with <inline-formula id="inf126">
<mml:math id="m135">
<mml:mi>&#x3c4;</mml:mi>
</mml:math>
</inline-formula> &#x3d; 0.01, <inline-formula id="inf127">
<mml:math id="m136">
<mml:mi>&#x3b7;</mml:mi>
</mml:math>
</inline-formula> &#x3d; 0.3, and <inline-formula id="inf128">
<mml:math id="m137">
<mml:mi>&#x3b5;</mml:mi>
</mml:math>
</inline-formula> &#x3d; 0.001, is selected. Choosing two groups of the different control objectives r &#x3d; 1.2 &#x2b; 0.1 sin(t) and r &#x3d; 1.2 controls the effect. The results are shown in <xref ref-type="fig" rid="F2">Figures 2</xref>, <xref ref-type="fig" rid="F3">3</xref>. It can be clearly seen that the power angle <inline-formula id="inf129">
<mml:math id="m138">
<mml:mi>&#x3b4;</mml:mi>
</mml:math>
</inline-formula> in the system is very unstable before the controller is connected, the fluctuations are large and irregular, and after adding the controller for 150&#xa0;s, the power angle <inline-formula id="inf130">
<mml:math id="m139">
<mml:mi>&#x3b4;</mml:mi>
</mml:math>
</inline-formula> tends to become stable, the fluctuation range becomes small and regular, and the system has almost no chattering.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Sliding mode control with target trajectory <italic>r</italic> &#x3d; 1.2 &#x2b; 0.1 sin(<italic>t</italic>).</p>
</caption>
<graphic xlink:href="fenrg-10-887561-g002.tif"/>
</fig>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Sliding mode control with target trajectory <italic>r</italic> &#x3d; 1.2.</p>
</caption>
<graphic xlink:href="fenrg-10-887561-g003.tif"/>
</fig>
<p>The symbolic function <inline-formula id="inf131">
<mml:math id="m140">
<mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>sgn</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> with the tracking control numerical simulation is carried out for the same control objectives r &#x3d; 1.2 &#x2b; 0.1 sin(<italic>t</italic>) and <italic>r</italic> &#x3d; 1.2, and the simulation results are shown in <xref ref-type="fig" rid="F4">Figure 4</xref> and <xref ref-type="fig" rid="F5">Figure 5</xref>, respectively. Obviously, the controller using the relay characteristic function <inline-formula id="inf132">
<mml:math id="m141">
<mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mi>s</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mi>&#x3b5;</mml:mi>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula> as the switching function has a better controlling effect and the system is smooth and stable without chattering. While using the conventional sign function as the switching function, although a good control effect can be obtained, the system chattering is obvious and the chattering frequency is high; at the same time, due to the severe chattering phenomenon, the numerical simulation takes a long time and the amount of data is large.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Sliding mode control of target trajectory <italic>r</italic> &#x3d; 1.2 &#x2b; 0.1 sin(<italic>t</italic>) using a sign function.</p>
</caption>
<graphic xlink:href="fenrg-10-887561-g004.tif"/>
</fig>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Sliding mode control with a sign function for the target trajectory <italic>r</italic> &#x3d; 1.2.</p>
</caption>
<graphic xlink:href="fenrg-10-887561-g005.tif"/>
</fig>
</sec>
<sec id="s5">
<title>5 Conclusion</title>
<p>This article proposes a chaotic oscillation control model of the power system under electromechanical power interference. The author proposes a study on the chaotic oscillation control model of a power system under electromechanical power disturbance. By establishing a fourth-order power system model with a power disturbance term, the Lyapunov exponent, bifurcation diagram, and spectral entropy are analyzed and the influence of the power disturbance term on the motion state of the power system is discussed. The simulation results show that the controller can quickly and smoothly suppress the chaotic oscillation of the system, and at the same time, it can effectively avoid the chattering problem and has strong robustness. Chaos control is one of the frontiers of nonlinear research, and the systematization and control methods of the chaos theory still need to be improved. Second, the majority of researchers should consider the engineering realization of chaos control. A large number of simulation experiments have theoretically proved the effectiveness of the control strategy but there is still a certain distance from being widely used in engineering applications. Therefore, researchers have to conduct in-depth research and engineering experiments and use interdisciplinary research methods for referencing to form practical research results.</p>
</sec>
</body>
<back>
<sec id="s6">
<title>Data Availability Statement</title>
<p>The original contributions presented in the study are included in the article/Supplementary Material, and further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s7">
<title>Author Contributions</title>
<p>The author confirms being the sole contributor to this work and has approved it for publication.</p>
</sec>
<sec sec-type="COI-statement" id="s8">
<title>Conflict of Interest</title>
<p>The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s9">
<title>Publisher&#x2019;s Note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Alam</surname>
<given-names>M. R.</given-names>
</name>
<name>
<surname>Bai</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Saha</surname>
<given-names>T. K.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Classification and Visualization of Power Quality Disturbance-Events Using Space Vector Ellipse in Complex Plane</article-title>. <source>IEEE Trans. Power Deliv.</source> <volume>36</volume> (<issue>99</issue>), <fpage>1</fpage>. <pub-id pub-id-type="doi">10.1109/TPWRD.2020.3008003</pub-id> </citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Banerjee</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Chaudhuri</surname>
<given-names>N. R.</given-names>
</name>
<name>
<surname>Kavasseri</surname>
<given-names>R. G.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>A Novel Explicit Disturbance Model-Based Robust Damping of Inter-area Oscillations through Mtdc Grids Embedded in Ac Systems</article-title>. <source>IEEE Trans. Power Deliv.</source> <volume>33</volume>, <fpage>1</fpage>. <pub-id pub-id-type="doi">10.1109/tpwrd.2018.2799170</pub-id> </citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Besselmann</surname>
<given-names>T. J.</given-names>
</name>
<name>
<surname>Alm&#xe9;r</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Ferreau</surname>
<given-names>H. J.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Model Predictive Control of Load-Commutated Inverter-Fed Synchronous Machines</article-title>. <source>IEEE Trans. Power Electro.</source> <volume>31</volume> (<issue>10</issue>), <fpage>7384</fpage>&#x2013;<lpage>7393</lpage>. <pub-id pub-id-type="doi">10.1109/TPEL.2015.2511095</pub-id> </citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bi</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Qin</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Martin</surname>
<given-names>K. E.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Disturbance Propagation Mechanism Based on the Electromechanical Wave Theory</article-title>. <source>IET Generation, Transm. Distribution</source> <volume>10</volume> (<issue>12</issue>), <fpage>2891</fpage>&#x2013;<lpage>2898</lpage>. <pub-id pub-id-type="doi">10.1049/iet-gtd.2015.1268</pub-id> </citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dong</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Shi</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Xiao</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Xia</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>A Finite Set Model Predictive Control Method for Quasi-Z Source Inverter-Permanent Magnet Synchronous Motor Drive System</article-title>. <source>IET Electric Power Appl.</source> <volume>13</volume> (<issue>3</issue>), <fpage>302</fpage>&#x2013;<lpage>309</lpage>. </citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Errouissi</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Al&#x2010;Durra</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Muyeen</surname>
<given-names>S. M.</given-names>
</name>
<name>
<surname>Leng</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Continuous&#x2010;time Model Predictive Control of a Permanent Magnet Synchronous Motor Drive with Disturbance Decoupling</article-title>. <source>Iet Electric Power Appl.</source> <volume>11</volume> (<issue>5</issue>), <fpage>697</fpage>&#x2013;<lpage>706</lpage>. <pub-id pub-id-type="doi">10.1049/iet-epa.2016.0499</pub-id> </citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hou</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Ke</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>A Flexible Constant Power Generation Scheme for Photovoltaic System by Error-Based Active Disturbance Rejection Control and Perturb &#x26; Observe</article-title>. <source>Energy</source> <volume>237</volume> (<issue>9</issue>), <fpage>121646</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2021.121646</pub-id> </citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Loh</surname>
<given-names>P. C.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Droop Control with Improved Disturbance Adaption for a Pv System with Two Power Conversion Stages</article-title>. <source>IEEE Trans. Ind. Electron.</source> <volume>63</volume> (<issue>10</issue>), <fpage>6073</fpage>&#x2013;<lpage>6085</lpage>. <pub-id pub-id-type="doi">10.1109/tie.2016.2580525</pub-id> </citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Miao</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Fan</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Han</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Markovian Switching Model and Non&#x2010;linear DC Modulation Control of AC/DC Power System</article-title>. <source>IET Generation, Transm. Distribution</source> <volume>11</volume> (<issue>10</issue>), <fpage>2654</fpage>&#x2013;<lpage>2663</lpage>. <pub-id pub-id-type="doi">10.1049/iet-gtd.2016.1862</pub-id> </citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Luo</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Lewis</surname>
<given-names>F. L.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Ouakad</surname>
<given-names>H. M.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Accelerated Adaptive Fuzzy Optimal Control of Three Coupled Fractional-Order Chaotic Electromechanical Transducers</article-title>. <source>IEEE Trans. Fuzzy Syst.</source> <volume>29</volume> (<issue>99</issue>), <fpage>1</fpage>. <pub-id pub-id-type="doi">10.1109/TFUZZ.2020.2984998</pub-id> </citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mishra</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Singh</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Panigrahi</surname>
<given-names>B. K.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Adaptive Current Control for a Bidirectional Interleaved Ev Charger with Disturbance Rejection</article-title>. <source>IEEE Trans. Industry Appl.</source> <volume>57</volume> (<issue>99</issue>), <fpage>1</fpage>. <pub-id pub-id-type="doi">10.1109/tia.2021.3074612</pub-id> </citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mousakazemi</surname>
<given-names>S. M. H.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Control of a Pwr Nuclear Reactor Core Power Using Scheduled Pid Controller with ga, Based on Two-point Kinetics Model and Adaptive Disturbance Rejection System</article-title>. <source>Ann. Nucl. Energ.</source> <volume>129</volume> (<issue>JUL</issue>), <fpage>487</fpage>&#x2013;<lpage>502</lpage>. <pub-id pub-id-type="doi">10.1016/j.anucene.2019.02.019</pub-id> </citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Overbye</surname>
<given-names>T. J.</given-names>
</name>
<name>
<surname>Klump</surname>
<given-names>R. P.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Determination of Emergency Power System Voltage Control Actions</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>13</volume> (<issue>1</issue>), <fpage>205</fpage>&#x2013;<lpage>210</lpage>. <pub-id pub-id-type="doi">10.1109/59.651637</pub-id> </citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Preece</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Milanovic</surname>
<given-names>J. V.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Efficient Estimation of the Probability of Small-Disturbance Instability of Large Uncertain Power Systems</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>31</volume> (<issue>2</issue>), <fpage>1063</fpage>&#x2013;<lpage>1072</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrs.2015.2417204</pub-id> </citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qin</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Bi</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Martin</surname>
<given-names>K. E.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Model Approach of Power Networks Based on Unsymmetrical Inertia Distribution for Disturbance Propagation Study</article-title>. <source>IET Generation, Transm. Distribution</source> <volume>12</volume> (<issue>9</issue>), <fpage>2055</fpage>&#x2013;<lpage>2064</lpage>. <pub-id pub-id-type="doi">10.1049/iet-gtd.2017.1243</pub-id> </citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ramanathan</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Vittal</surname>
<given-names>V.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Small-disturbance Angle Stability Enhancement through Direct Load Control Part Ii-Numerical Simulations and Results</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>21</volume> (<issue>2</issue>), <fpage>782</fpage>&#x2013;<lpage>790</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRS.2006.873025</pub-id> </citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shen</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>W. D.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Jin</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Frequency Response Model and its Closed-form Solution of Two-Machine Equivalent Power System</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>36</volume> (<issue>99</issue>), <fpage>1</fpage>. <pub-id pub-id-type="doi">10.1109/TPWRS.2020.3037695</pub-id> </citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Silva</surname>
<given-names>L. R. M.</given-names>
</name>
<name>
<surname>Kapisch</surname>
<given-names>E. B.</given-names>
</name>
<name>
<surname>Martins</surname>
<given-names>C. H. N.</given-names>
</name>
<name>
<surname>Filho</surname>
<given-names>L. M. A.</given-names>
</name>
<name>
<surname>Cerqueira</surname>
<given-names>A. S.</given-names>
</name>
<name>
<surname>Duque</surname>
<given-names>C. A.</given-names>
</name>
<etal/>
</person-group> (<year>2017</year>). <article-title>Gapless Power-Quality Disturbance Recorder</article-title>. <source>IEEE Trans. Power Deliv.</source> <volume>32</volume> (<issue>2</issue>), <fpage>862</fpage>&#x2013;<lpage>871</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrd.2016.2557280</pub-id> </citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Stankovic</surname>
<given-names>A. M.</given-names>
</name>
<name>
<surname>Dukic</surname>
<given-names>S. D.</given-names>
</name>
<name>
<surname>Saric</surname>
<given-names>A. T.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Approximate Bisimulation-Based Reduction of Power System Dynamic Models</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>30</volume> (<issue>3</issue>), <fpage>1252</fpage>&#x2013;<lpage>1260</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrs.2014.2342504</pub-id> </citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Characteristics of Electromechanical Disturbance Propagation in Non&#x2010;uniform Power Systems</article-title>. <source>IET Generation, Transm. Distribution</source> <volume>11</volume> (<issue>8</issue>), <fpage>1919</fpage>&#x2013;<lpage>1925</lpage>. <pub-id pub-id-type="doi">10.1049/iet-gtd.2016.1126</pub-id> </citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Rodriguez</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Fpga-based Continuous Control Set Model Predictive Current Control for Pmsm System Using Multistep Error Tracking Technique</article-title>. <source>IEEE Trans. Power Electron.</source> <volume>35</volume> (<issue>12</issue>), <fpage>13455</fpage>&#x2013;<lpage>13464</lpage>. <pub-id pub-id-type="doi">10.1109/tpel.2020.2984336</pub-id> </citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Liang</surname>
<given-names>H. R.</given-names>
</name>
<name>
<surname>Prokhorov</surname>
<given-names>A. V.</given-names>
</name>
<name>
<surname>Mokhlis</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Huat</surname>
<given-names>C. K.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Modal Control Design of Damping Controllers for Thyristor-Controlled Series Capacitor to Stabilize Common-Mode Torsional Oscillations of a Series-Capacitor Compensated Power System</article-title>. <source>IEEE Trans. Industry Appl.</source> <volume>55</volume>, <fpage>1</fpage>. <pub-id pub-id-type="doi">10.1109/tia.2019.2892343</pub-id> </citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Dong</surname>
<given-names>Z. Y.</given-names>
</name>
<name>
<surname>Wong</surname>
<given-names>K. P.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Post&#x2010;disturbance Transient Stability Assessment of Power Systems by a Self&#x2010;adaptive Intelligent System</article-title>. <source>IET Generation, Transm. Distribution</source> <volume>9</volume> (<issue>3</issue>), <fpage>296</fpage>&#x2013;<lpage>305</lpage>. <pub-id pub-id-type="doi">10.1049/iet-gtd.2014.0264</pub-id> </citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Jiao.</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Model-free Predictive Current Control of Dfig Based on an Extended State Observer under Unbalanced and Distorted Grid</article-title>. <source>IEEE Trans. Power Electron.</source> <volume>35</volume> (<issue>8</issue>), <fpage>8130</fpage>&#x2013;<lpage>8139</lpage>. <pub-id pub-id-type="doi">10.1109/tpel.2020.2967172</pub-id> </citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Kamwa</surname>
<given-names>I.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Guest Editorial: Next Generation of Synchrophasor-Based Power System Monitoring, Operation and Control</article-title>. <source>IET Generation Transm. Distribution</source> <volume>14</volume> (<issue>19</issue>), <fpage>3943</fpage>&#x2013;<lpage>3944</lpage>. <pub-id pub-id-type="doi">10.1049/iet-gtd.2020.1315</pub-id> </citation>
</ref>
</ref-list>
</back>
</article>