Research article Topical Sections

SVC control enhancement applying self-learning fuzzy algorithm for islanded microgrid

  • Maintaining voltage stability, within acceptable levels, for islanded Microgrids (MGs) is a challenge due to limited exchange power between generation and loads. This paper proposes an algorithm to enhance the dynamic performance of islanded MGs in presence of load disturbance using Static VAR Compensator (SVC) with Fuzzy Model Reference Learning Controller (FMRLC). The proposed algorithm compensates MG nonlinearity via fuzzy membership functions and inference mechanism imbedded in both controller and inverse model. Hence, MG keeps the desired performance as required at any operating condition. Furthermore, the self-learning capability of the proposed control algorithm compensates for grid parameter’s variation even with inadequate information about load dynamics. A reference model was designed to reject bus voltage disturbance with achievable performance by the proposed fuzzy controller. Three simulations scenarios have been presented to investigate effectiveness of proposed control algorithm in improving steady-state and transient performance of islanded MGs. The first scenario conducted without SVC, second conducted with SVC using PID controller and third conducted using FMRLC algorithm. A comparison for results shows ability of proposed control algorithm to enhance disturbance rejection due to learning process.

    Citation: Ahmed Eldessouky, Hossam Gabbar. SVC control enhancement applying self-learning fuzzy algorithm for islanded microgrid[J]. AIMS Energy, 2016, 4(2): 363-378. doi: 10.3934/energy.2016.2.363

    Related Papers:

    [1] Xiaoqiang Dai, Chao Yang, Shaobin Huang, Tao Yu, Yuanran Zhu . Finite time blow-up for a wave equation with dynamic boundary condition at critical and high energy levels in control systems. Electronic Research Archive, 2020, 28(1): 91-102. doi: 10.3934/era.2020006
    [2] Min Li . Long-wavelength limit for the Green–Naghdi equations. Electronic Research Archive, 2022, 30(7): 2700-2718. doi: 10.3934/era.2022138
    [3] Shasha Bian, Yitong Pei, Boling Guo . Numerical simulation of a generalized nonlinear derivative Schrödinger equation. Electronic Research Archive, 2022, 30(8): 3130-3152. doi: 10.3934/era.2022159
    [4] Xuerong Shi, Zuolei Wang, Lizhou Zhuang . Spatiotemporal pattern in a neural network with non-smooth memristor. Electronic Research Archive, 2022, 30(2): 715-731. doi: 10.3934/era.2022038
    [5] Yang Yang, Yun-Rui Yang, Xin-Jun Jiao . Traveling waves for a nonlocal dispersal SIR model equipped delay and generalized incidence. Electronic Research Archive, 2020, 28(1): 1-13. doi: 10.3934/era.2020001
    [6] Hyungyeong Jung, Sunghwan Moon . Reconstruction of the initial function from the solution of the fractional wave equation measured in two geometric settings. Electronic Research Archive, 2022, 30(12): 4436-4446. doi: 10.3934/era.2022225
    [7] José Luis Díaz Palencia, Abraham Otero . Oscillatory solutions and smoothing of a higher-order p-Laplacian operator. Electronic Research Archive, 2022, 30(9): 3527-3547. doi: 10.3934/era.2022180
    [8] Yang Liu, Wenke Li . A family of potential wells for a wave equation. Electronic Research Archive, 2020, 28(2): 807-820. doi: 10.3934/era.2020041
    [9] Hami Gündoğdu . RETRACTED ARTICLE: Impact of damping coefficients on nonlinear wave dynamics in shallow water with dual damping mechanisms. Electronic Research Archive, 2025, 33(4): 2567-2576. doi: 10.3934/era.2025114
    [10] Shan Zhao, Jun Feng . Qualitative analysis and traveling wave solutions of the (3+1)-dimensional Kadomtsev-Petviashvili-Sawada-Kotera-Ramani equation with conformable fractional derivative. Electronic Research Archive, 2025, 33(6): 3716-3732. doi: 10.3934/era.2025165
  • Maintaining voltage stability, within acceptable levels, for islanded Microgrids (MGs) is a challenge due to limited exchange power between generation and loads. This paper proposes an algorithm to enhance the dynamic performance of islanded MGs in presence of load disturbance using Static VAR Compensator (SVC) with Fuzzy Model Reference Learning Controller (FMRLC). The proposed algorithm compensates MG nonlinearity via fuzzy membership functions and inference mechanism imbedded in both controller and inverse model. Hence, MG keeps the desired performance as required at any operating condition. Furthermore, the self-learning capability of the proposed control algorithm compensates for grid parameter’s variation even with inadequate information about load dynamics. A reference model was designed to reject bus voltage disturbance with achievable performance by the proposed fuzzy controller. Three simulations scenarios have been presented to investigate effectiveness of proposed control algorithm in improving steady-state and transient performance of islanded MGs. The first scenario conducted without SVC, second conducted with SVC using PID controller and third conducted using FMRLC algorithm. A comparison for results shows ability of proposed control algorithm to enhance disturbance rejection due to learning process.


    First, I take this opportunity to express my great admiration toward Professor Goong Chen by dedicating this paper to him on the occasion of his 70th birthday.

    In this paper, we mainly consider the problem of vibrations governed by 1D wave equation $ w_{tt}-c^2w_{xx} = 0 $, where $ c $ denotes wave propagation speed, associated with a generalized boundary condition. Recall the definition of chaos for this kind of system, which is firstly introduced in [1], as below:

    Definition 1.1. Consider an initial-boundary problem $ (S) $ governed by 1D wave equation $ w_{tt}-c^2w_{xx} = 0 $ defined on a segment $ I $, where $ c > 0 $ denotes the propagation speed of the wave. The system is said to be chaotic if there exists a large class of initial data $ (w_0, w_1) $ such that

    (i) $ |w_x|+|w_t| $ is uniformly bounded,

    (ii) $ V(t)\overset{def}{ = }V_I(w_x(\cdot, t))+V_I(w_t(\cdot, t)) < +\infty $ for all $ t\geq0 $,

    (iii) $ \liminf\limits_{t\to+\infty}\frac{\ln V(t)}{t} > 0 $.

    Remark 1.1. When chaos occurs, the function $ x\mapsto \left(|w_x|+|w_t|\right)(x, t) $ is uniformly bounded, whereas the length of the curve $ \{(x, \; \left(|w_x|+|w_t|\right)(x, t)), \; x\in[0, 1]\} $ grows exponentially w.r.t. time $ t $. Therefore, the system $ (S) $ must undergo extremely complex oscillations as time $ t $ increasing.

    Remark 1.2. There are lots of work about chaos studies, see e.g., [2,3,4] and references therein. As we know, there is no a common mathematical definition for chaos, which is actually a challenge to give. However, Li-Yorke chaos is probably one of the most popular and acceptable notions of chaos. We will study the relationship between Li-Yorke chaos and chaos in the sense of Definition 1.1 in our future work. Particularly, if chaos happens and the solution $ (w_x, w_t) $ can be represented by two interval maps, denoted by $ K_1 $ and $ K_2 $, respectively, then $ K_1 $ and $ K_2 $ have positive entropy, which implies $ K_1 $ and $ K_2 $ are chaotic in the sense of Li-Yorke.

    Let us take a classical model to introduce the research background and our motivation, as below:

    $ {wttwxx=0,  x(0,1), t>0,wx(0,t)=ηwt(0,t), η1,  t>0,wx(1,t)=αwt(1,t)βw3t(1,t),  0<α<1, β0, t>0,w(x,0)=w0(x),wt(x,0)=w1(x),0x1,
    $
    (1.1)

    where $ \alpha, \beta $, and $ \eta $ are given constants. If $ \eta = 1 $, the system (1.1) is not well-posed. Thus throughout this paper we assume $ \eta\not = 1 $. The wave equation itself is linear and represents the infinite-dimensional harmonic oscillator. Let

    $ E(t) = \frac{1}{2}\int_0^1|\nabla w(x, t)|^2dx = \frac{1}{2}\int_0^1\left[w_x^2(x, t)+w_t^2(x, t)\right]dx, $

    be the energy function of this system. And assume that (1.1) admits a $ C^2 $ solution, then the boundary conditions show

    $ \frac{d}{dt}E(t) = \eta w_t^2(0, t)+w_t^2(1, t)[\alpha-\beta w_t^2(1, t)]. $

    The right-handed side boundary condition (at $ x = 1 $) is nonlinear when $ \beta\ne 0 $, which is usually called a van der Pol type boundary condition (see, e.g., [1,5,6,7,8,9]). The left-handed side boundary condition (at $ x = 0 $) is linear, where $ \eta > 0 $ indicates that energy is being injected into the system at $ x = 0 $. Thus if $ \eta > 0 $, the system (1.1) has a self-excited mechanism that supplies energy to the system itself, which induces irregular vibrations [1,6]. In particular, when $ \eta = 0 $, the free end at $ x = 0 $ has no effect to the energy, which is also called an energy-conserving boundary condition.

    The existence and uniqueness of the classical solution of (1.1) can be found in [5,6]. Furthermore, the system (1.1) has a smooth solution $ w\in \mathcal{C}^2 $ if the initial data satisfy

    $ w0C20([0,1]),w1C10([0,1]),
    $
    (1.2)

    where

    $ Ck0([0,1])={fCk([0,1]) | f(i)(0)=f(i)(1)=0, 0ik }, kN,
    $
    (1.3)

    see Theorem 6.1 in [6]. The weak solution as well as its numerical approximation are discussed in [9].

    The PDE system (1.1) has received considerable attention since it exhibits many interesting and complicated dynamical phenomena, such as limit cycles and chaotic behavior of $ (w_t, w_x) $ when the parameters $ \alpha, \beta $ and $ \eta $ assume certain values [1,6]. Different from dynamics of a system of ODEs, this is a simple and useful infinite-dimensional model for the study of spatiotemporal behaviors as time evolutes. For instance, the propagation of acoustic waves in a pipe satisfies the linear wave equation: $ {w}_{tt}-{w}_{xx} = 0. $ As we know, the solution of 1D wave equation describes a superposition of two traveling wave with arbitrary profiles, one propagating with unit speed to the left, the other with unit speed to the right. The boundary conditions appeared in (1.1) can create irregularly acoustical vibrations ([1,6,7]). This type of vibrations, for example, can be generated by noise signals radiated from underwater vehicles, and there are intensive research for the properties of acoustical vibrations in current literature (see e.g., [10] and references therein). Hence the study of this type of vibration is not only important but also may lead to a better understanding of the dynamics of acoustic systems.

    In this paper, we mainly consider the oscillation problems described by the following models:

    $ {wtt(x,t)c2wxx(x,t)=0,  x(0,1), t>0,wx(0,t)=ηwt(0,t), ηc1,  t>0,wx(1,t)=h(wt(1,t)), t>0,w0(x)=w(x,0),w1(x)=wt(x,0),0x1,
    $
    (1.4)
    $ {ztt(x,t)c2zxx(x,t)=0,  x(0,1), t>0,z(0,t)=0,  t>0,zx(1,t)=h(zt(1,t)), t>0,z0(x)=z(x,0),z1(x)=zt(x,0),0x1,
    $
    (1.5)

    where the function $ h\in\mathcal{C}^0(\mathbb{R}) $ satisfies

    (A1) $ \phi: t\mapsto \frac{1}{2}\left(h(t)-\frac{1}{c}t\right) $ strictly monotonically decreases on $ \mathbb{R} $,

    (A2) $ \phi(\mathbb{R}) = \mathbb{R} $.

    It is clear that $ \phi^{-1}:\mathbb{R}\to \mathbb{R} $ is well-defined and strictly decreases on $ \mathbb{R} $. In fact, the model (1.4) is a generalized case of the van del Pol type boundary condition, $ h(x) = \alpha x-\beta x^3 $, say. Denote the total energy of system (1.4) as

    $ E(t) = \frac{1}{2c^2}\int_0^1\left[c^2w_x^2(x, t)+w_t^2(x, t)\right]dx, $

    then the boundary conditions show

    $ \frac{d}{dt}E(t) = \eta w_t^2(0, t)+w_t^2(1, t)\cdot\frac{h(w_t(1,t))}{w_t(1,t)}. $

    Hence $ \eta > 0 $ can cause the energy of system to increase. Moreover, the system (1.1) has a generalized self-excited mechanism if we assume that, roughly speaking, $ \frac{h(y)}{y} > 0 $ if $ |y| $ is small and $ \frac{h(y)}{y} < 0 $ if $ |y| $ is large. Thus if $ \eta > 0 $, the system (1.4) has a self-excited mechanism that supplies energy to the system itself, which could induce irregular oscillations. The interesting part is that the system (1.5) may have chaos even though there is no energy supplier at the boundary.

    We can treat the initial-boundary problems (1.4) and (1.5) by using wave propagation method. Let's take (1.4) as a start. It is well-known that $ w $, a solution of 1D wave equation, has the following form:

    $ w(x,t)=L(x+ct)+R(xct),
    $
    (1.6)

    where $ L $ and $ R $ are two $ C^1 $ functions. Then the gradient of $ w $ can be represented as follows:

    $ wt(x,t)=cL(x+ct)cR(xct),wx(x,t)=L(x+ct)+R(xct),
    $
    (1.7)

    for $ x\in[0, 1], t\geq0 $. Introduce two new variables

    $ u(x,t)=L(x+ct),v(x,t)=R(xct),
    $
    (1.8)

    which are called the Riemann invariants. It is evident that $ u $ and $ v $ keeps constants along the lines $ x+ct = const. $ and $ x-ct = const. $, respectively, which are referred to as characteristics.

    When $ t > 0 $ and $ x = 0 $, from the boundary condition at the left end $ x = 0 $,

    $ v(0,t)=cwx(0,t)wt(0,t)2c=cηwt(0,t)wt(0,t)2c=1+cη2cwt(0,t)=1+cη2(u(0,t)v(0,t)),
    $

    that is

    $ v(0,t)=cη+1cη1u(0,t)def=γ(η)u(0,t).
    $
    (1.9)

    Without confusion, we also take $ \gamma(\eta) $ as a linear map.

    When $ t > 0 $ and $ x = 1 $, from the boundary condition at the right end $ x = 1 $, it follows:

    $ h\left(c(u(1,t)-v(1,t))\right)-\frac{1}{c}\cdot c\left(u(1,t)-v(1,t)\right)-2v(1,t) = 0, $

    that is

    $ v(1,t) = \phi\left(c(u(1,t)-v(1,t))\right), $

    which determines a reflection relationship between $ u(1, t) $ and $ v(1, t) $ as follow

    $ u(1,t)=1cϕ1(v(1,t))+v(1,t)def=φ(v(1,t)),
    $
    (1.10)

    where

    $ φ(x)=1cϕ1(x)+x.
    $
    (1.11)

    When $ t = 0 $, from the initial conditions,

    $ u0(x)def=u(x,0)=cw0(x)+w1(x)2c,v0(x)def=v(x,0)=cw0(x)w1(x)2c,
    $
    (1.12)

    which are referred to as the initial date of $ (u, v) $.

    For $ x\in[0, 1] $ and $ \tau\geq0 $, it follows from the boundary reflections (1.9) and (1.10) that

    $ u(x,τ+2c)=L(x+cτ+2)=u(1,1cx+τ+1c)=φ(v(1,1cx+τ+1c))=φ(R(xcτ))=φ(v(0,1cx+τ))=φ(γ(η)u(x,τ)),
    $

    and

    $ v\left(x,\tau+\frac{2}{c}\right) = \gamma(\eta)\circ \varphi\left(v(x,\tau)\right); $

    inductively,

    $ u(x,τ+2nc)=(φγ(η))n(u(x,τ)),v(x,τ+2nc)=(γ(η)φ)n(v(x,τ)),
    $
    (1.13)

    where the superscript $ n\in\mathbb{N} $ denotes the $ n $th iteration of a function. Analogously, for $ x\in[0, 1] $ and $ \tau\in[0, \frac{2}{c}) $,

    $ u(x,τ)={u0(x+cτ),0cτ1x,φ(v0(2xcτ)),1x<cτ2x,φγ(η)(u0(x+cτ2)),2x<cτ<2,
    $
    (1.14)

    and

    $ v(x,τ)={v0(xcτ),0cτx,γ(η)(u0(cτx)),x<cτx+1,γ(η)φ(v0(2+xcτ)),x+1<cτ<2.
    $
    (1.15)

    Equations (1.13)–(1.15) show that the system (1.4) is solvable and the dynamics of the solution $ (w, w_x, w_t) $ to the equation (1.4) can be uniquely determined by the initial data and the following two functions:

    $ ψη=γ(η)φ,g=φγ(η).
    $
    (1.16)

    Note that $ g = \gamma^{-1}(\eta)\circ \psi_\eta\circ\gamma(\eta) $, that is to say there is topological conjugacy between $ \psi_\eta $ and $ g $. Therefore one only needs to consider one of them, say $ \psi_\eta $.

    The paper will be organized as follows. In the next two sections, we will present the necessary and sufficient conditions to cause the onset of chaos, respectively. Section 4 shows, as a special case, in which the wave equation has a fixed end, that chaos can occur with the effects of self-regulations and energy-conservation. In Section 5, there are some applications of the theoretical outcomes. In the last section, it is the numerical simulations.

    In this section, we firstly give a necessary condition for the onset of chaos in the following system:

    $ {wtt(x,t)c2wxx(x,t)=0,  x(0,1), t>0,wx(0,t)=ηwt(0,t), ηc1,  t>0,wx(1,t)=h(wt(1,t)), t>0,w0(x)=w(x,0),w1(x)=wt(x,0),0x1,
    $
    (2.1)

    where $ h $ satisfies hypotheses $ (A1) $ and $ (A2) $. In addition, assume the function $ t\mapsto \left(h(t)+\frac{1}{c}t\right) $ is piecewise monotone.

    Theorem 2.1. Suppose the system $(2.1)$ is chaotic in the sense of Definition 1.1. Then $ h $ is not even and $ \eta\ne0 $.

    Proof. It is equivalent to prove that there is no chaos in the system (2.1) if $ \eta = 0 $ or $ h $ is an even function.

    Firstly, assume $ h $ is an even function. Let $ \eta\in\mathbb{R} $ and $ \eta\ne c^{-1} $. Recall the function $ \psi_\eta $ given by (1.16):

    $ \psi_\eta = \gamma(\eta)\circ\varphi,\; \varphi(x) = \frac{1}{c}\cdot\phi^{-1}\left(x\right)+x. $

    Introduce a new map $ Q $ from $ \mathbb{R} $ to $ \mathbb{R} $ as follow:

    $ Q(y)=12(h(y)+1cy).
    $
    (2.2)

    For $ x\in\mathbb{R} $, let $ y = \phi^{-1}(x) $, then

    $ ψη(x)=γ(η)(1cy+ϕ(y))=12γ(η)(1cy+h(y))=γ(η)Qϕ1(x).
    $

    Since $ \phi^{-1} $ is strictly decreasing and $ \gamma(\eta)\in\mathbb{R}^* $, $ \psi $ is monotonic if and only if $ Q $ is monotonic. Let $ y_1, \; y_2\in \mathbb{R} $ with $ y_1 < y_2 $. From the hypothesis $ h $ being an even function, it follows

    $ Q(y2)Q(y1)=12(h(y2)1c(y2))12(h(y1)1c(y1))=ϕ(y2)ϕ(y1)>0,
    $

    which implies $ Q $ strictly monotonically increases. Therefore $ \psi $ strictly monotonically increases (decreases) if $ \gamma(\eta) < 0 (\gamma(\eta) > 0) $. It is easily seen that

    $ wt(x,t)=c[γ1(η)ψn(γ(η)u(x,τ))ψn(v(x,τ))],wx(x,t)=γ1(η)ψn(γ(η)u(x,τ))+ψn(v(x,τ)).
    $

    Hence the dynamics of $ (w_x, w_t) $ is simple, and the chaos doesn't occur in the system.

    Next, assume $ \eta = 0 $, in other words, the system is free at the left end. In this case,

    $ \psi_0(x) = -x-\phi^{-1}(x). $

    We will show that there is no period point of $ \psi_0 $ with period 2. Let $ x_0\in\mathbb{R} $ satisfy $ x_0 < \psi_0(x_0) $. Define two functions as follows:

    $ q:x\mapsto-x+\psi_0(x_0)+x_0,\qquad k = \psi_0-q. $

    Since $ \phi^{-1}(\cdot) $ strictly monotonically decreases and

    $ k(x) = -x-\phi^{-1}(x)-q(x) = \phi^{-1}(x_0)-\phi^{-1}(x), $

    $ k $ strictly monotonically increases in $ [x_0, \; \psi_0(x_0)] $. Hence

    $ k(\psi_0(x_0)) = \psi_0^2(x_0)-x_0 > k(x_0) = 0. $

    That implies there are no period points of $ \psi_0 $ with period 2. By Sharkovsky's Theorem, there are no periods of $ \psi_0 $ with period lager than 2. By virtue of the Main Theorem 6 in [12], there is still no chaos in the system when $ \eta = 0 $.

    Next, we consider the following system governed by 1D wave equation with a fixed end:

    $ {ztt(x,t)c2zxx(x,t)=0,  x(0,1), t>0,z(0,t)=0,  t>0,zx(1,t)=h(zt(1,t)), t>0,z0(x)=z(x,0),z1(x)=zt(x,0),0x1,
    $
    (2.3)

    where $ h $ satisfies hypotheses $ (A1) $ and $ (A2) $.

    Theorem 2.2. Suppose the system $(2.3)$ is chaotic in the sense of Definition 1.1. Then$ h $ is neither an even function nor an odd function.

    Proof. Put

    $ \psi_{\infty}(x) = x+\phi^{-1}(x). $

    By the analysis in Section 1, it is clear that $ (z_x, z_t) $ can be represented by iterations of $ \psi_{\infty} $ and initial data. Hence, the dynamics of $ (z_x, z_t) $ is completely determined by $ \psi_{\infty} $.

    If $ h $ is even, $ \psi_{\infty} $ is monotonically monotone. If $ h $ is odd, $ -\psi_{\infty} $ has no periodic points of period larger than 2. Therefore, chaos never occurs in the system 2.3 if $ h $ is either an even function or an odd function.

    Remark 2.1. Chaos can definitely happen in the system 2.3 for a special kind of $ h $. One can find more details about that in the later section.

    In this section, we mainly try to determine some sufficient conditions to ensure the onset of chaotic oscillations in the following system:

    $ {wtt(x,t)c2wxx(x,t)=0,  x(0,1), t>0,wx(0,t)=ηwt(0,t), ηc1,  t>0,wx(1,t)=h(wt(1,t)), t>0,w0(x)=w(x,0),w1(x)=wt(x,0),0x1,
    $
    (3.1)

    where $ h\in\mathcal{C}^0(\mathbb{R}) $ satisfies $ h(0) = 0 $, hypotheses $ (A1) $ and $ (A2) $. When $ w_t(1, t)\equiv0 $, there should be no signals feedback to $ w_x(1, t) $, that's to say the right end should be free. Therefore it is reasonable to let $ h(0) = 0 $.

    According to the analysis in the first section, we have known that the function $ \psi_\eta $ given by (1.16) plays a vital role in studying the dynamics of system (3.1). We firstly give two lemmas that are useful in analyzing dynamics of $ \psi_\eta $, as follow:

    Lemma 3.1. Let $ I $ be a non-degenerate closed interval, $ J\subseteq\mathbb{R} $ and $ F:I\to J $, $ G:J\to \mathbb{R} $. Assume that $ V_JG < +\infty $ and $ F $ is piecewise monotone. Then

    $ V_{F(I)}G\leq V_IG\circ F < +\infty. $

    Proof. Let $ I_1 = [a, b] $ be a monotone interval of $ F $. Without loss of generality, let $ F $ monotonically increase in $ I_1 $. Let

    $ P:a = t_0\leq t_1\leq\cdots\leq t_n = b, $

    be a partition of $ I_1 $. Then

    $ V_{I_1}\left[G\circ F, P\right] = V_{F(I_1)}\left[G, P'\right]\leq V_{F(I_1)}G, $

    where

    $ P':F(a) = F(t_0)\leq F(t_1)\leq\cdots\leq F(t_n) = F(b), $

    is a partition of $ F(I_1) $. Consequently, $ V_{I_1}G\circ F\leq V_{F(I_1)}G $. Conversely, if

    $ P':F(a) = q_0\leq q_1\leq\cdots\leq q_n = F(b), $

    is a partition of of $ F(I_1) $, then

    $ P:a = \tilde{F}^{-1}(q_0)\leq \tilde{F}^{-1}(q_1)\leq\cdots\leq \tilde{F}^{-1}(q_n)\leq q_{n+1} = b, $

    where $ \tilde{F}^{-1}(q) = \min\{x\in I_1, F(x) = q\} $, is a partition of $ I_1 $. That implies

    $ V_{F(I_1)}\left[G, P'\right] = V_{I_1}\left[G\circ F, P\right]\leq V_{I_1}G\circ F. $

    Therefore, $ V_{I_1}G\circ F = V_{F(I_1)}G $.

    Let $ (I_j)_{1\leq j\leq n} $ be a finite sequence consisting of monotone intervals of $ F $. Assume $ \sqcup_{1\leq j\leq n} I_j = I $ and $ \#(I_i\cap I_j)\leq 1 $ provided $ i\ne j $. Then,

    $ VIGF=nj=1VIjGF=nj=1VF(Ij)GV1jnF(Ij)G=VF(I)G.
    $

    Lemma 3.2. Let $ I $ be a closed interval and $ F\in \mathcal{C}^0(I, I) $ be piecewise monotone. If there exist non-degenerate subintervals $ I_1, I_2\subseteq A $ with $ Card(I_1\cap I_2)\leq 1 $ such that $ I_2\subseteq F(I_1) $ and $ I_1\cup I_2\subseteq F(I_2) $, then

    $ lim infn+lnVI1I2Fnnln1+52.
    $
    (3.2)

    Proof. Let two subintervals $ I_1, I_2\subseteq A $ satisfy the hypothesis. Let $ n\in\mathbb{N} $. Put

    $ x_n = V_{I_2}F^n,\quad y_n = V_{I_1\cup I_2}F^n, $

    where $ F^0 = Id_{\mathbb{R}} $ if $ n = 0 $. Note that $ F(\cdot) $ is continuous and piecewise monotone, it follows from Lemma 3.1 that

    $ x_{n+1} = V_{I_2}F^n\circ F\geq V_{F(I_2)}F^n\geq V_{I_1\cup I_2}F^n = y_n, $

    and

    $ yn+2=VI1Fn+1F+VI2Fn+1FVF(I1)Fn+1+VF(I2)Fn+1VI2Fn+1+VI1I2Fn+1=xn+1+yn+1yn+yn+1.
    $

    Let $ (z_n)_{n\in\mathbb{N}} $ be Fibonacci sequence, i.e., $ z_{n+2} = z_{n+1}+z_n $, with the initial data $ z_0 = 0, z_1 = 1 $. It is well known that

    $ z_n = \frac{1}{\sqrt{5}}\left[\left(\frac{1+\sqrt{5}}{2}\right)^n-\left(\frac{1-\sqrt{5}}{2}\right)^n\right],\quad n\in\mathbb{N}. $

    It is clear that $ y_1 > 0 $ and $ \forall n\in\mathbb{N}, y_n\geq y_1z_n $. Therefore,

    $ \liminf\limits_{n\to+\infty}\frac{\ln V_{I_1\cup I_2}F^n}{n} = \liminf\limits_{n\to+\infty}\frac{\ln y_n}{n}\geq\lim\limits_{n\to+\infty}\frac{\ln y_1z_n}{n} = \ln\frac{1+\sqrt{5}}{2}. $

    This competes the proof.

    Proposition 3.1. Let $ h\in\mathcal{C}^0(\mathbb{R}) $ satisfy the hypothesises $ (A1)-(A2) $ and $ \psi_\eta $ be given by $(1.16)$. In addition, assume $ h(0) = 0 $ and

    (i) non-constant function $ Q:t\mapsto\frac12\left(h(t)+\frac{1}{c}t\right) $ is piecewise monotone,

    (ii) there is at least one solution to the equation $ h(y)+c^{-1}y = 0 $ in $ \mathbb{R}^* $.

    Then there exist $ A > 0 $ and a non-degenerate interval $ I\subseteq(0, +\infty)\setminus\{c^{-1}\} $ such that for any $ \eta\in I $, $ [0, A] $ is an invariant set of $ \psi_\eta $ and

    $ lim infn+lnV[0,A]ψnηnln1+52.
    $
    (3.3)

    Proof. We need to finish the following two steps:

    $ (1) $ Determine an invariant interval $ [0, A] $ of $ \psi_\eta $.

    $ (2) $ To make sure the hypothesises of Lemma 3.2 holds.

    It has been known that

    $ \psi_\eta = \gamma(\eta)\cdot\varphi = \gamma(\eta)\cdot Q\circ\phi^{-1}. $

    Put

    $ S^- = \{y < 0 \mid Q(y) = 0\ \text{and}\ \exists\ t\in(y,0), Q(t)\ne0\} $

    and

    $ S^+ = \{y > 0 \mid Q(y) = 0\ \text{and}\ \exists\ t\in(0,y), Q(t)\ne0\}. $

    By the hypothesises $ (i)-(ii) $, $ S^-\ne\emptyset $ or $ S^+\ne\emptyset $. Without loss of generality, assume $ S^-\ne\emptyset $. We take

    $ ˉy=maxS,A=c1ˉy,
    $
    (3.4)

    then

    $ ˉy<0, A=ϕ(ˉy)>0, ψη(A)=0.
    $
    (3.5)

    Moreover, let

    $ M = \max\{Q(t)\mid t\in[\bar{y},0]\},\quad m = \min\{Q(t)\mid t\in[\bar{y},0]\}. $

    By the definition of $ S^- $, it is evident that $ M = 0 $ or $ m = 0 $.

    Case 1: Assume $ M = 0 $. That implies $ m < 0 $. Put

    $ y_0 = \max\{ t\in[\bar{y},0]\mid Q(t) = m\},\quad x_0 = \phi(y_0). $

    For any $ \eta\in[0, c^{-1}) $, it is easily seen that

    $ \psi_\eta(x_0) = \max\{\psi_\eta(x)\mid x\in[0,A]\} = m\cdot\gamma(\eta). $

    and $ \psi_\eta(x) < \psi_\eta(x_0) $ provided that $ x\in[0, x_0) $. We need to prove that $ m\cdot\gamma(0) < A $. Proceed the proof by contradiction. Assume $ m\cdot\gamma(0)\geq A $, which implies

    $ \psi_0\left([0,x_0]\right)\cap\psi_0\left([x_0,A]\right)\supseteq [0,\psi_0(x_0)]\supseteq[0,A] = [0,x_0]\cup [x_0,A]. $

    Thus $ \psi_0 $ is turbulent. According to Lemma 3 in [11], $ \psi_0 $ has periodic points of all periods. But from the proof of Theorem 2.1, $ \psi_0 $ has no periodic points whose period is larger than 2. We have thus reached a contradiction. Since that $ m\cdot\gamma(\cdot) $ is strictly increasing in $ [0, c^{-1}) $ and $ m\cdot\gamma(\eta)\to+\infty $ as $ \eta\to(c^{-1})^- $, it is reasonable to put

    $ ˉη=γ1(mA)(0,c1).
    $
    (3.6)

    Then for $ \eta\in[0, \bar{\eta}] $, we have

    $ \forall x\in[0,A],\quad 0\leq \psi_\eta(x)\leq A, $

    that is to say $ [0, A] $ is an invariant set of $ \psi_\eta $ provided $ \eta\in[0, \bar{\eta}] $. Put

    $ η_={γ1(mx0), if mx0<1,0,  else,
    $
    (3.7)

    then $ \underline{\eta} < \bar{\eta} $ and $ \psi(x_0) > x_0 $ if $ \eta > \underline{\eta} $. For $ \eta\in[\underline{\eta}, c^{-1}) $, consider the following set:

    $ S(\eta,x_0) = \{x\in(0,x_0)\mid \gamma(\eta)\cdot\varphi(x) = x_0\}. $

    It is clear that $ S(\eta, x_0) $ is closed. Since $ x_0 $ is uniquely determined by $ h $ and $ h $ is independent of $ \eta $, the following function is well-defined:

    $ α(η)=maxS(η,x0),η[η_,c1).
    $
    (3.8)

    We first prove that $ \alpha(\cdot) $ is strictly monotonically decreasing in $ [\underline{\eta}, c^{-1}) $. Let $ t_1, t_2\in(\underline{\eta}, c^{-1}) $ with $ t_1 < t_2 $. From the definition of $ \alpha(\cdot) $, it follows

    $ γ(t1)φ(α(t1))=x0,γ(t1)φ(x)>x0,x(α(t1),x0],γ(t2)φ(α(t2))=x0,γ(t2)φ(x)>x0,x(α(t2),x0].
    $
    (3.9)

    Since that

    $ \gamma(t_1)\cdot\varphi(\alpha(t_2)) < \gamma(t_2)\cdot\varphi(\alpha(t_2)) = x_0,\ \gamma(t_1)\cdot\varphi(x_0) > x_0, $

    by applying the continuity of $ \gamma(t_1)\cdot\varphi(\cdot) $ there exists $ x\in(\alpha(t_2), x_0) $ such that $ \gamma(t_1)\cdot\varphi(x) = x_0 $. According to the definition of $ \alpha(t_1) $, we have

    $ \alpha(t_2) < x\leq\alpha(t_1). $

    Therefore, $ \alpha(\cdot) $ strictly monotonically decreases in $ [\underline{\eta}, c^{-1}) $.

    Next, prove that $ \alpha(\cdot) $ is continuous from the left. Let $ t_0\in(\underline{\eta}, x_0) $ be fixed. Since $ \alpha $ is monotone, we have

    $ α(t0)Ldef=limtt0α(t)<+.
    $
    (3.10)

    By the continuity of $ \eta\mapsto\gamma(\eta) $ and $ x\mapsto\varphi(x) $, we obtain

    $ limε0+γ(t0ε)=γ(t0),limε0+φ(α(t0ε))=φ(L),
    $
    (3.11)

    which gives $ \gamma(t_0)\cdot\varphi(L) = \gamma(t_0)\cdot\varphi(\alpha(t_0)) = x_0 $. From the definition of $ \alpha(\cdot) $ and (3.10), it follows

    $ L\leq\alpha(t_0)\leq L = \lim\limits_{t\to t^-_0}\alpha(t). $

    Therefore $ \alpha(\cdot) $ is continuous from the left.

    In particular, $ \psi_{\bar{\eta}}(\alpha(\underline{\eta})) = \gamma(\underline{\eta})\cdot\varphi\left(\alpha(\bar{\eta})\right) = x_0 $ and $ \psi_{\bar{\eta}}(x_0) = \gamma(\underline{\eta})\cdot\varphi(x_0) = A $. Define a function as follow:

    $ K(s)=α(η_s)ψ2ˉηs(x0).
    $
    (3.12)

    It is easily seen that $ K(\cdot) $ is continuous from right at $ s = 0 $. Note that

    $ K(0) = \alpha(\underline{\eta})-\psi_{\bar{\eta}}(\psi_{\bar{\eta}}(x_0)) = \alpha(\underline{\eta})-\psi_{\bar{\eta}}(A) = \alpha(\underline{\eta}) > 0, $

    hence there exists $ \rho_0\in(0, \bar{\eta}-\underline{\eta}) $ such that

    $ s[0,ρ0],K(s)=α(η_s)ψ2ˉηs(x0)>0.
    $
    (3.13)

    Put

    $ I=[ˉηρ0,ˉη]
    $
    (3.14)

    and

    $ J1=[α(η),x0], J2=[x0,ψη(x0)], ηI.
    $
    (3.15)

    Let $ \eta\in I $ be fixed. By (3.13),

    $ ψη(J1)[ψη(α(η)),ψη(x0)]=J2,ψη(J2)[ψ2η(x0),ψη(x0)][α(η),ψη(x0)]J1J2,
    $

    which is to say hypotheses of Lemma 3.2 holds. Therefore for any $ \eta\in I $, by applying Lemma 3.2 we have

    $ lim infn+lnV[0,A]ψnηnlim infn+lnVJ1J2ψnηnln1+52.
    $
    (3.16)

    For the case $ m = 0 $, one just needs to consider $ \eta > c^{-1} $. The proof is similar, we omit it.

    Theorem 3.1. Consider the system $(3.1)$ with the hypotheses $ (A1)-(A2) $. Suppose $ h(0) = 0 $ and

    (i) non-constant function $ Q:t\mapsto\frac12\left(h(t)+\frac{1}{c}t\right) $ is piecewise monotone,

    (ii) there is at least one solution to the equation $ h(y)+c^{-1}y = 0 $ in $ \mathbb{R}^* $.

    Then there exists a non-degenerate interval $ I $ such that for all $ \eta\in I $ the system $(3.1)$ is chaotic in the sense of Definition 1.1.

    Proof. Without loss of generality, assume the equation $ h(y)+c^{-1}y = 0 $ has at least one solution in $ (-\infty, 0) $. We still use the symbols given in the proof of Proposition 3.1, such as $ A $, $ x_0 $, $ \alpha(\eta) $ and so on. Let $ I $ and $ J_1, J_2 $ be given by (3.14) and (3.15), respectively. Take $ \eta\in I $.

    Let $ (w_0, w_1) $ be the initial data of system (3.1) and $ (w, w_x, w_t) $ be a solution of system (3.1). Recall the Riemann invariants as follow

    $ u(x,t) = \frac{cw_x(x,t)+w_t(x,t)}{2c},\ v(x,t) = \frac{cw_x(x,t)-w_t(x,t)}{2c},\ x\in[0,1],\ t\geq 0, $

    and $ u_0(\cdot) = u(\cdot, 0) $ and $ v_0(\cdot) = v(\cdot, 0) $ are the initial data. Let $ t > 2c^{-1} $, $ n = [\frac{ct}{2}] $ and $ \tau = t-2nc^{-1} $ be fixed. It is clear that $ 0\leq c\tau < 2 $. By using (1.13), (1.14) and (1.15), we obtain

    $ u(x,t)={γ1(η)ψnη(γ(η)u0(x+cτ)),0cτ1x,γ1(η)ψn+1η(v0(2xcτ)),1x<cτ2x,γ1(η)ψnη(γ(η)u0(x+cτ2)),2x<cτ<2,
    $
    (3.17)

    and

    $ v(x,t)={ψnη(v0(xcτ)),0cτx,ψnη(γ(η)u0(cτx)),x<cτx+1,ψn+1η(v0(2+xcτ)),x+1<cτ<2.
    $
    (3.18)

    Since $ [0, A] $ is an invariant set of $ \psi_\eta $, $ |u|+|v| $ is uniformly bounded if

    $ Range(γ(η)u0)Range(v0)[0,A].
    $
    (3.19)

    In addition, assume that $ u_0 $ and $ v_0 $ are piecewise monotone and

    $ [0,α(η)]Range(γ(η)u0)Range(v0).
    $
    (3.20)

    Consider the total variations of $ u(\cdot, t) $ and $ v(\cdot, t) $ on $ [0, 1] $. When $ 0\leq c\tau < 1 $, from (1.13)-(1.15) and Lemma 3.1, it follows

    $ V[0,1]u(,t)=V[0,1cτ]u(,t)+V[1cτ,1]u(,t)=V[cτ,1]γ1(η)ψnηγ(η)u0+V[1cτ,1]γ1(η)ψn+1ηv0=|γ1(η)|(V[cτ,1]ψnηγ(η)u0+V[1cτ,1]ψn+1ηv0),
    $
    $ V[0,1]v(,t)=V[0,cτ]v(,t)+V[cτ,1]v(,t)=V[0,cτ]ψnηγ(η)u0+V[0,1cτ]ψnηv0.
    $

    Note that $ 0 < \left|\gamma^{-1}(\eta)\right| < 1 $, by Lemma 3.1 and (3.20) we obtain

    $ V[0,1]u(,t)+V[0,1]v(,t)|γ1(η)|Vγ(η)u0([0,1])ψnη|γ1(η)|V[0,α(η)]ψnη|γ1(η)|Vψ2η[0,α(η)]ψn2η|γ1(η)|VJ1J2ψn2η.
    $

    Then Proposition 3.1 and (3.16) show that

    $ ln(V[0,1]u(,t)+V[0,1]v(,t))tln|γ1(η)|t+ln(VJ1J2ψn2η)tcln(VJ1J2ψn2η)2n+2+ln|γ1(η)|tc2ln1+52, as t+.
    $
    (3.21)

    When $ 1\leq c\tau < 2 $, in the same way we can obtain

    $ ln(V[0,1]u(,t)+V[0,1]v(,t))tcln(VJ1J2ψn1η)2n+2+ln|γ1(η)|tc2ln1+52, as t+.
    $
    (3.22)

    It is evident that

    $ V[0,1]u(,t)+V[0,1]v(,t)max{1,c1}(V[0,1]wx(,t)+V[0,1]wt(,t)).
    $
    (3.23)

    Hence

    $ \liminf\limits_{t\to+\infty}\frac{\ln\left(V_{[0,1]}w_x(\cdot,t)+V_{[0,1]}w_t(\cdot,t)\right)}{t}\geq\frac{c}{2}\ln\frac{1+\sqrt{5}}{2} > 0. $

    Therefore, the system (3.1) is chaotic in the sense of Definition 1.1.

    In this section, we mainly consider the oscillation problem governed by 1D wave equation with a fixed end. We will show that only the effect of self-regulation effect at one of the two ends can cause the onset of chaos. But if the fixed end is replaced by a free end, the system never has chaos. As we know, both a fixed end and a free end are called the energy-conserving boundary condition. Even though either a fixed end or a free end has the same effect to the energy of the system, the systems may show completely different dynamics: one with chaos, the other without chaos. Thus the relationship between chaos and energy of the system is much more complicated.

    Theorem 4.1. Consider an initial-boundary problem governed by 1D wave equation with a fixed end as follows:

    $ {ztt(x,t)c2zxx(x,t)=0,  x(0,1), t>0,z(0,t)=0,  t>0,zx(1,t)=h(zt(1,t)), t>0,z0(x)=z(x,0),z1(x)=zt(x,0),0x1.
    $
    (4.1)

    There exists a piecewise linear function $ h $ so that the system is chaotic in the sense of Definition 1.1.

    Proof. Without loss of generality, we take $ c = 1 $. Put

    $ hθ(x)={2x+2(θ23),  x[2(θ23),2),(1θ2)x, x[2,0],(θ1)x, x(0,1],2x+(θ+1), x(1,θ+1],x, else,
    $
    (4.2)

    where $ \theta\in(0, 1) $ is a fixed parameter. It is clear that $ h_\theta\in\mathcal{C}^0(\mathbb{R}) $ and

    $ \phi_{\theta}:x\mapsto\frac{1}{2}(h_\theta(x)-x) $

    is strictly monotonically decreasing and $ \phi_{\theta}(\mathbb{R}) = \mathbb{R} $. Then $ \phi_{\theta}^{-1} $ exists. Let

    $ Q_\theta(y) = \frac{1}{2}(h_\theta(y)+y), $

    then

    $ \psi_\theta(x) = x+\phi_{\theta}^{-1}(x) = Q_{\theta}\circ \phi_{\theta}^{-1}(x). $

    Let $ u = \frac{z_x+z_t}{2} $ and $ u = \frac{z_x-z_t}{2} $ be the Riemann invariants of system 4.1. For some initial date $ (u_0, v_0) $, suppose $ (u, v) $ is solved for $ t\leq2 $. Then for $ t = 2n+\tau $ with $ \tau\in[0, 2) $ and $ n\in\mathbb{N} $, we have

    $ u(x,t)=ψnθ(u(x,τ)),v(x,t)=ψnθ(v(x,τ)), x[0,1].
    $
    (4.3)

    Put

    $ I_1 = [-\theta-1,\frac12\theta-1],\ I_2 = [\frac12\theta-1,0],\ I_3 = [0,\theta^2]. $

    They are monotone intervals of $ \psi_\theta $, respectively. A simple calculation shows that

    $ \psi_\theta(-\theta-1) = 0,\ \psi_\theta(\frac12\theta-1) = \frac12\theta,\ \psi_\theta(0) = 0,\ \psi_\theta(\theta^2) = \theta^2-2. $

    Let $ \theta\in(0, 1) $ be sufficiently small such that $ \frac12\theta\geq\theta^2 $ and $ \theta^2-2\leq-\theta-1 $, which implies

    $ I_3\subseteq\psi_\theta(I_1)\cap\psi_\theta(I_2),\ I_1\cup I_2\subseteq\psi_\theta(I_3). $

    Let $ n\in\mathbb{N} $. Put

    $ x_n = V_{I_1}\psi^n_\theta,\ y_n = V_{I_2}\psi^n_\theta,\ z_n = V_{I_3}\psi^n_\theta. $

    Lemma 3.2 shows that

    $ x_{n+1}\geq V_{\psi_\theta(I_1)}\psi^n_\theta\geq V_{I_3}\psi^n_\theta = z_n,\ y_{n+1}\geq z_n,\ z_{n+2}\geq x_{n+1}+y_{n+1}, $

    which implies

    $ z_{n+2}\geq2z_{n}. $

    Inductively, we obtain

    $ z_{2n+2}\geq2^nz_2 > 0,\ z_{2n+1}\geq2^nz_1 > 0. $

    Therefore,

    $ lim infn+lnznn12ln2>0.
    $
    (4.4)

    A simple calculation shows that

    $ \psi_\theta(x) = 0,\ x\in(-\infty,-\theta-1]\cup[2(3-\theta^2),+\infty]. $

    By (4.3), $ |u|+|v| $ is always uniformly bounded for any initial date. Analogous to the proof of Theorem 3.1, there exists a large class of initial data such that

    $ \liminf\limits_{t\to+\infty}\frac{\ln\left(V_{[0,1]}z_x(\cdot,t)+V_{[0,1]}z_t(\cdot,t)\right)}{t}\geq\frac{1}{4}\ln2 > 0. $

    Therefore, the proof completes.

    We first consider a problem about perturbations at boundaries, described by the following model:

    $ {wtt(x,t)c2wxx(x,t)=0,  x(0,1), t>0,wx(0,t)=ηwt(0,t), ηc1,  t>0,wx(1,t)=c1wt(1,t)+εsin(wt(1,t)), t>0,w0(x)=w(x,0),w1(x)=wt(x,0),0x1.
    $
    (5.1)

    If $ \varepsilon = 0 $, for any initial data $ (w_0, w_1) $ and $ \eta\ne c^{-1} $, by using wave propagation method we can obtain

    $ \forall t > 2, x\in[0,1],\ w_x(x,t) = w_t(x,t) = 0. $

    That shows system (5.1) is of global asymptotical stability.

    Assume $ \varepsilon\ne0 $ and $ |\varepsilon| $ is sufficiently small. Put

    $ h_\varepsilon(y) = -c^{-1}y+\varepsilon\sin(y) $

    It is evident that $ i) \; y\mapsto h_\varepsilon(y)+c^{-1}y $ is piecewise monotone, $ ii) \; y\mapsto h_\varepsilon(y)-c^{-1}y $ strictly monotonically decreasing, $ iii) $ there are infinite solutions to the equation $ h_\varepsilon(y)+c^{-1}y = 0 $ in $ \mathbb{R}^* $. That shows the hypotheses of Theorem 3.1 hold. Therefore, system (5.1) is chaotic in the sense of Definition 1.1.

    The above analysis shows that this kind of system is no longer stable under small perturbations at one of the boundaries. That is to say this kind of system is not of structural stability.

    As we have stated in the previous sections, if the feedback control law is even chaos never happen. In particular, consider the system described by the following model:

    $ {ztt(x,t)zxx(x,t)=0,  x(0,1), t>0,zx(0,t)=ηzt(0,t), η1, or z(0,t)=0, t>0,zz(1,t)=cos(zt(1,t))1, t>0,z0(x)=z(x,0),z1(x)=zt(x,0),0x1.
    $
    (5.2)

    The feedback control law at the right end is periodic and even. Therefore there is no chaos in the system (5.2). Moreover, there should exists a compact attractor in the system 5.2.

    In this section, we will give some numerical simulations to validate the theoretical results of this paper. We first consider the system as follow:

    $ {ztt(x,t)zxx(x,t)=0,  x(0,1), t>0,z(0,t)=0,  t>0,zx(1,t)=h(zt(1,t)), t>0,z0(x)=z(x,0),z1(x)=zt(x,0),0x1,
    $
    (6.1)

    where

    $ h(x)={2x112,  x[112,2),34x, x[2,0],12x, x(0,1],2x+32, x(1,32],x, else.
    $
    (6.2)

    Choose

    $ w_0 = 0,\ w_1(x) = \frac14\sin^4(2\pi x),\ x\in[0,1], $

    as the initial data of system (6.1). Put

    $ ψ(x)={13x116,  x(14,112],7x, x(0,14],13x, x(34,0],13x+12, x[32,34],0, else.
    $
    (6.3)

    As we show in Section 4, the solution $ (z, z_x, z_t) $ can be represented by the iterations of this $ \psi_{\infty} $ and initial data $ (z_0, z_1) $. Note that $ h $ and $ \psi_\infty $ are piecewise linear, thus the feedback control at the boundaries is easy to implement and the solutions of system are much simpler to be solved and represented.

    We present the graphics in some detail, for $ z_x, \; z_t $ for $ 98\leq t\leq 100 $. Figure 1 shows that $ z_x, z_t $ are extremely oscillatory in every direction of space and time.

    Figure 1.  The profile of $ z_x $ (left) $ z_t $(right) for $ t \in [98,100] $.

    This work was supported by the National Natural Science Foundation of China under Grant 11671410 and the Natural Science Foundation of Guangdong Province, China (No.2022A1515012153).

    The authors declare there is no conflicts of interest.

    [1] Abdelazim T, Malik OP (2005) Power System Stabilizer Based on Model Reference Adaptive Fuzzy Control. Electr Pow Compo Sys 33: 985-998.
    [2] Abdelsalam AA, Gabbar HA, Sharaf AM (2014) Performance enhancement of hybrid AC/DC microgrid based D-FACTS. Int J Elec Power 63: 382-393.
    [3] Abdullah A, Zribi M (2009) Model reference control of LPV systems. J Franklin Inst 346: 854-871. doi: 10.1016/j.jfranklin.2009.04.006
    [4] Ali S, Qamar S, Khan L (2013) Hybrid Adaptive Recurrent Neurofuzzy Based SVC Control for Damping Inter-Area Oscillations. Middle-East Journal of Scientific Research 16: 536-547.
    [5] Azad ML, Singh SV, Khursheed A (2014) Improving Voltage Profile Of A Grid, Connected To Wind Farm Using Static Var Compensator. Int J Adv Engineer Technol 7: 1497-1506.
    [6] Bayat R, Ahmadi H (2013) Artificial Intelligence SVC Based Control of Two Machine Transmission System. I.J. Intelligent Systems and Applications 8: 1-8.
    [7] Bhole SS, Nigam P (2015) Improvement of Voltage Stability in Power System by Using SVC and STATCOM. International Journal of Advanced Research in Electrical, Electronics and Instrumentation Engineering 4: 749-755.
    [8] Biswas MM, Kamol KD (2011) Voltage Level Improving by Using Static VAR Compensator (SVC). GJRE: (J) General Engineering 11: 12-18.
    [9] Calderon J, Chamorro H, Ramos G (2012) Advanced SVC Intelligent Control to Improve Power Quality in Microgrids. Alternative Energies and Energy Quality (SIFAE), 2012 IEEE International Symposium on, 1-6.
    [10] Cerman O (2013) Fuzzy model reference control with adaptation mechanism. Expert Syst Appl 40: 5181-5187.
    [11] El-dessouky A, Tarbouchi M (2000) Model Reference Adaptive Fuzzy Controller For Induction Motor Using Auto-Attentive Approach. Industrial Electronics, 2000. ISIE 2000. Cholula, Puebla, Mexico.
    [12] El-dessouky A, Tarbouchi M (2001) Optimized Fuzzy Model Reference Learning Control for Induction Motor Using Genetic Algorithms. Industrial Electronics Society, 2001. IECON '01. The 27th Annual Conference of the IEEE 3: 2038-2043.
    [13] Fang DZ, Xiaodong Y, Chung TS, et al. (2004) Adaptive Fuzzy-Logic SVC Damping Controller Using Strategy of Oscillation Energy Descent. IEEE Transactions On Power Systems 19: 1414-1421.
    [14] Farsangi MM, Nezamabadi-pour H, Yong-Hua S, et al. (2007) Placement of SVCs and Selection of Stabilizing Signals in Power Systems. IEEE Transactions On Power Systems 22: 1061-1071. doi: 10.1109/TPWRS.2007.901285
    [15] Fernando I, Kwasnicki W, Gole A (1996). Modelling of conventional and advanced Static var compensators in an electromagnetic transient simulation program. Proceedings of the International Symposium on Modern Electric Power Systems 1: 60-70.
    [16] GENC I, USTA ö (2005) Impacts of Distributed Generators on the Oscillatory Stability of Interconnected Power Systems. Turk J Electr Eng Co 13: 149-161.
    [17] Goléa N, Goléa A, Benmahammed K (2002). Fuzzy Model Reference Adaptive Control. IEEE Transactions On Fuzzy Systems 10: 436-444.
    [18] Guang Z, Lijuan Z, Quanhai W, et al. (2007) The Research of Model Reference Adaptive Control of Static Var Compensator (SVC). 7th Internatonal Conference on Power Electronics, 300-304.
    [19] Hatziargyriou ND, Meliopoulos AP (2002) Distributed energy sources: Technical challenges. IEEE Power Eng Soc Winter Meeting 2: 1017-1022.
    [20] Hušek P, Cerman O (2013) Fuzzy Model Reference Control with Adaptation of Input Fuzzy Sets. Knowl-Based Syst 49: 116-122.
    [21] Ian A, Hiskens DJ (1992). Incorporation of SVCS Into Energy Function Methods. IEEE Transaction on Power Systems 7: 133-140.
    [22] Katiraei F, Iravani R, Hatziargyriou N, et al. (2008) Microgrids management. IEEE Power Energy Mag 6: 54-65.
    [23] Khan UN, Yan L (2008) Power Swing Phenomena and its Detection and Prevention. 7th EEEIC International Workshop on Environment and Electrical Engineering, 69-72.
    [24] Kroposki B, Lasseter R, Ise T, et al. (2008) Making microgrids work. IEEE Power and Energy Mag 6: 40-53.
    [25] Kumar R, Choubey A (2014) Voltage Stability Improvement by using SVC with Fuzzy Logic Controller in Multi-Machine Power System. International Journal of Electrical and Electronics Research 2: 61-66.
    [26] Layne J, K Passino (1993) Fuzzy model reference learning control for cargo ship steering. IEEE Control Systems Magazine 13: 23-34.
    [27] Mishra Y, Mishra S, Dong Z (2008) Rough Fuzzy Control of SVC for Power System Stability Enhancement. J Electr Eng Technol 3: 337-345.
    [28] Noroozian M, Andersson G (2001) A Robust Control Strategy for Shunt and Series Reactive Compensators to Damp Electromechanical Oscillations . IEEE Transactions on Power Delivery 16: 812-817.
    [29] Procyk T, Mamdani E (1979) A Linguistic Self-Organizing Process Controller. Automatica 15: 15-30.
    [30] Rogers KM (2009) Power System Control with Distributed Flexible AC Transmission System Devices. Master thesis, University of Illinois at Urbana-Champaign.
    [31] Sebastiana R, Quesadab J (2006) Distributed control system for frequency control in a isolated wind system. Renewable Energy 31: 285-305.
    [32] Wang J, Fu C, Zhang Y (2008) SVC Control System Based on Instantaneous Reactive Power Theory and Fuzzy PID. IEEE Transactions On Industrial Electronics 55: 1658-1665.
    [33] Xiongfei W, Guerrero J, Zhe C (2010) Control of grid interactive AC microgrids. IEEE International Industrial Electronics (ISIE), IEEE International Symposium on : 2211-2216.
    [34] Yen J, Langari R (1998) Fuzzy Logic: Intelligence, Control, and Information. Prentice Hall.
  • This article has been cited by:

    1. Xu Zhang, Iram Hussan, Manyu Zhao, Qigui Yang, Li–Yorke Chaos in Partial Differential Equations with or Without Energy Injection, 2025, 35, 0218-1274, 10.1142/S0218127425500245
  • Reader Comments
  • © 2016 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(6164) PDF downloads(1233) Cited by(0)

Figures and Tables

Figures(11)  /  Tables(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog