lin_iirmse.m

来自「基于Matlab」· M 代码 · 共 20 行

function mse = iirmse(c,spacing,SNR)

dw=2*pi/256;
w=0:dw:2*pi;
if (spacing==1),
   C=freqz(c,1,w);
   T=abs(C).^2;
else
   ca=c(1:2:length(c));
   cb=c(2:2:length(c));
   Ca=freqz(ca,1,w);
   Cb=freqz(cb,1,w);
   T=abs(Ca).^2+abs(Cb).^2;
end;

for i=1:length(SNR),
   lambda=10^(-SNR(i)/10);
   mse(i)=real(1/(2*pi)*sum(lambda./(T+lambda)*dw));
end;