1 二进制双极性通信系统

1、二进制双极性通信系统的蒙特卡罗仿真蒙特卡罗计算机仿真用于估算数字通信系统的误码率，特别适用于难于对判决器的性能进行分析的情况。利用蒙特卡罗仿真估算二进制双极性通信系统的误码率具体方法可以用下图表示：二进制双极性通信系统的蒙特卡罗仿真模型如图，首先仿真产生判决器的输入随机变量。使用均匀分布随机数发生器产生来自于二进制数据源的二进制01信息序列，该01序列被映射为土E的序列，E代表信号能量，且可归一化为1。使用高斯噪声发生器产生均值为0,方差为的高斯随机数序列n。判决器将随机变量和判决电平0相比较。如果大于0，判决发送比特是0，否则判决发送比特是1。该判决输出与所发送的信息比特序列相比较，并计算误

2、比特数和误码率。进行编码仿真时，该模型应在D和D处加上相应的编码器和译码器。下面绘制了蒙特卡罗计算机仿真的基本流程图。蒙特卡罗计算机仿真的基本流程图绘制蒙特卡罗仿真曲线的一般流程图MATLAB源程序functionp=smldPe55(snr_in_dB)%smldPe55.m二进制双极性通信系统的蒙特卡罗仿真%snr_in_dB信噪比%p误码率E=1;SNR=exp(snr_in_dB*log(10)/10);%Eb/N0sgma=E/sqrt(2*SNR);N=10A5;%仿真序列长度10A5，运行时间约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二进制双极性通信系统的蒙特卡罗仿真%ce55.m仿真绘图语句%仿真序列长度10八5,运行时间约15分钟echoonSNRindB551=0:16;

4、SNRindB552=0:0.1:16;smld_err_prb55=zeros(1,length(SNRindB551);%计算仿真误码率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(二进制双极性通信系统的蒙特卡罗仿真结果与理论值比较);holdonsemilogy(SNRindB552,theo_err_prb55,b:);%ce55.m仿真绘图语句运行结果极性通信系统的蒙特卡罗仿