Hilbert Results

function [ia,iph,ifq] = HilbertRes(s,fs)
%HilbertRes Obtains useful results from the Hilbert Transform
%   Calculates the instant amplitude, phase and frequency from
%   the Hilbert Transform of the signal provided

% Input Arguments
%   s   :  signal

% Output Arguments
%   ia   :  instant amplitude
%   iph  :  instant phase [Hz]
%   ifq  :  instant frequency

H = hilbert(s); HT = imag(H);

ia = sqrt(s.^2 + HT.^2);

ifq = angle(H(2:end).*conj(H(1:end-1)))*fs/(2*pi);

iph = 2*pi*fs*cumtrapz(ifq);

end