function?[PSNR?MSE]?=?psnr（X?Y）
%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%?計算峰值信噪比PSNR
%?將RGB轉成YCbCr格式進行計算
%?如果直接計算會比轉后計算值要小2dB左右（當然是個別測試）
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%
?if?size（X3）~=1???%判斷圖像時不是彩色圖，如果是，結果為3，否則為1
???org=rgb2ycbcr（X）;
???test=rgb2ycbcr（Y）;
???Y1=org（::1）;
???Y2=test（::1）;
???Y1=double（Y1）;??%計算平方時候需要轉成double類型，否則uchar類型會丟失數據
???Y2=double（Y2）;
?else??????????????%灰度圖像，不用轉換
?????Y1=double（X）;
?????Y2=double（Y）;
?end
?
if?nargin<2????
???D?=?Y1;
else
??if?any（size（Y1）~=size（Y2））
????error（‘The?input?size?is?not?equal?to?each?other!‘）;
??end
?D?=?Y1?-?Y2;?
end
MSE?=?sum（D（:）.*D（:））?/?numel（Y1）;?
PSNR?=?10*log10（255^2?/?MSE）;

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

?????文件????????764??2015-08-28?15:47??PSNR和SSIM—matlab\psnr.m

?????文件???????7325??2015-08-28?16:32??PSNR和SSIM—matlab\ssim.m

?????文件?????????22??2015-08-28?16:34??PSNR和SSIM—matlab\程序說明.txt

?????目錄??????????0??2015-08-28?16:37??PSNR和SSIM—matlab

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

?????????????????8111????????????????????4