matlab - 平均二乗誤差 (MSE) の問題

Question

異なるコンピューター (MATLAB) で平均二乗誤差を見つけるために次のコードを実行すると、異なる答えが返されます。なぜこれが当てはまるのか、誰かが考えを持っていますか?

function [mse,psnr] = mse_psnr(im1,im2)

N = size(im1);
acc = 0;

for k1=1:N(1)
  for k2=1:N(2)
     acc = acc+ ( im1(k1,k2) - im2(k1,k2) )^2;   
  end
end

mse  = acc/(N(1)*N(2));
psnr = 10*log10((255^2)/mse);

Answer 1

コードをベクトル化して、ネストされたループを削除してください。

function [mse psnr] = mse_psnr( im1, im2 )
% make sure input data type is double
im1 = double(im1);
im2 = double(im2);
acc = sum( (im1(:)-im2(:)).^2 ); % instead of nested loop...
mse = acc / numel(im1);
psnr = 10*log10( (255^2)/mse );