A?=?double（imread（‘1.bmp‘））;%讀入圖片
B?=?double（imread（‘2.bmp‘））;%同上
[mn]?=?size（A）;
[m1n1]?=?size（B）;
ps=0;
for?i?=?1:m
????for?j?=?1:n
????????ps?=?ps+A（ij）^2;
????end
end
pn=0;
for?i?=?1:m1
????for?j?=?1:n1
????????pn?=?pn+（A（ij）-B（ij））^2;
????????%pn?=?pn+（（A（ij））^2-（B（ij））^2）;
????end
end
snr?=?10*log10（ps/pn）

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????文件?????????343??2019-07-20?19:45??\wudixiaofeixia-9544697-SNR.m
?????文件??????????36??2019-07-20?19:45??\no.txt