home *** CD-ROM | disk | FTP | other *** search
- function PHI=phase(G)
- %PHASE Computes the phase of a complex vector
- %
- % PHI=phase(G)
- %
- % G is a complex-valued row vector and PHI is returned as its
- % phase (in radians), with an effort made to keep it continuous
- % over the pi-borders.
-
- % L. Ljung 10-2-86
- % Copyright (c) 1986-90 by the MathWorks, Inc.
- % All Rights Reserved.
-
- PHI=atan2(imag(G),real(G));
- N=length(PHI);
- DF=PHI(1:N-1)-PHI(2:N);
- I=find(abs(DF)>3.5);
- for i=I
- if i~=0,
- PHI=PHI+2*pi*sign(DF(i))*[zeros(1,i) ones(1,N-i)];
- end
- end
-
-