SNRScript

function [x, calcSNR] = ourSNR(signal, Tsymbol, fs)

y = signal(1:Tsymbol*fs);           % First symbol of the signal - test sample
N = length(y);
Y = fft(y);
Yabs = abs(Y);
YPower = Yabs.^2;
YdB = 10*log10(YPower);             % Converting to dB

peak = max(YdB);                    % Find peak in dB
noise = mean(YdB);                  % Finding noise in dB

x = peak - noise - 10*log10(N/2);   % Calculate SNR

calcSNR = snr(y);                   % To test SNR function

end