motor.m

根据《物理学报》2006年第10期上的4篇论文

function dy = Motor(t,y)
dy = zeros(7,1);
mu=1;alfa=1.9;
dy(1) = -mu*y(1) + y(2)*(y(3)+alfa);
dy(2) = -mu*y(2) + y(1)*(y(3)-alfa);
% --------------------Choas characteristic---------------------------
% ******************** Note the size of dy **************************
% dy(3) = 1-y(1)*y(2);
% -----------------------Choas Control--------------------------------
% ******************** Note the size of dy **************************
% k=2;z0=3;
% dy(3) = 1-y(1)*y(2) - k*(y(3)-z0);
% ----------------------Synchronization-------------------------------
% ******************** Note the size of dy **************************
belta=5;
dy(3) = 1-y(1)*y(2);
dy(4) = -mu*y(4) + y(5)*(y(6)+alfa) - y(7)*(y(4)-y(1));
dy(5) = -mu*y(5) + y(4)*(y(6)-alfa);
dy(6) = 1-y(4)*y(5);
dy(7)=belta*(y(4)-y(1))*(y(4)-y(1));