1 二進(jìn)制雙極性通信系統(tǒng)

1、二進(jìn)制雙極性通信系統(tǒng)的蒙特卡羅仿真蒙特卡羅計(jì)算機(jī)仿真用于估算數(shù)字通信系統(tǒng)的誤碼率，特別適用于難于對判決器的性能進(jìn)行分析的情況。利用蒙特卡羅仿真估算二進(jìn)制雙極性通信系統(tǒng)的誤碼率具體方法可以用下圖表示：二進(jìn)制雙極性通信系統(tǒng)的蒙特卡羅仿真模型如圖，首先仿真產(chǎn)生判決器的輸入隨機(jī)變量。使用均勻分布隨機(jī)數(shù)發(fā)生器產(chǎn)生來自于二進(jìn)制數(shù)據(jù)源的二進(jìn)制01信息序列，該01序列被映射為土E的序列，E代表信號能量，且可歸一化為1。使用高斯噪聲發(fā)生器產(chǎn)生均值為0,方差為的高斯隨機(jī)數(shù)序列n。判決器將隨機(jī)變量和判決電平0相比較。如果大于0，判決發(fā)送比特是0，否則判決發(fā)送比特是1。該判決輸出與所發(fā)送的信息比特序列相比較，并計(jì)算誤

2、比特?cái)?shù)和誤碼率。進(jìn)行編碼仿真時，該模型應(yīng)在D和D處加上相應(yīng)的編碼器和譯碼器。下面繪制了蒙特卡羅計(jì)算機(jī)仿真的基本流程圖。蒙特卡羅計(jì)算機(jī)仿真的基本流程圖繪制蒙特卡羅仿真曲線的一般流程圖MATLAB源程序functionp=smldPe55(snr_in_dB)%smldPe55.m二進(jìn)制雙極性通信系統(tǒng)的蒙特卡羅仿真%snr_in_dB信噪比%p誤碼率E=1;SNR=exp(snr_in_dB*log(10)/10);%Eb/N0sgma=E/sqrt(2*SNR);N=10A5;%仿真序列長度10A5，運(yùn)行時間約1分鐘temp=0;dsource=0;decis=0;numoferr=0;p=nu

3、moferr/N;fori=1:Ntemp=rand;if(temp0.5)dsource=0;elsedsource=1;endif(dsource=0)r=-E+gngauss(sgma);elser=E+gngauss(sgma);endif(r0)decis=0;elsedecis=1;endif(decis=dsource)numoferr=numoferr+1;endendnumoferr,p=numoferr/N;%smldPe55.m二進(jìn)制雙極性通信系統(tǒng)的蒙特卡羅仿真%ce55.m仿真繪圖語句%仿真序列長度10八5,運(yùn)行時間約15分鐘echoonSNRindB551=0:16;

4、SNRindB552=0:0.1:16;smld_err_prb55=zeros(1,length(SNRindB551);%計(jì)算仿真誤碼率theo_err_prb55=zeros(1,length(SNRindB552);fori=1:length(SNRindB551)smld_err_prb55(i)=smldPe55(SNRindB551(i);endfori=1:length(SNRindB552)SNR55=exp(SNRindB552(i)*log(10)/10);theo_err_prb55(i)=Qfunct(sqrt(2*SNR55);%Qfuncty=(1/2)*erfc(x/sqrt(2);理論誤碼率公式endfigure;semilogy(SNRindB551,smld_err_prb55,r*-);axis(0,16,1e-5,1);xlabel(Eb/N0indB);ylabel(Pe);title(二進(jìn)制雙極性通信系統(tǒng)的蒙特卡羅仿真結(jié)果與理論值比較);holdonsemilogy(SNRindB552,theo_err_prb55,b:);%ce55.m仿真繪圖語句運(yùn)行結(jié)果極性通信系統(tǒng)的蒙特卡羅仿