output = awgn((input*2)-1, 20*log10(1/(3*sigma^2)));