声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1966|回复: 10

[分形与混沌] 混沌同步编程

[复制链接]
发表于 2009-11-29 11:57 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
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 编辑 ]
回复
分享到:

使用道具 举报

发表于 2009-11-30 19:24 | 显示全部楼层
难道你计算出了什么问题?
这样看没有什么问题好像
 楼主| 发表于 2009-12-8 13:19 | 显示全部楼层

谢谢

我的同步不能算啊,不知道哪里问题,x(5)^2+x(6)^2与x(9)^2+x(10)^2同步怎么算啊
 楼主| 发表于 2009-12-10 11:06 | 显示全部楼层
就是不能同步,为什么啊
 楼主| 发表于 2009-12-10 11:07 | 显示全部楼层
请各位帮帮忙啊
 楼主| 发表于 2009-12-13 10:58 | 显示全部楼层

回复 沙发 无水1324 的帖子

x(5)^2+x(6)^2与x(9)^2+x(10)^2同步怎么编写啊
能帮我看看吗?
我编的就是不能同步。
大概是k=0.3时同步
谢谢啦。
 楼主| 发表于 2009-12-14 08:16 | 显示全部楼层

回复 沙发 无水1324 的帖子

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)')
 楼主| 发表于 2009-12-14 08:17 | 显示全部楼层

回复 沙发 无水1324 的帖子

这是我做的,就是不能同步,帮我看看,谢谢
发表于 2009-12-16 16:50 | 显示全部楼层

回复 6楼 aachangshuai 的帖子

真的不好意思,我对这个同步,我都不了解
发表于 2009-12-24 18:08 | 显示全部楼层

怀疑方程

你的驱动方式没错吧
 楼主| 发表于 2010-3-14 13:58 | 显示全部楼层

回复 10楼 hehy350901 的帖子

没错啊,就是调解k的值啊,就是不对啊,可是不知道错哪里啦
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-25 03:57 , Processed in 0.090749 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表