switch8.m

来自「混沌跳转系统」· M 代码 · 共 29 行

function xdot=Switch8(t,x) ;
a=3; b=20; c=-20;d=10;e=-10;k=4; m=4; 
if x(3)>0 & (x(3)+sqrt(x(1)^2+x(2)^2))>k
xdot=[(a-k)*x(1)+b*x(2);
      -b*x(1)+(a-k)*x(2);
      c*x(3)+k*d];
else if x(3)<0 & (x(3)-sqrt(x(1)^2+x(2)^2))<-m
xdot=[(a-m)*x(1)+b*x(2);
      -b*x(1)+(a-m)*x(2);
      c*x(3)+m*e];
else 
xdot=[a*x(1)+b*x(2);
      -b*x(1)+a*x(2);
      c*x(3)];
end
end

% function f=Switch(t,x) ;
% a=3; b=20; c=-20;k=4; d=10; 
% if x(3)+sqrt(x(1).^2+x(2).^2)>k
% f(1)= (a-k)*x(1)+b*x(2);
% f(2)= -b*x(1)+(a-k)*x(2);
% f(3)= c*x(3)+k*d;
% else
% f(1)= (a-k)*x(1)+b*x(2);
% f(2)= -b*x(1)+(a-k)*x(2);
% f(3)= c*x(3)+k*d;
% end
% f=f(:);