loading
I=inv(Hls)* Y; for k=1:64
if(real(I(k))>0)%Putting it through a slicer I(k)=1; else I(k)=-1; end end for k=1:64
if(I(k)~=d(k))
error_count_ls=error_count_ls+1; end end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % I:MMSE ESTIMATOR BASED RECEIVER:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %I(k) represents the decision matrix.. I=inv(Hmmse)* Y; for k=1:64
if(real(I(k))>0)%Putting it through a slicer I(k)=1; else I(k)=-1; end end for k=1:64
if(I(k)~=d(k))
error_count_mmse=error_count_mmse+1; end end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% End %End of the 1000 run simulation.. ser_ls(n)=error_count_ls/64000; ser_mmse(n)=error_count_mmse/64000; ser_ls ser_mmse
SNR(n)=SNR_send; end; %作图：
semilogy(SNR,ser_ls,'-r*'); hold on;
semilogy(SNR,ser_mmse,'-b*'); grid on;
xlabel('SNR ( dB )'); ylabel('SER');
legend('SER-ls','SER-mmse');
title('plot of SNR and SER for an OFDM system with MMSE/LS estimator');