getHermitianData.m

Go to the documentation of this file.

function [hermitianFreq, hermitianData] = getHermitianData(posFreqs,posData)
%GETNEGATIVEFREQDATA Obtains the corresponding negative frequencies for an analytic
%signal
%   Obtains the negative conjugate symmetric frequency domain data
%   when supplied with positive freq data only.
%   Note the positive frequency data should not contain the last points
%   i.e. data(1) ~= data(end), otherwise periodicity will be broken.
 
%Determine frequency resolution
deltaF = diff(posFreqs(1:2));
 
%Create negative frequency and data vectors
negFreqs = -posFreqs(end:-1:2); %Don't include DC
negData = conj(posData(end:-1:2)); %negData is conjugate of posData
 
%Remove highest frequency from positive freq and data vectors
posFreqs(end) = [];
posData(end) = [];
 
hermitianFreq = [negFreqs(:); posFreqs(:)];
hermitianData = [negData(:); posData(:)];
end