clear;
close?all;
I=imread（‘lena‘‘bmp‘）;%I為原圖像?
%?I1=wiener2（I[3?3]）;
%?L=medfilt2（I1[3?3]）;
%?I2=imadjust（L[0.1?0.9][0?1]1.2）;%I2為變換后圖像
I2=imread（‘lenaq0‘‘bmp‘）
%以上I是原圖像，I2是變換后的圖像，具體變換可以自定義
%下面是均方誤差mse峰值信噪比psnr和熵entropy的計算
X=I;
Y=I2;
mse_m=double（zeros（256256））;?%均方誤差
for?i=1:256
for?j=1:256
mse_m（ij）=（X（ij）-Y（ij））^2;
end
end
mse=sum（mse_m（1:256））/256*256;
psnr=10*log（double（255*255/mse））/log（10）;%峰值信噪比
%?%?下面計算熵
%?t（1:256）=0;?
%?%給每個灰度值計數
%?for?i=1:180
%?for?j=1:180
%?t（Y（ij）+1）=t（Y（ij）+1）+1;
%?end
%?end
%?%計算各個灰度出現頻率
%?for?i=1:256
%?t（i）=t（i）/（256*256）;
%?end
%?entropy=0;?%計算熵
%?for?i=1:256
%?if?t（i）>0
%?entropy=entropy-t（i）*log（t（i））/log（2）;
%?end
%?end

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

????I.A....???????802??2013-03-31?17:15??水印評價代碼\pingjia.m

?????文件????????884??2011-04-07?07:55??水印評價代碼\PSNR&NC.txt

????I.A....???????419??2013-03-31?20:42??水印評價代碼\PSNR和歸一化相關系數NC值.txt

?????目錄??????????0??2013-03-31?20:51??水印評價代碼

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

?????????????????2105????????????????????4