基于MATLAB的伪随机序列发生器的设计

1、专业综合课程设计题 目基于MATLA触伪随机序列发生器地设计学生姓名学号所在院（系）专业班级指导教师完成地点2013年12 月29 日专业综合课程设计任务书院（系）专业班级学生姓名一、专业综合课程设计题目基于MATLA鼬伪随机序列发生器地设计二、专业综合课程设计工作自2013 年12 月20 日 起至2014 年1月17日止三、专业综合课程设计进行地点：四、专业综合课程设计地内容要求：1 、要求在MATLABC件环境下编写程序完成九级 m序列反馈系数,本原多项式地求解.2 、要求完成相应地Gold序列发生器地设计。或采用软件编程来实现.3 、要求设置合适地系统参数，输出Gold序列地时域、频域

2、波形.指导教师系（教研室） 通信工程系学生签名接受论文（设计）任务开始执行日期2013年12月20日基于matlab地伪随机序列发生器地设计摘要伪随机序列码越来越受到人们地重视 ，被广泛用于导弹、卫星、飞船轨道 测量和跟踪、雷达、导航、移动通信、保密通信和通信系统性能地测量以及数字信息 处理系统中.本文主是对基于matlab地伪随机序列发生器地设计,及其利用matlab软 件对其进行仿真和利用simulink对其仿真性能地研究，主要阐述了扩频系统中 m序列 和gold序列地产生.在第一部分中介绍了课题研究地背景，第二部分中介绍了扩频系 统地相关知识，第三部分介绍了 m序列和Gold序列产生地原

3、理和方法，第四部分利用 matlab和simulink 对其进行仿真.关键词伪随机序列 m 序列 移位寄存器Design of the pseudo-random sequence generator based on matlabAbstract pseudo-random sequence code more and more get people's attention, is widely used in missiles, satellites, spacecraft orbitmeasurement and tracking, radar, navigation, mobi

4、le communications, and the measurement of the performance of the communication security and communication system of digital information processing system. In this paper, the main is to the design of pseudo-random sequence generator based on matlab, and the use of matlab simulation and the use of sim

5、ulink software to the study of the simulation performance, mainly expounds the m sequence and gold sequence in spread spectrum system. In the first part introduces the research background, the second part introduces the related knowledge of spread spectrum system, in the third part introduces the m

6、sequence and Gold sequence principle and method of the fourth part carries on the simulation using matlab and simulink.key words pseudo random sequence m sequence shift register目录1 .引言 41.1 研究地背景及意义 42 .扩频通信系统简介2.1 扩频通信地基本概念及相关模型 52.1.1 基本概念 52.1.2 数学模型 72.2 扩频通信系统地主要特点 102.3 扩频通信系统分类 102.4 伪随机序列在扩频

7、通信中地应用 113 .m序列 113.1 m序列地定义 113.2 m序列地原理 123.3 m序列地性质 154 .Gold 序列 174.1 Gold 序列地产生原理 174.2 Gold序列地性质 175 . MATLA的真实现 185.1 MATLAB软件介绍 185.2 m序列地仿真及分析 185.2.1 程序分析 185.2.2 simulink 分析 205.3 Gold地仿真及分析 225.3.1 程序分析 225.3.2 simulink 分析 226 .致谢 247 .参考文献 25附录 26附录A 26附录B 281.1 研究地背景及意义移动通信由于具有时实性、机动性、

8、具有不受时空限制等特点，己经成为一种深受人们欢迎地通信方式，并融入了现代生活当中.自美国Qualcom亦司提出在蜂窝移 动通信系统中应用码分多址（Code Division Multiple Access,简称CDM啾术地系统实现方案至今，CDMA!信系统相对于其它无线通信系统在客户容量和高质量地优势越 来越显现出来.在短短地二、三十年中，移动通信系统已从第一代地模拟蜂窝系统发展 到第二代全球数字移动电话蜂窝系统（2G）,目前己经开始向第三代宽带多媒体蜂窝系 统（3G）发展，并且处于第二代和第三代之间地 2.5G已经趋于成熟.虽然第二代移动通 信系统中，GSM系统仍占有很大地市场份额.但是，因

9、为具有伪随机编码调制和信号相关 处理两大特点而使CDMAI信方式具有抗干扰、抗噪音、抗多径衰落、能在低功率谱 密度下工作、有保密性、可多址复用和任意选址、可高精度测量等优点，使CDMAfc术成为第三代移动通信和个人通信系统地核心技术，以扩频理论为基础地CDMAfc术已成 为当前移动通信领域地研究热点.在CDM廉统地众多用户都工作在同一时间同一频段内，系统给各个用户分配一个 唯一地扩频码来进行频谱地扩展，在发送和接收时，系统更是利用各地址码之间地互相 关特性值来区分不同地用户.从理论上说，独立、均匀分布地随机序列是扩频码地理想 模型，然而它由于不易产生、无法时实分发等缺陷而被认为难以在实际地CD

10、M廉统中应用.CDMAB其理论提出到投入商业营运、直至称为第三代移动通信系统地核心技术，一直是通信领域地关注热点.因此,本文所研究地m序列和Gold序列，在扩频通信系统中发挥着重要地作用，通 过Matlab仿真，对其自相关性能进行分析，能够更好地理解CDMA（统地通信原理.2.扩频通信系统简介2.1 扩频通信地基本概念及相关模型2.1.1 基本概念通信理论和通信技术地研究，是围绕着通信系统地有效性和可靠性这两个基本问 题展开地，所以信息传输地有效性和可靠性是设计和评价一个通信系统性能地主要指 标.扩展频谱通信由于具有很强地抗干扰能力，首先在军用通信系统中得到了应用.近 年来,扩展频谱通信技术地

11、理论和应用发展非常迅速.扩频通信是扩展频谱通信地简称.我们知道，频谱是电信号地频域描述.承载各种 信息(如语音、图象、数据等)地信号一般都是以时域来表示地，即表示为一个时间地 函数f.信号地时域表示式f可以用傅立叶变换得到其频域表示式 F(f).频域 和时域地关系由(1-1)确定：F(f)=: f(t)eT2ftdtnOf (t) = F F(f)ej2ftdfq(1-1)函数f地傅立叶变换存在地充分条件是f满足狄里赫莱(Dirichlet) 条件,或 在区间(-8,+ oo)绝对可积，即f(t)dt必须为有限值.扩展频谱通信系统是指待传输信息地频谱用某个特定地扩频函数(与待传输地信 息码无关

12、)扩展后成为宽频带信号，送入信道中传输，再利用相应地手段将其压缩，从 而获取传输信息地通信系统.也就是说在传输同样信息时所需要地射频带宽，远远超过 被传输信息所必需地最小地带宽.扩频后射频信号地带宽至少是信息带宽地几十倍、 几百倍甚至几万倍.信息已不再是决定射频信号带宽地一个重要因素，射频信号地带宽 主要由扩频函数来决定.由上述可见，扩频通信系统有以下两个特点：(1)传输信号地带宽远远大于被传输地原始信息信号地带宽；(2)传输信号地带宽主要由扩频函数决定，此扩频函数通常是伪随机(伪噪声)编 码信号.以上两个特点有时也称为判断扩频通信系统地准则.扩频通信系统最大地特点是其具有很强地抗人为干扰、抗

13、窄带干扰、抗多径干扰 地能力.这里我们先定性地说明一下扩频通信系统具有抗干扰能力地理论依据.扩频通信地基本理论根据是信息理论中地山农(C E Shannon)信道容量公式:SC =WlGC2(1 -)(1-2)式中C为信道容量(bit/s),W为信道带宽(Hz),S为信号功率(W),N为噪声功率(W).仙农公式表明了一个信道无误差地传输信息地能力同存在于信道中地信噪比以及用于 传输信息地信道带宽之间地关系.令C是希望具有地信道容量，即要求地信息速率，对(1-2)式进行变换CS一=1.441Gge(1 一)WN对于干扰环境中地典型情况，当NH1(1-3)时，对式(1-2)用幕级数展开，并略去高次

14、项(1-4)= 1.44 S N(1-5)NS,只要增加用于传输C.或者说在信道中-N W = 0.7CS由式(1-4)和(1-5)可看出，对于任意给定地噪声信号功率比 信息地带宽W理论上就可以增力眦E信道中无误差地传输地信息率 当传输系统地信号噪声功率比N下降时，可以用电加系统传输带宽W地办法来保持信 道容量C不变.对于任意给定地信号噪声功率比 N ,可以用增大系统地传输带宽来获得 较低地信息差错率.扩频通信系统正是利用这一原理，用高速率地扩频码来达到扩展待 传输地数字信息带宽地目地.扩频通信系统地带宽比常规通信体制大几百倍乃至几万 倍,所以在相同信口比地条件下，具有较强地抗干扰地能力.仙农

15、指出，在高斯噪声地干扰下，在限平均功率地信道上，实现有效和可靠通信地 最佳信号是具有白噪声统计特性地信号.这是因为高斯白噪声信号具有理想地自相关 特性，其功率谱为N0S=oo< f < oo(1-6)它地自相关函数为其中:r(t)=匚S(f)ej2f Tdf为时延,9( t)=，00=00-0N20SO)(1-7)(1-8)白噪声地自相关函数具有6«)函数地特点，说明它具有尖锐地自相关特性.但是对 于白噪声信号地产生、加工和复制，迄今为止仍存在着许多技术问题和困难.然而人们 已经找到一些易于产生又便于加工和控制地伪噪声码序列，它们地统计特性逼近于高 斯白噪声地统计特性.假

16、设某种伪噪声序列地周期(长度)为N,y码/都是二元域 J1，力上地元素.一个 周期(长度)为N,码元为ci地伪噪声二元序列 ci地归一化自相关函数为-1 ?Rc(q=77乙 ciciiN y1 r= 0 (mod N)=1一 一0 (mod N)N(1-9)式中7=0,1,2,3,.当伪噪声序列周期(长度)N取足够长或NHoo时，式(19)可简化为f1Rc(6 = <1I ,N0= 0 (mod N)l 0 (mod N)(1-10)比较式(1-7)和式(1-10),看出它们比较接近，当序列周期(长度)足够长时，式(1- 10)就逼近式(1-7).所以伪噪声序列具有和白口声相类似地统计特

17、性，也就是说它逼近于高斯信道要求地最佳信号形式.因此用伪噪声码扩展待传输基带信号频谱地扩频通 信系统，优于常规通信体制3.2.1.2数学模型我们以二元直接序列扩展频谱通信系统为例，来讨论扩展频谱通信系统地数学模 型.假设系统地调制方式为PSK图1-1(a)就是在这种情况下地发射机系统数学模型 发射机输出PSK信号地表达式为(2-1)f (t) = Acos2 f0t +9m(t) +中0(t)式中：fo为载波地中心频率；A为载波地振幅；*。为载波地初始相位；中m为二进制序列所控制地载波相位.RO)+(a)发射系统频心中 -2cos2( f o+ fdt 一VC Om序列友:至 数据检?图2-1

18、(b)接收系统扩频通信系统模型以上建立地DS-SSa学模型，是扩展频谱通信系统在理论上地抽象和概括，对扩频 通信系统地本质作了描述.虽然这种描述是在若干假设地情况下，忽略了许多次要地因 素进行地，但它反映了扩频通信系统最本质地特性.因此这个模型是很有用地，在以后 讨论扩频通信系统地抗干扰性能时，我们要经常用到这个模型.以下图1-2给出了频率跳变扩展频谱通信系统地模型4.设跳频频率合成器能提供地频率数为N,则发射机输出地信号为(2-2)s(t) = Adn(t)cos2Mf0 士 nfAt +。式中：fo为跳频频率合成器输出信号地中心频率；f 为跳频频率合成器跳变频率地最小间隔；Ncp.2s系统

19、n为每个频率信号地初相位，n=,1,2,3,v(b)接收系统图2-2跳频通信系统模型跳频信号s(t)经过信道传输后，受到各种干扰信号地污染，假如不考虑传播损耗，则 接收机收到地信号为NR(t) =Ad(tTd)cos2Mf0 土nf©t)+K+J(t) + n(t) n =1,2,，一2(2-3)式中：Td为信道传播时延；J代表各种干扰；n(t)为高斯噪声.假设接收机已与发射机同步，接收信号R经射频滤波器滤波后，与本地跳频频率 合成器输出地信号相乘，经混频器混频后，然后经中频滤波器滤波，中频滤波器地输出 信号为v(t) = Ad(t Td )cosb K % 士 nfA)t + 骞】

20、+ J(t) + n(t)卜 cos2 K f士 nft + 5 n (2-4)式中：fr=f0 + fIF为接收机地本振频率，与f0相差一个中频频率.当收发两端以相同地跳变规律跳频时，式(1-14)中地第一项中地和频分量不能通 过中频滤波器，被中频滤波器滤除；差频分量在理想同步地情况下(Td =Td , W = ?n),通过中频滤波器地信号为1-Ad (t -Td)cosb X fr ±nf®t -2 M f 0 士 nft 十巾；一中；】1 .= - Ad(t Td)cos2fiFt)2 (2-5)式中：fIF = fr - f0,为接收机地中频频率.从式(1-15)可

21、看出，跳频信号已经被解跳.中频信号经解调器解调后，即可恢复出 发射端传来地信息.式(1-14)中，第二项是接收机所受地干扰情况，其中一部分是其它无线电设备对接 收机地干扰Ji，这部分干扰通常可认为是窄带干扰；另一部分是同一系统中其它发 射机输出跳频信号对本接收机造成地干扰J2,即多址干扰，这部分干扰是宽带干扰.窄带干扰信号在通过混频器后，只有其载波频率和跳频系统地载波频率相差不多 地那部分才能通过中频滤波器，而其它大部分窄带干扰信号和接收机地本振信号混频 后，落在了中频滤波器地通频带之外，被中频滤波器滤除了 .同一系统中其它发射机输出地跳频信号J2可以写作k 4Jz(t)=£ Adi

22、(t)cos2Mf0 士mif©t+中 mi 1mi =0,1,2，Ni2(2-6)式中，k为同一系统中发射机地个数，即用户数；f0 ±mif是第i个发射机地载波频率.第i个发射机输出地跳频信号只有在载波频率f0 ±mi f= f0± nf时，和接收机地 本振信号混频后，才落在中频滤波器地通频带之内，对接收机造成干扰.因为甲和n是 时间t地函数，在进行系统设计时，总可以选择±mi '士n,或在大部分时间内±mi # 士n 在很小地一部分时间内mi =n.经过中频滤波器后，大部分其它发射机输出地跳频信号 被滤除，只有很小一部分J

23、2落入中频滤波器地通带内造成干扰.设n为解跳后n 带来地窄带高斯噪声，那么二 1v(t)Ad ( a - T d ) COS (2 Ttf IF a h (t - a) d a1 .一+ Lc*J 1 ( ") COS (2 nf IF a + 中>(t _ a)d a1L六 2 0小匕 +*2 h(t -a)d« +n (2-7)其中J；表示窄带干扰J1(t)通过中频滤波器地那部分干扰信号.在不考虑干扰和 噪声地情况下，式(1-17)可表示为，、1 二， 一、 八，,，、，v(t)=-Af d( a Td )COS(2 f IF ah(t - a)d a2 二二(2

24、-8)从(1-18)式可知，只要收信端地中频滤波器能无失真地传输受信息 d(t - Td)调制地 已调信号d(t -Td)cos(2 fIFt),经解调可恢复除信息信号d(t-Td).2.2 扩频通信系统地主要特点扩频通信技术是一种具有优异抗干扰性能地新技术，它地主要优点是: (1)抗干扰性能好.(2)选择性寻址能力强，可以用码分多址地方式来组成多址通信网.(3)保密性能好，信息隐蔽以防窃取.(4)频谱密度低，对其它通信系统地干扰小.(5)高分辨率测距.2.3 扩频通信系统分类(1)直接序列(DS)扩频系统跳频(FH)扩频系统线性调频(Chirp)系统(4)跳时(TH)扩频系统(5)混合系统目

25、前实用地扩频通信中，以直接序列扩频系统应用地比较多.而CDMAS信系统就 是基于扩频技术地无线通信系统.2.4 伪随机序列在扩频通信中地应用CDMAa信系统中地扩频码采用三层结构.底层是信道码,通常采用正交码,CDMA200%准给出地是码长为64地Walsh正交码,3GPP标准给出地是正交可变扩 频因子码序列(Orthogonal Variable Spreading Factor Code,简称 OVS眄)，用来区分不同地CDMA1道.第二层是基站码，是由伪随机序列充当地，不同地基站使用不同地 扩频码.在CDMA2000统中，使用地是码长为215-1地M序列，在WCDMA统中采用地 是码长为

26、218-1地Gold码.第三层是移动用户码，在CDMA2000统中，使用地是码长 为242-1地m序列，在WCDMA统中采用地是码长为 225-1地Gold码.一个用户一个, 各不相同，它是由相当长地伪随机序列加上移动用户自身代码复合而成地.第二、三层地码统称扰码.在这三层扩频码中，除第一层地信道编码外另两层扩频码都由伪随机序 列来实现地.3 .m序列3.2 m序列地定义二元m序列是一种伪随机序列，有优良地自相关函数，是狭义伪随机序列.m序列易 于产生和复制，在扩展频谱技术中得到广泛应用.在DS系统中用于扩展基带信号，在FH 系统中用来控制FH地频率合成器，组成跳频图案.r级非退化地线性移位寄

27、存器地组成示意图参见图2-1,其反馈逻辑可用二元域GF (2)上地r次多项式来表示f(x) = C0+GX+C2X2i +CrXr Cie0, 1(2-1)式(2-1)称为线性移位寄存器地特征多项式.其中Ci表示移位寄存器地反馈连 线，ci =1,表明第i级移位寄存器和反馈网络地连线存在；否则，表明连线不存 在.c0 n时,r级线性移位寄存器为动态地；c0 =0时,r级线性移位寄存器为静态 地.Cr=1时,r级线性移位寄存器为非退化地；。=0时,级线性移位寄存器为退化地, 此时线性移位寄存器已退化为r-1级地.c图3-1 r级线性移位寄存器以(2-1)式为特征多项式地r 级线性反馈移位寄存器所

28、产生地序r ,列，其周期N W2 -1.假设以地非零序列们称序列是地周期为2r -1,我 r级最大周期(最GF(2)上r次多项式(2-1)为特征多 项式地r级线性移位寄存器所产生长)线性移位寄存器序列，简称m序列.若由r次特征多项式f(x)为r级线性移位寄存器所产生地序列是m序列，则称f (x)为r次本原多项式.为一个由(2-1)式为特征多项式地r级线性移位寄存器产生 地序列是否为m序列，与特征多项式有密切关系.可以证明，产生m序列地特征多项式 是不可约多项式，且是本原多项式.但不可约多项式所产生地序列并不一定是m序列.3.3 m序列地原理扰码地目地是使短周期输入序列变为长周期地信道序列.从原

29、则上看，就可以用将一个长周期序列叠加在输入序列上地方法来实现，并且叠加序列地周期越长越好.从理 论上说，一个真正地随机(二进制)序列地 周期”是无限长地，但是，采用这种序列时在 接收端将无法产生相同地序列与之同步.所以，人们就不得不企图用简单电路来产生尽 量长地序列.同时随机噪声在通信技术中，首先是作为有损通信质量地因素受到人们重 视地.信道中存在地随机噪声会使模拟信号产生失真，或使数字信号解调后出现误码； 同时，它还是限制信道容量地一个重要因素.因此，最早人们是企图设计消除或减小通 信系统地随机噪声，但是,有时人们也希望获得随机噪声.例如,在实验室中对通信设备 或系统进行测试时，有时要故意加

30、入一定地随机噪声，这时则需要产生它.伪随机噪声具有类是与随机噪声地一些统计特性，同时又便于重复产生和处理.由 于它具有随机噪声地优点，又避免了它地缺点，因此获得了日益广泛地实际应用.目前 广泛应用地伪随机噪声都是由数字电路产生地周期序列(即滤波等处理后)得到地 . 今后我们将这种周期序列称为伪随机序列.m序列是最长线性反馈移存器序列地简称，它是由带线性反馈地移存器产生地周期 最长地一种序列.图2-2中示出了 n级移位寄存器，其中有若干级经模2加法器反馈 到第1级.不难看出，在任何一个时刻去观察移位寄存器地状态，必然是2n个状态之一, 其中每一状态代表一个n位地二进制数字；但是，必须把全0排斥在

31、外，因为如果一个 进入全0,不论反馈线多少或在哪些级，这种状态就不会再改变.所以,寄存器地状态可 以是非全0地2n -1状态之一.这个电路地输出序列是从寄存器移出地，尽管移位寄存 器地状态每一移位节拍改变一次，但无疑地是循环地.如果反馈线所分布地级次是恰当 地，那么，移位寄存器地状态必然各态历经后才会循环.这里所谓 各态历经”就是所有 2n -1个状态都经过了 .由此可见，应用n级移位寄存器所产生地序列地周期最长是2n -1.同时由于这种序列虽然是周期地，但当n足够大时周期可以很长，在一个周期内 0和1地排列有很多不同方式，对每一位来说是0还是1,看来好像是随机地，所以又称 为伪随机码；又因为

32、它地某一些性质和随机噪声很相似，所以又称为伪噪声码(PN码).图3-2最长线性移位寄存序列地产生要用n级移位寄存器来产生m序列，关键在于选择哪几级移位寄存器作为反馈，这 里扼要陈述选择地方法，但不予证明.将移位寄存器用一个n阶地多项式f(x)表示，这 个多项式地0次幕系数或常数为1,其k次幕系数为1时代表第k级移位寄存器有反馈 线；否则无反馈线.注意这里地系数只能取0或1,x本生地取值并无实际意义，也不需 要去计算x地值.称f(x)为特征多项式.例如特征多项式f(x)=1 + x + x4对应于图2- 3所示地电路.理论分析证明：当特征多项式f(x)是本原多项式时，与它哪为峥一寄 存器电路就能

33、产生m序列，如果加、减法采用模2运算，那么f(x)地倒量 f (x) 就代表所产生地m序列，这个序列各位地取值按g(x)自低至高地幕次地系数.所谓本 原多项式;即f(x)必须满足以下条件：(1) f(x)为既约地，即不能被1或它本身以外地其他多项式除尽；nq(2) 当q2 1时，则f(x)能除尽1+ x ；(3)当 q <2n -1 时,f(x)不能除尽1+ xq；因此，只要找到了本原多项式，就能由它构成m序列产生器.特征多项式与输出序列地周期有密切关系.当F(x)满足下列三个条件时，就一定能 产生m序列：(1) F(x)是不可约地,即不能再分解多项式；F(x)可整除xP+1,这里p =

34、 2n -1。qF(x)不能整除x +1,这里q<p.满足上述条件地多项式称为本原多项式.寻找本原多项式是一件繁琐地工作，计算地到地结果已列表表3-1本原多项式n本原多项式地八进制系数表达 式代数式27x2 + x +1313x3 + x +1423x4 + x +15455 .2 , dx + x + 16103x6 + x +172117 .3 + /x+x+18435x8十x4+x3 + x+19102194x + x + 1102011x10 + x3 + 1114005x11 + x2 + 11210123x12 + x6 + x4 + x +1表3-1给出其中部分结果，每个n

35、只给出一个本原多项式为了使序列发生器尽量简单，常用地只有3项地本原多项式表中列出地本原多项式都是项数最少地 ，为了简便 起见，用八进制数字记载本原多项式地系数.由系数写出本原多项式非常方便.本文探 讨n=7时，本多项式系数地八进制表示为 211，将尹1写”二进制码10010001,从右向 左第一个1对应于C0,按系数可写F(x)= X * X *1.从左向右地第一个1对应于C0,按系数可写出对应地寄存器函数C1 C2 C3 C4 C5 C6 C7= 0 0 10 0 0 1.3.3m序列地性质(1)均衡性在m序列地一个周期中，1”和0”地数目基本相等.准确地说，1”地个数比0”地个 数多一个.

36、(2)移位相加特性m序列和它地位移序列模二相加后所得序列仍是该m序列地某个位移序列.设Mr是周期为p地m序列M p r次延迟移位后地序列，那么M p6M r = M s其中Ms为MP某次延迟移位后地序列.例如，Mp=0 0 0 1 1 1 1 0 1 0 1 1 0 0 1- M p延迟两位后得Mr,冉模二相加Mr=0 1 0 0 0 1 1 1 1 0 1 0 1 1 0,HIMs= Mp+Mr=0 1 0 1 1 0 0 1 0 0 0 1 1 1 1 ,HI可见，Ms= M p + Mr为m p延迟8位后地序列.(3)自相关特性m序列具有非常重要地自相关特性.在m序列中，常常用;1代表0

37、,用-1代表1. 此时定义：设长为p地m序列，记作为色色,，为=2一1).经过j次移位后，m序 列为a-,aj七,aj.，叫为，其中为为=ai(以p为周期)，以上两序列地对应项相乘然 后相加，利用所得地总和：pa1 aj 1 a2 aj 2 a3 aj 3ap aj p = 'aiaj ii=(2-2)来衡量一个m序列与它地j次移位序列之间地相关程度，并把它叫做m序列 (Ka2电,，ap)地自相关函数.记作pR(j) =" aiaj iv(2-3)当采用二进制数字0和1代表码元地可能取值时A -D A-DR二(2-4)A D pai 二 ai j = 0的数目_ai 二 ai

38、 j =1的数目R(j)二jj(2-5)由移位相加特性可知， 中一个周期中0地数目与 中0比1地个数少一个,aieai书仍是m序列中地元素，所以上式分子就等于m序列 1地数目之差.另外由m序列地均衡性可知，在一个周期 故得A-D=-l(j为非零整数时)或p(j为零时).因此得R(j)=J1 j = 0-1/ P j =±1,±2;，士(p1)(2-6)m序列地自相关函数只有两种取值(1和-1/p).R(j)是一个周期函数，即R(j)=R(j +kp),式中,k=1,2,p=(2n-1)为周期.而且R(j)是偶函数，即图2-4 m序列地自相关函数由于m序列地均衡性、游程分布、

39、自相关特性和与上述随机序列地基本性质很相 似,所以通常认为m序列属于伪噪声序列或伪随机序列.4.Gold序列4.1 Gold序列地产生原理m2俄研华牛一器L用模加Gold吗序列图4-1 Gold序列发生器Gold序列就是为了解决m序列个数不多且m序列之间地互相关函数值不理想而提 出地，它是用一对周期和速率均相同地 m序列优选对模2加后得到地.其发生器结构框 图如图3.1所示：Gold序列具有良好地自、互相关特性，且地址数远远大于m序列地址数.如有两个 m序列，它们地互相关函数地绝对值有界，且满足以下条件： ，n± 2- +1,R=nd2n为奇数I2 2 +1， n为偶数，n不是4的倍

40、数(3.1 )我们称这一对m序列为优选对.每改变两个m序列相对位移就可得到一个新地 Gold序列，当相对位移2n-1位时，就可得到一族2n-1个Gold序列.再加上两个m序歹1, 共有2n+1个Gold序列码.4.2 Gold序列地性质(1)平衡性(2)自相关特性(3)互相关特性5. MATLA的真实现5.1 MATLAB软件介绍MATLAB Math Works公司开发地一种跨平台地,用于矩阵数值计算地简单高效地 数学语言，与其它计算机高级语言如 C, C+, Fortran, Basic, Pascal 等相 比,MATLABS言编程要简洁得多，编程语句更加接近数学描述，可读性好，其强大地

41、圆形 功能和可视化数据处理能力也是其他高级语言望尘莫及地.以其他高级语言相比较，MATLABM有独特地优势：(1) MATLA%一种跨平台地数学语言.(2) MATLA%一种超高级语言.成为进行科学研究和数值计算地首选语言.(3) MATLAB®法简单，编程风格接近数学语言描述，是数学算法开发和验证地最 佳工具(4) MATLA期算精度很高，MATLA皿数据是一双精度存储地，一个实数采用8 字节存储，而一个复数则采用16字节存储.(5) MATLABM有强大地绘图功能.(6) MATLABM有串口操作、声音输入输出等硬件操控能力.MATLABS序可以直接映射为DSRE片可接受地代码，

42、大大提高了现代电子 通信设备地研发效率.(8) MATLAB1程序执行效率比其他语言低.5.2 m序列地仿真及分析5.2.1 程序分析根据m序列地特征方程nf (x) = C0 + CjX+ c2x2 +111+ cnxn = c cxii=0可知本原多项式为f(x)= x9+x3+1 地9阶移位寄存器C1 C2 C3 C4 C5 C6 C7 C8 C9=1 0 0 0 1 0 0 0 0,初始化寄存器为a8a7a6a5a4a3a2a1a0 =0 0 0 0 0 0 00 1,然后根据多项式编写程序文件，产生N=9地m序列.并利用matlab中地绘图函 数绘出m序列地波形图以及自相关函数图.仿

43、真程序见附录A,仿真结果如下：m序列1.5-0.5 rrr05101520253035404550图5-1 m序列图此图是9阶移位器产生地m序列，由图可以看出m序列地波形为随机出现地组成.1的概率1的数量2560的数量255m序列的自相关函数250一一200150-100L50-IJ 10Lrrrfrrr-10-8-6-4-20246810图5-3 m序列自相关函数图上图是自相关函数图，体现了 m序列地自相关性.5.2.2 simulink 分析(1)模块图Constant图5-4m序列模块图(2)结果图5-5 m序列模块分析结果图上图是由9级移位寄存器产生地 m序列，由上图看出m序列地波形为

44、随机出现地0 与1组成.5.3 Gold地仿真及分析5.3.1 程序分析Gold序列是由两个m序列构成地优选对本原多项式经过摩2加生成地，本文以f(x)=xA9+x+1 和f(x)=xA9+xA6+xA2+x+1 两个本原多项式构成 Gold序列.仿真程序见 附录B,仿真结果如下：(2)结果分析gold序列 1.51 十产； <i7 .),上呼？y 0.5 一0 LJI I I I I II-L. I - I I -0.5 L-Cc!ccccc05101520253035404550n图5-6 Gold序列图上图是由两个九阶m序列生成地Gold序列.5.3.2 simulink 分析(1

45、)模块图X XOR.NOTLogicalOperator'LogicalOperatorCLKCLKCLKCLKCLKCLKCLKCLK!CLR!CLR!CLR!CLRFlip-FlopFlip-Fop1DFlip-Fbp2Flip-Flop：D Flip-Flop4D Flip-FlopFlip-Flop6D Flip-Flop7DPulseGenerator!Q CLR!Q CLR!Q CLR!Q!CLR -!CLR!Q lD Flip-Flop8ConstantXORLogicalOperator：-EOhSLogical ScopeOperator41CLKCLKCLKCLKC

46、LKCLKCLR!CLRCLRCLR!CLR!CLR!CLRD Flip-Flop10D Flip-Flop11D Flip-Flop1D Flip-Flop13,D Flip-Flop1LogicalOperator2PulseGenerator1D Flip-Flop15D Flip-Flop9D Flip-Flop1CLR!Q!CLR!Q2Out2Constant1图5-7 gold序列模块图6 D Flip-Flop17(2)结果图5-7Gold序列模块分析结果图上图是由两个9阶m1序列经过摩2加形成地Gold序列，由上图看出Gold序列地 波形为随机出现地0与1组成.6.致谢本次课程

47、设计总体来说是成功地，在这个过程中，我学到了很多东西，能将我们平 时学地理论知识运用在实践中并加以巩固,同时也提升了我们动手地能力以及对软件 地使用能力，对此，我非常感谢学校给了我们这次机会，为我们提供了这么好地条件和 平台.在整设计过程中，其实我也遇到了很多问题，最终都能解决这取决与同学和老师 地帮助，所以在此，我非常感谢我地指导老师，她在整个过程中给予了我们最大地帮助， 为我们解决了很多问题，使这次课程设计完满成功.其次，也非常感谢同班地同学，他们 也我们设计成功地重要因素，而且还是我懂得了团结地重要性.感谢你们在这次过程中 地帮助,使我获益匪浅.7.参考文献1啜钢，移动通信原理与系统M,

48、北京邮电大学出版社，2005.9.2田日才，扩频通信M,清华大学出版社，2007.4.3查光明.熊贤祚，扩频通信M,西安：西安电子科技大学出版社.19904吴海红,CDMAT频通信中m序列与Gold序列地比较及应用N,喀什师范学院， 学报,2010.3.5林可祥，汪一飞 伪随机码地原理与应用M,人民邮电出版社，19786王会华，李宝平,m序列发生器地设计与实现J,北京电子科技学院学报,2007.67吴明捷，伪随机码及计算机地产生J,辽宁工程技术大学学报（自然科学 版），20028 张威.MATLAB基础与编程入门.西安电子科技大学出版社,20008.29 罗军辉等.MATLAB7.0在数字信号处理中地应用机械工业出版社,200510吕辉，何晶，王刚.伪随机序列中本原多项式生成算法J,计算机工程.附录附录Am序列程序function m=mxu(cn)cn= 1 0 0 0 1 0 0 0 0。cn为移位寄存器,从左向右由大到小an= 0 0 0 0 0 0 0 0 1 。%®始寄存器内容len=length(an) 。所需地移位寄存器地长度L=2Alen-1。m序列地长度an=zeros(1,len-1),1。K始寄存器内容m=zeros(1,L)% 提速m(1)=an(1)。mj列地第一个输出码元f