CN104202155B - 一种异分数阶时滞混沌系统的延时同步控制电路设计方法 - Google Patents

Abstract

一种异分数阶时滞混沌系统的延时同步控制电路设计方法，包括如下步骤：构建一个三维异分数阶混沌系统；引入时滞变量τ_i＞0(i＝1,2,3)，构建其含时滞量的动力学方程；利用分数阶链型、树型、混合型和新型4种电路单元的电路，设计(S1)和(S2)中的系统方程的组合电路；引入时间变量ξ_i(i＝1,2,3)，构造一个延时响应控制系统；设计(S4)中的延时响应控制系统的电路原理图，并通过驱动‑响应系统电路实现异分数阶时滞混沌系统间的同步控制。本发明设计了该异分数阶混沌系统的含时滞量τ_i(i＝1,2,3)的混沌系统，并对其进行了延时驱动‑响应同步控制，设计的电路实现了异分数阶时滞混沌系统的延时同步控制。

Description

一种异分数阶时滞混沌系统的延时同步控制电路设计方法

技术领域

本发明属于非线性动力学和自动控制领域，特别涉及分数阶混沌系统电路构建和延时同步控制技术。

背景技术

混沌现象是自然界中普遍存在的一种宏观无序、微观有序的非线性现象，近年来混沌理论获得了巨大而深远的发展，各种新混沌系统的提出层出不穷。分数阶微积分是研究任意阶次的微分、积分算子特性及应用的数学问题，是整数阶微积分概念的延伸和推广。由于分数阶混沌系统模型本身的复杂性，微分动力学方程中相同分数阶的混沌系统研究较多，相应的仿真电路也集中在整数阶或同分数阶设定下的实现。关于微分方程中各变量不同阶数(异阶)及其电路单元交叉的异元组合电路的仿真鲜有研究报道。

混沌系统对初始条件的极端敏感性，导致许多学者曾经一度认为两个混沌系统不可能实现同步。自从Pecora和Carroll在1990年首次用电子电路实现了两个混沌系统的同步，混沌系统的同步问题就引起了许多学者的强烈关注，各种同步控制的方法相继被提出，例如完全同步、时滞同步、相位同步、反相同步、部分同步、广义同步、脉冲同步、投影同步等。直到如今，混沌系统的同步问题依然是混沌领域的一个研究热点问题。

严格来说，任何实际系统的当前状态不可避免的受到过去状态的影响，即当前状态变化率不仅依靠当前时刻的状态，而且也依赖于过去某时刻或者某段时间的状态，系统的这种特性称之为时滞，具有时滞的系统称为时滞系统。时滞广泛存在于多种物理系统中，如振荡电路、激光、核反应、神经网络和通信网络等。时滞系统为无穷维状态空间，能够产生多于维数的正的Lyapunov指数，因此结构简单的时滞系统也可具有非常复杂的动力学行为。由于信号传输速度的限制，任何信号的传输都需要一定的时间，因而每路信号的传输时间也会不尽相同，所以响应系统的状态一定延时于驱动系统的状态，因此研究系统状态的延时同步具有极其重要的价值。

发明内容

本发明的目的是提出一种异分数阶时滞混沌系统的延时同步控制电路设计方法。

本发明包括以下步骤：

(S1)、基于分数阶理论构建一个三维异分数阶混沌系统(动力学方程)；

(S2)、在(S1)所述三维异分数阶混沌系统的动力学方程的变量中，引入时滞变量τ_i＞0(i＝1,2,3)，构建其含时滞量的动力学方程；

(S3)、异分数阶组合电路设计。采用Multisim模拟电路仿真平台，利用分数阶链型、树型、混合型和新型4种电路单元的电路，对(S1)和(S2)系统进行不同分数阶值组合电路的仿真实验，设计(S1)和(S2)中的系统方程的组合电路原理图并仿真；

(S4)、延时响应控制系统设计。以(S2)系统动力学方程为驱动系统，其中x,y,z为驱动系统的状态变量，构造一个新的延时响应控制系统，其中x',y',z'为驱动系统的状态变量，ξ_i(i＝1,2,3)为延时时间；

(S5)、延时驱动-响应同步控制分析及电路设计。分析变量x′由变量y和z驱动，变量y′由变量x和z驱动，变量z′由变量x和y驱动，对其进行数值仿真，分析设计的延时驱动-响应同步控制器能否对(S2)中的系统方程进行同步控制；

延时驱动-响应同步控制电路设计。采用Multisim电路仿真平台，设计(S4)中的延时响应控制系统的电路原理图，并通过驱动-响应系统电路实现异分数阶时滞混沌系统间的同步控制。

进一步说，本发明的具体步骤如下：

步骤1：分数阶混沌系统模型的构建。

基于分数阶理论构建一个分数阶混沌系统，其动力学方程为：

其中，0＜q_i≤1(i＝1,2,3)为系统(1)的阶数，x,y,z为状态变量，a,b,c,d,e是系统参数，并且均为实数。当a＝2,b＝3.65,c＝8,d＝3,e＝2时，系统(1)的三个Lyapunov指数分别为L₁＝0.7730，L₂＝0.00008，L₃＝-7.1231。由于该系统的三个Lyapunov指数中一个为正，一个趋近于零，一个为负，并且其和小于零，系统的Lyapunov维数：

所以系统(1)存在一个典型的混沌吸引子如图1所示。

步骤2：异分数阶时滞混沌系统的构建。

对于系统(1)，引入时滞变量τ_i(i＝1,2,3)，构建其时滞动力学方程为：

其中τ_i＞0(i＝1,2,3)为系统的时滞常数。

为了不失一般性，本发明中τ_i(i＝1,2,3)取值全不相同即异时滞系统。当q₁＝q₂＝0.95，q₃＝0.9，a＝2,b＝3.65,c＝8,d＝3,e＝2,τ₁＝0.03,τ₂＝0.05,τ₃＝0.01时，计算可知其最大Lyapunov指数为L_max＝1.2635，因此系统(3)此时处于混沌状态。

步骤3：异分数阶组合电路设计与仿真。

(1)异分数阶组合电路设计。

对于一个特定的三维分数阶混沌系统，当阶次q₁、q₂和q₃取不同情况组合时，排列的结果参见表1，一共有6+18+3＝27种排列组合方式。由于对于每一个q_i(i＝1,2,3)值均有链型、树型、混合型和新型4种电路单元选择，这样任何一种组合就有4³＝64种电路单元设计方式，从而对于任意三维的分数阶q₁、q₂和q₃取值，分数阶混沌系统组合电路方式共有64×27＝1728种。本发明中，q₁、q₂和q₃的取值可以为0.9和0.95，因此对应的选择种取值方式分别如下：

方式一：q₁＝q₂＝q₃＝0.9；q₁＝q₂＝q₃＝0.95(q_i全部相同)；

方式二：q₁＝0.9,q₂＝q₃＝0.95；q₁＝q₂＝0.9,q₃＝0.95；q₁＝0.95,q₂＝q₃＝0.9；q₁＝q₂＝0.95,q₃＝0.9；q₁＝q₃＝0.95,q₂＝0.9；q₁＝q₃＝0.9,q₂＝0.95(q_i不全相同)。

对于上述两种取值方式中的每一组异元电路组合数量均为C¹ ₄C¹ ₃C¹ ₂＝4×3×2＝24种。

表1不同q₁,q₂,q₃排列的组合数

(2)时滞电路仿真。

为了简化电路设计而又不失一般性，本发明从q_i(i＝1,2,3)不全相同组合方式中任意选择了一种q₁＝q₂＝0.95，q₃＝0.9组合方式进行多元电路仿真实验。由于各电子元件允许电压的有限性，因此，为了可靠的进行电路实验，需要先将系统的输出信号减小为原来的1/2，取a＝2,b＝3.65,c＝8,d＝3,e＝2，q₁＝q₂＝0.95，q₃＝0.9并且对应的单元电路分别为树型、新型、链型时，设计系统(1)的电路原理图并进行仿真实验，如图2所示，各变量的电路仿真相图如图3～图5所示，与图1的数值计算结果有极其相似的吻合性，从而验证电路设计的正确性；再取时滞τ₁＝0.03,τ₂＝0.05,τ₃＝0.01，根据延时单元电路如图6所示，其时滞τ的近似表达式如下：

其中n为LCL滤波器数目且n≥1，设计系统(3)的电路原理图进行仿真实验，如图7所示，各变量的电路仿真相图如图8～图10所示。

步骤4：延时驱动-响应同步控制。

(1)延时驱动-响应控制分析与设计。

以系统(3)作为驱动系统，讨论分析分析变量x′由变量y和z驱动，变量y′由变量x和z驱动，变量z′由变量x和y驱动；响应系统依次为：

其中x,y,z分别为驱动系统(3)的状态变量，x′,y′,z′依次为响应系统(5)-(7)的状态变量，ξ_i(i＝1,2,3)为延时时间。

通过对上述三组驱动-响应系统进行误差计算，结果表明只有使用变量z′由变量x和y驱动这种方法，系统才能达到同步，即选用延时响应系统(7)，才能与驱动系统(3)趋于同步，误差仿真如图11所示。

(2)延时驱动-响应同步控制电路仿真。

设计的延时响应系统(7)的电路图如图12所示，其中x，y和z′是状态变量，延时单元lag5即ξ₃的延时为0.04，延时驱动-响应同步控制电路原理图及同步仿真结果如图13和图14所示。

本发明以新构建的三维分数阶混沌系统为基础，对其进行了基本的动力学特性分析，证实了系统的混沌特性并确定其混沌吸引子的存在。利用分数阶链型、树型、混合型和新型4种电路单元的电路进行了不同分数阶值组合电路的仿真实验，组合电路数共有1728种。在新建的系统中，阶数q_i(i＝1,2,3)的取值为0.9和0.95，本发明中从q_i不全相同这类组合方式中，任意选择一种异分数阶电路(电路单元互不相同)组合进行了分析和仿真实验。实验结果表明，异分数阶电路仿真与计算机数值计算具有极高的吻合度，证实了异分数阶电路设计的有效性和灵活性，同时验证了该分数阶混沌系统在物理上的可实现性。本发明设计了该异分数阶混沌系统的含时滞量τ_i(i＝1,2,3)的混沌系统，并对其进行了延时驱动-响应同步控制，设计的电路实现了异分数阶时滞混沌系统的延时同步控制。

附图说明

图1为本发明异分数阶混沌系统的各空间混沌吸引子相图。

图2为本发明异分数阶混沌系统的组合电路原理图。

图3为本发明异分数阶混沌系统的组合电路仿真x-y相图。

图4为本发明异分数阶混沌系统的组合电路仿真x-z相图。

图5为本发明异分数阶混沌系统的组合电路仿真y-z相图。

图6为本发明时滞单元电路图。

图7为本发明异分数阶时滞混沌系统的组合电路原理图。

图8为本发明异分数阶时滞混沌系统组合电路仿真x-y相图。

图9为本发明异分数阶时滞混沌系统组合电路仿真x-z相图。

图10为本发明异分数阶时滞混沌系统组合电路仿真y-z相图。

图11为本发明延时驱动-响应同步控制误差曲线图。

图12为本发明延时驱动-响应控制系统电路图。

图13为本发明延时驱动-响应同步控制电路原理图。

图14为本发明延时同步控制电路原理图状态变量z-z′投影曲线。

具体实施方式

以下将结合附图对本发明作进一步的详细描述。将通过以下实施例作进一步说明。

实施例1。设计实现异分数阶混沌系统(1)的组合电路。

为了简化电路设计而又不失一般性，本发明从q_i(i＝1,2,3)不全相同组合方式中任意选择了一种q₁＝q₂＝0.95，q₃＝0.9组合方式进行多元电路仿真实验。由于各电子元件允许电压的有限性，因此，为了可靠的进行电路实验，需要先将系统的输出信号减小为原来的1/2，取q₁＝q₂＝0.95，q₃＝0.9并且对应的单元电路分别为树型、新型、链型时，其Multisim电路原理图如图2所示。根据系统电路原理图及电路基本理论，可得系统的数学方程如(8)式所示。

将方程(8)和(1)进行比较，可得：

令C₁＝C₂＝C₃＝33nF，R_f1＝R_f2＝R_f3＝100kΩ，R₃＝R₆＝R₉＝50kΩ，R₁₂＝166.7kΩ，R₁＝R₂＝R₄＝R₅＝R₇＝R₈＝10kΩ，R₁₁＝R₁₃＝R₂₂＝R₂₃＝R₃₂＝50kΩ，R₂₁＝27.4kΩ，R₃₁＝12.5kΩ时，用Multisim对该电路方程进行了仿真实验，各变量相图仿真结果如图3～图5所示。与图1比较可以看出，电路仿真实验结果与数值计算结果十分吻合，故该异分数阶组合混沌系统电路是可以物理实现的。为了确保异元电路设计的有效性，对方式二中的其余5种取值分别也进行了异元电路仿真，仿真结果与使用基于广义的Adams-Bashforth-Moulton算法的Matlab数值计算高度吻合，进一步说明了该设计思想的有效性。

实施例2。设计实现异分数阶时滞混沌系统(3)的组合电路。

本发明所采用时滞单元电路如图6所示，时滞τ的近似表达式如(4)式所示。由于低通滤波器网络会受到信号频率的限制，而该时滞单元电路在截止频率f_c＝1kHz以下时具有平稳特性。当噪声频率和有用信号频率接近时，单级滤波器无法达到预期效果。此时需要选用多级滤波器来避免噪声干扰，在输入与输出之间设n＝10组T型滤波器，端口处配置匹配电阻R₂₈＝R₃₀＝1kΩ，且通带内的特性阻抗为常数。取R₂₆＝R₂₇＝R₂₉＝10kΩ，R₃₃＝22kΩ。

对于系统(3)，当延时单元lagi(i＝1,2,3)的对应时滞τ_i(i＝1,2,3)的数值分别为τ₁＝0.03,τ₂＝0.05,τ₃＝0.01时，由(4)式可以算出各自时滞单元电路中LC的数值分别为1mH，4.5nF；1mH，12.5nF；1mH，0.5nF。当q₁＝q₂＝0.95，q₃＝0.9，a＝2,b＝3.65,c＝8,d＝3,e＝2,τ₁＝0.03,τ₂＝0.05,τ₃＝0.01时，系统(3)的电路原理图如图7所示，各变量相图仿真结果如图8～图10所示，电路仿真结果表明此异分数阶时滞混沌系统在物理上是可以实现的。

实施例3。设计实现延时驱动-响应同步控制。

对于驱动系统(3)和响应系统(7)，当初值选为[x,y,z,z']＝[2,-6,0,10]，延时时间δ₃＝0.04时，matlab数值仿真同步误差曲线如图11所示。延时响应系统(7)的电路如图12所示，其中x，y和z′是状态变量；延时单元lag5即δ₃的延时为0.04，由(4)式可以算出时滞单元电路中LC的数值为1mH，16nF；R₅₀＝12.5KΩ，R₅₁＝R₆₀＝50KΩ，R₆₁＝R₆₂＝10KΩ，R_f4＝100KΩ。延时驱动-响应同步控制电路原理图如图13所示，z-z'同步仿真结果如图14所示，仿真结果表明此延时驱动-响应同步控制电路设计的有效性。

Claims (1)

1.一种异分数阶时滞混沌系统的延时同步控制电路设计方法，其特征是包括如下步骤：

S1：基于分数阶理论构建一个三维异分数阶混沌系统，其动力学方程为：

$\left\{\begin{array}{c}\frac{d^{q_1}x}{{\mathrm{dt}}^{q_1}}=-ax+dyz+y^2\\ \frac{d^{q_2}y}{{\mathrm{dt}}^{q_2}}=by-ez-xz\\ \frac{d^{q_3}z}{{\mathrm{dt}}^{q_3}}=-cz+xy\end{array}\right.$

其中，0＜q_i≤1,i＝1,2,3为系统的阶数，x,y,z为状态变量，a,b,c,d,e是系统实数参数；

S2：在S1所述三维异分数阶混沌系统的动力学方程的变量中，引入时滞变量τ_i＞0,i＝1,2,3，构建其含时滞量的动力学方程：

$\left\{\begin{array}{c}\frac{d^{q_1}x}{{\mathrm{dt}}^{q_1}}=-ax(t-{\mathrm{\τ}}_1)+dy\left(t\right)z\left(t\right)+y^2\left(t\right)\\ \frac{d^{q_2}y}{{\mathrm{dt}}^{q_2}}=by(t-{\mathrm{\τ}}_2)-ez\left(t\right)-x\left(t\right)z\left(t\right)\\ \frac{d^{q_3}z}{{\mathrm{dt}}^{q_3}}=-cz(t-{\mathrm{\τ}}_3)+x\left(t\right)y\left(t\right)\end{array}\right.$

其中τ_i＞0,i＝1,2,3为系统的时滞常数；

S3：采用Multisim模拟电路仿真平台，利用分数阶链型、树型、混合型和新型4种电路单元的电路，对S1和S2中系统进行不同分数阶值组合电路的仿真实验，设计S1和S2中的系统方程的组合电路原理图；

S4：以S2中系统动力学方程为驱动系统，引入时间变量ξ_i,i＝1,2,3，构造一个延时响应控制系统：

$\frac{d^{q_1}{x}^{\′}\left(t\right)}{{\mathrm{dt}}^{q_1}}=-{\mathrm{ax}}^{\′}(t-{\mathrm{\ξ}}_1)+dy\left(t\right)z\left(t\right)+y^2\left(t\right)$

$\frac{d^{q_2}{y}^{\′}\left(t\right)}{{\mathrm{dt}}^{q_2}}={\mathrm{by}}^{\′}(t-{\mathrm{\ξ}}_2)-ez\left(t\right)-x\left(t\right)z\left(t\right)$

$\frac{d^{q_3}{z}^{\′}\left(t\right)}{{\mathrm{dt}}^{q_3}}=-{\mathrm{cz}}^{\′}(t-{\mathrm{\ξ}}_3)+x\left(t\right)y\left(t\right)$

其中x,y,z分别为S2中驱动系统的状态变量，x′,y′,z′依次为响应系统的状态变量，ξ_i,i＝1,2,3为延时时间；

S5：采用Multisim电路仿真平台，设计S4中的延时响应控制系统的电路原理图，并通过驱动-响应系统电路实现异分数阶时滞混沌系统间的同步控制。