daniele@162: function x = snr(s,r)
daniele@162: % SNR calculates the signal-to-noise ratio between the signal s and its
daniele@162: % representation r defined as:
daniele@162: % SNR = 20*log10(||s||_2/||s-r||_2)
daniele@162: err = s-r;
daniele@162: x = mag2db(norm(s(:))/norm(err(:)));