基于卷积编码的扩频通信系统软件平台方案设计书.docx
《基于卷积编码的扩频通信系统软件平台方案设计书.docx》由会员分享，可在线阅读，更多相关《基于卷积编码的扩频通信系统软件平台方案设计书.docx（26页珍藏版）》请在冰豆网上搜索。
基于卷积编码的扩频通信系统软件平台方案设计书
摘要1
前言2
1设计内容及要求3
1.1设计目的3
1.2设计任务3
1.3设计内容3
2设计原理4
2.1数据源4
2.2卷积编码4
2.3m序列发生器4
2.4扩频与解扩4
2.4.1扩频4
2.4.2解扩4
2.5Viterbi译码5
2.6可靠度评估6
3软件仿真及结果分析7
4心得体会8
5附录9
参考文献25
摘要
MATLAB是一种用于科学计算的高性能语言。
它可将计算、可视化和程序设计集成在一个非常容易使用的环境中。
MATLAB是一个交互系统，它的基本数据元素是数组，尤其适合解决用矩阵和向量组织数据的科学技术计算问题。
此次课设就是利用MATLAB软件编写程序实现卷积编码和Viterbi译码的算法，并且实现扩频通信系统的扩频与解扩，从而得出信噪比与误码率之间的关系。
前言
扩频通信是现代通信系统中新的通信方式，它具有较强的抗干扰、抗衰落和抗多径性能，频谱利用率高。
扩频信号是不可预测的、伪随机的宽带信号，其带宽远大于要传输的数据（信息）带宽，同时接收机中必须有与宽带载波同步的副本。
在发端输人的信息先调制形成数字信号，然后由扩频码发生器产生的扩频码序列去调制数字信号以展宽信号的频谱，展宽后的信号再调制到射频发送出去。
在接收端收到的宽带射频信号，变频至中频，然后由本地产生的与发端相同的扩频码序列去相关解扩，再经信息解调，恢复成原始信息输出。
通信系统都要进行3次调制和相应的解调。
一次调制为信息调制。
与一般通信系统比较，多了扩频调制和解扩部分。
扩频通信具备如下特征：
（1）数字传输方式；
（2）传输信号的带宽远大于被传信息带宽；（3）带宽的展宽，是利用与被传信息无关的函数（扩频函数）对被传信息的信元重新进行调制实现的；（4）接收端用相同的扩频函数进行相关解调（解扩），求解出被传信息的数据。
用扩频函数（也称伪随机码）调制和对信号相关处理是扩频通信有别于其他通信的两大特点。
1设计内容及要求
1.1设计目的
锻炼学生综合运用所学知识设计实际通信系统的能力，为学生走向工作岗位处理本专业方向的各种实际问题打好基础，做到学以致用。
具体到本次设计：
（1）学习通信中的纠错编码技术及其主要应用；
（2）掌握扩频通信原理技术及其特点，了解其在通信系统中特别是第三代移动通信系统所采用的码分多址（CDMA）技术中的应用；
（3）学会用C或MATLAB软件方法设计一般通信系统；
（4）学习CDMA手机的信息处理过程，掌握其基带电路主要组成部分和工作原理及功能；模拟CDMA手机的信息处理过程。
完成基带电路主要组成部分的软件平台设计并评估系统的可靠性。
1.2设计任务
（1）建立CDMA手机处理信息的卷积编码扩频通信系统
（2）软件实现卷积编码和Viterbi译码的算法
（3）软件实现扩频通信系统的扩频与解扩
（4）评估卷积编码扩频通信系统的可靠性
1.3设计内容
本次课程设计主要实现的任务有：
建立卷积编码扩频通信系统；软件实现卷积编码和Viterbi译码的算法；软件实现扩频通信系统的扩频与解扩；评估卷积编码扩频通信系统的可靠性。
系统设计的总体框图如图1-1：
m序列发生器
卷积编码
解扩
判决数据
Viterbi
译码
信道
扩频
数据源
图1-1
2设计原理
2.1数据源
随机产生b（t）={b（0）,b
（1）,b
（2）,……}的二进制数据
2.2卷积编码
卷积码（又称连环码），是由伊莱亚斯（P.Elis）提出的一种非分组码。
它把k比特信息段编成n比特的码组，该码组不仅同当前的k比特信息段有关，而且还同前面的（N-1）个信息段有关联（N为大于1的整数）。
通常，把卷积码记作（n,k,N），其中k为输入码元数，n为输出码元数，N为约束长度，表示编码器的存储器级数。
卷积编码属于信道编码，主要用来纠正码元的随机差错，它是以牺牲效率来换取可靠性，利用增加监督位，进行检错和纠错。
卷积码编码器是一个由k个输入端、n个输出端，且具有（N-1）或m节移位寄存器构成的有限状态记忆系统，通常称为时序网络，其原理如图2-1所示。
图2-1
2.3m序列发生器
二进制的m序列是一种重要的伪随机序列，有良好的自相关特性，有时称为伪噪声序列。
m序列是最长线性移位寄存器序列的简称。
顾名思义，m序列是由多级移位寄存器或其延迟元件通过线性反馈产生的最长的码序列。
它能产生的的最大长度的码序列为2^n-1位。
2.4扩频与解扩
将编码输出与扩频码相乘，即完成扩频；将扩频后的信号（叠加有噪声）与扩频码相乘，即为解扩。
2.4.1扩频
扩频通信的基本特点是传输信号所占用的频带宽度（W）远大于原始信息本身实际所需的最小（有效）带宽（DF），其比值称为处理增益Gp，其中Gp=W/DF。
任何信息的有效传输都需要一定的频率宽度，为了充分利用有限的频率资源，增加通路数目，人们广泛选择不同调制方式，采用宽频信道（同轴电缆、微波和光纤等）和压缩频带等措施，同时力求使传输的媒介中传输的信号占用尽量窄的带宽。
因现今使用的电话、广播系统中，无论是采用调幅、调频或脉冲编码调制制式，Gp值一般都在十多倍范围内，统称为“窄带通信”，而扩频通信的Gp值，高达数百、上千，称为 “宽带通信”。
如图图2-2为扩频系统能够一般原理图：
图2-2
2.4.2解扩
正如在一般的窄带通信中，已调信号在接收端都要进行解调来恢复所传的信息。
在扩频通信中接收端则用与发送端相同的扩频码序列与收到的扩频信号进行相关解调，恢复所传的信息。
换句话说，这种相关解调起到解扩的作用。
即把扩展以后的信号又恢复成原来所传的信息。
这种在发端把窄带信息扩展成宽带信号，而在收端又将其解扩成窄带信息的处理过程，会带来一系列好处。
弄清楚扩频和解扩处理过程的机制，是理解扩频通信本质的关键所在。
基本的扩频过程就是在收端产生与发端完全相同的PN码，对收到的扩频信号，在平衡调制器中再一次进行二相相移键控调制。
可以看出发端相移键控调制后的信号在收端又被恢复成原来的载波信号。
当然一个必要的条件是本地的PN码信号的相位必须和收到的相移后的信号在相移点对准，才能正确地将相移后的信号再翻转过来。
2.5Viterbi译码
Viterbi译码算法是一种卷积码的解码算法。
卷积码的Viterbi译码是根据接收码字序列寻找编码时通过网格图最佳路径的过程，找到最佳路径即完成了译码过程，并可以纠正接收码字中的错误比特。
所谓“最佳”，是指最大后验条件概率：
P（C/R）=maxj[P（Cj/R）]，其译码流程为：
（1）根据接收码符号计算出相应的分支量度值；
（2）进入某一状态的分支量度与其前状态路径量度累加求和；
（3）比较到达当前状态的新的路径量度的大小，选择最大者作为新的状态路径量度存储起来，并保存与此路径对应的码字；
（4）对所有状态都实施上述加、比、选运算；
（5）在每一译码时刻，满足延时就从256条留存路径中，选择路径量度最大的一条路径作为译码数输出；
（6）进入下一译码时刻，重复以上步骤，直至译码结束。
2.6可靠度评估
用误码率（BER）指标，计算公式
……
比较两个数据
3软件仿真及结果分析
本次课设需要仿真的内容如下：
（1）直接进行信号传输,即没有进行卷积和扩频的程序一.。
本程序仿真得到的图形见附录中图5-1；
（2）没有进行卷积，进行扩频，扩频周期L=15的程序二。
本程序仿真得到的图形见附录中图5-2；
（3）没有进行卷积，进行扩频，扩频周期L=31的程序三。
本程序仿真得到的图形见附录中图5-3；
（4）进行卷积，进行扩频，扩频周期L=15的程序四。
本程序仿真得到的图形见附录中图5-4；
（5）改变卷积码的值，扩频周期仍为15的程序五。
本程序仿真得到的图形见附录中图5-5；
通过以上图形可以得出的结论如下：
在加入卷积后，系统的抗干扰能力比没有加卷积的时候强；加入扩频比不加扩频的时候抗干扰能力强，并且扩频因子越大，系统的抗干扰能力越强。
在相同条件下，码率越小系统的抗干扰能力越好。
4心得体会
转眼间为期三周的专业方向课程设计就结束了。
这次课程设计也是我们大学学习中的最后一次课程设计了，因此对于这次课程设计我也更加珍惜，比以前更加认真对待。
回想以前的多次课程设计，似乎也觉得这次比以前学到的东西要多。
毕竟这次课程设计是我们专业方向的课设，对专业只是的要求也比较高，不仅让我们能巩固到进入大学以来所学的专业知识，更让我们提前体会到进入社会后从事专业技术工作的状态。
这次课程设计要求我们用到MATLAB软件工具或者是C语言进行软件的编程，来实现课设中所要求的卷积编码和Viterbi译码的算法以及实现扩频通信系统的扩频与解扩。
这就要求我们对MATLAB和C语言能够熟练应用，面对复杂的程序编写刚开始我和我的小组成员有点儿手足无措，不知从哪里入手，但经过老师对课设原理的基本讲解后，我们对课设实现的基本模块也有了初步了解。
在这个基础上我们有从图书馆借来MATLAB的相关书籍来作为参考，也把和《通信原理》课本相关的知识内容进行了巩固。
最后我们针对课设原理的模块各个击破，最后终于把程序完成了，病进行了相关的仿真和调试。
通过这次课设又让我进一步学习到了小组合作的重要性，也学到了在面对困难及问题的时候不能慌张，要仔细认真的从中一步步找错误，急于求成只能让问题更复杂。
在以后毕业的工作中更需要这种合作精神以及严谨认真的态度。
这些对我们以后的发展都很重要。
5附录
程序一：
本程序未加扩频码，未加扩频
EN=input（'请输入信噪比（用分贝表示）:
'）。
L=input（'请输入数据帧长度（=信息比特数+尾比特数,默认值:
400）:
'）。
ifisempty（L）
L=400。
end
ferrlim=input（'请输入程序终止的帧错误数,默认值:
15:
'）。
ifisempty（ferrlim）
ferrlim=15。
end
c=[-1-1-11111-11-111-1-11]。
N=length（c）。
g=[101。
111]。
k0=1。
forj=1:
length（EN）。
errs（j）=0。
nferr（j）=0。
nframe=0。
whilenferr（j）nframe=nframe+1。
tem=rand（1,L）。
fori=1:
L
iftem（i）<=0.5
d（i）=0。
else
d（i）=1。
end
end
n=（1/10^（EN（j）/10））*randn（1,L）。
b=2*d-1。
fori=1:
L
y（i）=b（i）+n（i）。
ify（i）>=0
Db（i）=1。
else
Db（i）=0。
end
end
err=length（find（Db~=d））。
iferr>0
nferr（j）=nferr（j）+1。
end
errs（j）=errs（j）+err。
ifrem（nframe,3）==0|nferr（j）==ferrlim
%Biterrorrate
ber（j）=errs（j）/nframe/（L）。
%Frameerrorrate
fer（j）=nferr（j）/nframe。
%Displayintermediateresultsinproce
升级会员 