dfe_iirmse.m

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

function mse = dfe_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)=exp(real(1/2/pi*sum(log(lambda./(T+lambda))*dw)));
end;