function [out,freqout]=AWGN_PSD(sig_square,freq)
%sig_square is the variance of the white noise process (No/2)
%freq are the frequencies to be evaluated
out=sig_square*ones(1,length(freq));
freqout=freq;