%文件名：PSNR.m
%編??寫：郭林庚
%編寫時間：2005.11.01
%函數(shù)功能：本函數(shù)將完成對輸入圖像的峰值信噪比計算
%輸入格式舉例：psnr=PSNR（‘lena.jpg‘‘lenawater.jpg‘）;
%參數(shù)說明：
%original為原始圖像
%test為加有水印的圖像
%psnrvalue為兩者峰值信噪比
%實際上也計算了MSE，SNR值

function?psnrvalue=PSNR（originaltest）;

%計算原始圖像的信號功率
A=imread（original）;
%A=rgb2gray（A）;
A=double（A）;
B=imread（test）;
%B=rgb2gray（B）;
B=double（B）;

%計算MSE
%判斷輸入圖像是否有效
[mn]=size（A）;
[m2n2]=size（B）;
if?m2~=m||n2~=n
????err