matlab 分数阶混沌系统的完全同步控制

1、内容简介


625-可以交流、咨询、答疑

2、内容说明

分数阶微积分这一重要的数学分支,其诞生在1695年,几乎和经典微积分同时出现。那一年,德国数学家Leibniz 和法国数学家L'Hopital 通信,探讨当导数的阶变为1/2时,其意义是什么?当时Leibniz也不知道定义与意义,只是回复道:“”这会导致悖论,终有一天将会是一个很有用的结果”。分数阶微积分狭义上主要包括分数阶微分与分数阶积分,广义上同时包括分数阶差分 [1]  与分数阶和商。由于近一些年分数阶微积分的理论成功应用到各大领域中,人们逐渐发现分数阶微积分能够刻画自然科学以及工程应用领域一些非经典现象。分数阶微积分比较热门领域包括:分数阶数值算法,分数阶同步等问题。

混沌系统是指在一个确定性系统中,存在着貌似随机的不规则运动,其行为表现为不确定性、不可重复、不可预测,这就是混沌现象。混沌是非线性动力系统的固有特性,是非线性系统普遍存在的现象。按照动力学系统的性质,混沌可以分成四种类型:时间混沌、 空间混沌、时空混沌、功能混沌。

3、仿真分析

clear,close all
q1=1.1;q2=1.1;q3=1.1;
h=0.01;N=4000;
a=0;
t(N+1)=[0];
x0=1;y0=2;z0=3;
%x0=-3.5;y0=4.2;z0=25
M1=0;M2=0;M3=0;
x(N+1)=[0];y(N+1)=[0];z(N+1)=[0];
x1(N+1)=[0];y1(N+1)=[0];z1(N+1)=[0];

x1(1)=x0+h^q1*(25*a+10)*(y0-x0)/(gamma(q1)*q1);
y1(1)=y0+h^q2*((28-35*a)*x0-x0*z0+(29*a-1)*y0)/(gamma(q2)*q2);
z1(1)=z0+h^q3*(x0*y0-(8+a)*z0/3)/(gamma(q3)*q3);

x(1)=x0+h^q1*(25*a+10)*(y1(1)-x1(1))+q1*(25*a+10)*(y0-x0)/gamma(q1+2);
y(1)=y0+h^q2*((28-35*a)*x1(1)-x1(1)*z1(1)+(29*a-1)*y1(1)+q2*((28-35*a)*x0-x0*z0+(29*a-1)*y0))/gamma(q2+2);
z(1)=z0+h^q3*(x1(1)*y1(1)-(8+a)*z1(1)/3+q3*(x0*y0-(8+a)*z0/3))/gamma(q3+2);

 

 

 

4、参考论文

  • 0
    点赞
  • 22
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
分数混沌系统是指动力系统方程中存在分数导数的混沌系统。其中,分数导数是连续时间的非整数导数。同步是指两个或多个系统的状态变量在某种条件下实现完全或部分的一致。下面是一个用MATLAB编写的两个分数混沌系统同步的程序示例。 首先,定义分数混沌系统的方程: 系统1: 带有分数导数的Lorenz系统方程 dx1/dt = a1(y1 - x1^α) dy1/dt = x1^(1-α) - x1z1 - y1 dz1/dt = x1y1 - bz1 系统2: 带有分数导数的Rössler系统方程 dx2/dt = -(y2 + z2) dy2/dt = x2 + ay2^α dz2/dt = b2 + z2(x2 - c2) 其中α为分数指数,a1,a2,b,和 c为系统参数。 下面是MATLAB程序示例: alpha = 0.8; % 分数指数 a1 = 10; a2 = 0.2; b = 8/3; c1 = 28; c2 = 14; tspan = [0 100]; x0 = [1 2 3]; % 初始状态变量 % 定义分数混沌系统函数 f1 = @(t, x) [a1*(x(2) - x(1)^alpha); x(1)^(1-alpha) - x(1)*x(3) - x(2); x(1)*x(2) - b*x(3)]; f2 = @(t, x) [-x(2) - x(3); x(1) + a2*x(2)^alpha; b + x(3)*(x(1) - c2)]; % 求解系统1的状态变量随时间的演化 [t1, x1] = ode23(f1, tspan, x0); % 求解系统2的状态变量随时间的演化 [t2, x2] = ode23(f2, tspan, x0); % 打印系统1和系统2的状态变量 disp("系统1状态变量x1的值:") disp(x1) disp("系统2状态变量x2的值:") disp(x2) % 绘制系统1和系统2的状态变量随时间的变化曲线 figure plot(t1, x1(:,1), 'r', t2, x2(:,1), 'b') xlabel('时间') ylabel('状态变量x') legend('系统1', '系统2')

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值