QPSK和16QAM信號通過AWGN信道的誤比特率

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%本程序用于仿真QPSK，16QAM，通過AWGN信道的BER隨SNR的變化規(guī)律。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clearall;clcM_4psk=4;%4psk的進制數(shù)k_4psk=log2(M_4psk);%4psk的一個符號包含的比特數(shù)M_16qam=16;%16qam的進制數(shù)k_16qam=log2(M_16qam);%16qam一個符號包含的比特數(shù)NumSymbols=10^6;%定義要發(fā)送的比特數(shù)rngdefault%使用默認的隨機數(shù)生成器dataIn=randi([01],NumSymbols,1);%生成二進制發(fā)送數(shù)據(jù)dataInMatrix_4psk=reshape(dataIn,length(dataIn)/2,2);%將發(fā)送二進制數(shù)據(jù)整理成兩個一行，準備進行進制轉換dataSymbolsIn_4psk=bi2de(dataInMatrix_4psk);%轉換成十進制dataInMatrix_16qam=reshape(dataIn,length(dataIn)/4,4);%將發(fā)送二進制數(shù)據(jù)整理成四個一行，準備進行進制轉換dataSymbolsIn_16qam=bi2de(dataInMatrix_16qam);%轉換成十進制dataModG_4psk=pskmod(dataSymbolsIn_4psk,M_4psk,0,'gray');%4psk調制以及格雷碼編碼dataModG_16qam=qammod(dataSymbolsIn_16qam,M_16qam,0,'gray');%16qam調制以及格雷碼編碼EbNo=[-2:2:15];%設置SNR范圍berG_4psk=zeros(1,length(EbNo));%定義4psk誤碼率數(shù)組berG_16qam=zeros(1,length(EbNo));%定義16qam誤碼率數(shù)組forindex=1:length(EbNo)%計算以db表示的SNRsnr_4psk=EbNo(index)+10*log10(k_4psk);snr_16qam=EbNo(index)+10*log10(k_16qam);%信號通過AWGN信道receivedSignalG_4psk=awgn(dataModG_4psk,snr_4psk,'measured');receivedSignalG_16qam=awgn(dataModG_16qam,snr_16qam,'measured');%解調以及格雷碼解碼dataSymbolsOutG_4psk=pskdemod(receivedSignalG_4psk,M_4psk,0,'gray');dataSymbolsOutG_16qam=qamdemod(receivedSignalG_16qam,M_16qam,0,'gray');dataOutMatrixG_4psk=de2bi(dataSymbolsOutG_4psk,k_4psk);%十進制轉二進制，準備計算BERdataOutG_4psk=dataOutMatrixG_4psk(:);%將數(shù)據(jù)轉化成一維向量dataOutMatrixG_16qam=de2bi(dataSymbolsOutG_16qam,k_16qam);%十進制轉二進制，準備計算BERdataOutG_16qam=dataOutMatrixG_16qam(:);%將數(shù)據(jù)轉化成一維向量[numErrorsG_4psk,berG_4psk(index)]=biterr(dataIn,dataOutG_4psk);%計算4psk誤碼率[numErrorsG_16qam,berG_16qam(index)]=biterr(dataIn,dataOutG_16qam);%計算16qam誤碼率end%計算信號通過AWGN信道的理論誤比特率==========================================TheoryBER_4psk=berawgn(EbNo,'psk',4,'nondiff');%計算4psk通過AWGN理論誤比特率TheoryBER_16qam=berawgn(EbNo,'qam',16);%計算16qam通過AWGN理論誤比特率%==========================================================================%畫圖semilogy(EbNo,TheoryBER_4psk,'b-',EbNo,TheoryBER_16qam,'r-');%畫出理論BER曲線holdon;semilogy(EbNo,berG_4psk,'b*',EbNo,berG_16qam,'r*');%畫出仿真BER曲線gridon;xlabel('Eb/N0(dB)');y