function PHI= phaseM(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 1986-2004 The MathWorks, Inc.

%PHI = unwrap(angle(G));
[nr,nc] = size(G);
if min(nr,nc) > 1
    error(sprintf(['PHASE applies only to row or column vectors.'...
        '\nFor matrices you have to decide along which dimension the'...
        '\nphase should be continuous.']))
end
if nr>nc
    G = G.';
end
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 > N/6  % the first sixth without changes
        if i~=0,
            PHI=PHI+2*pi*sign(DF(i))*[zeros(1,i) ones(1,N-i)];
        end
    end
end
if nr>nc
    PHI = PHI.';
end