驱动和响应系统实现chen氏混沌同步

1、主函数

文件名：chen_main.m

function chen_main

% 耦合系数对同步的影响

global m n;

format long;

tspan=0:0.001:5;

Y0=[3 4 20 4 5 21];

hold on

m=0.5;n=0.5;

[t,y]=ode45(@chen,tspan,Y0);

plot(t,y(:,1)-y(:,4),'r')

legend('m=n=0.5')

2、微分函数

函数名：

代码: chen.m

function dy=chen(t,y)

format long

a=35;b=3;c=28;

% dy=zeros(3,1);

% 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);

% 同步

global m n;

u=5;

dy=zeros(6,1);

D1=funD(y(1),y(2),y(3));

D2=funD(y(4),y(5),y(6));

% 驱动系统

dy(1)=a*(y(2)-y(1))+m*0;

dy(2)=(c-a)*y(1)-y(1)*y(3)+c*y(2)+m*(D1(2,:)-D2(2,:)); dy(3)=y(1)*y(2)-b*y(3)+m*(D1(3,:)-D2(3,:));

% 响应系统

dy(4)=a*(y(5)-y(4))+n*0;

dy(5)=(c-a)*y(4)-y(4)*y(6)+c*y(5)+n*(D2(2,:)-D1(2,:)); dy(6)=y(4)*y(5)-b*y(6)+n*(D2(3,:)-D1(3,:));

3、非线性部分子函数

函数名：funD.m

代码：

function out=funD(x,y,z) c=28;u=5;

out=[0;(u-c)*y+x*z;-x*y];