funp893.m

来自「matlab算法集 matlab算法集」· M 代码 · 共 12 行

function dx = funp893 (t,x)
%-----------------------------------------------------------------------
% Description: Phase-locked loop, Example 8.9.3
%-----------------------------------------------------------------------
   global wc w0 alpha
   tau = 3*pi/wc;
   u = funu893(t);
   dx = [ (alpha*u*cos(x(2))-x(1))/tau;
          x(1) + wc ];