求助：谁有二阶锁相环相图程序
帮忙提供以下，谢谢！相图形式见附件。我编的程序如下：子程序
function dy=ODEfun2(x,y,k,tau1,tau2)
dy = (-k*tau2/tau1)*cos(x)-(k/tau1)*sin(x)/y;
retur
主程序
tau1=4;tau2=10;k=12;
w【-4，4];
for i=1:17
[T,Y]=ode45('ODEfun22',[-pi,pi],w(i));
plot(T,Y);hold on;grid on;
end
出错，请哪位高手指点一下，感激不尽。
相图.rar
(52.42 KB)
