直接序列扩频课程设计文档格式.doc
《直接序列扩频课程设计文档格式.doc》由会员分享,可在线阅读,更多相关《直接序列扩频课程设计文档格式.doc(35页珍藏版)》请在冰点文库上搜索。
![直接序列扩频课程设计文档格式.doc](https://file1.bingdoc.com/fileroot1/2023-4/29/731505b1-3886-4ca4-91df-4eb3b4487aa6/731505b1-3886-4ca4-91df-4eb3b4487aa61.gif)
3.4MATLAB仿真关键程序及实现 15
3.4.1工具设计 15
3.4.2M序列 18
3.4.3系统主程序实现 19
3.5结果分析 22
3.4.1无白噪声通信系统 22
3.4.1高斯白噪声信道系统分析 26
第4章结束语 30
参考文献 31
附件:
32
II
成都学院(成都大学)课程设计报告
第1章绪论
1.1背景
直接序列扩频(DSSS—DirectSequenceSpreadSpectrum)技术是当今人们所熟知的扩频技术之一。
这种技术是将要发送的信息用伪随机码(PN码)扩展到一个很宽的频带上去,在接收端,用与发端扩展用的相同的伪随机码对接收到的扩频信号进行相关处理,恢复出发送的信息。
它是二战期间开发的,最初的用途是为军事通信提供安全保障,是美军重要的无线保密通信技术。
这种技术使敌人很难探测到信号。
即便探测到信号,如果不知道正确的编码,也不可能将噪声信号重新汇编成原始的信号。
有关扩频通信技术的观点是在1941年由好莱坞女演员HedyLamarr和钢琴家GeorgeAntheil提出的。
基于对鱼雷控制的安全无线通信的思路,他们申请了美国专利#2.292.387。
不幸的是,当时该技术并没有引起美国军方的重视,直到十九世纪八十年代才引起关注,将它用于敌对环境中的无线通信系统。
直序扩频解决了短距离数据收发信机、如:
卫星定位系统(GPS)、3G移动通信系统、WLAN(IEEE802.11a,IEEE802.11b,IEE802.11g)和蓝牙技术等应用的关键问题。
扩频技术也为提高无线电频率的利用率(无线电频谱是有限的因此也是一种昂贵的资源)提供帮助。
直序扩频通信系统在发端输入的数字信号信息,先由扩频码发生器产生的扩频码序列去调制数字信号以展宽信号的频谱,扩频码序列一般采用PN码。
展宽后的信号再调制到射频发送出去。
调制多采用BPSK、DPSK、MPSK等调制方式。
在接收端收到的信号进行解调(一般采用相干解调)。
然后由本地产生的与发端相同的扩频码序列去相关解扩。
恢复成原输入的信息输出。
由此可见,—般的扩频通信系统都要进行两次调制和相应的解调。
一次调制为扩频调制,二次调制为射频调制,以及相应的解扩和射频解调。
与一般通信系统比较,扩频通信就是多了扩频调制和解扩部分。
1.2目的和意义
本文所选课题为直接序列扩频通信系统的设计和仿真,直接扩频通信系统具有以下优点和不足。
1)抗干扰性强
抗干扰是扩频通信主要特性之一,比如信号扩频宽度为100倍,窄带干扰基本上不起作用,而宽带干扰的强度降低了100倍,如要保持原干扰强度,则需加大100倍总功率,这实质上是难以实现的。
因信号接收需要扩频编码进行相关解扩处理才能得到,所以即使以同类型信号进行干扰,在不知道信号的扩频码的情况下,由于不同扩频编码之间的不同的相关性,干扰也不起作用。
正因为扩频技术抗干扰性强,美国军方在海湾战争等处广泛采用扩频技术的无线网桥来连接分布在不同区域的计算机网络。
2)隐蔽性好
因为信号在很宽的频带上被扩展,单位带宽上的功率很小,即信号功率谱密度很低,信号淹没在白噪声之中,别人难以发现信号的存在,加之不知扩频编码,很难拾取有用信号,而极低的功率谱密度,也很少对于其他电讯设备构成干扰。
3)直扩通信速率高
直扩通信速率可达2M,8M,11M,无须申请频率资源,建网简单,网络性能好。
4)有很强的保密性能。
对于直扩系统而言,射频带宽很宽,谱密度很低,甚至淹没在噪音中,就很难检查到信号的存在。
由于直扩信号的频谱密度很低,直扩系统对其它系统的影响就很小。
本设计为掌握利用计算机来加深对所学知识的理解和掌握,通过MATLAB/Simulink仿真平台,运用所学的理论和方法进行仿真、解决问题。
探究直接扩频通信系统的性能及痛惜各部分信号的波形和频率特征。
1.3主要内容
(1)对于直接扩频及PSK调制解调基本原理进行研究,实际包括抽样定理、快速傅里叶变换FFT、数字滤波器原理等。
(2)对系统进行抽样建模,并完成各部分功能描述及实现算法研究。
(3)程序编写,包括主系统函数,各部分子系统函数及工具设计。
(4)设置参数,运行程序,对于实验结果进行分析。
第2章直接扩频通信系统的理论基础
2.1直接序列扩频的概念及理论基础
2.1.1直接序列扩频的概念
所谓直接序列(DS:
DirectSequence)扩频,就是直接用具有高码率的扩频码序列在发送端去扩展信号的频谱。
而在接收端,用相同的扩频码序列去进行解扩,把展宽的扩频信号还原成原始的信息。
2.1.2扩频通信的理论基础
长期以来,人们总是想方设法使信号所占的频谱尽量窄,以充分提高十分宝贵的频率资源利用率。
但扩频通信在发送端用扩频码调制,使信号所占的频带宽度远大于所传信息必须的带宽,在接收端采用相同的扩频码进行相关解扩以恢复出所传信息数据。
为什么要用宽频带信号来传输窄带信息呢?
主要是为了通信的安全可靠性。
这可用信息论和抗干扰理论的基本观点来说明:
根据香农(C.E.Shannon)在信息论研究中总结出的信道容量公式,香农公式:
(2-1)
式中:
C--信息的传输速率(信道容量)单位b/s;
S--信号平均功率单位W;
B--频带宽度单位Hz;
N--噪声平均功率单位W。
由式中可以看出:
为了提高信息的传输速率C,可以从两种途径实现,既加大带宽B或提高信噪比S/N。
换句话说,当信号的传输速率C一定时,信号带宽B和信噪比S/N是可以互换的,即增加信号带宽可以降低对信噪比的要求,当带宽增加到一定程度,允许信噪比进一步降低,有用信号功率接近噪声功率甚至淹没在噪声之下也是可能的。
扩频通信就是用宽带传输技术来换取信噪比上的好处。
柯捷尔尼可夫在其潜在抗干扰性理论中得到如下关于信息传输差错概率的公式
(2-2)
此公式指出,差错概率Pe是信号能量E与噪声功率谱密度之比的函数。
设信息持续时间为T,或数字信息的码元宽度为T,则信息的带宽Bm为
(2-3)
信号功率S为
(2-4)
已调(或已扩频)信号的宽度为B,则噪声功率为
(2-5)
将式(4-3)~(4-5)代入式(4-2),可得
(2-6)
上面公式指出,差错概率Pe是输入信号与噪声功率之比(S/N)和信号带宽与信息带宽之比(B/Bm)二者乘积的函数,信噪比与带宽是可以互换的。
它同样指出了用增加带宽的方法可以换取信噪比上的好处。
综上所述:
将信息带宽扩展100倍,甚至用1000倍以上的带宽信号来传输信息,就是为了提高通信的抗干扰能力,即在强干扰条件下保证可靠安全的通信。
这就是扩频通信的基本思想和理论依据。
2.1.3直接序列扩频的基本原理
所谓直接序列扩频(DS),就是直接用具有高速率的扩频码序列在发送端去扩展信号的频谱。
而接收端,用相同的扩频码序列进行解扩,把展宽的扩频信号还原成原始信息。
图2-1示出了直接扩频通信系统的原理方框图
图2-1直扩通信系统的组成框图
图2-2直扩通信系统的主要相位或波形
如图2-2在发送端输入信息码元m(t),它是二进制数据,图中为0、1两个码元,其码元宽度为。
加入扩频解调器,图中为模2加法器,扩频码为一个伪随机码(PN码),记作p(t)。
伪码的波形如图2-2中的第
(2)个波形,其码元宽度为,且取。
通常在DS系统中,伪码的速率远远大于信码速率,即,也就是说,伪码的宽度远远小于信码的宽度,即,这样才能展宽频谱。
模2加法器的运算规则可用下式表示
(2-7)
当m(t)与p(t)符号相同时,c(t)为0;
而当m(t)与p(t)不同时,则为1。
c(t)的波形如图4-2所示中的第(3)个波形。
由图可见,当信码m(t)为0时,c(t)与p(t)相同;
而当信码m(t)为1时,则c(t)为p(t)取反既是。
显然,包含信码的c(t)其码元宽度已变成了,即已进行了频谱扩展。
其扩展处理增益也可用下式表示
(2-8)
在一定的情况下,若伪码速率越高,即伪码宽度(码片宽度)越窄,则扩频处理增益越大。
经过扩频,还要载频调制,以便信号在信道上有效的传输。
图中采用二相相移键控方式。
调相器可由环行调制器完成,即将c(t)与载频相乘,输出为。
即
(2-9)
式中,
(2-10)
因此,经过扩频和相位调制后的信号为
(2-11)
由上面讨论可知,经过扩频调制信号c(t)可看作只取1的二进制波形,然后对载频进行调制,这里是采用调相(BPSK)。
所谓调制,就是指相乘过程,可采用相乘器,环行调制器(或平衡调制器),最后得到的是抑制载波双边带振幅调制信号。
这里假定平衡调制器是理想对称,码序列取+1、1的概率相同,即调制信号无直流分量,这样平衡调制器输出的已调波中,无载波分量。
通过发射机中推动级、功放和输出电路加至天线发射出去。
通常载波频率较高,或者说载波周期较小,它远小于伪码的周期,即满足。
但图2-2中(4)示出的载波波形的宽度为,这是为了便于看清楚一些,否则要在一个期间内画几十个甚至几百个正弦波。
对于PSK来说,主要是看清楚已调波与调制信号之间的相位关系。
图4-2中的第(5)个图为已调波的波形。
这里,当c(t)为一码时,已调波与载波取反相;
而当c(t)为0码时,取同相。
已调波与载波的相位关系如图4-2中的第(6)个图所示。
接收端的工作原理:
假设发射的信号经过信道传输,不出现差错,经过接收机前端电路(包括输入电路、高频放大器等),输出仍为。
这里不考虑信道衰减问题,因为对PSK调制信号而言,重要的是相位问题,这里的假定对分析工作原理是不受影响的。
相关器完成相干解调和解扩。
接收机中的本振信号频率与载波相差为一个固定的中频。
假定收端的伪码(PN)与发端的PN码相同。
接收端本地调相情况与发端相似,这里的调制信号是p(t),即调相器输出信号的相位仅取决于p(t),当p(t)=1时,的相位为;
当p(t)=0时,的相位为0。
信号的相位如图4-2中(7)所示。
相关器的作用在这里可等效为对输入相关器的、相位进行模2加。
对二元制的0、而言,同号模2加为0,异号模2加为。
因此相关器的输出的中频相位如图4-2中的(8)所示。
然后通过中频滤波器。
滤除不相关的各种干扰,经解调恢复出原始信息。
这一过程说明了直扩系统的基本原理和它是怎样通过对信号进行扩频与解扩处理从而获得提高输出信噪比的好处的。
它体现了直扩系统的抗干扰能力。
2.1.4m序列
m序列是最长线性移位寄存器序列的简称。
二进制的m序列是一种重要的伪随机序列,有优良的自相关特性。
容易产生、规律性强,但其随机性接近于噪声和随机序列。
m序列在扩展频谱及码分多址技术中有着广泛的应用,并在m序列基础上还能够成其它码序列,因此无论从m序列直接应用还是从掌握伪随机序列基本理论而言,应该熟悉m序列的产生及其主要特性。
顾名思义,m序列是由多级移位寄存器或其他延迟元件通过线性反馈产生的最长的码序列。
在二进制移位寄存器发生器中,若n为级数,则所能产生的最大长度的码序列为2n-1位。
现在来看看如何由多级移位寄存器经线性反馈产生周期性的m序列。
图2-3为一最简单的三级移位寄存器构成的m序列发生器。
图2-3三级移位寄存器构成的m序列发生器
图中Dl、D2、D3为三级移位寄存器,为模二加法器。
移位寄存器的作用为在时钟脉冲驱动下,能将所暂存的“1”或“0”逐级向右移。
模二加法器的作用为图中(b)所示的运算,即0十0=0,0十1=1,1十0=l,l十1=0。
图(a)中D2、D3输出的模二和反馈为Dl的输入。
在图(c)中示出,在时钟脉冲驱动下,三级移位寄存器的暂存数据按列改变。
D3的变化即输出序列。
如移位寄存器各级的初始状态为111时,输出序列为1110010。
在输出周期为-1=7的码序列后,D1、D2、D3又回到111状态。
在时钟脉冲的驱动下,输出序列作周期性的重复。
因7位为所能产生的最长的码序列,1110010则为m序列。
这一简单的例子说明:
m序列的最大长度决定于移位寄存器的级数,而码的结构决定于反馈抽头的位置和数量。
不同的抽头组合可以产生不同长度和不同结构的码序列。
有的抽头组合并不能产生最长周期的序列。
对于何种抽头能产生何种长度和结构的码序列,已经进行了大量的研究工作。
现在已经得到3-100级m序列发生器的连接图和所产生的m序列的结构。
例如4级移位寄存器产生的15位的m序列之一为111101011001000。
同理我们不难得到31、63、127、255、511、l023…位的m序列。
一个码序列的随机性由以下三点来表征:
1.一个周期内“l”和“0”的位数仅相差1位。
2.一个周期内长度为l的游程(连续为“0”或连续为“l”)占1/2,长度为2的游程占l/4,长度3的游程占l/8。
只有一个包含n个“l”的游程,也只有一个包含(n—1)个“0”的游程。
“l”和“0”的游程数相等。
3.一个周期长的序列与其循环移位序列比较,相同码的位数与不相同码的位数相差l位。
在m序列中一个周期内“1”的数目比“0”的数目多l位。
例如上述7位码中有4个“1”和3个“0”。
在15位码中有8个“l”和7个“0”。
在表3-2中列出长为15位的游程分布。
表2-1111101011001000游程分布
游程长度(比特)
“1”的游程数
“0”的游程数
所包含的比特数
1
2
4
3
游程总数8
合计15
一般说来,m序列中长为R的游程数占游程总数的l/2k。
m序列的自相关函数由下式计算:
(2-12)
A为“0”的位数;
D为“1”的位数
令p=A+D=2n-1则:
(2-13)
设n=3,p=23-1=7,则:
(2-14)
它正是图3-1(d)中所示的二值自相关函数。
m序列和其移位后的序列逐位模二相加,所得的序列还是m序列,只是相移不同而已。
例如1110100与向右移三位后的序列1001110逐位模二相加后的序列为0111010,相当于原序列向右移一位后的序列,仍是m序列。
m序列发生器中移位寄存器的各种状态,除全0状态外,其他状态只在m序列中出现一次。
如7位m序列中顺序出现的状态为111,110,101,010,100,00l和011,然后再回到初始状态111。
m序列发生器中,并不是任何抽头组合都能产生m序列。
理论分析指出,产生的m序列数由下式决定
(2-15)
其中由F(X)为欧拉数(即包括1在内的小于X并与它互质的正整数的个数)。
例如5级移位寄存器产生的31位m序列只有6个。
2.2直扩系统的性能
2.2.1直扩系统的抗干扰性
直扩系统最早应用是在军事通信中作为很强抗干扰性的通信手段。
直扩系统对窄带干扰、宽带干扰等,都具有抗干扰能力,其抗干扰能力大小就是前面提出的扩频处理增益,越大,抗干扰能力就越强。
下面就来分析直扩系统抗宽带干扰和抗窄带干扰的原理
这里的宽带干扰是泛指的与扩频信号不相关的,在CDMA通信网中,其它用户的信号就是一种宽带干扰。
相关处理前,信号频谱是很宽的,经相关处理后,有用信息被解扩,其功率谱集中于信息带宽内,而宽带干扰通过相关器,其功率谱密度基本不变。
由于解扩后必然连接窄带滤波器,保证信号能顺利通过,对信号频带之外的各种干扰起到很大的抑制作用,从而提高了输出的信噪比。
对单频或窄带干扰,直扩系统有很强的抗干扰能力。
图4-4(a)为解扩前的功率谱,窄带干扰功率很大,由于干扰与本地扩频码(PN码)是不相关的。
对干扰来说,相关器起到扩展频谱的目的,功率谱密度就大大下降,其中对信号有害的干扰分量只有落入信息带宽部分,从而抑制了大部分干扰。
由于有用信号能顺利通过窄带滤波器,因此提高了输出的信噪比。
2.2.2直扩信号的抗截获性
截获敌方信号的目的在于:
1.发现敌方信号的存在,
2.确定敌方信号的频率,
3.确定敌方发射机的方向。
论分析表明,信号的检测概率与信号能量与噪声功率谱密度之比成正比,与信号的频带宽度成反比。
直扩信号正好具有这两方面的优势,它的功率谱密度很低,单位时间内的能量就很小,同时它的频带很宽。
因此,它具有很强的抗截获性。
如果满足直扩信号在接收机输入端的功率低于或与外来噪声及接收机本身的热噪声功率相比拟的条件、则一般接收机发现不了直扩信号的存在。
另外,由于直扩信号的宽频带特性,截获时需要在很宽的频率范围进行搜索和监测,也是困难之一。
因此,直扩信号可以用来进行隐藏通信。
至于如何发现敌方直扩信号的存在,和弄清楚其参数,即直扩信号的检测与估值问题。
第3章直接序列扩频通信系统的MATLAB仿真
3.1matlab简介
MATLAB原意为“矩阵实验室—MATrixLABoratory”,它是目前控制界国际上最流行的软件,它除了传统的交互式编程之外,还提供了丰富可靠的矩阵运算、图形绘制、数据和图象处理、Windows编程等便利工具。
MATLAB还配备了大量工具箱,特别是还提供了仿真工具软件SIMULINK。
MATLAB在80年代一出现,首先是在控制界得到研究人员的瞩目。
随着MATLAB软件的不断完善,特别是仿真工具SIMULINK的出现,使MATLAB的应用范围越来越广。
随着MATLAB库函数和仿真工具箱的不断扩充,使其在系统仿真与分析、信号处理、图像处理等方面的应用越来越广泛。
MATLAB具有3大特点:
(1)功能强大,包括数值计算和符号计算,计算结果和编程可视化,数学和文字统一处理,离线和在线皆可计算;
(2)面友好,语言自然。
MATLAB以复数矩阵为计算单元,指令表达与标准教科书的数学表达式相近;
(3)开放性强。
MATLAB有很好的可扩充性,可以把它当作一种更高级的语言去使用,可容易地编写各种通用或专用应用程序。
正是由于MATLAB的这些特点,使它获得了对应用学科(特别是边缘学科和交叉学科)的极强适应力,并很快成为应用学科计算机辅助分析设计、仿真、教学乃至科技文字处理不可缺少的基础软件,成为欧美高等院校、科研机构教学与科研必备的基本工具。
本文就是利用编制MATLAB仿真m文件进行DS扩频通信系统的性能仿真及其分析。
3.2设计思路及模型抽象
图3-1设计思路流程图
本文在进行系统设计时,首先参考移动通信教科书相关理论知识对直接扩频通信系统进行仔细研究,划分出PN码生成子系统,PSK调制子系统等五大子系统,然后对各子系统的功能及实现原理进行细致研究选出合适的算法完成各系统模块的编程。
接着利用各子系统完成系统搭建,运行仿真系统,对结果进行分析,对各项参数和指标进行评估。
以判断系统是否符合要求。
图3-2直接序列扩频通信系统模型
经过研究,为了完成本设计要求,根据通信系统基本模型和本设计要求完成系统模型抽象,主要包含信源(随机数字信号发生器)、信道编码器、调制器、信道(无噪声信道和高斯白噪声信道两种)、解调器、信道解码器、信宿等八个子系统。
为了观察仿真结果需要数字信号示波器、模拟信号示波器和频谱仪几大反馈子系统。
3.3系统仿真流程图
根据以上分析现作出通信系统主要流程图:
图3-3直接扩频通信系统流程图
本系统中信源为二进制随机信号生成器,信号频率为fc。
PN码为m码生成器,采用寄存器移位模二加运算生成,信号频率为fm=100*fc。
PSK载波为fp=fm*100的余弦信号,通过与双极性基带信号进行乘法运算实现PSK调制。
系统中抽样频率为fs=64*fp。
3.3.1发射机
发射机部分主要包含信道编码器和调制器两个部件。
信道编码器包含PN序列生成,源码抽样和编码三个步骤。
调制器包含载波生成,基带信号极性变换和调制三个步骤。
其模型如下:
图3-4发射机模型
3.3.1接收机
接收机机部分主要包含解调器和信道解码器两个部件。
解调器包含带通滤波,载波生成,解调,低通滤波以及抽样判决等几个步骤。
信道解码器包含PN序列生成,解码和抽样判决三个步骤。
接收机最终输出的信号应该与源信号完全一致。
图3-5接收机模型
3.4matlab仿真关键程序及实现
3.4.1工具设计
1)模拟信号示波器
本设计中模拟示波器较为简单主要是根据序列长度求出及序列频率求出其横坐标时间值。
然后利用plot函数绘出其波形。
图3-5模拟信号示波器程序流程图
程序如下:
functionplot_double(m,s,fs)%m为输入序列,s为图片名,fs为输入序列频率
N=length(m);
%计算序列长度
n=0:
1:
N-1;
t=n/fs;
%求出横坐标
%绘波形图
plot(x,m,'
k'
);
xlabel('
时间t'
ylabel('
幅度h'
title(s);
2)数字示波器
本数字示波器主要利用matlab中plot函数完成波形绘制,但是需要对横纵坐标进行坐标变换方可画出其波形。
其流程图如下。
图3-6数字示波器流程图
坐标变换原理如下:
常规plot函数横纵坐标序列如下。