一个Chen系统的激活控制同步——MATLAB实现

最新推荐文章于 2023-10-08 09:40:58 发布

MATLAB机器学习

最新推荐文章于 2023-10-08 09:40:58 发布

阅读量624

点赞数 5

文章标签： matlab 矩阵线性代数控制器同步

本文链接：https://blog.csdn.net/weixin_44457137/article/details/124529949

版权

对于 $C h e n$ 系统，驱动系统定义为：
$\left\{ \begin{array}{l} \dot{x}=a\left( y-x \right)\\ \dot{y}=\left( c-a \right) x-xz+cy\\ \dot{z}=xy-bz\\ \end{array} \right.$
响应系统为：
$\left\{ \begin{array}{l} \dot{x}_1=a\left( y_1-x_1 \right) +\mu _1\left( t \right)\\ \dot{y}_1=\left( c-a \right) x_1-x_1z_1+cy_1+\mu _2\left( t \right)\\ \dot{z}_1=x_1y_1-bz_1+\mu _3\left( t \right)\\ \end{array} \right.$
其中参数 $a = 35, b = 3, c = 28$ ,且 $\mu _i\left( t \right) \left( i=1,2,3 \right)$
为控制函数。
这里我们定义误差变量为： $e_1=x_1-x,e_2=y_1-y,e_3=z_1-z$
则可以得到如下误差系统：
$\left\{ \begin{array}{l} \dot{e}_1=a\left( e_2-e_1 \right) +\mu _1\left( t \right)\\ \dot{e}_2=\left( c-a \right) e_1-x_1e_3-z_1e_1+e_1e_3+\mu _2\left( t \right)\\ \dot{e}_3=x_1e_2+y_1e_1-e_1e_2-be_3+\mu _3\left( t \right)\\ \end{array} \right.$
定义激活控制函数如下：
$\left\{ \begin{array}{l} \mu _1\left( t \right) =V_1\left( t \right)\\ \mu _2\left( t \right) =x_1e_3+z_1e_1-e_1e_3+V_2\left( t \right)\\ \mu _3\left( t \right) =-x_1e_2-y_1e_1+e_1e_2+V_3\left( t \right)\\ \end{array} \right.$
这里的控制输入 $V_i\left( t \right) \left( i=1,2,3 \right)$
选择如下：
$\left( \begin{array}{c} V_1\left( t \right)\\ V_2\left( t \right)\\ V_3\left( t \right)\\ \end{array} \right) =A\left( \begin{array}{c} e_1\\ e_2\\ e_3\\ \end{array} \right)$
其中：
$A=\left( \begin{matrix} a-1& -a& 0\\ a-c& -\left( 1+c \right)& 0\\ 0& 0& b-1\\ \end{matrix} \right)$
由于闭环系统的特征根为 $- 1, - 1, - 1$ ，因此当 $t\rightarrow \infty$ 时，误差变量 $e_1,e_2,e_3$
收敛于零，这就意味着两个 $C h e n$ 系统能够实现混沌同步。
其混沌吸引子可表现如下：
在这里插入图片描述实现上述过程的MATLAB代码为：

clear,clc;
X0 = [0.5,0.3,0.78];
tspan = [0,100];
[T,Y]=ode45(@Chen,tspan,X0);
plot3(Y(100:end,1),Y(100:end,2),Y(100:end,3),'color',[0.3,0.2,0.7]);
xlabel('x','FontName','Times New Roman','FontSize',24);
ylabel('y','FontName','Times New Roman','FontSize',24);
zlabel('z','FontName','Times New Roman','FontSize',24);
function dy = Chen(t,y)
a = 35;b = 3;c = 28;
dy = zeros(3,1);   
 % a column vector
dy(1) = a*(y(2) - y(1));
dy(2) = (c-a)*y(1)-y(1)*y(3)+c*y(2);
dy(3) = y(1)*y(2)-b*y(3);
end

其一维时间序列图如下图所示：
在这里插入图片描述

实现上述过程的MATLAB代码如下：

function ChenSynchronization
clear,clc;
Y0 = [0.5,0.3,0.78,0.9,-1,3];
tspan = [0,100];
[t,Y]=ode45(@Chen,tspan,Y0);
             figure(1)
   plot(t,Y(:,4)-Y(:,1),'color',[0.3,0.2,0.7])
    axis([0 50 -4 4])
    xlabel('t','FontName','Times New Roman','FontSize',20)
    ylabel('e_1','FontName','Times New Roman','FontSize',20)
    set(gca,'FontName','Times New Roman','FontSize',20,'LineWidth',2)
               figure(2)
    plot(t,Y(:,5)-Y(:,2),'color',[0.3,0.2,0.7])
     axis([0 50 -4 4])
    xlabel('t','FontName','Times New Roman','FontSize',20)
    ylabel('e_2','FontName','Times New Roman','FontSize',20)
    set(gca,'FontName','Times New Roman','FontSize',20,'LineWidth',2)
               figure(3)
    plot(t,Y(:,6)-Y(:,3),'color',[0.3,0.2,0.7])
    axis([0 50 -4 4])
    xlabel('t','FontName','Times New Roman','FontSize',20)
    ylabel('e_3','FontName','Times New Roman','FontSize',20)
    set(gca,'FontName','Times New Roman','FontSize',20,'LineWidth',2)
end
function dy = Chen(~,y)
a = 35;b = 3;c = 28;
dy = zeros(6,1); 

E = [(y(4)-y(1));(y(5)-y(2));(y(6)-y(3))];
A = [(a-1),-a,0;
    (a-c),-(1+c),0;
    0,0,b-1];
V = A*E;
U =zeros(3,1);
U(1) = V(1);
U(2) = y(4)*(y(6)-y(3))+y(6)*(y(4)-y(1))-(y(4)-y(1))*(y(6)-y(3))+V(2);
U(3) = -y(4)*(y(5)-y(2))-y(5)*(y(4)-y(1))+(y(4)-y(1))*(y(5)-y(2))+V(3);

dy(1) = a*(y(2) - y(1));
dy(2) = (c-a)*y(1)-y(1)*y(3)+c*y(2);
dy(3) = y(1)*y(2)-b*y(3);

dy(4) = a*(y(5)-y(4))+U(1);
dy(5) = (c-a)*y(4)-y(4)*y(6)+c*y(5)+U(2);
dy(6) = y(4)*y(5)-b*y(6)+U(3);
end