求两幅图像的均方差（MATLAB实现）

求两幅图像的均方差

I=imread('moban.jpg');

I1=double(I);
B=imread('h.jpg');
B1=double(B);
s1=0;
s2=0;
[m n]=size(I1);
for i=1:1:m;
    for j=1:1:n;
       a=[I1(i,j)-B1(i,j)].^2;
       b=I1(i,j).^2;
    s1=s1+a;
    s2=s2+b;   
  end
end
s=s1/s2;
disp(s');