混沌同步编程

aachangshuai

function dx=chaos(t,x)
dx=zeros(12,1);
E=7.5;
dx(1)=-x(1)+x(2)+x(1)*x(3)+x(4)*x(2)+E;
dx(2)=-x(1)-x(2)-x(2)*x(3)+x(4)*x(1);
dx(3)=-2*x(1)^2+2*x(2)^2-x(3)+x(4);
dx(4)=-4*x(1)*x(2)-x(3)-x(4);
k=0.05;
dx(5)=-x(5)+x(6)+x(5)*x(7)+x(8)*x(6)+E-k*x(1)+k*x(2);
dx(6)=-x(5)-x(6)-x(6)*x(7)+x(8)*x(5)-k*x(1)-k*x(2);
dx(7)=-2*x(5)^2+2*x(6)^2-x(7)+x(8);
dx(8)=-4*x(5)*x(6)-x(7)-x(8);
dx(9)=-x(9)+x(10)+x(9)*x(11)+x(12)*x(10)+E-k*x(1)+k*x(2);
dx(10)=-x(9)-x(10)-x(9)*x(11)+x(12)*x(9)-k*x(1)-k*x(2);
dx(11)=-2*x(9)^2+2*x(10)^2-x(11)+x(12);
dx(12)=-4*x(9)*x(10)-x(11)-x(12);
各位好！
请问各位高手我的混沌同步方程对吗？
谢谢啦，请指教

[ 本帖最后由 无水1324 于 2009-11-30 19:24 编辑 ]

无水1324

难道你计算出了什么问题？
这样看没有什么问题好像

aachangshuai

我的同步不能算啊，不知道哪里问题，x(5)^2+x(6)^2与x(9)^2+x(10)^2同步怎么算啊

aachangshuai

就是不能同步，为什么啊

aachangshuai

请各位帮帮忙啊

aachangshuai

x(5)^2+x(6)^2与x(9)^2+x(10)^2同步怎么编写啊
能帮我看看吗？
我编的就是不能同步。
大概是k=0.3时同步
谢谢啦。

aachangshuai

function dx=chaos5(t,x)
dx=zeros(8,1);
E=7.5;
dx(1)=-x(1)+x(2)+x(1)*x(3)+x(4)*x(2)+E;
dx(2)=-x(1)-x(2)-x(2)*x(3)+x(4)*x(1);
dx(3)=-2*x(1)^2+2*x(2)^2-x(3)+x(4);
dx(4)=-4*x(1)*x(2)-x(3)-x(4);
k=0.3;
dx(5)=-x(5)+x(6)+x(5)*x(7)+x(8)*x(6)+E-k*x(1)+k*x(2);
dx(6)=-x(5)-x(6)-x(6)*x(7)+x(8)*x(5)-k*x(1)-k*x(2);
dx(7)=-2*x(5)^2+2*x(6)^2-x(7)+x(8);
dx(8)=-4*x(5)*x(6)-x(7)-x(8);

function chaos6
[T,x]=ode45('chaos5',[0 100],[0.1;0.1;0.1;0.1;0.2;0.2;0.1;0.1]);
subplot(2,1,1)
plot(T,x(:,5).^2+x(:,6).^2)
t=0:0.1:100;
a=spline(T,x(:,5).^2+x(:,6).^2,t)
for m=1:1000
    l(m,1)=a(m)
end
[T,x]=ode45('chaos9',[0 100],[0.1;0.1;0.1;0.1;0.1;0.1;0.2;0.2]);
subplot(2,1,2)
plot(T,x(:,5).^2+x(:,6).^2)
t=0:0.1:100;
b=spline(T,x(:,5).^2+x(:,6).^2,t)
for m=1:1000
    k(m,1)=b(m)
end
figure
i=1:1000
plot(0.1*i,k(i,1)-l(i,1))
xlabel('t'),ylabel('I(A_1)-I(B_1)')

aachangshuai

这是我做的，就是不能同步，帮我看看，谢谢

无水1324

真的不好意思，我对这个同步，我都不了解

hehy350901

你的驱动方式没错吧

aachangshuai

没错啊，就是调解k的值啊，就是不对啊，可是不知道错哪里啦