基于Matlab的CDMA通信系统仿真

基于Matlab的CDMA通信系统仿真-1-通信系统综合设计与实践1绪论1.1课题背景及目的20世纪60年代以来，随着民用通信事业的发展，频带拥挤问题日益突出。CDMA（CodeDiveiionMultipleAcce,码分多址）通信，在使用相同频率资源的情况下，理论上CDMA移动网比模拟网容量大20倍，实际使用中比模拟网大10倍，比GSM要大4〜5倍，所以在通信领域中起着非常重要的作用。CDMA的基本原理是利用互相正交（或尽可能正交）的不同编码，分配给不同用户调制信号，实现多用户同时使用同一频率接入系统和网络的通信。由于利用互相正交（或尽可能正交）的编码去调制信号，会将原信号的频谱带宽扩展，因此，这种通信方式，又称为扩频通信。本论文所完成的CDMA通信仿真系统，是结合CDMA的实际通信情况，利用MATLAB的通信工具箱一SIMULINK组建出完整的CDMA通信系统，完成整体设计方案，实现完整的发送到接收的端到端的CDMA无线通信系统的建模、仿真和分析。教学实践表明，该系统的完成使得比较抽象的概念得以直接表示，烦琐的计算得以大大简化，提高上机效率，在通信原理课程教学中起到良好的辅助作用。1.2课题研究方法为了研究CDMA通信系统的通信方式，我们对两种扩频码（m序列和正交gold序列）经过衰落信道后再解扩，通过比较两种扩频码的误比特率与信噪比的关系得出用来扩频的PN码哪种更好。使其更符合CDMA通信的抗干扰能力强的要求和实现多用户同时在同一频率互不干扰进行通信而误比特率性能不随着用户数的增加而恶化这样的目的进行仿真实验。2-通信系统综合设计与实践2CDMA基础及原理口CDMA多址技术的原理是基于扩频技术，即将需传送的具有一定信号带宽的信息数据，用一个带宽远大于信号带宽的高速伪随机码进行调制，使原数据信号的带宽被扩展，再经载波调制并发送出去。接收端由使用完全相同的伪随机码，与接收的带宽信号作相关处理，把宽带信号换成原信息数据的窄带信号即解扩，以实现信息通信。2.1扩频通信扩频通信技术是一种信息传输方式，其信号所占有的频带宽度远大于所传信息所需的最小带宽；频带的扩展是通过一个独立的码序列来完成，用编码及调制的方法来实现的，与所传信息数据无关；在接收端则用同样的码进行相关同步接受、解扩及恢复所传信息数据。2.1.1扩频通信理论基础香农公式：C=Wlog2（1+S/N）口1、在给定的传输速率C不变的条件下，频带宽度W和信噪比S/N是可以互换的。即可通过增加频带宽度的方法，在较低的信噪比情况下，传输信息。2、扩展频谱换取信噪比要求的降低，正是扩频通信的重要特点，并由此为扩频通信的应用奠定了基础。2.1.2扩频通信系统的分类（1）直接序列扩频（DS）口⑵跳频扩频（FH）（3）跳时扩频（TH）□（4）混合方式（以上三种基本方式的不同组合）3-通信系统综合设计与实践在实际的CDMA系统中，直接序列扩频得到了广泛的认可和应用，所以，在本次实验中主要研究直接序列扩频技术。.2直接序列扩频DS-SS是直接用具有高码率的扩频码序列在发送端去扩展信号的频谱。而在收端，用相同的扩频码序列去进行解扩，把展宽的扩频信号还原成原始的信息。扩频码发生器载波发生器载波发生器解扩码发生器二进制信息扩频调制数据调制信道数据解调解扩恢复信息图2-1直接序列扩频的发射机和接收机框图.3扩频码序列在扩频系统中，信号频谱的扩展是通过扩频码实现的，扩频系统的性能与扩频码的性能有很大的关系，对扩频通常提出下列要求：易于产生；具有随机性；扩频码应该具有尽可能长的周期，使干扰者难以从扩频码的一小段中重建整个码序列；扩频码应该具有良好的自相关和互相关特性，以利于接收时的捕获和跟踪，以及多用户检测等。扩频码中应用最多的是M序列，又称最大长度序列，还有GOLD序列、WALSH码序列等，本次试验中我们主要运用了M序列和正交GOLD序列。2.3.1PN码口CDMA信道的区分是靠PN码来进行的，因而要求PN码自相关性要好，互相关性要弱，实现和编码方案简单等。目前的CDMA系统就是采用一种基本的PN序列m序列作为地址码，利用它的不同相位来区分不同用户。口M序列是最长线性移位寄存器序列的简称，它通常是由反馈移位寄存器产生-4-通信系统综合设计与实践的具有像随机噪声波形的周期性二元序列。Gold码序列是一种基于m序列的码序列。Gold码是m序列的组合码，由同步时钟控制的两个m序列逐位模2加得到。这两个码发生器的周期相同，速率也相同，因而两者保持一定的相位关系，这样产生的组合码与这两个子码序列的周期也相同。正交Gold码是在优选对产生的Gold码末尾加0，使序列的长度为偶数。3.1DS-CDMA通信系统原理图口图3-1DS-CDMA通信系统原理图口在实验中，我们采用；两种扩频码进行仿真。M序列和正交Gold序列，扩频后的数据通过脉冲成型滤波器后通过信道同时到达接收端，在接收端分别对不同用户信息数据进行解扩，恢复各个用户的原始信息。-5-通信系统综合设计与实践-6--6-3.2仿真流程设计与介绍3.2.1仿真设计步骤1、m序列的DS-CDMA在AWGN下的性能的仿真。2、正交Gold序列的DS-CDMA在AWGN下的性能的仿真。口3、对比两个扩频码在AWGN信道下的性能曲线，分析哪种扩频码更适合在AWGN信道中传输，传输衰减最小，对比m序列和正交Gold序列的抗干扰能力。4、为了证明这一结论，把两种PN码再经过Rayleigh衰落信道下的性能进行验证。（代码见附录）3.2.2程序编程设计最外层循环的时候是信噪比，根据每种信噪比下，首先产生各个用户的发射数据，根据用户数，每一行是每个用户的发射数据，然后进行Gray编码、4-QAM调制，调制完后，对信号进行扩频，扩频部分是通过pread函数完成的，是完成扩频的功能。扩频时，看输入参数是否满足要求，满足要求后，用扩频码与每一个数据符号相乘，对原始符号进行扩展，完成扩频。扩频完成后，再通过脉冲成形滤波器计算每个用户信号功率。当用户数大于1时，所有用户数据相加。相加完成后，如果需要通过瑞利衰落信道，再让信道系数与发射数据相乘。接收端加入高斯白噪声，首先根据信噪比计算高斯白噪声标准差，根据标准差加入白噪声。如果通过瑞利衰落信道后，我们还需要进行一个信道补偿。最后通过脉冲成形滤波器进行降采样，之后进行数据解扩，解扩部分是通过depread函数完成的，解扩完成后，对信号进行4-QAM解调、Gray编码逆映射，分别对用户原始数据和解扩后数据进行比较，最后统计误比特率。以上就是完成扩频的主程序。M序列和正交Gold序列脚本程序的介绍：用户数统一定义的是1、4、7,扩频码分别选择的是M序列和正交Gold序列，然后在每种用户数下仿真扩频系统的性能，最后画出仿真得到的误比特率。通信系统综合设计与实践4仿真结果分析实验仿真运行结果100m序列DS-CDMA在AWGN信道下的性能uer=1uer=4uer=710误比特率(BER)-^0-2100-3123456信噪比EbNo(dB)78910图4TM序列DS_CDMA在AWGN下的性能口从图4-1中可以看出，由于m序列即使在完全同步时，之间的互相换值也不为0，因此随着用户数的增加(解调其中一个用户的数据时，其他用户会对该用户的解调产生干扰)，干扰越来越大，导致系统的误码率性能下降。通信系统综合设计与实践100正交Gold序列DS-CDMA在AWGN信道下的性能uer=1uer=4uer=710-1误比特率(BER)10-210-3100-412345信噪比EbNo(dB)6789图4-2正交Gold序列在AWGN下的性能口从图4-2中可以看出，由于正交Gold序列在完全同步时，它们的互相关值为0，因此其BER性能并不随着用户数的增加而恶化。在信噪比为2、4、6时误比特率基本相同。由图可以看出来当误比特率仿真时候小于10时结果就会有一定的偏差。但还是可以看出正交Gold序列比M序列的性能要好。所以，初步得出结论：正交Gold序列在AWGN信道下的抗干扰能力更强，衰减比,序列要小。为证明这一结论进行了下面的实验(图4-3、图4-4)。-3通信系统综合设计与实践100m序列DS-CDMA在Rayleigh信道下的性能uer=1uer=4uer=7误比特率(BER)10-^00-2123456信噪比EbNo(dB)78910图4-3M序列在Rayleigh衰落信道下的性能口对比图4-1和图4-3,显然M序列在瑞利衰落信道下的性能要比能过AWGN信道下的性能要差。口通信系统综合设计与实践100正交Gold序列DS-CDMA在Reyleigh信道下的性能uer=1uer=4uer=7误比特率(BER)10-1100-2123456信噪比EbNo(dB)78910图4-4正交Gold序列在Rayleigh衰落信道下的性能口正交Gold序列在的瑞利衰落信道下的性能也要比AWGN信道下的性能要差，可以看出，DS-CDMA系统在Rayeligh衰落信道下的性能要比AWGN信道下的性能差。在图4-3上可以看出，随着用户数的增加，其误比特率还是会增高，但是由图4-4可以看出，信号的误比特率甚至由于用户的增多而降低，性能变的更优，此次实验，可能存在些许误差，而且为了使程序更快运行出结果我们给他加了信道补偿，也可能会由于描点过少，造成图像与真实有些不符，但从大体方向上仍然能够得出结论，即使在瑞利衰落信道下正交Gold序列的性能依然优于m序列。口通信系统综合设计与实践总结对于本仿真而言，已基本完成了设计任务书的设计和研究目的，通过用MATLAB对DS-CDMA系统的仿真调试、结果分析，让我组熟悉YDS-CDMA的工作原理，加深了对扩频通信的认识，并深刻的了解PN码在不同信道的衰落情况。通过仿真结果中波形的直观方式，更让我们了解到了系统衰减的规律。但是在本设计中也存在着缺陷与不足1、在本设计中过程中可能由于程序过于复杂，信息本身所占用的带宽偏大，用来传输信息的带宽相对不够大，造成结果并不明显，信息衰减严重。2、在设计中只考虑到加性高斯信道所带来的干扰，在实际通信信道却是复杂多变，存在着各种各样的情况，所以最后的接收信号是在很简单的干扰下得出。要想应用于实际中，必须加入各种噪声来考虑，以实现真实系统的设计。3、实验中，我们选择加大传输带宽，来弥补其他干扰带来的影响，可能会造成消耗过大，浪费带宽的结果，实际应用中是不实用的。只能作为仿真，进行研究。即使如此，在本次设计的整个过程中，以上的结果已经令我们受益匪浅。通信系统的性能分析和仿真，随着通信技术、信息技术和计算机技术的发展以及网络系统的大量应用，显得越来越重要。通信系统综合设计与实践致谢在本次仿真设计中，感谢指导老师在论文选题、设计任务书中对我的无私帮助，加深了我对专业基础理论的理解，拓宽了我的专业知识面，实在是受益匪浅。在此，谨对所有在论文写作中帮助过我的老师、同学表示衷心的感谢和由衷的敬意！谢谢大家！通信系统综合设计与实践参考文献［1］刘学勇,编著.通信系统建模与仿真（电子工业出版社）［2］邓薇，编著^人1人18函数速查手册（人民邮电出版社）口［3］【美】亚〃廿@山^丫」00,著.移动通信工程理论和应用（第二版）（人民邮电出版社）［4］王华奎，李艳萍等编著.移动通信原理与技术（清华大学出版社）通信系统综合设计与实践附录(各部分编程)：1、M序列口function[mout]=meq(n,tap,inidata,num)0/甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘时时时时时时时时时时时时时时时时时时时时时时时时时时时时时时时时某某某%n：m序列的阶数n%tap：反馈寄存器的连接位置％inidata:寄存器的初始值序列%num:输出的m序列的个数口%mout:输出的m序列，如果num>1,则每一行为一个m序列口0/甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘时时时时时时时时时时时时时时时时时时时时时时时时时时时时时时时时某某某ifnargin<4num=1;endmout=zero(num,2'T);fpo=zero(n,1);口fpo(tap)=1;forii=1:2八n-11mout(1,ii)二inidata(n);%寄存器的输出值temp=mod(inidata某fpo,2);%计算反馈数据口inidata(2:n)=inidata(1:n—1);%寄存器移位一次inidata(1)=temp;%更新第1个寄存器的值endifnum>1%如果要输出多个m序列，生成其他m序列口forii=2:nummout(ii,:)=hift(mout(ii-1,:),1);通信系统综合设计与实践endend关于hift函数的代码：口function[outregi]=hift(inregi,hiftr)0/甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘/O^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘某某某％inrege：输入序列口%hiftr:循环右移的位数口%outregi:输出序列口0/甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘/O^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘某某某v=length(inregi);outregi=inregi;hiftr=rem(hiftr,v);ifhiftr>0outregi(:,1:hiftr)=inregi(:,v-hiftr+1:v);%循环移位outregi(:,1+hiftr:v)=inregi(:,1:v-hiftr);eleifhiftr<0outregi(:,1:v+hiftr)=inregi(:,1-hiftr:v);outregi(:,v+hiftr+1:v)=inregi(:,1:-hiftr);end%某某某某某某某某某某某某某某某某某某某某某endoffile某某某甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘-15-通信系统综合设计与实践2、Gold序列：口function[gout]=goldeq(m1,m2,num)%某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某某%由1:田序列1%m2:m序列2口%num:生成的Gold序列个数%8。口已生成的Gold序列输出口0/甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘时时时时时时时时时时时时时时时时时时时时时时时时时时时时时时时时某某某ifnargin<3%如果没有指定生成的Gold序列个数，默认为1num=1;endgout=zero(num,length(m1));forii=1:num%根据Gold序列生成方法生成Gold序列gout(ii,:)二某or(m1,m2);m2=hift(m2,1);end%某某某某某某某某某某某某某某某某某某某某某某endoffile某某甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘其hift函数代码同M序列口-16-通信系统综合设计与实践3、仿真部分主程序：%直接序列扩频主程序代码function[ber]=dcdma(uer,eq)%uer:同时进行扩频通信的用户数口%eq：扩频码1：M-序列2:Gold序列3:正交Gold序列%ber：该用户数下的误码率%某某某某某某某某某某某某某某某某某某某某某某某某某某初始化次公甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘r=25600.0;%符号速率口nSymbol=1000;%每种信噪比下发送的符号数M=4;%4-QAM调制br=r某log2(M);%比特速率graycode=[0132];%Gray编码规则口%某某某某某某某某某某某某某某某某某某某某某某脉冲成形滤波器余料甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘delay=10;%升余弦滤波器时延F=8;%滤波器过采样数rolloff=0.5;%升余弦滤波器滚降因子rrcfilter=rcoine(1,F,'fir/qrt',rolloff,delay);%设计根升余弦滤波器%某某某某某某某某某某某某某某某某某某某某某某扩频码产生参数甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘%uer二uerl;%用户数tage=3;%m序列的阶数ptapl=[13];%m序列1」1的寄存器连接方式ptap2=[23];%m序列2的寄存器连接方式regi1=[111];%m序列1的寄存器初始值regi2=[111];%m序列2的寄存器初始值%某某某某某某某某某某某某某某某某某某某某扩频码的生成某某某甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘witcheqcae1%M-序列code=meq(tage,ptap1,regi1,uer);口cae2%Gold序列m1=meq(tage,ptap1,regi1);m2=meq(tage,ptap2,regi2);code=goldeq(m1,m2,uer);EbNo=0:2:10;%Eb/No变化范围口-17-通信系统综合设计与实践cae3%正交Gold序列m1=meq(tage,ptap1,regi1);m2=meq(tage,ptap2,regi2);code=[goldeq(m1,m2,uer),zero(uer,1)];endcode=code某2T;clen=length(code);口%某某某某某某某某某某某某某某某某某某某某某某某某某某衰落信道参数某某某某某某某某某某某某某某某某某某某某某某某某某某t=1/F/r/clen;%信道采样时间间隔口t=(0:nSymbol某F某clen-1+2某delay某F)某t;%每种信噪比下的符号传输时间y=160;%多普勒频移[Hz]%h=rayleigh(fd,t);口0/甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘代/O^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^KI/J吉XC心甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘真开始某某某某某某某某某某某某某某某某某某某某某某某某某某某某forind某=1:length(EbNo)ind某口%某某某某某某某某某某某某某某某某某某某某某某某某某某发射端甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘data二randrc(uer,nSymbol,[0:3]);%产生各个用户的发射数据data1=graycode(data+1);%Gray编码data1=qammod(data1,M);%4-QAM调制卜口土]=pread(data1,code);%扩频口out1=rcoflt(out.',r,F某r,'filter',rrcfilter);口%通过脉冲成形滤波器pow=um(ab((out1))」2)/nSymbol;%计算每个用户信号功率ifuer>1%用户数大于1时，所有用户数据相加outl二um(out1.');eleout1=out1.';end%某某某某某某某某某某某某某某某某某某某某通过瑞利衰落信道某甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘某%0口土1=h.某outl;口%某某某某某某某某某某某某某某某某某某某某某某某某某接收端某甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘igma=qrt(0.5某pow某r/br某10X—EbNo(ind某)/10));%根据信噪比计算高斯白噪声方差y二口；口-18-通信系统综合设计与实践forii=1:uery(ii,:)=out1+igma(ii).某(randn(1,length(out1))+j某randn(1,length(out1)));%加入高斯白噪声(AWGN)口%y(ii,:)=y(ii,:)./h;%假设理想信道估计end口y二rcoflt(y.',r,F某r,'F/filter',rrcfilter);口%通过脉冲成形滤波器进行滤波y二downample(y,F);%降采样forii=1:uery1(:,ii)=y(2某delay+1:end-2某delay,ii);end口yd二depread(y1.',code);%数据解扩demodata=qamdemod(yd,M);%4-QAM解调demodata=graycode(demodata+1);%Gray编码逆映射口[err,ber(ind某)]二biterr(data,demodata,log2(M));%统计误比特率end口主程序中pread函数的代码：口%扩频函数function[out]=pread(data,code)0/甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘某某某%data:输入数据序列%。^©:扩频码序列口%out:扩频后的输出数据序列口0/甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘/0^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K7rC甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘某某某witchnargincae{0,1}%如果输入参数个数不对，提示错误error('缺少输入参数');end-19-通信系统综合设计与实践[hn,vn]=ize(data);[hc,vc]=ize(code);ifhn>hc%如果扩频码数小于输入的待扩频的数据序列，提示错误error('缺少扩频码序列')；end口out=zero(hn,vn某vc);口forii=1:hnout(ii,:)=rehape(code(ii,:).'某data(ii,:),1,vn某vc);end口0/甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘/0^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K7rCendoffile某某某某某某某某某某某某某某某某某某某某某某某某某口解扩函数depread的代码：口%信号解扩functionout=depread(data,code)0/甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘/O^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘7K7K7K7K7K7K7K7K7K7K7K7K7K7K7K7K7K7K7K7K7K7K7K7K7K7K7K7K7K7K7K7K某某某%data:输入数据序列口%code:解扩使用的扩频码序列口%out:解扩后的输出数据序列口0/甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘/O^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K^K甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘某某某witchnargin%如果输入参数个数不对，提示错误cae{0,1}口error('缺少输入参数')；end口[hn,vn]=ize(data);[hc,vc]=ize(code);out=zero(hc,vn/vc);forii=1:hc某某二rehape(data(ii,:),vc,vn/vc);out(ii,:)=code(ii,:)某某某/vc;end0/甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘甘/0^K^K^K^K^K^K^K^K^K