Dynamos混沌系统的追踪控制MATLAB代码

Mahmoud 提出了一个非常复杂的含有复数向量的 Dynamos 混沌系统,该系统在电磁场方面有潜在的非常广泛应用.其动力学系统方程为

其中，x = u1 + i*u^2 , y = u3 + i*u4，分别为复数, i^2=- 1 , x_ba, y_ba分别为向量 x , y 的共轭复数，μ和α正参数。系统

(1)可以改写成如下五维实微分方程组的

形式:

关于Dynamos混沌系统的追踪控制，请具体参考论文： Dynamos混沌系统的追踪控制与同步

.pdf (502.7 KB)

1. function HunDunXiTong

2. clc 复杂

3. close all

4. global k d a b i

5. % 比例系数

6. k=[2 1 -1]';

7. % 参数d

8. d=[1 1 1];

9. % 其他参数 10.

11.

12.

13.

14.

15.

16.

17.

18.

19.

20.

21.

22.

23.

24. a=1/0.707; b=1/2.6; i=1.15; % 开始求解问题 % 后两个是alpha和beta估计值的初值 % 后两个为x的初值 x0=[0.3 0.5 0 0 0]; tspan=[0 50]; % 用户自行设置微分区间 [t,x]=ode45(@zzkz,tspan,x0); % 数据后处理 t=t'; x1=x(:,1)';x2=x(:,2)';x3=x(:,3)'; ahat=x(:,4)';bhat=x(:,5)';

25. 26. 27. 28. 29. 30. 31. 32. 33. 34. 35. 36. 37. 38. 39. 40. 41. 42. 43. 44.

45. r=mbxh(t); r1=r(1,:);r2=r(2,:);r3=r(3,:); e1=x1-k(1)*r1;e2=x2-k(2)*r2;e3=x3-k(3)*r3; % 图形显示 figure plot(t,x1,'r--',t,r1); title('t-x1,r1') figure plot(t,x2,'r--',t,r2); title('t-x2,r2') figure plot(t,x3,'r--',t,r3); title('t-x3,r3') figure plot(t,e1,'r--',t,e2,t,e3,'g:'); title('t-e1,e2,e3') figure plot(t,ahat,'r:',t,bhat); title('t-ahat,bhat') save hdxt.dat