flipReal.m

Go to the documentation of this file.

function [stable_p] = flipReal(p)
%FLIPREAL Flips real part of complex number if it is positive
%  This function is for use with vector fitting. If a pole becomes unstable
%  (moves into the RHP), we make it stable by flipping it back into the LHP
 
%Create negative mask if real part positive
mask = -1*sign(real(p));
 
%Flip all positive real parts
zR = mask.*real(p);
 
%Create new stable complex poles
stable_p = zR + 1j*imag(p);
 
end