chua系统matlab代码

本文介绍了使用MATLAB编程实现Chua混沌系统的详细过程,包括两个关键文件:chua.m和main.m。运行代码后,将得到Chua系统混沌行为的模拟结果。
摘要由CSDN通过智能技术生成

chua系统matlab代码

chua.m

function out = chen(t,X)
global a b A B E1 ;
x=X(1);
y=X(2);
z=X(3);
if  (x<-E1)
       out=[A*(y-b*x+0.5*(a-b)*2*E1);
            x-y+z;
以下是一个简单的混合级联混沌同步的Matlab示例代码,其中包括两个Chua混沌振荡器的级联同步: ``` % Chua混沌同步示例 clear all clc % 设置参数 a = 15.6; b = 28; c = 10; h = 0.001; t = 0:h:100; % 初始化变量 x1 = zeros(1,length(t)); y1 = zeros(1,length(t)); z1 = zeros(1,length(t)); x2 = zeros(1,length(t)); y2 = zeros(1,length(t)); z2 = zeros(1,length(t)); x1(1) = -10; y1(1) = -10; z1(1) = -10; x2(1) = -8; y2(1) = -8; z2(1) = -8; % Chua混沌系统1 for i = 1:length(t)-1 x1(i+1) = x1(i) + h*a*(y1(i)-x1(i)-c*z1(i)); y1(i+1) = y1(i) + h*(x1(i)-y1(i)+z1(i)); z1(i+1) = z1(i) + h*(-b*z1(i)+(x1(i)+y1(i))/3); end % Chua混沌系统2 for i = 1:length(t)-1 x2(i+1) = x2(i) + h*a*(y2(i)-x2(i)-c*z2(i)); y2(i+1) = y2(i) + h*(x2(i)-y2(i)+z2(i)); z2(i+1) = z2(i) + h*(-b*z2(i)+(x2(i)+y2(i))/3); end % 建立级联同步 K = 0.1; % 耦合强度 for i = 1:length(t)-1 x_diff = x1(i) - x2(i); y_diff = y1(i) - y2(i); z_diff = z1(i) - z2(i); x2(i+1) = x2(i) + h*(a*(y2(i)-x2(i)-c*z2(i)) + K*x_diff); y2(i+1) = y2(i) + h*((x2(i)-y2(i)+z2(i)) + K*y_diff); z2(i+1) = z2(i) + h*(-b*z2(i)+(x2(i)+y2(i))/3 + K*z_diff); end % 画图 figure(1) plot(t,x1,'b',t,x2,'r') title('Chua混沌系统1和系统2') xlabel('时间') ylabel('x') legend('系统1','系统2') figure(2) plot3(x1,y1,z1,'b',x2,y2,z2,'r') title('Chua混沌系统1和系统2的相空间轨迹') xlabel('x') ylabel('y') zlabel('z') legend('系统1','系统2') ``` 此代码将生成两个Chua混沌系统的轨迹以及它们的相空间轨迹,其中一个系统被作为“驱动系统”来同步另一个系统。你可以通过更改参数和耦合强度来尝试不同的级联混沌同步方案。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值