16qam_test.m

基于MIMO通信系统中的各种空时编译码算法

clear all;

tx_bits=randint(100000,1);
M=16;
EbNo=0:2:18;
P_qam=10;
temp=reshape(tx_bits,log2(M),[]).';
tx_num=bi2de(temp,'left-msb');
tx_sig=qammod(tx_num,M);
%% AWGN
SNR=EbNo+10*log10(log(M));
for ii=1:length(SNR) 
    rx_sig=tx_sig+sqrt((0.5*P_qam)/(10^(SNR(ii)/10)))*(randn(size(tx_sig))+j*randn(size(tx_sig)));
    rx_num=qamdemod(rx_sig,M);
    rx_bits=reshape(de2bi(rx_num,log2(M),'left-msb').',1,[]).';
    [num_err ber(ii)]=biterr(tx_bits,rx_bits);
end
%% 
semilogy(EbNo,ber,'r-o');
xlabel('EbNo');
ylabel('BER');
title([num2str(M) ' QAM']);
grid on;
clear all;