getPColFromNumeric_v1.m

Go to the documentation of this file.

function [PCol] = getPColFromNumeric_v1(waveform, varargin)

%Input parser
nSignals = size(waveform, 2);
if nargin == 1
    SNR = inf(1, nSignals);
elseif nargin == 2
    SNR = varargin{1};
    if isscalar(SNR)
        SNR = repmat(SNR, 1, nSignals);
    elseif (length(SNR) ~= nSignals)
        robolog('SNR must be scalar or match number of waveforms', 'ERR')
    end
else
    robolog('Too many input arguments', 'ERR')
end

%main calculation
avpow = pwr.meanpwr(waveform);
PCol = pwr(SNR(1), {avpow(1),'W'});
for jj=2:nSignals
    PCol(jj) = pwr(SNR(jj), {avpow(jj),'W'});
end

end