psnr.m

来自「QMF滤波器在语音编码中的应用,通过matlab 实验程序设计和实现二通道镜像滤」· M 代码 · 共 16 行

% function for finding MSE and PSNR
function p = psnr(x,xc,g)
x=double(x);
xc=double(xc);
t=(x-xc).^2;
disp('MSE:');
s=sum(sum(t));
%disp(s);
%p=g*g*(size(x,1)*size(x,2))/(s);
n=size(x,1)*size(x,2);

disp('PSNR:');
p = (g*g)/(s/n);
p = 10*log10(p);
disp(p);