% <Liu混沌系統的線性反饋同步控制及電路實驗的研究>,物理學報,2006年10月
clear all,clc
%options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[T,Y] = ode45(@Liu_fun,[0 12],[0.01,0.01,0.01,0.001,0.001,0.001]);
% -------------------Choas characteristic---------------------------
figure
plot3(Y(:,1),Y(:,2),Y(:,3))
figure
plot(T,Y(:,1),'r-',T,Y(:,2),'b-.',T,Y(:,3),'g.')
% --------------------Synchronization---------------------------
figure
subplot(3,1,1)
plot(T,Y(:,1)-Y(:,4))
subplot(3,1,2)
plot(T,Y(:,2)-Y(:,5))
subplot(3,1,3)
plot(T,Y(:,3)-Y(:,6))