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);