snr.m

signal procesing toolbox

function v = snr(x,y)

% snr - signal to noise ratio
%
%   v = snr(x,y);
%
% v = 20*log10( norm(x(:)) / norm(x(:)-y(:)) )
%
%   x is the original clean signal (reference).
%   y is the denoised signal.
% 
%   Copyright (c) 2008 Gabriel Peyre

v = 20*log10(norm(x(:))/norm(x(:)-y(:)));