% <變形耦合發電機系統中的混沌控制>,物理學報,2006年10月
clear all,clc
%options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[T,Y] = ode45(@Motor2,[0 20],[1,2,3]);
% ------------------- Choas characteristic ----------------------
% figure
% plot3(Y(:,1),Y(:,2),Y(:,3))
% figure
% plot(T,Y(:,1),'r-',T,Y(:,2),'b-.',T,Y(:,3),'k.')
% ---------------- Choas Control Using Single Variable ----------
% ---------------- Choas Control Using Double Variable ----------
% ---------------- Choas Control Using Non feedback Method 1 ---- X
% ---------------- Choas Control Using Non feedback Method 2 ---- X
figure
plot3(Y(:,1),Y(:,2),Y(:,3))
figure
plot(T,Y(:,1),'r-',T,Y(:,2),'b-.',T,Y(:,3),'g--')