phaseshifting.m

open source not mine

function phase = PhaseShifting(phasein,teta)
%
% phase = PhaseShifting(phasein,teta),
% Phase shifter.
%
% inputs:
% phasein: calculated ECG phase
% teta: desired phase shift. teta>0 and teta<0 corresponds with phase leads
% and phase lags, respectively.
%
% output:
% phase: the shifted phase.
%
%
% Open Source ECG Toolbox, version 1.0, November 2006
% Released under the GNU General Public License
% Copyright (C) 2006  Reza Sameni
% Sharif University of Technology, Tehran, Iran -- LIS-INPG, Grenoble, France
% [e-mail]

% This program is free software; you can redistribute it and/or modify it
% under the terms of the GNU General Public License as published by the
% Free Software Foundation; either version 2 of the License, or (at your
% option) any later version.
% This program is distributed in the hope that it will be useful, but
% WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General
% Public License for more details. You should have received a copy of the
% GNU General Public License along with this program; if not, write to the
% Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,
% MA  02110-1301, USA.

phase = phasein + teta;

phase = mod(phase + pi, 2*pi) - pi;
% I = find(phase>pi);
% phase(I) = phase(I) - 2*pi;
% 
% I = find(phase<-pi);
% phase(I) = phase(I) + 2*pi;