function srr=SRRA(A00,A11) %performs signal to noise ratio (18) % \SNRA(\A1)=10\log_{10}\frac{\nn{\A0}}{\nn{ \A1-\A0}}. % % Input % A0 reference filter % A1 estimated filter % Output % srr signal to noise ratio % Written by Alexis Benichoux, 2011 [M N useless]=size(A00); for m=1:M for n=1:N A0=A00(m,n,:); A1=A11(m,n,:); r=(A0(:)-A1(:)); s=sum(r.*r); s0=sum(A0(:).*A0(:)); srra(m,n)=10*log10(s0/s); end end srr=mean(srra(:)); end