erms的m函数

function rmse=compare(f1,f2,scale)
%COMPARE Computers and displays the error between two matrices.
%RMSE = COMPARE(F1, F2, SCALE) 返回输入F1和F2之间的均方根误差。
%返回输入F1和F2之间的均方根误差，显示差值的直方图 。
%并显示一个缩放的差值图像。当SCALE被省略时，使用一个为1的缩放系数。

%检查输入参数并设置默认值
error(nargchk(2,3,nargin));
if nargin<3
    scale=1;
end

%计算均方根误差
e=double(f1)-double(f2);
[m,n]=size(e);
rmse=sqrt(sum(e(:).^2)/(m*n));

%如果出现错误，输出错误图像和柱状图（即rmse~=0）。
if rmse
    %来自错误的柱状图。
    emax=max(abs(e(:)));
    [h,x]=hist(e(:),emax);
    if length(h)>=1
        figure;bar(x,h,"k");
        
        %对称地缩放错误图像并显示
        emax=emax/scale;
        e=mat2gray(e,[-emax,emax]);
        figure;imshow(e);
    end
end