function?snr=SNR（IIn）
%?計算信號噪聲比函數
%?by?Qulei
%?I?:original?signal
%?In:noisy?signal（ie.?Original?signal?+?noise?signal）
%?snr=10*log10（sigma2（I2）/sigma2（I2-I1））

[rowcolnchannel]=size（I）;

snr=0;
if?nchannel==1%gray?image
Ps=sum（sum（（I-mean（mean（I）））.^2））;%signal?power
Pn=sum（sum（（I-In）.^2））;%noise?power
snr=10*log10（Ps/Pn）;
elseif?nchannel==3%color?image
for?i=1:3
Ps=sum（sum（（I（::i）-mean（mean（I（::i））））.^2））;%signal?power
Pn=sum（sum（（I（::i）-In（::i））.^2））;%noise?power
snr=snr+10*log10（Ps/Pn）;
end
snr=snr/3;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1283??2012-11-23?10:25??snr_practice.m

?????文件????????558??2012-11-21?10:25??SNR.m

-----------?---------??----------?-----??----

?????????????????1841????????????????????2