function x = discreteAWGN(data, N0)
  x = data + normrnd(0, sqrt(N0 / 2), 1, length(data));
end