文章快速检索 高级检索
  力学学报  2016, Vol. 48 Issue (4): 953-962  DOI: 10.6052/0459-1879-16-044
0

页岩气专题论文

引用本文 [复制中英文]

陈章耀, 王亚茗, 张春, 毕勤胜. 双状态切换下BVP振子的复杂行为分析[J]. 力学学报, 2016, 48(4): 953-962. DOI: 10.6052/0459-1879-16-044.
[复制中文]
Chen Zhangyao , Wang Yaming , Zhang Chun , Bi Qinsheng . COMPLICATED BEHAVIORS AS WELL AS THE MECHANISM IN BVP OSCILLATOR WITH SWITCHES RELATED TO TWO STATES[J]. Chinese Journal of Ship Research, 2016, 48(4): 953-962. DOI: 10.6052/0459-1879-16-044.
[复制英文]

基金项目

国家自然科学基金(11472115, 11572141, 11502091) 和镇江市科技攻关基金(GY2013032, GY2013052) 资助项目.

通讯作者

毕勤胜, 教授, 主要研究方向:非线性动力学.E-mail:qbi@ujs.edu.cn

文章历史

2016-02-02 收稿
2016-04-21 录用
2016-04-27网络版发表
双状态切换下BVP振子的复杂行为分析
陈章耀, 王亚茗, 张春, 毕勤胜     
江苏大学土木工程与力学学院, 镇江 212013
摘要: 非线性切换系统具有广泛的工程背景,而传统的非线性理论不能直接用来解决其中的问题,因而成为当前国内外热点和前沿课题之一. 目前相关工作大都是围绕固定时间或单状态切换开展的,而实际工程系统大都属于多状态切换问题,同时多状态切换涉及到更为丰富的动力学行为. 本文基于两广义BVP 振子,通过引入双向切换开关,构建了双状态切换下的非线性动力学模型,进而研究状态切换导致的各种运动模式及其相应的产生机制. 应用非光滑系统的Poincaré映射理论,推导了双状态切换下的Lyapunov 指数的计算公式,结合子系统的分岔分析,得到了切换系统随分岔参数变化的动力学演化过程及其相应的最大Lyapunov 指数的变化情况. 得到了双状态切换条件下系统存在着各种形式的振荡行为,分析了诸如周期突变等现象及通往混沌的倍周期分岔道路,揭示了不同运动模式的产生机制及倍周期序列的本质. 与固定时间切换和单状态切换系统不同,双临界状态切换系统存在着更为丰富的非线性现象,其主要原因在于双状态切换会产生更多的切换点,且切换点的位置更加多变. 同时切换系统的倍周期分岔序列与光滑系统中的倍周期分岔序列不同,切换系统的倍周期分岔序列只对应于切换点数目的成倍增加,而其相应的周期一般不对应于严格的周期倍化过程.
关键词: 切换系统    广义BVP 振子    双状态切换    Lyapunov 指数    
0 引言

非线性切换系统具有广泛的工程背景,其复杂动力学行为及其机理研究是当前国内外动力学与控制领域的热点课题之一[1-7],如电路系统中继电、饱和、滞环等环节以及发电机组或大型设备进出电网[8],汽车行驶过程中的换档[9]等都有非线性切换系统的应用.目前关于非线线性切换系统的工作大都都局限于系统的稳定性及稳定条件分析,如Zhang等[10]给出了一类非线性切换系统的稳定性条件,Xiang等[11]基于状态更新探讨了一类不确定非线性切换系统的鲁棒性,Jin等[12]提出了一种控制非线性切换序列稳定性的新方法.但是,上述相关工作包括其他一些相关类型的研究[13-17]基本上都是围绕平衡点的稳定性开展的,即使对于控制,其目标也是尽快实现平衡点的稳定[18-24],而实际系统存在着各种各样的平衡态,可能是周期振荡、概周期振荡、有时甚至是混沌振荡.同时,不同参数的变化会影响切换系统的特性,因此需要深入探讨切换系统的动力学演化过程,进而揭示其复杂性本质.

非线性切换系统有不同的切换模式,主要有时间切换、状态切换和混合切换等模式.对于时间切换模式如周期切换等非线性切换系统,国内外学者开展了研究工作,取得了一定成果[25-28].而切换模式依赖于状态变量的非线性切换系统的研究相对较少,如Kousaka等[29]介绍了一种二维非线性切换系统在状态切换模式下的动力学特性,给出了一种切换系统分岔的一种验证方法,并且发现切分岔和全局分岔是产生各种周期解和混沌的主要原因;另外张春等[30]在Duffing振子和vanderPol振子之间引入混合切换条件,给出了相关Lyapunov指数的算法,并研究了其动力学特性.但是Kousaka给出的系统维数较低,而张春所做的混合切换没有具体的物理模型进行验证,还存在可优化之处,因此相关工作尚处在初步阶段,许多问题尤其是各种状态切换下的振荡行为及机理分析还有待进一步深入.

本文为了分析双状态切换面导致的各种复杂行为及其演化过程,引入典型的广义BVP振子[31],构建双状态切换下的非线性系统模型,以期揭示其复杂本质.这对于深入理解非线性状态切换系统的各种新的复杂现象,进而发展非线性理论具有一定的科学意义,同时对于实际切换系统中的新现象解释、参数识别、模型论证以及控制技术等等也具有一定的实际应用价值.

1 切换系统的数学模型

考察引入双向开关切换的广义BVP振子,如图 1所示.该系统由一个线性电阻$r$,2个电容$C$,一个电感$L$和2个不同的非线性电阻${R_{N1}}$$R_{\rm N2} $组成. 开关S随电压$v_1 $的不同会在${\rm SW}_1$${\rm SW}_2 $两点间切换. 其切换规律为,开关首先位于${\rm SW}_1$,在系统振荡过程中电压$v_1 $的值达到$v_1 = v_{\rm S2}$时,开关切换到${\rm SW}_2 $,当系统振荡过程中$v_1 $的值达到$v_1 =v_{\rm S1} $时,开关切回到${\rm SW}_1 $. 当开关固定于${\rm SW}_1$${\rm SW}_2 $时,均为一广义BVP振子.系统状态可以由如下方程描述

图 1 广义BVP切换电路 Fig. 1 The extended BVP switching circuit
$\begin{equation}\label{eq1} \left. {\begin{array}{l} C\dfrac{{\rm d}v_1 }{{\rm d}t} =-i_1-g_\sigma (v_1 ),\\ C\dfrac{{\rm d}v_2 }{{\rm d}t} = i_1-\dfrac{v_2 }{r},\\ L\dfrac{{\rm d}i_1 }{{\rm d}t} = v_1-v_2 . \\ \end{array}}~~\sigma = 1,2\right\}\end{equation}$ (1)

其中$\sigma = 1$为开关S位于${\rm SW}_1$,系统受子系统1控制;$\sigma = 2$则开关位于${\rm SW}_2$,系统受子系统2控制. $g_\sigma (v_1 ) =-a_\sigma v_1-b_\sigma\tanh (c_\sigma v_1 )$表示非线性电阻的电压和电流特性,$a_\sigma $,$b_\sigma $$c_\sigma $为常数($\sigma = 1,\;2$).

引入$\tau = \dfrac{t}{\sqrt {LC} }$$A_\sigma =a_\sigma \sqrt {\dfrac{L}{C}} $$B_\sigma = b_\sigma c_\sigma\sqrt {\dfrac{L}{C}} $$\delta = \dfrac{1}{r}\sqrt{\dfrac{L}{C}} $$x = \dfrac{v_1 }{b_\sigma }\sqrt {\dfrac{C}{L}}$$y = \dfrac{v_2 }{b_\sigma }\sqrt {\dfrac{C}{L}}$$z = \dfrac{i_1 }{b_\sigma }$($\sigma =1,\;2$),式(1)可表示为如下无量纲形式

$\begin{equation}\label{eq2} \left. {\begin{array}{l} \dot x=-z + A_\sigma x + \tanh B_\sigma x \\ \dot y = z-\delta y \\ \dot z= x-y \\ \end{array}},~~~ \sigma = 1,2\right\}\end{equation}$ (2)

切换系统方程(2)可简写为

$\begin{equation}\label{eq3} \dot { X}= { F}_\sigma ({ X}),\quad\sigma =1,2\end{equation}$ (3)

其中${ X} = \left\{ {x,y,z} \right\}^{\rm T}$${ F}_\sigma$为相应的向量场,$\sigma $表示上述切换律,数学定义如下

$\sigma = \left\{ {\matrix{ {1,\quad \sigma = 2} \hfill \cr {2,\quad \sigma = 1} \hfill \cr } {\rm{and}}\matrix{ {x = {x_{01}}} \cr {x = {x_{02}}} \cr } } \right.$

其相应的无量纲数学模型以一定的规律在两个子系统之间切换,显然该切换系统存在两个不同的状态切换面,即$\varSigma_1 $:${ X} = { X}_1 $$\varSigma _2 $:${ X} = { X}_2$. 双状态切换可能导致更为丰富的动力学行为.

2 子系统分岔分析

由于切换系统(3)涉及到两个向量场,其动力学行为与两个向量场的动力学特性密切相关.为揭示该切换系统的复杂行为,首先讨论两子系统的分岔特性.

为方便起见,在此假设参数$\delta = 1.2$$A_\sigma = B_\sigma =q$,($\sigma = 1,2)$.显然,随$q$的变化,两子系统具有相同的动力特性(见图 2,图中实线表示稳定的平衡点,虚线表示不稳定的平衡点).

图 2可以看出,随着参数$q$的增加,当$q <0.425$时,子系统具有稳定的焦点,记为$E_0 $. 在$q = 0.425$$H_1$点处,其相应的特征值出现一对纯虚根,进一步分析可以发现,此时会产生超临界Hopf分岔,也即从平衡点$E_0$处出现稳定的极限环,记为$LC_0 $图 2(a)中用实心点表示).当$q$增大到$q = 0.62$时,平衡点数量发生变化,出现叉式分岔.不稳定的平衡点$E_0 $分岔为三个不稳定的平衡点,记为$E_0 ,E_\pm $.其中$E_\pm $相互对称. 同时$E_\pm $$H_2$相应的特征值又出现一对纯虚根,也即从$E_\pm$由超临界Hopf分岔导致两对称的极限环$LC_\pm $.

图 2 子系统随参数变化的分岔图 Fig. 2 The bifurcation diagram of the subsystem

数值仿真可以发现与上述平衡点分岔特性相吻合的动力学演化过程,即在$H_1$处出现围绕平衡点$E_0 $的稳定的极限环. 当$q$增大到$q =0.62$时有三个稳定的极限环共存,其中围绕$E_0 $的极限环在$q =0.833$时折回消失. 而围绕平衡点$E_\pm$的极限环具有相互对称的演化特征,所以图 2(b)仅给出了围绕$E_ +$的极限环$LC_ + $的演化过程. 稳定的极限环$LC_ + $$q =0.802$会失稳导致混沌运动. 在$q = 0.863$,由$LC_\pm$分别产生的两相互对称的混沌运动会相互作用,导致混沌激变现象,产生扩大了的混沌吸引子.$LC_\pm $$q = 1.17$处折回消失.

3 切换系统的李雅普诺夫指数计算

由于切换系统是一个非光滑系统,其诸如周期运动和混沌运动的动力学行为必须通过一定的统计规律加以验证,本文用Lyapunov指数加以验证.

对于平衡点和周期运动,很多学者常用系统在平衡点的Jacobi矩阵特征值实部以及周期运动扰动方程的特征指数实部来判断其稳定性.Lyapunov指数是对上述特征值和特征指数的一种推广,给出了对系统任意相邻相轨线平均发散程度或平均收敛程度的一种度量.因此,它可以度量系统吸引子的稳定性,还可以判定吸引子的类型,如平衡点、周期运动、概周期运动及混沌运动等,是目前判断各种动力学行为的一种可靠定量方法.但对于非光滑系统(尤其是状态切换系统),李雅普诺夫指数的计算不能直接使用光滑系统的方法,需要采用复合系统的微分包含理论和Poincaré映射的数学方法,对于切换部分需要使用复合映射及复合映射的链式法则进行计算,并进行必要的处理. 下面分2节介绍其具体计算原理.

3.1 Lyapunov指数

在计算非线性切换动力系统Lyapunov指数时,采用Poincaré映射的数学方法,即应用Poincaré映射$P$将三维切换系统(3)转化为三维离散动力系统

$\begin{equation}\label{eq4} { X}_{(k)} = P({ X}_{(k-1)} ),{ X} \in R^3,k = 1,2,\cdots ,N\end{equation}$ (4)

注意到在$(k-1)$时刻,若基准轨道${ X}_{^{(k-1)}} ({X}_{^0} )$和邻近轨道${ X}_{^{(k-1)}} ({ X}_{^0} + \delta{ X}_{^0} )$的距离

$\delta { X}_{^{(k-1)}} = {X}_{^{(k-1)}} ({ X}_{^0} + \delta { X}_{^0} )-{X}_{^{(k-1)}} ({ X}_{^0} )$

充分小,式(4)在${X}_{^{(k-1)}} $处的线性化方程为

$\delta { X}_{^{(k)}} = DP({ X}_{^{(k-1)}} )\delta {X}_{^{(k-1)}}$

由上式递推得到

$\delta { X}_{^{(k)}} = DP_{^k} ({ X}_{^{(0)}} )\delta {X}_{^{(0)}} ,\quad {{ X} \in R^3},$k = 1,2,\cdots ,N$

式中

$\begin{equation}\label{eq5}\begin{array}{c} DP_k ({ X}_{(0)} ) = DP({ X}_{(k-1)} \times \\ DP({ X}_{(k-2)} )\times \cdots \times DP({ X}_{(0)} ) \\ \end{array}\end{equation}$ (5)

由此离散系统可得整个切换系统的Lyapunov指数

$\begin{equation}\label{eq6} \lambda = \mathop {\lim }\limits_{N \to \infty }\frac{1}{N}\sum\limits_{k = 0}^{N-1} {\ln \left| {DP_{^k} ({X}_{^{(0)}} )} \right|}\end{equation}$ (6)

对于非光滑系统如切换系统必须使用复合映射及复合映射求导的链式法则求出对应的Poincaré映射$P({X})$$DP({ X})$.

3.2 三维状态切换系统的复合映射

$\varphi _1 ({ X},\tau )$,$\varphi _2 ({ X},\tau)$分别为切换系统(3)中两子系统的解,由于整个系统是由两个不同子系统通过状态切换面连接的切换系统,因此存在两个切换面即局部超曲面$\varSigma _1 $$\varSigma_2$,另外由于状态切换,这里还存在一个投影面$\varSigma$(如图 3),3个超曲面分别表示为

$\varSigma{}_1=\{({ X},t) = (x,y,z,t) \in R^3\times R^+|x =x_{01}\}$

其中,$x_{01} $为固定值,即切换曲面方程为 $q_1 ({ X}) = x-x_{01} = 0$

$\varSigma{}_2 { = \left\{ {({ X},t) = (x,y,z,t) \in R^3\times R^ + \left| {x = x_{02} } \right.} \right\}}$

其中,$x_{02} $为固定值,即切换曲面方程为 $q_2 ({ X}) = x-x_{02} = 0$

$\varSigma{} { = \left\{ {{ X} = (x,y,z) \in R^3\left| {x =x_{\rm ref} } \right.} \right\}}$

这样在$\varSigma_1 $$\varSigma_2$之间就存在两个局部映射$P_1$$P_2 $. 而在$\varSigma_1 $$\varSigma$之间存在两个投影映射$h$$h^{-1}$.Poincaré映射$P$将被分成4个部分

$\varSigma \buildrel {h^{-1}} \over \longrightarrow \varSigma{}_1\buildrel {P_1 } \over \longrightarrow \varSigma{}_2 \buildrel {P_2 } \over \longrightarrow \varSigma{}_1 \buildrel h \over \longrightarrow \varSigma $

其中

$\begin{array}{lll}P_1:\varSigma _1 \to \varSigma _2 \\{ X}_{2i-1} \mapsto { X}_{2i} = \varphi _1 (T_1 ,{X}_{2i-1} )\\ P_2 : \varSigma _2 \to \varSigma _1\\ { X}_{2i + 1} \mapsto \varphi _2 (T_2 ,{ X}_{2i} ) = \varphi _2 (T_2 ,\varphi (T_1 ,{ X}_{2i-1} ))=\\ \quad~~ \varphi _2 \circ \varphi _1 (T_1 + T_2 ,{ X}_{2i-1} )\\ h:\varSigma _1 \to \varSigma ,~~(x,y,z) \mapsto (y,z)~~\mbox{即}\\ h_1 (x,y,z) = y,~~ h_2 (x,y,z) = z\\ h^{-1}:\varSigma \to \varSigma _1,~~(y,z) \mapsto (x,y,z)~~\mbox{即}\\ h_1^{-1} (y,z) = x,~~h_2^{-1} (y,z) = y\\ h_3^{-1} (y,z) = z\end{array}$
图 3 Poincaré截面和Poincaré映射示意图 Fig. 3 Poincaré section and Poincaré mapping

于是$P$的复合映射为$P = h \circ P_2 \circ P_1 \circ h^{-1}$.由复合映射求导的链式法则得出

$\begin{equation}\label{eq7} { D}P = \nabla h \circ { D}P_2 \circ { D}P_1\circ \nabla h^{-1}\end{equation}$

设映射$P_\sigma $

$\frac{{\rm d}{ X}}{{\rm d}t} = { F}_\sigma ({ X}),\quad \sigma = 1,2$ (7)

由于$P_1 $为状态切换[20],所以

${ D}P_1 = \left( {{ I}-\frac{{ F}_1 ({ X}_{2i})\times \nabla q_1 }{\nabla q_1 \cdot { F}_1 ({ X}_{2i} )}}\right)\times \left( {\frac{\partial \varphi _1 }{\partial { X}_{2i-1} }} \right)$

同理$P_2 $也为状态切换,所以

${ D}P_2 = \left( {{ I}-\frac{{ F}_2 ({ X}_{2i + 1})\times \nabla q_2 }{\nabla q_2 \cdot { F}_2 ({ X}_{2i + 1})}} \right)\times \left( {\frac{\partial \varphi _2 }{\partial { X}_{2i} }} \right)$

其中,$\nabla q_\sigma = [1,0,0],(\sigma =1,2)$为两切换曲面的法向量. 而$\dfrac{\partial \varphi _1}{\partial { X}_{2i-1} }$$\dfrac{\partial \varphi _2}{\partial { X}_{2i} }$

$\left. {\begin{array}{ll} \dfrac{{\rm d}}{{\rm d}t}\left[{\dfrac{\partial \varphi _\sigma }{\partial { u}}} \right] =\dfrac{\partial { F}_\sigma }{\partial x}\dfrac{\partial x}{\partial { u}} \\ \dfrac{\partial { X}}{\partial { u}} = { I} \\ \end{array}},\quad \sigma = 1,2\right\}$

方程组求出,这里${ I}$为单位矩阵,${ u}$${ X}_{2i-1}$${ X}_{2i} $,$\sigma $分别等于1或2. 而

$\nabla h = \left( {{\begin{array}{*{20}c} {\dfrac{\partial h_1 }{\partial x}} & {\dfrac{\partial h_1 }{\partial y}} & {\dfrac{\partial h_1 }{\partial z}} \\ {\dfrac{\partial h_2 }{\partial x}} & {\dfrac{\partial h_2 }{\partial y}} & {\dfrac{\partial h_2 }{\partial z}} \\ \end{array} }} \right) =\left[{{\begin{array}{*{20}c} 0 & 1 & 0 \\ 0 & 0 & 1 \\ \end{array} }} \right] $
$\nabla h^{-1} = \left( {{\begin{array}{*{20}c} {\dfrac{\partial h_1^{-1} }{\partial y}} & {\dfrac{\partial h_1^{-1} }{\partial z}} \\ {\dfrac{\partial h_2^{-1} }{\partial y}} & {\dfrac{\partial h_2^{-1} }{\partial z}} \\ {\dfrac{\partial h_3^{-1} }{\partial y}} & {\dfrac{\partial h_3^{-1} }{\partial z}} \\ \end{array} }} \right) =\left[{{\begin{array}{*{20}c} 0 & 0 \\ 1 & 0 \\ 0 & 1 \\ \end{array} }} \right]~~~~~~~~~~~$

注意到,由于此处切换面均与$x$有关,由上式可知$\nabla h,\nabla h^{-1}$两个投影映射的第一列和第一行元素均为零,根据式(7),${D}P$变成了$2\times 2$的矩阵,故Lyapunov指数只有两个.通过式(4)$\sim$式(7)推导,作者编写了状态切换系统随参数变化的Lyapunov指数的计算程序,对状态切换系统进行了数值模拟,以便于分析其动力学行为.

4 切换系统的动力学行为分析

由于状态切换条件的引入,整个系统轨迹会依次受到不同向量场的控制,从而引起复杂的动力学行为. 首先考察两子系统的运动特征.

在子系统1中选取参数为$\delta = 1,2$,$A_1 = B_1 =0.7$,由第2节分析可知,子系统1在该参数下存在3个稳定周期运动,如图 4(a)所示.由于选取切换系统的初始条件为(0.1,0.1,0.1),故该子系统表现为围绕平衡点$E_1 $(1.629,1.629,1.955)的周期振荡,如图 4(b)所示.

图 4 系统的稳定吸引子($\delta = 1.2,A_1 = B_1 = 0.7)$ Fig. 4 The stable attractors of subsystem ($\delta = 1.2,A_1 =B_1 =0.7$)

在子系统二中固定参数$\delta = 1.2$,取参数$A_2 = B_2 =q$为分岔参数,两切换面为

$\eqalign{ & \sum {_1} = \left\{ {(X,t) = (x,y,z,t) \in {R^3} \times {R^ + }\left| {x = 0.0} \right.} \right\} \cr & \sum {_2} = \left\{ {(X,t) = (x,y,z,t) \in {R^3} \times {R^ + }\left| {x = 2.0} \right.} \right\} \cr} $

确定两子系统参数及切换面之后,利用数值仿真画出切换系统随参数$q$变化的分岔图和最大Lyapunov指数图(图 5).从图 5(a)中可以看出,当混沌在$q = 0.631$消失后,在$q \in[0.631,0.780]$区域会产生各种类型的稳定的周期振荡,这也可以从其相应的最大Lyapunov指数得到证实(见图 5(b)).

图 5 切换系统的分岔图和最大Lyapunov指数图 Fig. 5 The bifurcation characteristic of switched system

稳定的周期振荡在$q =0.775$会进入混沌运动,也即由周期振荡直接进入混沌,从图 5(b)中也可以发现,其最大Lyapunov指数在$q= 0.775$时穿越零值,而在$q =0.894$时,其最大Lyapunov指数再次穿越零值,导致混沌现象消失,产生周期振荡.随着$q$的变化,切换系统会产生周期突变的特殊现象并且存在特殊的通往混沌的路径,下面分2小节对这2种行为进行分析.

4.1 周期突变现象

值得注意的是,当参数$q = 0.666,0.682,0.698,0.734$$0.746$时,切换系统虽然均表现为周期振荡,但在这些参数值处,其相应的分岔图上对应的$z$值却出现了不连续现象,这说明在这些参数的相邻参数$q$两区间中,其对应的周期振荡的形式不一样,也即在这些参数处周期振荡会产生突变现象.现以$q = 0.666$$q = 0.698$为例,说明其周期振荡的突变行为.

图 6分别给出了$q = 0.650$$q = 0.670$,即突变点$q =0.666$两边的不同周期振荡行为. 其中轨迹从$A_1$点出发,受子系统1控制.显然两周期运动的轨迹可以分为4个部分,亦即$A_1 B_1 $,$B_1 A_2 $,$A_2 B_2 $$B_2 A_1 $,其中$A_1 B_1 $$A_2 B_2$受子系统1控制,而$B_1 A_2 $$B_2 A_1$受子系统2控制,其振荡轨迹的运动可以详细描述如下:当轨迹运行到位于分界面$\varSigma_1 $上的$A_1$点时,系统受子系统1控制,此时子系统1存在着稳定的周期振荡,也即围绕平衡点$E_1$的周期吸引子(见图 4(b)). 显然,轨迹表现为以$A_1$点为初值,向稳定的极限环逼近的暂态过程.在这过程中,其振荡幅值不断增加,当运行到$B_1$点时,轨迹穿越分界面$\varSigma_2$,使得系统转向受子系统2控制.由图 2(a)中的分岔图可知,此时子系统2也存在着围绕平衡点$E_+$的稳定极限环,因而,轨迹将以$B_1$为初值,向该极限环逐渐逼近.当运动到$A_2$点时,轨迹再次穿越分界面$\varSigma_1$,使得其转向受子系统1控制,当轨迹运行到$B_2$点处,穿越分界面$\varSigma_2$,导致系统受子系统2控制.当轨迹回到出发点$A$时,完成一个周期的运动.

图 6 突变点$q = 0.666$两相邻参数域内的不同周期运动 Fig. 6 The suddenly changed different periodic movements for $q =0.666$

图 6可以看出,两周期振荡的轨迹存在着很大的差异,导致这种差异的原因在于参数$q$的不同,使得轨迹在两子系统之间切换的具体位置会存在着一定的区别.虽然从$A_1$点出发的轨迹受子系统1的控制,到达$B_1$点的过程不受$q$变化的影响(图 6$A_1B_1$一段),但当轨迹从$B_1$点出发受子系统2控制时,由于$q$值的不同,子系统2存在的稳定极限环也不同,$q=0.650$时的稳定极限环振动幅值要略小于$q=0.670$时的稳定极限环,从而导致从相同的$B_1$出发的轨迹在子系统2控制下到达截面$\varSigma_1$ 时的$A_2$点位置不同,进一步计算可以发现,在$q=0.650$时,$A_2$点到其相应的极限环的距离为$d_1=0.342$,而在$q=0.670$ 时,$A_2$点到其相应的极限环的距离为$d_1=0.0431$${d_1} \gg {d_2}$,导致在各自的子系统2控制下,在$q=0.650$时,轨迹趋于其相应的极限环存在较长的暂态过程,才能到达截面$\varSigma_2$,而在$q=0.670$时,由于$A_2$点离其相应的稳定极限环很近,轨迹很快就到达截面$\varSigma_2$,从而造成两周期振荡存在着周期突变现象. 进一步计算可得,$q=0.650$时其周期为$T_1=41.077$,而$q=0.670$时其相应的周期为$T_1=32.725$.

同样在突变点$q=0.698$ 两侧,也会存在着不同的周期振荡.图 7分别给出了$q=0.695$$q=0.710$ 时的周期振荡.显然图 7(a)中的周期轨迹可以分为4部分,而图 7(c)中的周期轨迹则可以分为2部分,分别受子系统1和2交替控制.

依然从相同的$A_1$点出发,系统受子系统1控制,在到达位于截面$\varSigma_2$ 上的$B_1$点时,轨迹经历了相同的过程(图 7(b)图 7(d)中的$A_1B_1$ 段). 从$B_1$点出发,轨迹将分别受$q=0.695$$q=0.710$ 时相应的子系统2控制.由于在不同的$q$ 值时子系统2对应于不同的稳定极限环,在$q=0.710$$B_1$出发的轨迹将很快到达截面$\varSigma_1$上的$A_1$点,完成了一个周期的振荡,而在$q=0.695$ 时,$B_1$出发的轨迹也会较快到达截面$\varSigma_1$上的与出发点$A_1$不同的$A_2$点,然后受子系统1控制,经历较长的振荡过程,到达截面$\varSigma_2$上的$B_2$点,再受子系统2控制,到达截面$\varSigma_1$上的$A_1$点,完成一个周期的振荡.造成这一差别的原因依然在于切换点的具体位置和相应极限环的不同.在$q=0.710$时,两子系统的稳定极限环非常接近,从而导致只存在两个切换点,而在$q=0.695$时,其相应的稳定极限环与子系统1的稳定极限环存在较大的差距,同时,经过子系统2控制后的$A_2$点与子系统1的稳定极限环之间的距离较大,${d_{\rm{s}}}$大致等于 $0.126$,导致在子系统1的控制下,轨迹需要很长的暂态过程才能到达截面$\varSigma_2$上的$B_2$点,从而产生周期突变现象,进一步计算可得,$q=0.695$时振荡的周期为$T_3=24.790$,而$q=0.710$时振荡的周期为$T_4=8.044$.

图 7 突变点$q = 0.695$两相邻参数域内的不同周期运动 Fig. 7 The suddenly changed different periodic movements\for $q =0.695$
4.2 混沌道路

稳定的周期振荡在$q=0.775$会由倍周期分岔进入混沌,图 8(a)为切换系统$q=0.775$附近的分岔图,图 8(b)$\sim$图 8(d)分别给出了混沌道路上的周期二和周期四窗口以及混沌振荡,随着$q$的增加,在两截面上的切换点的数目会成倍地增加,从而导致倍周期分岔序列.

图 8 周期直接分岔为混沌 Fig. 8 Period-doubling bifurcation to chaos

必须指出的是,切换系统的倍周期分岔序列与光滑系统中的倍周期分岔序列不同,因为在光滑系统中只存在单个系统,而在切换系统中,轨迹的运动涉及到两个不同的子系统,而且,其中的一个子系统随参数$q$的变化而变化,从而导致其倍周期分岔序列不是严格的周期成倍增加,而是切换点数目的成倍增加,如图 8中周期二窗口周期$T_5=22.497$,而周期四窗口的周期$T_6=60.539$.

5 结论

与固定时间切换和单状态切换系统不同,双临界状态切换系统存在着更为复杂的非线性现象,如周期突变现象,即在突变点参数的两侧,虽然同样是周期运动,但其运动模式和相应的周期均会存在很大的差异,导致这一现象的主要原因在于切换点的位置及相应的子系统的吸引子存在明显的区别.同时,系统也存在着倍周期分岔序列导致混沌的道路,然而,与光滑系统中的倍周期分岔序列不同,切换系统的倍周期分岔序列只对应于切换点数目的成倍增加,而其相应的周期一般不对应于严格的周期倍化过程.

参考文献
[1] Pérez C, Bensítez F. Switched convergence of second-order switch nonlinear systems[J]. Journal of Control, Automation, and Systems,2012, 10 (5) : 920-930. DOI: 10.1007/s12555-012-0508-0. (0)
[2] 王仁明, 关治洪, 刘新芝. 一类非线性切换系统的稳定性分析[J]. 系统工程与电子技术,2004, 26 (1) : 68-71. ( Wang Renming, Guan Zhihong, Liu Xinzhi. Stability analysis for a class of nonlinear switched systems[J]. Systems Engineering and Electronics,2004, 26 (1) : 68-71. (in Chinese) ) (0)
[3] Davila J, Pisano A, Usai E. Continuous and discrete state reconstruction for nonlinear switched system via high-order sliding-mode observers[J]. International Journal of System Science,2011, 42 (5) : 725-735. DOI: 10.1080/00207721.2010.518254. (0)
[4] 刘正凡, 蔡晨晓, 段文勇, 等. 不确定切换时滞非线性系统状态切换的指数稳定性[J]. 控制与决策,2014, 29 (12) : 2247-2252. ( Liu Zhengfan, Cai Chenxiao, Duan Wenyong, et al. Exponential stability of uncertain switched nonlinear time-delay systems with statedependent switching[J]. Control and Decision,2014, 29 (12) : 2247-2252. (in Chinese) ) (0)
[5] Oishi M, Mitchell I, Bayen AM, et al. Hybrid verification of an interface for an automatic landing[J]. Proceeding of the IEEE Conference on Decision and Control,2002, 2 : 1607-1613. DOI: 10.1109/CDC.2002.1184749. (0)
[6] Bishop BE, Spong MW. Control of redundant manipulators using logic-based switching[J]. Proceedings of the 36th IEEE Conference on Decision and Control,1998, 2 : 16-18. (0)
[7] Horowitz R, Varaiya P. Control design of an automated highway system[J]. Proceedings of IEEE,2000, 88 (7) : 913-925. DOI: 10.1109/5.871301. (0)
[8] 方志明. 切换系统稳定性分析与优化控制若干问题研究. [博士论文]. 南京: 南京理工大学, 2012 ( Fang Zhiming. Study on several problems of stability and optimal control of switched systems. [PhD Thesis]. Nanjing: Nanjing University of Science & Technology, 2012 (in Chinese) ) (0)
[9] Goncalves JM. Constructive global analysis of hybrid systems. [PhD thesis]. MA, USA:Massachusetts Institute of Technology, 2000 (0)
[10] Zhang JF, Han ZZ, Zhu FB, et al. Absolute exponential stability and stabilization of switched nonlinear systems[J]. Systems & Control Letters,2014, 66 (1) : 51-57. (0)
[11] Xiang WM, Xiao J, Iqba MN. Robust fault detection for a class of uncertain switched nonlinear systems via the state updating approach[J]. Nonlinear Analysis: Hybrid Systems,2014, 12 (2) : 132-146. (0)
[12] Jin Y, Fu J, Zhang YM, et al. Reliable control of a class of switched cascade nonlinear systems with its application to flight control[J]. Nonlinear Analysis: Hybrid Systems,2014, 11 (1) : 11-21. (0)
[13] Yu Aleksandrov A, Chen Y, Platonov AV, et al. Stability analysis for a class of switched nonlinear systems[J]. Automatica,2011, 47 (10) : 2286-2291. DOI: 10.1016/j.automatica.2011.08.016. (0)
[14] Moulay E, Bourdais R, Perruquetti W. Stabilization of nonlinear switched systems using control Lyapunov functions[J]. Nonlinear Analysis: Hybrid Systems,2007, 1 (4) : 482-490. DOI: 10.1016/j.nahs.2005.12.001. (0)
[15] Liu YY, Zhao J. Stabilization of switched nonlinear systems with passive and non-passive subsystems[J]. Nonlinear Dyn,2012, 67 (3) : 1709-1716. DOI: 10.1007/s11071-011-0098-z. (0)
[16] Colaneria P, Geromelb JC, Astolfic A. Stabilization of continuoustime switched nonlinear systems[J]. Systems & Control,2008, 57 (1) : 95-103. (0)
[17] Alessandria A, Sanguineti M. Connections between Lp stability and asymptotic stability of nonlinear switched systems[J]. Nonlinear Analysis: Hybrid Systems,2007, 1 (1) : 501-509. (0)
[18] Hajiahmadi M, De Schutter B, Hellendoorn H. Stabilization and robust H1 control for sector-bounded switched nonlinear systems[J]. Automatica,2014, 50 (10) : 2726-2731. DOI: 10.1016/j.automatica.2014.08.015. (0)
[19] Han TT, Ge SS, Lee TH. Adaptive neural control for a class of switched nonlinear systems[J]. Systems & Control Letters,2009, 58 (2) : 109-118. (0)
[20] Elfarra NH, Mhaskar P, Christofider PD. Output feedback control of switched nonlinear systems using multiple Lyapunov functions[J]. Systems & Control Letters,2005, 54 (12) : 1163-1182. (0)
[21] Chiang ML, Fu LC. Rubust output feedback stabilization of switched nonlinear systems with average dwell time[J]. Asian Journal of Control,2014, 16 (1) : 264-276. DOI: 10.1002/asjc.2014.16.issue-1. (0)
[22] Jouili K, Braiek NB. Stabilization of non-minimum phase switched nonlinear systems with the concept of multi-diffeomorphism[J]. Commun Nonlinear Sci Numer Simulat,2015, 23 (1-3) : 282-293. DOI: 10.1016/j.cnsns.2014.11.018. (0)
[23] Mojica-Nava E, Quijano N, Rakoto-Ravalontsalama N. A polynomial approach for optimal control of switched nonlinear systems[J]. Robust Nonlinear Control,2014, 24 (12) : 1797-1808. DOI: 10.1002/rnc.v24.12. (0)
[24] Lin Q, Loxton R, Teo KL. Optimal control of nonlinear switched systems: computational methods and applications[J]. Journal of the Operations Research Society of China,2013, 1 (3) : 275-311. DOI: 10.1007/s40305-013-0021-z. (0)
[25] 陈章耀, 雪增红, 张春, 等. 周期切换下Rayleigh 振子的振荡行为及机理[J]. 物理学报,2014, 63 (1) : 1-8. ( Chen Zhangyao, Xue Zenghong, Zhang Chun, et al. Oscillation behaviors and mechanism of Rayleigh oscillator with periodic switches[J]. Acta Phys. Sin,2014, 63 (1) : 1-8. (in Chinese) ) (0)
[26] Jerome JW. Time dependent closed quantum systems: nonlinear Kohn–Sham potential operators and weak solutions[J]. Journal of Mathematical Analysis and Applications,2015, 429 (2) : 995-1006. DOI: 10.1016/j.jmaa.2015.04.047. (0)
[27] Yang JQ, Chen YT, Zhu FL, et al. Synchronous switching observer for nonlinear switched systems with minimum dwell time constraint[J]. Journal of the Franklin Institute,2015, 352 (11) : 4665-4681. DOI: 10.1016/j.jfranklin.2015.07.010. (0)
[28] Ma RC, Zhao SZ, Wang M. Global robust stabilisation of a class of uncertain switched nonlinear systems with dwell time specifications[J]. International Journal of Control,2014, 87 (3) : 589-599. DOI: 10.1080/00207179.2013.852247. (0)
[29] Kousaka T, Ueta T, Kawakami H. Bifurcation of switched nonlinear dynamical systems[J]. IEEE Trans. Circuits Syst,1999, 46 (7) : 878-885. DOI: 10.1109/82.775383. (0)
[30] Zhang C, Bi QS, Han XJ, et al. On two-parameter bifurcation analysis of switched system composed of Duffing and van der Pol oscillators[J]. Communications in Nonlinear Science and Numerical Simulation,2014, 19 (3) : 750-757. DOI: 10.1016/j.cnsns.2013.06.028. (0)
[31] Nishiuchi Y, Ueta T, Kawakami H. Stable torus and its bifurcation phenomena in a simple three-dimensional autonomous circuit[J]. Chaos, Solitons and Fractals,2006, 27 (4) : 941-951. DOI: 10.1016/j.chaos.2005.04.092. (0)
COMPLICATED BEHAVIORS AS WELL AS THE MECHANISM IN BVP OSCILLATOR WITH SWITCHES RELATED TO TWO STATES
Chen Zhangyao, Wang Yaming, Zhang Chun, Bi Qinsheng     
Faculty of Civil Engineering and Mechanics, Jiangsu University, Zhenjiang 212013, China
Abstract: The dynamics of nonlinear switched systems which possess wide engineering background and cannot be explored directly by traditional nonlinear theory, become one of hot and frontier tasks at home and abroad for the time being. The complicated behaviors as well as the mechanism of the vector field alternated between two subsystems by two different critical states are investigated in this paper. Upon employing the typical generalized BVP oscillator as an example, by introducing bilateral switch, the nonlinear dynamical model alternated between two subsystems related two states is established, the different movement forms as well as the dynamical evolution of which caused by switches are explored in details. Based on the Poincaré theory of nonlinear system, the computational equation of Lyapunov exponents of switched system is derived. Combined with the bifurcation analysis of subsystems, different oscillations of the system are discussed, upon which the nonlinear behaviors such as sudden changes of period in periodic oscillations and the route to chaos with period-doubling bifurcations as well as the related essence are presented. Different from the systems with fixed time or single state switch, much more nonlinear phenomena may be observed in the dynamic systems with two state switches in which there may exist more switch points with changeable positions. Furthermore, different from the cascading of period-doubling bifurcations in smooth systems, the period-doubling bifurcations in switched systems correspond to the doubling of the number of switch points, which usually does not correspond to the doubling of the real periodic length of the movements.
Key words: switched system    generalized BVP oscillator    switch related to two states    Lyapunov exponents