单片机matlab方式产生伪随机m序列

信息科学与技术学通信原理课程设I课题名称： 伪随机序列发生器的设计学生姓名：张昕灏学院： 信息科学与技术学院专业年级：电子信息项目201级指导教师：田敏副教授完成日期： 二O一三年七月十二日目录前言1第一章设计内容及要求2.设1计内容21.设2计要求21.方3案选择2第二章序列的特性分析序列的原理.均2衡特性52.游3程分布5.线4性叠加性62.自5相关特性6第三章序列设计.设1计流程图83.特2征多项式确定83.本3原多项式确定10序列的最终产生以五阶移位寄存器举例）第四章设计成果分析及总结13.仿1真结果分析13.设2计总结14心得体会15参考文献16附录程序附录51单片机实现方法18电路图18设计说明18结果验证18代码及与对应 代码数模转换输出代码：20反馈链接状态及波形输出控制代码22使用器前言扩展频谱通信是一种不同于常规通信系统的新调制理论和技术，简称扩频通信［1。］其设计思想是将待传输的信息信号用特定的扩频码扩展频谱后成为宽带信号进行传输，接收时再采用相应的技术手段将频谱压缩，恢复原来待传信息信号的带宽，从而实现通信。扩频通信具有两个特点：传输信号的带宽远大于原始信息信号的带宽；传输信号的带宽主要有扩频码决定，此扩频码通常是伪随机码。伪随机码 简称码，可以人为产生与复制，具有类似白噪声的性质，相关函数具有尖锐的特性，功率谱占据很宽的频带，易于从其他信号或干扰中分离出来，具有优良的抗干扰特性，其特点是：具有尖锐的自相关函数；互相关函数值应足够小；有足够长的码周期，以确保抗侦破与抗干扰的要求；码的数量足够多，以实现码分多址的要求；平衡性好，以满足抗干扰的要求；项目上易于产生、加工、复制与控制［2。］扩频通信的优势主要来自于伪随机码具有白噪声的统计特性。而随着扩频速率的不断提高，扩频码的长度急剧增加，利用计算机设计并验证扩频码的各项指标能大大提高效率。通过对伪随机码中常用的序列的结构和性质进行了分析，给出了基于平台的序列生成算法及代码伪随机序列分析软件平台选择 ，是美国 公司推出的一种以矩阵为基本编程单位的高效数值计算语言。 在编程效率、可读性、可移植性与可扩充性上，远远优于其他高级编程语言，是公认的最优秀的科技应用软件。第一章设计内容及要求1.设1计内容题目：伪随机序列发生器的设计1、建立通信系统的基本数学模型要根据通信系统的基本原理，确定总的系统功能，将各部分功能模块化，并找出各部分之间的关系，画出系统框图。、按设计的系统框图组建系统。、根据系统性能指标，设置和调整各模块参数。、实现系统运行仿真，观察分析结果（分析窗口、动态探针、实时显示>。1.设2计要求拿到题目后首先根据通信系统原理框图进行建模，然后在微机上进行设计和软件仿真，如仿真结果不满足要求，则修改设计，直到满足要求为止。<1）按题目要求的功能进行设计，电路各个组成部分须有设计说明；<2）设计软件可以是systemview、MaxPlusn<QuartusII）、MatLab、LabView、Protel、Multisim等软件中的一种或几种。1.方3案选择方案一：用 单片机设计，观察波形优点：直接接触硬件，设计结果更真实可靠。缺点：成本高；不能直接得到序列；当寄存器数目发生变化则相应的程序和硬件都要发生变化，不利于扩展。方案二：采用VHDL硬件描述语言设计优点：可以模拟设计出任何需要的逻辑器件，以实现相应的功能。缺点：语言描述复杂，实际生活中有些逻辑器件是很难甚至不可能实现的。方案三：使用 编程，直接输出序列基于上述原因考虑，使用纯软件方式设计是最合适的。首先，不需要依赖硬件环境；其次，程序只需一次编写，即可处处运行；再次，无论是波形或是序列码均可直接得到；最后，在各种设计语言中， 是最精通于科学数值计算和图像生成。第二章 序列的特性分析序列的原理伪随机信号具有类似于随机噪声的一些统计特性,同时又便于重复产生和处理。目前广泛使用的伪随机信号都是由数字电路产生的周期序列得到的。产生伪随机序列的电路包括线性反馈的移位寄存器：序列和非线性反馈移存器；序列［3。］序列是最长线性反馈移位寄存器序列的简称，序列是由带线性反馈的移位寄存器产生的由级串联的移位寄存器和和反馈逻辑线路可组成动态移位寄存器，如果反馈逻辑线路只由模2和构成，则称为线性反馈移位寄存器。带线性反馈逻辑的移位寄存器设定初始状态后，在时钟触发下，每次移位后各级寄存器会发生变化。其中任何一级寄存器的输出，随着时钟节拍的推移都会产生一个序列，该序列称为移位寄存器序列。序列码发生器是一种反馈移位型结构的电路，它由位移位寄存器加异或反馈网络组成，其序列长度=n只有一个多余状态即全状态，所以称为最大线性序列码发生器。因为其结构已定型，且反馈函数和连接形式都有一定的规律，因此利用查表的方式就设计出序列码。列出部分序列码的反馈函数和移存器位数的对应关系。如果给定一个序列信号长度，则根据=求出，由查表 便可以得到相应的反馈函数。表反馈函数r^i反馈函数i ।山1=—g=I均2衡特性伪随机序列的平衡性是指序列中"”的数目只比"”的数目多。码的平衡性由码序列中的直流分量决定。平衡性好，则载波抑制度大，从而有利扩频通信的抗干扰能力以及保密和抗侦破能力。其物理意义是不平衡码会使扩频后的信号中出现一些稳定的信号，从而易于被检测而导致保密能力的降低。2.游3程分布游程是一个序列中取值（1或0＞相同连在一起的元素的统称，所谓的游程长度就是一个游程中元素的个数。而序列的分布特性为1序列的一个周期 由，游程总数为n。2当1游程长度kWn游程数目占总数ok当1WkWn连“1”和连“”游程各占一半。2.线4性叠加性序列和其移位后的序列逐位模相加所得的序列还是序列只是相移不同而已.例如1110与1向0右0移3位后的序列1001逐1位1模02相加后的序列为 相当于原序列向右移位后的序列仍是序列。2.自5相关特性周期为的序列的自相关函数为:W-0 _W-口A+D p其中：—该序列与其次移位序列一个周期中对应元素相同数目—该序列与其次移位序列一个周期中对应元素不同数目P－序列周期上式可改写为：R5_［/㊉.+/=0］的数目—［%的的数目］由移位相加特性和均衡特性，可知序列的自相关函数为：1 j=0+J氏（/）=r_1炉.P 」=12……P-1J从序列的自相关系数可以看出序列是一个狭义伪随机码。当序列的移位值为其周期的整数倍时，其自相关值取得最大值为 ,移位值取其他值时，其自相关值恒为-1。其自相关函数如图2-所1示［4。］图2-1m序列的自相关函数第三章序列设计设设1计流程图级线性移位寄存器的如图3-所2示：图 级线性移位寄存器级线性移位寄存器的如图3-所2示：图 级线性移位寄存器图中」表示反馈线的两种可能连接方式，表示连线接通，第 级输出加入反馈中；出加入反馈中；回表示连接线断开，第 级输出未参加反馈。因此，一般形式的线性反馈逻辑表达式为：将上式左面的H移至右面，并将代入上式，则上式可改写定义一个与上式相对应的多项式根据上式可以确定序列的特征X线性移位信描述线性E因此，一般形式的线性反馈逻辑表达式为：将上式左面的H移至右面，并将代入上式，则上式可改写定义一个与上式相对应的多项式根据上式可以确定序列的特征X线性移位信描述线性E其中多项式。本原三特征多;一定能产件全接状态。式来存器的特征条件时，就满足上述条件的多项式称为本原多项式这样产生序列的充要条件就变成了如何寻找本原多项式。寻找本原多项式是一件繁琐的工作，计算得到的结果已列表。表3-本1原多项式系数

本原多项式的八进制系数表达式代数式1X11X1i^i1 11X11X11X|■ 1表 给出其中部分结果，每个只给出一个本原多项式为了使序列发生器尽量简单，常用的只有3项的本原多项式表中列出的本原多项式都是项数最少的，为了简便起见，用八进制数字记载本原多项式的系数。由系数写出本原多项式非常方便。本文探讨时，本多项式系数的八进制表示为，将写为二进制码 ，从右向左第一个对应于，按系数可写出。从左向右的第一个对应于」，按系数可写出对应的寄存器函数下4HHH 。序列的最终产生以五阶移位寄存器举例）根据序列的特征方程/⑴+5*+%++…+疗=z=0可知本原多项式为的5阶移位寄存器为回HHHH，移位寄存器吉构为初始化寄存器为 HWWW回，寄0存器0首先1左］移位。可知区,这时依据特征式得知反馈。故要循环次，得到所矗|程序代码如下：。因为为序列。为移位寄存器阶寄存器，码长所＞需的移位寄存器的长度序列的长度1初始寄存器内容序列的第一个输出码元S寄存器与反馈的模和移1位后的寄存器a新的寄存器输出%存对寄＞序列绘图将以上代码命名为以上可根据阶位不同而做出不同的寄序列第四章设计成果分析及总结4.仿1真结果分析里输入:

得到图4-所得到图4-所1示：1 0Columns23through311 1 0 1 01 0 0 11 1 0 1 0通过运行结果可以看出是由0、1组成的阶梯形图形，stairs函数功能是画阶梯形，生成一系列的m序列，成功实现了要求。因为m序列全为01的集合，直接的stairs绘制使得0，1的位置不清楚，使用stem(mxu＞重新运行得到图4-2。本文设计的题目是序列发生器的设计，文中介绍的扩频伪随机码序列的生成源程序是基于线性反馈移位寄存器结构而编写的，更清楚的描述了序列的生成过程中的数字逻辑。本方法应用移位寄存器理论，从序列的本原多项式出发，其算法核心是找到序列本原多项式与线性序列移位寄存器反馈逻辑式之间的关系，然后采用语言编程，并借助开发平台实现序列。序列可以软件实现，也可以硬件实现，但是通过本次设计可以看到软件设计的许多优点。在课程设计的过程中，查询了大量的资料，通过相关资料的查阅，还掌握了通信领域的有关知识，扩大了知识面。心得体会通信原理是一门很有趣的课程，任何一个系统都是一个复杂的整体，学习通信原理是要涉及到整体的每一部分。讨论某一部分原理时又要涉及到其它部分的工作原理。这样一来，不仅不能在短时间内较深入理解通信系统的工作原理，而且也很难孤立地理解某一部分的工作原理。所以，在循序渐进的课堂教案过程中，我总是处于“学会了一些新知识，弄清了一些原来保留的问题，又出现了一些新问题”的循环中，直到课程结束时，才把保留的问题勉强搞清楚。实验中，我不仅复习和巩固了通信原理中的很多知识，而且学到了很多在书本上所没有学到过的知识。且充分了解到画程序流程图的必要性。通过程序流程图，在做设计的过程中，我们每一步要做什么，每一步要完成什么任务都有一个很清楚的思路，在程序测试的过程中也有利于查错。参考文献贾怀义，毕红军，宫剑基于序列扩频码的研究北方交通大学学报，2001[2杨]睿.论伪随机序列及其应用.沈阳项目学院学报＜自然科学版），2009[3樊]昌信，曹丽娜.通信原理（第6版＞.北京:国防工业出版社，2006廖振明 序列的互相关特性研究国防科学技术大学，张志军基于 的伪随机序列的实现新乡高等师范学校学报，附录程序为移位寄存器。所需的移位寄存器的长度l 序列的长度1初始寄存器内容。 序列的第一个输出码元寄存器与反馈的模和移1位后的寄存器。新的寄存器输出对序列绘图附录51单片机实现方法电路图附录图单片机连接电路产生序列的电路图设计说明因为使用硬件设计的序列发生器，与硬件和程序都紧密相关，如果序列长度发生变化，则相应的程序电路都会变化，所以在此仅以四级设计方式举例。附录图中的开关从左至右表示 2开关合上表示连接。下面的开关5，是输出波形开关，开关合上一次输出一个完整序列波形，如果一直合上则连续输出。在本设计中，最组要的模块是数模转换，因为使用单片机得到的仅仅是一

组数字。要应用到实际系统中，就必须转换为电压或电流值，才有实际意义。结果验证因为篇幅有限，仅以cn=[1 0 01]为例演示验证。因为单片机不能输出序列，故而将程序略加修改成标位语言，在计算机上得出的结果所示：附录图标准语言方式验证得到的序列该结果与实际计算和 所得结果完全一致。接下来验证单片机上的序列波形，波形图如下：附录图标准语言方式验证得到的序列该结果与实际计算和 所得结果完全一致。接下来验证单片机上的序列波形，波形图如下：附录图3单片机仿真输出的波形从波形图上可以看到，序列的起点处是在第一下降沿处，从波形宽度来看，大致满足1这一0变化1规1律。

看，大致满足代1码及与对应代码//an=[zeros(1,len-1),1];代1码及与对应代码b[01=a[O] Nm序列的第一个输出码元数模转换输出代码：#include＜intrins.h＞#definenops(＞。do{_nop_(＞。_nop_(＞。_nop_(＞。_nop_(＞。_nop_(＞。}while(0〉#definePCF85910x90sbitADDA_SCL=P3A7。sbitADDA_SDA=P3A6。bitack。voidStart_I2c(＞。voidStop_I2c(＞。voidSendByte(unsignedcharc＞。voidtransformData(unsignedchara＞。bitDACconversion(unsignedcharsla,unsignedcharc,unsignedcharval＞。voidStart_I2c(＞//{ADDA_SDA=1。nops(＞。ADDA_SCL=1。nops(＞。ADDA_SDA=0。nops(＞。ADDA_SCL=0。nops(>。}/***•近仃2C%脆0话%b*/voidStop_I2c(>//{ADDA_SCL=0。nops(>。ADDA_SDA=0。nops(>。ADDA_SCL=1。nops(>。ADDA_SDA=1。nops(>。}/****/voidSendByte(unsignedcharc>//{unsignedchari。for(i=0。i<8。i++>{if((c<<i>&0x80>ADDA_SDA=1。elseADDA_SDA=0。nops(>。ADDA_SCL=1。nops(>。ADDA_SCL=0。}ADDA_SDA=1。nops(>。ADDA_SCL=1。nops(>。if(ADDA_SDA==1>ack=0。elseack=1。ADDA_SCL=0。}bitDACconversion(unsignedcharsla,unsignedcharc,unsignedcharval>//Start_I2c(＞。SendBy