A?=?imread（‘3.bmp‘）;%讀入圖片
A?=?double（A）;??%把矩陣A定義為double型
B?=?imread（‘333.bmp‘）;%同上
B?=?double（B）;

%%?計算MSE

[mn]?=?size（A）;%計算A的像素（長和寬）
[m1n1]?=?size（B）;
if?m~=m1||n~=n1?????%如果A?B長寬不同就顯示“圖像大小不一致”
????error（‘圖像大小不一致‘）;
end

msevalue?=?0;??%定義初始值
for?i?=?1:m
????for?j?=?1:n????%循環用于求和
????????msevalue?=?msevalue+（A（ij）-B（ij））^2;????%每循環一次就加一次
????end
end
msevalue?=?msevalue/（m*n）;??%用求的和除以m*n
if?msevalue?==?0??????%表示A（ij）=B（ij）
????error（‘圖像完全相同‘）;
end

%%?計算峰值信噪比
psnrvalue?=?255^2/msevalue;??%這個你懂得?????第二個程序類似的
psnrvalue?=?10*log10（psnrvalue）

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

?????文件????????693??2012-05-27?17:13??信噪比和峰值信噪比計算\test1.m

?????文件????????735??2012-05-30?16:13??信噪比和峰值信噪比計算\test2.m

?????目錄??????????0??2012-06-06?13:13??信噪比和峰值信噪比計算

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

?????????????????1428????????????????????3