基于matlab产生gold序列课程设计报告.docx
《基于matlab产生gold序列课程设计报告.docx》由会员分享，可在线阅读，更多相关《基于matlab产生gold序列课程设计报告.docx（28页珍藏版）》请在冰豆网上搜索。
基于matlab产生gold序列课程设计报告
专业课程设计报告
题目：
基于matlab产生gold序列
姓名：
朱维娜
专业：
通信工程
班级学号：
09042106
同组人：
尧国振
指导教师：
夏思满
南昌航空大学信息工程学院
2021年6月29日
专业课程设计任务书
2010－2011学年第2学期 第17周－20周
题目
基于Matlab产生Gold序列
内容及要求
前提：
掌握PN序列的相关知识，掌握Gold序列的产生原理
设计要求：
（1）使用Matlab的m文件，通过编程生成任意长度的Gold序列（提示：
Gold序列由两个不同m序列模二加得到）；
（2）对以上特定的Gold序列，设计m文件，分析该序列的相关性。
。
进度安排
17周：
查找资料，进行系统软件方案设计；
18周：
软件的分模块调试；
19周：
系统联调；
20周：
设计结果验收，报告初稿的撰写。
学生姓名：
指导地点：
E楼610室
任务下达
2011年6月13日
任务完成
2011年7月8日
考核方式
1.评阅□ 2.答辩□□□
指导教师
夏思满
系（部）主任
摘要
Gold序列是R·Gold提出的一种基于m序列的码序列，这种序列有较优良的自相关和相互关特性，构造简单，产生的序列数多，因此取得了普遍的应用。
本文第一介绍了扩频通信中伪随机序列性质，在介绍伪随机中经常使用的m序列和Gold序列码产生的方式原理和性质，先用matlab编程产生m序列，在用模二加产生gold序列，再运用Matlab对Gold自相关和相互关进行了仿真分析。
且Matlab在编程效率、可读性、可移植性与可扩充性上，远远优于其他高级编程语言，是公认的最优秀的科技应用软件。
Gold序列能够用软件也能够用硬件二种方式实现，可是通过本次设计能够看见软件设计的许多优势
关键词：
伪随机序列；Gold序列；m序列；Matlab仿真
第一章设计要求和引言
1.设计内容
前提：
把握PN序列的相关知识，把握Gold序列的产生原理
设计要求：
（1）利用Matlab的m文件，通过编程生成任意长度的Gold序列（提示：
Gold序列由两个不同m序列模二加取得）；
（2）对以上特定的Gold序列，设计m文件，分析该序列的相关性。
1.1.2设计流程框图
图1-1实验设计流
1.2研究的背景及意义
移动通信由于具有时实性、机动性、具有不受时空限制等特点，己经成为一种深受人们欢迎的通信方式，并融入了现代生活当中。
自美国Qualcomm公司提出在蜂窝移动通信系统中应用码分多址（CodeDivisionMultipleAccess，简称CDMA）技术的系统实现方案至今，CDMA通信系统相关于其它无线通信系统在客户容量和高质量的优势愈来愈显现出来。
在短短的二、三十年中，移动通信系统已从第一代的模拟蜂窝系统进展到第二代全世界数字移动蜂窝系统（2G），目前己经开始向第三代宽带多媒体蜂窝系统（3G）进展，而且处于第二代和第三代之间的已经趋于成熟。
尽管第二代移动通信系统中，GSM系统仍占有专门大的市场份额。
可是，因为具有伪随机编码调制和信号相关处置两大特点而使CDMA通信方式具有抗干扰、抗噪音、抗多径衰落、能在低功率谱密度下工作、有保密性、可多址复用和任意选址、可高精度测量等优势，使CDMA技术成为第三代移动通信和个人通信系统的核心技术，以扩频理论为基础的CDMA技术已成为当前移动通信领域的研究热点。
在CDMA系统的众多用户都工作在同一时刻同一频段内，系统给各个用户分派一个唯一的扩频码来进行频谱的扩展，在发送和接收时，系统更是利用各地址码之间的相互关特性值来区分不同的用户。
因此，扩频码的特性直接阻碍到CDMA系统的捕捉同步性能、抗干扰性能和多址能力。
从理论上说，独立、均匀散布的随机序列是扩频码的理想模型，但是它由于不易产生、无法时实分发等缺点而被以为难以在实际的CDMA系统中应用。
CDMA自其理论提出到投入商业营运、直至称为第三代移动通信系统的核心技术，一直是通信领域的关注热点。
作为CDMA的基础技术之一的PN码的选择和产生也是倍受业内人士关注的，如何找到易生成且相关特性好的PN码成为研究人员追求的目标之一。
为此，人们设计了各类确信性的伪随机序列来代替随机序列作为扩频码。
迄今为止，世界各国的学者在伪随机序列的设计与选择方面己做了大量的工作，例如，由m序列优选对生成的Gold序列己被用作第三代移动通信系统中WCDMA的扩频码；和通过对m序列添加一个全“0”状态取得的M序列和m序列也已被用作第三代移动通信系统中CDMA2000的扩频码。
m序列、Gold序列等线性序列多由线性移位寄放器所产生，有易于实现、具有较好的相关特性等优势。
实际应用的CDMA通信系统采纳复合扩频技术，即用正交码（Walsh函数序列，OVSF码族）作为信道化码来区分小区、用Gold序列或M序列作为扰码来区分用户。
因此，本文所研究的Gold序列，在扩频通信系统中发挥着重要的作用，通过Matlab仿真，对其自相关性能进行分析，能够更好的明白得CDMA系统的通信原理。
1.3CDMA通信技术简介
1.3.1扩频的理论基础
在信息论中，关于持续信道，若是信道带宽为B，且受到加性高斯白噪声干扰，那么其信道容量的理论公式（香农公式）为：
（1-2）
其中C—信道容量，单位bit/s;B一信道带宽，单位Hz；S一信号平均功率，单位w；N一噪声平均功率，单位w。
从香农公式可知：
（1）要增大信息传输速度，就必需增大信道带宽B或信噪比S。
由于公式中对数部份转变得比较缓慢，因此增加B比增加S加更有效，也确实是说若是传输信号的带宽变窄，将致使信号功率的大幅提高。
而若是通过增加带宽去换取信号功率的减小，就能够节省较大的信号功率能源。
即B增加时，信道容量增加较快。
（2）当信道容量为常量时，信道带宽与信噪比存在互换关系。
在C恒定的情形下，能够通过减少发送功率，增加信道带宽的方式维持信道容量不变的目标。
也能够通过减小带宽，增强信号功率的方式。
信道容量能够通过带宽与信噪比的互换而维持不变。
（3）当带宽增加到必然程度时，信道容量也不能无穷增加。
这是因为噪声功率N=n0B，当信道带宽B增加时，N也随着增加，因此C有一个极限值。
扩频通信是指系统所传输的信号（带宽为Bm）被扩展至一个很宽的频带Bc。
用来传输信息的信号带宽远远大于信息本身带宽的一种通信方式。
它利用高速度的扩频码来达到扩展传输信号的带宽，从而减小了发送功率。
对扩频通信来讲Bc/Bm的值一样为100～1000。
1.2.2扩频通信的分类
（1）直接序列（DS）扩频系统：
用一组高速数字编码序列直接扩展频谱，由于编码序列的带宽远远大于原始信号的带宽，从而扩展了发射信号的频谱。
（2）跳频（FH）扩频系统：
使发射机频率在一组预先制定的频率上依照编码序列所规定的顺序离散的跳变，从而扩展发射波的频谱。
一样来讲，跳频图案由伪随机码操纵，从而使载频的跳变具有均匀散布的性质。
（3）线性调频（Chirp）系统：
在这种系统中，载频在一给定的脉冲时刻距离内线性的扫过一个宽的频带，从而扩展发射波的频谱。
（4）跳时（TH）扩频系统：
这种系统与跳频系统类似，区别在于一个是操纵频率，而另一个是操纵时刻，即TH系统是用伪随机码操纵发射时刻和时刻的长短。
（5）混合系统：
前述几种方式的某种形式的组合，如DS/FH系统、DS/TH系统、FH/TH系统、DS/TH/FH系统等。
目前有效的扩频通信中，以直接序列扩频系统应用的比较多。
而CDMA通信系统确实是基于扩频技术的无线通信系统。
1.2.3CDMA扩频通信系统的组成
CDMA通信系统是最具代表性的扩频通信技术应用，它的大体工作方式有直接序列扩频（DirectSequenceSpreadSpectrum，简称DS）方式、跳变频率（FrequencyHopping，简称FH）方式和跳变时刻（TimeHopping，简称TH）方式三种。
其中，直扩（DS）方式同另外两种方式比较，实现频谱扩展方便，不管对通信、测距应用仍是其它应用都很适合，因此在目前利用的最多，也是最典型的一种扩频通信方式。
CDMA扩频通信系统包括两个大体技术：
一个是码分技术，其基础是扩频技术；另一个是多址技术。
目前的CDMA系统确实是采纳m序列及由其产生的其它PN序列作为地址码，利用它们的不同相位来区分不同用户。
在第三代移动通信系统中别离采纳了m序列、Gold序列及M序列作为地址编码，用Walsh序列作为信道编码。
因此，PN码的选择直接阻碍到CDMA系统的容量、抗干扰能力、接入和切换速度等性能。
CDMA信道的区分也是靠PN序列来进行的，PN序列较好的相关特性—自相关特性尖锐，相互关特性较弱，加上实现和编码方案简单等特点，使其在以后的移动通信系统中处于相当重要的位置。
1.2.4伪随机序列在CDMA通信系统中的应用
CDMA通信系统中的扩频码采纳三层结构。
底层是信道码，通常采纳正交码，CDMA2000标准给出的是码长为64的Walsh正交码，3GPP标准给出的是正交可变扩频因子码序列（OrthogonalVariableSpreadingFactorCode，简称OVSF码），用来区分不同的CDMA信道。
第二层是基站码，是由伪随机序列充当的，不同的基站利用不同的扩频码。
在CDMA2000系统中在WCDMA系统中采纳的是码长为218-1的Gold码。
第三层是移动用户码，在CDMA2000系统中，利用的是码长为242-1的m序列，在WCDMA系统中采纳的是码长为225-1的Gold码。
一个用户一个，各不相同，它是由相当长的伪随机序列加上移动用户自身代码复合而成的。
第二、三层的码统称扰码。
在这三层扩频码中，除第一层的信道编码外另两层扩频码都由伪随机序列来实现的。
通信的码分系统和超短波战术通信的码分系统。
民用通信方面，也接踵显现一些具体的方案。
第二章伪随机序列
2.1伪随机序列相关概念
伪随机序列作为扩频通信系统中的一部份是十分关键的，它关系到扩频系统的性能。
四十年代末，信息论的奠基人香农（C.E.Shannon）提出的编码定理指出：
只要信息速度Rb小于信道容量C，那么总能够找到某种编码方式，在码周期相当长的条件下，能够几乎无过失的从收到高斯噪声干扰的信号中复制出原发信息。
那个地址有两个条件，一是Rb<=C，二是编码的码周期足够长。
同时香农在证明编码定理的时候，提出用具有白噪声统计特性的信号来编码。
白噪声是一种随机进程，它的瞬时值服从正态散布，功率谱在很宽频带内都是均匀的。
可是至今无法实现对白噪声放大、调制、检测、同步及操纵等，而只能用具有类似于限带白噪声统计特性的伪随机序列信号来逼近它，并作为扩频系统的扩频码。
六十年代末，一些易于产生、加工和复制且具有白噪声性质的“伪噪声编码技术”日趋成熟，因此高效抗干扰编码通信变得蓬勃进展起来。
同时用各类不同波形的正交码来实现波形分割的码分多址通信也接踵显现，实现了无线用户的随意呼唤通信。
这种技术在地面多址
2.1.1伪随机序列的数学概念
伪随机序列（伪随机码）的一样概念是：
若是一个序列，一方面它的结构（或形式）是能够预先确信的，而且是能够重复地产生和复制的；另一方面它又有某种随机序列的随机特性（即统计特性），咱们称这种序列为伪随机序列（伪随机码）。
伪随机序列尽管只有两个电平，但却具有类似白噪声的相关特性，只是幅度概率散布再也不服从高斯散布。
它应具有如下特性：
（l）每一周期内0和1显现的次数近似相等。
（2）每一周期内，长度为n比特的游程显现的次数比长度为n+1比特游程次数多一倍（游程是指相同码元的码元串）。
（3）关于狭义伪随机序列，将给定随机序列位移任何一个非零数量个元素，所得的序列将和原序列有一半的元素相同，一半的元素不同。
白噪声是一种随机进程，瞬时值服从正态散布，自相关函数和功率谱密度有极好的相关性，伪随机序列是针对白噪声演化而来的，只有“0”和“1”两种电平，因此伪随机编码概率散布不具有正态散布形式。
但当序列足够长时，由中心极限定理可知，它趋近于正态散布，由此，伪随机序列概念如下：
（1）凡自相关函数具有
（2-1）
式的序列称为狭义伪随机序列。
（2）凡自相关函数具有
（2-2）
形式的序列，成为第一类广义伪随机序列。
（3）凡相互关系数具有
或
（2-3）
形式的序列，称为第二类广义伪随机序列。
（4）凡相关函数知足
（1）、
（2）、（3）三者之一的序列，统称为伪随机序列。
由上面的四种概念能够看出，狭义伪随机序列是第一类广义伪随机序列的一种特例。
2.1.2随机序列的相关特性
扩频系统中，对伪随机序列而言，最关切的问题确实是其相关特性，包括自相关性、相互关性及部份相关性。
下面别离给出这些相关函数的概念。
设有两条长为N的序列{a}和{b}，序列中的元素别离为ai，bi，（i=1,2,3,…,N）。
那么序列的自相关函数概念为：
（2-4）
由于{a}是周期为P的序列，故有ai+p=ai,其归一化自相关函数ρa（τ）概念为：
（2-5）
序列{a}和{b}的相互关函数概念为：
（2-6）
归一化相互关函数概念为：
（2-7）
关于二进制序列，能够表示为：
（2-8）
其中，A为序列{a}和{b}对应码元相同的数量，D为不相同的数量。
若ρab（τ）=0，那么序列{a}和序列{b}正交，概念{a}的部份相关函数和归一化部份相关函数为（式中t为某一整数）：
（2-9）
概念序列{a}和序列{b}的部份相互关函数和归一化部份相互关函数别离为：
（2-10）
第三章m序列的产生
3.1m序列产生原理
伪随机信号具有类似于随机噪声的一些统计特性,同时又便于重复产生和处置。
目前普遍利用的伪随机信号都是由数字电路产生的周期序列取得的。
产生伪随机序列的电路包括线性反馈的移位寄放器：
m序列和非线性反馈移存器。
m序列是最长线性移位寄放器序列的简称。
它是由多级移位寄放器或其他延迟元件通过线性反馈产生的最长的码序列。
由于m序列容易产生、规律性强、有许多优良的性能，在扩频通信中最先取得普遍的应用。
如图2.1所示，m序列可由二进制线性反馈移位寄放器产生。
它要紧由n个串联的寄放器、移位脉冲产生器和模2加法器组成。
图中第i级移存器的状态ai表示，ai=0或ai=1，i=整数。
反馈线的连接状态用ci表示，ci=1表示此线接通（参加反馈），ci=0表示此线断开。
由于反馈的存在，移存器的输入端受控地输入信号。
不难看出，假设初始状态为全“0”,那么移位后取得的仍为全“0”，因此应幸免显现全“0”状态，又因为n级移存器共有2n-1种可能的不同状态，除全“0”状态外，剩下2n-1种状态可用。
每移位一次，就显现一种状态，在移位假设干次后，必然能重复显现前某一状态，其后的进程便周而复始了。
反馈线位置不同将显现不同周期的不同序列，咱们希望找到线性反馈的位置，能使移存器产生的序列最长，即达到周期P=2n-1。
按图中线路连接关系，能够写为：
（模2）（3-1）
该式称为递推方程。
图3-2线性反馈移位寄放器
上面曾经指出，ci的取值决定了移位寄放器的反馈连接和序列的结构。
此刻将它用以下方程表示：
（3-3）
这一方程称为特点多项式。
式中xi仅指明其系数ci的值（1或0），x本身的取值并无实际意义，也不需要去计算x的值。
例如，假设特点方程为f（x）=1+x+x4那么它仅表示x0,x1和x4的系数c0=c1=c4=1,其余为零。
经严格证明：
假设反馈移位寄放器的特点多项式为本原多项式，那么移位寄放器能产生m序列。
只要找到本原多项式，就可组成m系列发生器。
特点多项式与输出序列的周期有紧密关系.当F（x）知足以下三个条件时，就必然能产生m序列：
（1）F（x）是不可约的，即不能再分解多项式；
（2）F（x）可整除
那个地址
;
（3）F（x）不能整除
，那个地址q。
寻觅本原多项式是一件繁琐的工作，计算的到的结果已列表。
n
本原多项式的八进制系数表达式
代数式
2
7
3
13
4
23
5
45
6
103
7
211
8
435
9
1021
10
2011
11
4005
12
10123
表3-4本原多项式系数
表3-2给出其中部份结果，每一个n只给出一个本原多项式为了使序列发生器尽可能简单，经常使用的只有3项的本原多项式表中列出的本原多项式都是项数最少的，为了简便起见，用八进制数字记载本原多项式的系数。
由系数写出本原多项式非
常方便。
本文探讨n=5时，本多项式系数的八进制表示为45，将45写为二进制码100101，从右向左第一个1对应于C0，按系数可写出F（x）=X5+x2+1从左向右的第一个1对应于C0，按系数可写出对应的寄放器函数[
]=[00101]。
m序列的大体性质如下：
（1）周期性：
m序列的周期p取决于它的移位寄放器的级数,p=2n-1
（2）平稳特性：
m序列中0和1的个数接近相等；m序列中一个周期内“1”的数量比“0”的数量多1个。
（3）游程特性：
m序列中长度为1的游程约占游程总数的1/2，长度为2的游程约占游程总数的1/22,长度为3的游程约占游程总数的1/23…
（4）线性叠加性：
m序列和其移位后的序列逐位模2相加，所得的序列仍是m序列，只是相移不同罢了。
例如1110100与向右移3位后的序列1001110逐位模2相加后的序列为0111010，相当于原序列向右移1位后的序列,仍是m序列。
用公式表示为:
U（i）（模2加）Up（i）=Up（i）（3-5）
其中:
u（i）、up（i）、uq（i）别离为原序列、平移p个元素后的序列及平移相加后取得的序列中的第i个元素。
（5）二值自相关特性：
码位数越长越接近于随机噪声的自相关特性。
m序列的自相关函数计算式为
（3-6）
其中:
为码序列的最大长度,亦即m序列的周期；
Tc为m序列码的码元宽度。
可见,相关函数是个周期函数。
（6）m序列发生器中，并非是任何抽头组合都能产生m序列。
理论分析指出，产生的m序列数由下式决定：
（3-4）（3-7）
其中Φ（x）为欧拉数（即包括1在内的小于x并与它互质的正整数的个数）。
例如5级移位寄放器产生的31位m序列只有6个。
3.3序列产生流程图
图3-7产生m序列流程图
第四章Gold序列
m序列尽管性能优良，但一样长度的m序列个数不多，且序列之间的相互关值并非都好。
R·Gold提出了一种基于m序列的码序列，称为Gold码序列。
随着级数n的增加，Gold码序列的数量远超过同级数的m序列的数量，且Gold码序列具有良好的自相关特性和相互关特性，取得了普遍的应用。
4.1Gold序列的产生原理
Gold序列确实是为了解决m序列个数不多且m序列之间的相互关函数值不睬想而提出的，它是用一对周期和速度均相同的m序列优选对模2加后取得的。
其发生器结构框图如图3.1所示：
图4-1Gold序列发生器
Gold序列具有良好的自、相互关特性，且地址数远远大于m序列地址数。
如有两个m序列,它们的相互关函数的绝对值有界，且知足以下条件：
（4-2）
咱们称这一对m序列为优选对。
每改变两个m序列相对位移就可取得一个新的Gold序列，当相对位移2n-1位时，就可取得一族2n-1个Gold序列。
再加上两个m序列，共有2n+1个Gold序列码。
4.2Gold序列的大体性质
（1）平稳性:
Gold码序列分为平稳码和非平稳码。
Gold序列的平稳特性有3种,也确实是Gold序列有3种“0”和“1”情形:
①“1”码元数量仅比“0”码元数量多一个,这确实是平稳Gold序列。
②“1”码元过量。
③“1”码元过少。
后两种序列是不平衡Gold序列。
当n为奇数时,在周期N的N+2个Gold序列中，有2n-1个序列是平稳的。
即平稳码数量占50%,非平稳码数量占50%。
当n为偶数，但不能被4整除时,在周期N=2n-1的N+2个Gold序列中,平稳码占75%,非平稳码占25%。
相比较而言,m序列是平稳的,“1”码和“0”码的个数大体相等。
（2）自相关特性：
Gold证明了Gold码序列的自相关函数的所有非最顶峰的取值是三值。
其自相关函数值所有非最顶峰取值R如下式。
其中p=2n-1,p为Gold码序列的周期。
（4-3）
在位移k=0，R取得最顶峰，即R=1,现在同m序列一样，具有尖锐的自相关峰值。
因此Gold码序列应具有四个值的自相关函数值。
（3.2）式中t的取值是，当n为奇数时，
。
当n为偶数,且不是4的整倍数时，
。
（3）相互关特性：
Gold码序列具有较好的相互关特性，Gold码序列的相互关函数值的最大值不超过其m序列优选对的相互关值。
Gold码也具有三值相互关函数值，其取值同（3.2）式。
当n为奇数时，序列族中约50%的码序列的相互关函数值为-1/p；而n为偶数时，有75%的码序列的相互关函数值为-1/p。
第五章基于matlab产生gold序列
图5-1产生gold序列流程图
5.2.gold序列的产生和相关系分析
Gold序列是有两个m序列模2加取得的，因此第一步确实是要产生二个m序列，且m序列的产生要用到本来多项式，先确信移位寄放器的位数，在查表取得本来多项式，在编写程序产生m序列，检查是不是正确，，在编程进行摸2加取得gold序列，关于matlab的函数用stem不用flot，stem能专门好显示0,1。
程序调试产生gold序列及结果分析
Matlab语言是现今国际上科学界最具阻碍力，也是最有活力的软件。
它起源于矩阵运算,并已经进展成一种高度集成的运算机语言。
Matlab具有壮大的数学运算能力，方便有效的画图功能及语言的高度集成。
Matlab是矩阵实验室（MatrixLaboratory）之意。
Matlab除具有卓越的数值计算能力外，它还提供了专业水平的符号计算、文字处置、可视化建仿照真和实时操纵等功能。
取得gold序列值为：
通过运行结果能够看出是由0、1组成的阶梯形图形，stem函数使结果明显，其中r=6时，l=2^r-1=63位。
通过计算两个m序列摩尔加取得的gold序列与理论计算值一致，达到了实验要求。
5.2.2Gold序列自相关性其自相关性
自相关函数反映一个信号在不同时刻上去相关联程度，且峰值越平稳相关性越好。
对已产生的gold序列进行自相关分析，先要变成双极性，在编程进行相关性分析。
运行结果分析：
上图可知gold序列应具有四个值的自相关函数值，在位移k=0，R取得最顶峰，即R=1,现在同m序列一样，具有尖锐的自相关峰值，gold序列自相关性良好。
自相关的运行结果与理论分析相符，达到了实验的要求。
升级会员 