用MATLAB求RMSE

function f=RMSE(h1,h2)%RMSE return RMSE(均方根误差) 求两图像的均方根误差%input must be a imagehandle 输入图像句柄%image fusion evaluate parameter 图像融合评价参数% example% 标准图像 h1% 融合后图像 h2% f=RMSE(h1,h2);%融合图像与标准图像差异程度,差异越小说明融合图像与标准图像越接近s=size(size(h1));%判断是灰度图还是RGBif s(2)==2f1=h1;f2=h2;elsef1=rgb2gray(h1);f2=rgb2gray(h2);end G1=double(f1);G2=double(f2);[m1,n1]=size(G1);[m2,n2]=size(G2);m=min(m1,m2);n=min(n1,n2);c=0;for i=1:mfor j=1:nw=G1(i,j)-G2(i,j);c=c+w^2;endendf=sqrt(c/(m*n));