数字通信系统误码率仿真分析

3G移动通信实验报告实验名称：数字通信系统误码率仿真分析学生姓名：学生学号：学生班级:所学专业：实验日期:1.实验目的1.掌握几种典型数字通信系统误码率分析方法。2.掌握误码率对数字通信系统的影响及改进方法。2.实验原理1、数字通信系统的主要性能指标通信的任务是传递信息，因此信息传输的有效性和可靠性是通信系统的最主要的质量指标。有效性是指在给定信道内能传输的信息内容的多少，而可靠性是指接收信息的准确程度。为了提高有效性，需要提高传输速率，但是可靠性随之降低。因此有效性和可靠性是相互矛盾的，又是可交换的。可以用降低有效性的办法提高可靠性，也可以用降低可靠性的办法提高有效性。数字通信系统的有效性通常用信息传输速率来衡量。当信道一定时，传输速率越高，有效性就越好。传输速率有三种定义：码元速率（sR）：单位时间内传输的码元数目，单位是波特（Baud），因此又称为波特率；信息速率（bR）：单位时间内传输的信息量（比特数），单位是比特/秒（b/s），因此又称为比特率；消息速率（MR）：单位时间内传输的消息数目。对于M进制通信系统，码元速率与信息速率的关系为：sbMRRsb/log2baudMRRbs2log特别说明的是，在二进制数字通信系统中信源的各种可能消息的出现概率相等时，码元速率和信息速率相等。在实际应用中，通常都默认这两个速率相等，所以常常简单地把一个二进制码元称为一个比特。数字通信系统的可靠性的衡量指标是错误率。它也有三种不同定义：误码率（eP）：指错误接收码元数目在传输码元总数中所占的比例，即传输总码元数错误接收码元数eP误比特率（bP）：指错误接收比特数目在传输比特总数中所占的比例，即传输总比特数错误接收比特数bP误字率（WP）：指错误接收字数在传输总字数中所占的比例。若一个字由k比特组成，每比特用一码元传输，则误字率等于keWPP11对于二进制系统而言，误码率和误比特率显然相等。而M进制信号的每个码元含有Mn2log比特，并且一个特定的错误码元可以有1M种不同的错误样式。其中恰好有i比特错误的错误样式有inC个。假设这些错去样式以等概出现，则当一个码元发生错误时，在n个比特种错误比特所占比例的数学期望值等于12121111MMMCniMEnEninin一个码元中的比特数错误比特数所以，当M比较大时，误比特率等于eeebPPMMPnEP2112实际应用中，在计算误比特率时，0/nEb是个十分重要的参数。bE为单位比特的平均信号能量，0n为噪声的单边功率谱密度。但人们在实际系统中能够直接测量到的是平均信号功率S和噪声功率N，并且由此可计算出信噪比NS/。若一个二进制通信系统的信息传输速率为bR，接收机带宽为B，则信噪比可表示为BRnENSbb0这里，BRb/为单位频带的比特率，它表示特定调制方案下的频带利用率，又称频带效率。在二进制数字调制时，接收端可以采用相干解调，也可以采用非相干解调，它们的抗噪声能力不同，误码性能也不同。下面就分相干解调与非相干解调来分析数字通信系统的误码率。2、相干解调时的误码率相干解调需要在接收端产生用于相干的参考载波信号，其一般模型为：LPFttCcdcostSptSdtS图10-1相干解调的一般模型设接收信号为ttSttStScQcIsincos与同频同相相干载波相乘后得ttSttStSttStScQcIIcp2sin212cos2121cos经低通滤波后得tftStSId21在实际应用中，经常把相干解调与最佳接收混为一谈，确切地说，只有当相干解调中的滤波器严格按照匹配滤波器的要求来设计，才是真正的最佳接收。由通信原理知识可知，在匹配滤波器条件下求得的二进制调制的最小误比特率为10nEQPbb其中Q函数的定义为dyeQy2/221。为相关系数，取值范围为1,1，取决于发送的信号tStS21和的相似程度，即21012ssTEEdttStSs1sE、2sE分别为tStS21和在sTt0内的能量。下面分析几种典型的二进制数字调制系统在相干解调时的误码率：对于2ASK信号，sCTttAtStS0,cos,021，由于sTssdttSEE02221,0，则0。而2121ssbEEE，因此有02,/nEQPbASKb对于2PSK信号，有sCsbCsbTttTEtStTEtS0,cos2,cos221，可求出1，因此02,/2nEQPbPSKb对于2FSK信号，有ssbsbTttTEtStTEtS0,cos2,cos22211，其相关系数为ssTsTsbbdtttTtdttTEE01212012coscos1coscos21当sTnff2/12时tStS21和相互正交，0。当12很大时，也接近于0。对于0的2FSK信号，其误比特率为02,/nEQPbFSKb3、非相干解调时的误码率非相干解调的最大优点是不需要在接收端产生用于相干的参考载波。因此其设备相对要简单便宜一些，但性能相对与相干解调来说有点差。对于2FSK信号，有0,2/exp21nEPbNCFSKb（2）2ASK与2FSK相同，即0,,2/exp21nEPPbNCFSKbNCASKb（3）对于2DPSK信号，有02,/exp21nEPbDPSKb下图显示了不同二进制数字通信系统在相干解调和非相干解调时的误码率比较。图10-2二进制调制的误比特率曲线由上图可以看出，使用相干解调时的误码率整体上比使用非相干解调时要低，这也说明使用相干解调时的抗噪声能力比较强。4、多进制数字调制的误码率在多进制数字调制中，每个符号间隔sTt0内，可能发送的符号有M种：tS1，tS2，…，tSM。在实际应用中，通常取nM2，n为大于1的正整数。当携带信息的参数分别为载波的幅度、频率或相位时，对应的有M进制幅度键控（MASK）、M进制频移键控（MFSK）、M进制相移键控（MPSK）。也可以把其中的两个参数组合起来调制，例如M进制正交幅度键控（QAM）。多进制数字调制的调制解调方法与二进制的类似，主要差别在基带成型阶段。解调多采用相干解调，下面就列出几种典型多进制数字调制在相干解调时的误码率公式：M进制幅度键控（MASK）BRnEMMQMMNSMQMMPsbMASKS0222,1log3121312M进制频移键控（MFSK）MnEQMPbMFSKs20,log1M进制相移键控（MPSK）MnEQMnEerfcPssMPSKs2020,sin22sin（4）多进制正交幅度调制（MQAM）022,1log612nELLQLLPbMQAMs这里2LM。3.实验内容1、编写MATLAB程序实现Rake接收机。2、修改信噪比，观察三种合并算法的误码率。4.实验代码clearall;dn=100;d=randint(1,dn);d4_qpsk=d(1,2:2:dn)+2*d(1,1:2:dn-1);fori=1:dn/2ifd4_qpsk(1,i)==0I(1,i)=1;Q(1,i)=1;elseifd4_qpsk(1,i)==1I(1,i)=-1;Q(1,i)=1;elseifd4_qpsk(1,i)==2I(1,i)=1;Q(1,i)=-1;elseifd4_qpsk(1,i)==3I(1,i)=-1;Q(1,i)=-1;endendendendenddin_2psk=2*d-1;en_min=-10;en_max=0;en_d=0.5;cn=1000;s_n=1;foren=en_min:en_d:en_maxern_2psk=0;ern_2ask=0;ern_qpsk=0;snr=en+10*log(2)/log(10);fory=1:cndout_2psk=awgn(din_2psk,snr);dout_2ask=awgn(d,snr);Iout=awgn(I,snr);Qout=awgn(Q,snr);dout_2psk=sign(dout_2psk);fori=1:dnifdout_2ask(1,i)0.5dout_2ask(1,i)=1;elsedout_2ask(1,i)=0;endifdout_2ask(1,i)~=d(1,i)ern_2ask=ern_2ask+1;endenddout_2psk=(dout_2psk+1)*0.5;fori=1:dnifdout_2psk(1,i)~=d(1,i)ern_2psk=ern_2psk+1;endendfori=1:dn/2ifIout(1,i)0ifQout(1,i)0dout_qpsk(1,2*i)=0;dout_qpsk(1,2*i-1)=0;elsedout_qpsk(1,2*i)=0;dout_qpsk(1,2*i-1)=1;endelseifQout(1,i)0dout_qpsk(1,2*i)=1;dout_qpsk(1,2*i-1)=0;elsedout_qpsk(1,2*i)=1;dout_qpsk(1,2*i-1)=1;endendifdout_qpsk(1,2*i)~=d(1,2*i)ern_qpsk=ern_qpsk+1;endifdout_qpsk(1,2*i-1)~=d(1,2*i-1)ern_qpsk=ern_qpsk+1;endendenderror_rate_2psk(1,s_n)=ern_2psk/dn/cn;error_rate_2ask(1,s_n)=ern_2ask/dn/cn;error_rate_qpsk(1,s_n)=ern_qpsk/dn/cn;Pbpsk(1,s_n)=qfunc(sqrt(10^(snr/10)));P2ask(1,s_n)=qfunc(sqrt(0.25*10^(snr/10)));Pqpsk(1,s_n)=0.5*2*qfunc(sqrt(2*2*10^((en)/10)*(sin(pi/4))^2));s_n=s_n+1;endx=en_min:en_d:en_max;subplot(1,2,1);semilogy(x,error_rate_2psk,'o-',x,Pbpsk,'*-',x,error_rate_2ask,'o-',x,P2ask,'*-');legend('bpsk仿真值','bpsk理论值','2ask仿真值','2ask理论值')title('2PSK+2ASK误码率分析');xlabel('Eb/n0dB');ylabel('error_ratedB');gridon;subplot(1,2,2);semilogy(x,error_rate_qpsk,'o-',x,Pqpsk,'*-');legend('qpsk仿真值','qpsk理论值')title('QPSK误码率分析');xlabel('Eb/n0dB');ylabel('error_ratedB');gridon;;5.实验结果