PSNR用于matlab中的音频信号

Question

我正在尝试找出项目中音频文件的MSE和PSNR的值。我到目前为止编写的代码如下：

[y1,fs1, nbits1,opts1]=wavread('one.wav');
[y2,fs2, nbits2,opts2]=wavread('newOne.wav');
[c1x,c1y]=size(y1);
[c2x,c2y]=size(y1);
if c1x ~= c2x
    disp('dimeonsions do not agree');
 else
 R=c1x;
 C=c1y;
 err = (((y1-y2).^2)/(R*C));
 MSE=sqrt(err);
 MAXVAL=65535;
  PSNR = 20*log10(MAXVAL/MSE); 
  disp(['mse=' num2str(MSE) ' PSNR=' num2str(PSNR)]);
end

但我收到的错误如下：

连接的矩阵的维度不一致。

我做错了什么？

Answer 1

您需要求和平方误差才能计算MSE - 更改：

err = (((y1-y2).^2)/(R*C));

为：

err = sum((y1-y2).^2)/(R*C);