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

3G移动通信实验报告实验名称：数字通信系统误码率仿真分析学生：学生学号：学生班级:所学专业：实验日期:1. 实验目的1. 掌握几种典型数字通信系统误码率分析方法2. 掌握误码率对数字通信系统的影响及改良方法2. 实验原理1、数字通信系统的主要性能指标通信的任务是传递信息，因此信息传输的有效性和可靠性是通信系统的最主要的质量指标有效性是指在给定信道能传输的信息容的多少，而可靠性是指接收信息的准确程度为了提高有效性，需要提高传输速率，但是可靠性随之降低因此有效性和可靠性是相互矛盾的，又是可交换的可以用降低有效性的方法提高可靠性，也可以用降低可靠性的方法提高有效性数字通信系统的有效性通常用信息传输速率来衡量当信道一定时，传输速率越高，有效性就越好传输速率有三种定义：码元速率〔〕：单位时间传输的码元数目，单位是波特〔Baud〕，因此又称为波特率；信息速率〔〕：单位时间传输的信息量〔比特数〕，单位是比特/秒〔b/s〕，因此又称为比特率；消息速率〔〕：单位时间传输的消息数目对于M进制通信系统，码元速率与信息速率的关系为：特别说明的是，在二进制数字通信系统源的各种可能消息的出现概率相等时，码元速率和信息速率相等。

在实际应用中，通常都默认这两个速率相等，所以常常简单地把一个二进制码元称为一个比特数字通信系统的可靠性的衡量指标是错误率它也有三种不同定义：误码率〔〕：指错误接收码元数目在传输码元总数中所占的比例，即误比特率〔〕：指错误接收比特数目在传输比特总数中所占的比例，即误字率〔〕：指错误接收字数在传输总字数中所占的比例假设一个字由k比特组成，每比特用一码元传输，那么误字率等于对于二进制系统而言，误码率和误比特率显然相等而M进制信号的每个码元含有比特，并且一个特定的错误码元可以有种不同的错误样式其中恰好有i比特错误的错误样式有个假设这些错去样式以等概出现，那么当一个码元发生错误时，在n个比特种错误比特所占比例的数学期望值等于所以，当M比拟大时，误比特率等于实际应用中，在计算误比特率时，是个十分重要的参数为单位比特的平均信号能量，为噪声的单边功率谱密度但人们在实际系统中能够直接测量到的是平均信号功率S和噪声功率N，并且由此可计算出信噪比假设一个二进制通信系统的信息传输速率为，接收机带宽为B，那么信噪比可表示为这里，为单位频带的比特率，它表示特定调制方案下的频带利用率，又称频带效率在二进制数字调制时，接收端可以采用相干解调，也可以采用非相干解调，它们的抗噪声能力不同，误码性能也不同。

下面就分相干解调与非相干解调来分析数字通信系统的误码率2、相干解调时的误码率相干解调需要在接收端产生用于相干的参考载波信号，其一般模型为：图10-1相干解调的一般模型设接收信号为与同频同相相干载波相乘后得经低通滤波后得在实际应用中，经常把相干解调与最正确接收混为一谈，确切地说，只有当相干解调中的滤波器严格按照匹配滤波器的要求来设计，才是真正的最正确接收由通信原理知识可知，在匹配滤波器条件下求得的二进制调制的最小误比特率为其中Q函数的定义为为相关系数，取值围为，取决于发送的信号的相似程度，即、分别为在的能量下面分析几种典型的二进制数字调制系统在相干解调时的误码率：对于2ASK信号，，由于，那么而，因此有对于2PSK信号，有，可求出，因此对于2FSK信号，有，其相关系数为当时相互正交，当很大时，也接近于0对于的2FSK信号，其误比特率为3、非相干解调时的误码率非相干解调的最大优点是不需要在接收端产生用于相干的参考载波因此其设备相对要简单廉价一些，但性能相对与相干解调来说有点差对于2FSK信号，有〔2〕2ASK与2FSK一样，即〔3〕对于2DPSK信号，有下列图显示了不同二进制数字通信系统在相干解调和非相干解调时的误码率比拟。

图10-2二进制调制的误比特率曲线由上图可以看出，使用相干解调时的误码率整体上比使用非相干解调时要低，这也说明使用相干解调时的抗噪声能力比拟强4、多进制数字调制的误码率在多进制数字调制中，每个符号间隔，可能发送的符号有M种：，，…，在实际应用中，通常取，n为大于1的正整数当携带信息的参数分别为载波的幅度、频率或相位时，对应的有M进制幅度键控〔MASK〕、M进制频移键控〔MFSK〕、M进制相移键控〔MPSK〕也可以把其中的两个参数组合起来调制，例如M进制正交幅度键控〔QAM〕多进制数字调制的调制解调方法与二进制的类似，主要差异在基带成型阶段解调多采用相干解调，下面就列出几种典型多进制数字调制在相干解调时的误码率公式：M进制幅度键控〔MASK〕M进制频移键控〔MFSK〕M进制相移键控〔MPSK〕〔4〕多进制正交幅度调制〔MQAM〕这里3. 实验容1、编写MATLAB程序实现Rake接收机2 、修改信噪比，观察三种合并算法的误码率4. 实验代码clear all;dn=100;d=randint(1,dn);d4_qpsk=d(1,2:2:dn)+2*d(1,1:2:dn-1);fori=1:dn/2if d4_qpsk(1,i)==0I(1,i)=1;Q(1,i)=1;elseif d4_qpsk(1,i)==1I(1,i)=-1;Q(1,i)=1;elseif d4_qpsk(1,i)==2I(1,i)=1;Q(1,i)=-1;elseif d4_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;=1000;s_n=1;for en=en_min:en_d:en_max ern_2psk=0; ern_2ask=0;ern_qpsk=0;snr=en+10*log(2)/log(10);for y=1: dout_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:dnif dout_2ask(1,i)>0.5 dout_2ask(1,i)=1;else dout_2ask(1,i)=0;endif dout_2ask(1,i)~=d(1,i) ern_2ask=ern_2ask+1;endend dout_2psk=(dout_2psk+1)*0.5;fori=1:dnif dout_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;endendend error_rate_2psk(1,s_n)=ern_2psk/dn/; error_rate_2ask(1,s_n)=ern_2ask/dn/;error_rate_qpsk(1,s_n)=ern_qpsk/dn/;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/n0 dB');ylabel('error_rate dB');grid on;subplot(1,2,2);semilogy(x,error_rate_qpsk,'o-',x,Pqpsk,'*-');legend('qpsk仿真值','qpsk理论值')title('QPSK误码率分析');xlabel('Eb/n0 dB');ylabel('error_rate dB');grid on;;5. 实验结果 / 。