function?[PSNRmse]=psnr（XY）
%1．求圖像峰值信噪比PSNR的代碼
%?function?[PSNRmse]=psnr（XY）
%?Peak?signal?to?noise?ratio?of?the?difference?between?images?and?the?mean?square?error
%?If?the?second?input?Y?is?missing?then?the?PSNR?and?MSE?of?X?itself?becomes?the?output?（as?if?Y=0）.

if?nargin<2?D=X;
else
if?any（size（X）~=size（Y））?error（‘The?input?size?is?not?equal?to?each?other!‘）;
end
D=X-Y;
end

mse=sum（D（:）.*D（:））/numel（X）;
PSNR=10*log10（255^2/mse）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????472??2013-09-25?22:02??psnr.m
?????文件?????????271??2013-09-24?22:41??SNR.m