snr.m

此程序主要是通过使用小波达到图像文件压缩的目的

function[r]=snr(X,Y)
[u,v]=size(Y);
%SNR=10*log10(Ps/Pn).单位dB(decibel,分贝），Ps、Pn分别为信号与噪音的功率，如果是振幅相比的话，前面乘20。
N=0;S=0;
for i=1:u
        for j=1:v
          N=N+(abs(X(i,j)-Y(i,j)))^2;
          S=S+X(i,j)^2; 
        end
end

r=10*log10(S/N);