equation.m

该程序是二阶琐相环的防震实例 仿真结果给出了二阶琐相环的相平面

function yp=equation(x,y);
%
%BBI 99
global K fi fo
RC=25/(2*pi*fo);   %滤波器的时间常数
yp(1)=y(2);
yp(2)=K/RC*sin(2*pi*fi*x).*cos(2*pi*po*x+y(1))-y(2)/RC;
v=version;v=v(1);
if strcmp(v,'5');yp=yp';end;