直接序列扩频通信系统的误码率仿真.doc

《直接序列扩频通信系统的误码率仿真.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、直接序列扩频通信系统的误码率仿真1．引言扩展频谱通信系统是将基带信号的频谱扩展至很宽的频带上，然后再进行传输的一种通信系统,即将待传送的信息数据用伪随机编码调制，实现频谱扩展后再传输，接收端则采用同样的编码进行解调及相关处理，恢复原始信息数据。扩频通信的基础理论根据信息论中的shannon公式式中，C是系统的信道容量，B是系统信道带宽,N是噪声功率,S为信号的功率,S/N即为信噪比。Shannon公式表明了一个系统信道无误差的传输信息的能力与存在于信道中的信噪比以及用于传输信息的系统信道带宽之间的关系。该公式说明了两个极为重要

2、的概念：一是在一定的信道容量条件下，可以用减少发送信号功率、增加带宽的方法来达到信道容量的要求；另一个是可以采用减少带宽而增加信号功率的方法来达到信道容量的要求。这也就说明了信道容量可以通过带宽与信噪比的互换来保持不变。在实际的工程应用中，改变信号的功率并不容易，相比较而言，扩展信号的带宽更容易操作，所以，要提高信道容量，采用增加信号的带宽比提高信号功率的方法要有效的多。由于扩频通信系统可以在信号功率远低于噪声功率的环境中工作，因此扩频通信系统具有抗干扰能力强，保密性强等优点，在现在通信领域内的应用越来越广泛。2．系统概述本次

3、仿真实验是以MATLAB为仿真平台，信号是8位双极性二进制信号,由1和-1组成。随后对产生的双极性信号进行时域抽样，得到基带信号s，是一组1024位的信息码。伪随机序列由mgen函数产生，共有1024个码元。对已得到的基带信号进行扩频调制，直接把基带信号S与产生的伪随机序列相乘，得到扩频信号。然后对已作扩频处理的信号作BPSK载波调制，得到发射信号。发射信号通过存在高斯白噪声的信道，到达接到端，接收端首先对信号进行解调，将信号解调到中频段，然后将解调的信号与伪随机信号码作乘法运算，对解调信号进行解扩。解扩信号中由于存在高斯白噪

4、声，所有要对噪声进行滤波处理，由于信号为低频基带信号，所以用由汉明窗设计成的FIR低通数字滤波器对解扩信号作低通滤波，产生的信号在经过检验判决，最后还原为原始信号。误码率的计算是由一个比较程序完成的，对判决得到的信号与原始信号进行位对位的比较，计算判决得到的信号中相比原始信号不同的码元的个数，然后将错误码元个数除以总的码元个数，得到在此信噪比下的误码率。3.系统参数信息码元数8位二进制码调制方式BPSK信息码速率128bps扩频码速率8bps采样率128Hz信道参数高斯白噪声信道，信噪比为自定义值4.系统结构扩频通信系统的结构

5、框图如图1所示:信息码产生信息码恢复信道滤波扩频解扩白噪声扩频码产生图1扩频通信系统的系统结构原始信息码是8位双极性二进制信号，然后再通过时域上的抽样，最后变成信息码。扩频码是由函数mgen产生的十六位伪随机信号，然后对随机信号进行抽样，然后再复制加长，最后产生扩频码。mgen函数的程序代码如下:.function[out]=mgen(g,state,N)gen=dec2bin(g)-48;%二进制数转换成相应的数字向量M=length(gen);curState=dec2bin(state,M-1)-48;fork=1:No

6、ut(k)=curState(M-1);a=rem(sum(gen(2:end).*curState),2);curState=[acurState(1:M-2)];end信息码与扩频码进行乘法运算，完成扩频过程。对扩频信号进行载波调制，应用BPSK（二进制相位键控）把扩频信号调制到发射频率上，形成发射信号。发射信号加高斯白噪声，通过函数awgn(signal,snr,’measured’)对扩频信号加入噪声干扰。这个函数的好处是可以通过对snr的改变，实现信号的信噪比的变化。加入噪声后的信号进入接收端，成为接受信号，对接受信

7、号进行BPSK解调处理，然后对调后的信号进行解扩。解扩的方法是用与发射端进行扩频处理时使用的伪随机码完全一样的伪随机序列与解调的信号相乘，将信号解扩。解扩后得到的信号由于仍存在噪声的干扰，因此进行低通滤波，滤波器选用比较常用的FIR数字滤波器。滤波后的信号进行检验判决，最后恢复原始信号。最后，通过对恢复信号码与原始信号码逐位的比较，得到误码的个数，最后经过计算，得出误码率。5．系统工作原理（1）程序代码流程图（2）扩频通信系统的误码率仿真程序如下：code_length=8;x=[11-11-11-1-1];%产生信号码for

8、i=1:8%对信号进行抽样，抽样率为128s(1+(i-1)*128:i*128)=x(i);endn=1:1:1024;subplot(3,2,1)plot(n,s);axis([11024-1.51.5])ylabel('信息码');figure(1)PN_length=