% snr (Signal-to-noise ratio function) by Dan Piehl (2006)
%
% This will determine the signal-to-noise ratio given both
% a signal contaminated with noise, and the noiseless signal
%

function b = snr(noisya,noiselessa)

  sq_noise=sum(sum(sum(abs(noisya-noiselessa).^2)));
  sq_signal=sum(sum(sum(abs(noiselessa).^2)));
  b=sqrt(sq_signal/sq_noise);