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;