与this question类似,我努力让代码工作以计算PSNR。
我得到的错误是" Matrix尺寸必须同意"。
以下是代码:
[y1,fs1]=audioread('track.mp3');
[y2,fs2]=audioread('track.mp3');
[c1x,c1y]=size(y1);
[c2x,c2y]=size(y1);
R=c1x;
C=c1y;
%err = sum((y1-y2).^2)/(R*C);
err = sum((y1-y2)/(R*C));
MSE=sqrt(err);
MAXVAL=65535;
PSNR = 20*log10(MAXVAL/MSE);
disp(['mse=' num2str(MSE) ' PSNR=' num2str(PSNR)]);