test.m

来自「本代码是基于OFDM系统的一种转换域估计算法」· M 代码 · 共 20 行

% ls mse
figure(1);

SNR_dB=[0 5 10 15 20 25 30 35];
snr=zeros(1,length(SNR_dB));
for i=1:length(SNR_dB)
    snr(i)=10^(SNR_dB(i)/10);
end
semilogy(SNR_dB,1./snr,'b-*');

%lmmse mse
sum_lamiga=sum(dlamda);
deta=dlamda./(dlamda+((17/9)/snr));
mse=0;
for i=1:128
    mse=mse+(1/128)*{(dlamda.*((1-deta).^2))+(17/(9*snr))*deta.^2};
end