motor_sync.m

来自「根据《物理学报》2006年第10期上的4篇论文」· M 代码 · 共 25 行

% <耦合发电机系统的自适应控制与同步>,物理学报,2006年10月
clear all,clc
%options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[T,Y] = ode45(@Motor,[0 20],[1,2,3,2,3,1,0]);
% -------------------Choas characteristic----------------------
% figure
% plot3(Y(:,1),Y(:,2),Y(:,3))
% figure
% plot(T,Y(:,1),'r-',T,Y(:,2),'b-.',T,Y(:,3),'g.')
% -------------------Choas Control-----------------------------
% figure
% % plot3(Y(:,1),Y(:,2),Y(:,3))
% figure
% plot(T,Y(:,1),'r-',T,Y(:,2),'b-.',T,Y(:,3),'g--')
% --------------------Synchronization---------------------------- X
% ****** Note: [1,2,3] instead by [1,2,3,2,3,1,0] in ode45() ******
% figure
% plot3(Y(:,4),Y(:,5),Y(:,6))
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))