哈工大扩频通信实验报告

1、苗需麻y至sfc季HarbinInstituteofTechnology扩频通信实验报告课程名称：扩频通信实验题目：Gold码特性研究院系：电子与信息工程学院班级：通信一班姓名：学号：指导教师：迟永钢时间：2012年5月哈尔滨工业大学第1章绪论错误!未定义书签。引言错误!未定义书签。实验内容错误!未定义书签。第2章m序列实验错误！未定义书签m序列相关概念错误！未定义书签m码序列的定义错误！未定义书签m码序列的自相关函数错误！未定义书签m序列抽取结果分析错误！未定义书签m码序列的抽取及反多项式错误！未定义书签m码序列线性移位寄存器结构及相应序列.错误！未定义书签第3章m序列优选对实验错误！未定义

2、书签m序列优选对的查找错误！未定义书签m序列优选对的定义错误！未定义书签m序列优选对的查找结果错误！未定义书签m序列优选对的自相关及互相关函数错误！未定义书签第4章Gold序列实验错误！未定义书签。Gold序列的定义错误！未定义书签。Gold序列的生成及特性错误！未定义书签。Gold序列生成结果错误！未定义书签。Gold码自相关及互相关特性错误！未定义书签。平衡Gold码相关实验错误！未定义书签。平衡Gold码的定义错误！未定义书签。平衡Gold码的判定错误！未定义书签。第5章总结错误！未定义书签。实验小结错误!未定义书签。实验心得错误!未定义书签。附录错误!未定义书签。第1章绪论引言伪随机信

3、号既有随机信号所具有的优良的相关性，又有随机信号所不具备的规律性。因此，伪随机信号既易于从干扰信号中被识别和分离出来，又可以方便地产生和重复，其相关函数接近白噪声的相关函数。m序列是目前广泛应用的一种伪随机序列，其在通信领域有着广泛的应用，如扩频通信，卫星通信的码分多址，数字数据中的加密、加扰、同步、误码率测量等领域。而Gold序列是m序列的复合码序列，由两个码长相等、码时钟速率相同的m序列优选对的模2和序列构成.每改变两个m序列相对位移就可得到一个新的Gold序列。本实验研究针对r=5级得m序列、Gold序列生成过程以及它们的自相关函数和互相关函数的特性，同时可挑出平衡Gold序列，从而分析

4、以上集中序列的性质和应用。实验内容1 .以r=5145E为基础，抽取出其他的m序列，请详细说明抽取过程；2 .画出r=5的全部m序列移位寄存器结构，并明确哪些序列彼此是互反多项式；3 .在生成的m序列集中，寻找出m序列优选对，请确定优选对的数量，并画出它们的自相关和互相关函数图形；4 .依据所选取的m序列优选对生成所有Gold序列族，确定产生Gold序列族的数量，标出每个Gold序列族中的所有序列，并实例验证族内序列彼此的自相关和互相关特性；5 .在生成的每个Gold序列族内，明确标出平衡序列和非平衡序列，并验证其分布关系。第2章m序列实验2.1 m序列相关概念2.1.1 m码序列的定义r级非

5、退化的线性移位寄存器的组成如图2-1所示，r级线性移位寄存器的反馈逻辑可用二元域GF（2）上的r次多项式来表示2rfxc0c1xc2xLcrx,ci0,1（2-1）式（2-1）称为线性移位寄存器的特征多项式。对于动态线性寄存器，其反馈逻辑也可以用线性移位寄存器的递归线性关系式来表示aic1ai 1 c2ai 1 Lcrai r ,ci0,12-2)C0= 1输出akana0An ai图2-1 r级线性移位寄存器以式(2-1)为特征多项式的r级线性移位寄存器所产生的序列。假设以GF(2)域上r次多项式(2-1)为特征多项式的r级线性移位寄存器所产生的非零序列ai的周期为N2r1，称序列ai是最长

6、周期的r级线性移位寄存器序列，简称m序列。并且此r次特征多项式f(x)为r次本原多项式。m序列的最大长度决定于移位寄存器的级数，而码的结构决定于反馈抽头的位置和数量。不同的抽头组合可以产生不同长度和不同结构的码序列。有的抽头组合并不能产生最长周期的序列。2.1.2 m码序列的自相关函数自相关函数定义为Rn(T)二+T)dt,其中f(t)为捕获序列，也常用c(t)表示。互相关函数定义为Rn(T)=-I)dt,其中f(t)、g(t)为两个码序列。对二进制时间离散码序列，自相关函数和互相关函数的计算可简化如下：把两个码序列进行逐对和逐比特比较(模2加)，则自相关(或互相关)值为一致比特数减不一致比特

7、数，逐次改变T从0-(2门-1)Tc,则可得到自相关(或互相关)函数。如把相关值除以(2n-1),称为归一化相关函数。显然，自相关函数的最大值为1。根据m序列的性质，得到其自相关函数：HT二mNR(t)二，_；工羊出/二0,±±2|(2-3)m序列自相关函数曲线，如图2-2所示：2.2 m序列抽取结果分析本次实验是应用MATLABS言进行的仿真。2.2.1 m码序列的抽取及反多项式(一)m序列的抽取以r=5145E为基础，N=31,所以Nc=31,gcd(N,q)=1，并且q<31。此外，u(q)=u(2iq),因为m序列相同，只是相位不同，依据此方法，可以抽出所有的

8、m序列。例如：q=1,2,4,8,1632(mo。=1所产生的序列均和q=1时是同样的。则下一个m序列，从q=3开始计算；直至所有的q<31都取遍了，就可得出所有的m序列。按上述原理经编程可得：(1)r=5时：1248153122417&102015714282511221326211530292723135711可知每一行都是同一m序列在不同相位的情况。当q=1,3,5,7,11,15时是不同的m序列，则r=5时，可产生6个不同的m序列。(2)同理，r=6时，按上述原理经编程可得:124sIS325102口40IF34n2244255037132g5241193g2346295

9、8534331626159E5471511132331则q=1,5,11,13,23,31时是不同的m序列，则r=6时，可产生6个m序列。(3)同理，r=4时，按上述原理经编程可得：124S714131117则当q=1,7时是不同的m序列，则r=4时，可产生2个m1序列。通过以上实验结果，验证了图2-3中第二列，即：r=4时，有2个m序列优选对；r=5时，有6个m序列优选对；r=6时，有6个m序列优选对。N口何枚，of¥R值的iiirvqiienfnid幻十Mn322420563后6y图2-3m序列优选对及最大连通集个数(二)判断m序列的互反多项式以r=5145E为例，WJN=31。

10、令k=(N-1)/2=15,易知u(q)=u(q+31*i),i=1,2,3,因为u(q)也是以N=31为周期的序列。如果q+31*i为k=15的倍数,设为x倍，则q此时的取值与q=x时的本原多项式为互反多项式。例如：q=7,7+31=38,7+62=69;q=14,14+31=45=15*3,则，q=3和q=7的多项式为互反多项式。以下数字均表示q的值，经实验结果显示:r =5时，有1-15375-117-311-5 lb-4,即1-15 , 3-7和5-11 （q=4与q=1是同一 m序列）这三组是分别互为反多项式的;r =6时，有1-315-2311-1313-1123-531-4,即1

11、-31, 5-23, 11-13这三组是分别互为反多项式的;r =4时，有1-77-4,即1和7互为反多项式。2.2.2m码序列线性移位寄存器结构及相应序列经查表得出的本原多项式145E、375G、567H。q=1时，45E可表示为100101,对应的本原多项式是f(x)x5x21,其移位寄存器结构如下：图2-445E对应线性移位寄存器结构其反多项式为f(x)x5x31,即q=i5时，其移位寄存器结构如下:图2-545E反多项式线性移位寄存器结构如图2-6图所示，依次是以100000为起始状态放入移位寄存器中，当q=1和q=15时本原多项式产生的m序列。如果以别的起始状态放入移位寄存器中，产生

12、的m序列其实是一样的，只是可能差了几个序列的位置。图2-6第1组和第2组m序列q=3时，75G可表示为111101,对应的本原多项式是5432f(X)xxxx1,其移位寄存器结构如下：图2-7 75G对应线性移位寄存器结构如图2-9图所示，依次是以100000为起始状态放入移位寄存器中，当q=3和q=7时本原多项式产生的m序列17n7工?工ArJ&CCGJK0.21015202530图2-9第3组和第4组m序列q=5 时67H可表示为 110111对应的本原多项式是43x xx 1 ,即q=11时，其移位寄存器结构如下:图2-11 67H反多项式线性移位寄存器结构其

13、反多项式为f(x) x5如图2-12图所示，依次是以100000为起始状态放入移位寄存器中，当q=1和q=15时本原多项式产生的m序列图2-12 第5组和第6组m序列综上所述，以100000为起始状态放入以上6种位寄存器中，会产生以下的6种m序列:NO. L is口孙2 is1 0*3 is1 0IK).4 is1 0 国.5 i?1 0DO.fi is l 0ODO0 L 1 l 1 L 0 0 0 1 10 1110 10 (10 3 0 I 1 I fl j 1 110 0 10 0第3章m序列优选对实验3.1 m序列优选对的查找3.1.1 m序列优选对的定义m序列对的相关值可能是三值的

14、、四值的或者多值的。一些特殊的m序列对的11,t（n）,t（n）2互相关是三值的，此三值为：N，其中:n 11 2t(n)n 21 2r,n为奇数,n为偶数,被成为m序列的优选对。此外，m序列优选对也可以指在m序列集中，其互相关函数绝对值的最大值（称为峰值互相关函数）R' 最接近或达到互相关值下限（最小值）的一对 m序列。设m和bj分别是由本原多项式F1（x）和F2（x）产生的m序歹I，若和bj的峰值互相关函数（非归一化）满足为奇数T寸J22三十1为偶数且不是4的倍数则F1(x)和F2(x)所产生的m序列GJ和bj构成m序列优选对3.1.2 m序列优选对的查找结果(1)在本实验中，当r

15、=5时，既可以得出Et)Ie=9。将不同多项式产生的m序列依次求互相关函数，然后找到归一化后最大的数值乘以31,则是所需找的互相关值最大的数。结果显示如下：ini=10010;&&&&&45E10100;&&&&&45E反11110;&&&&&75G10111;&&&&&75G反11011;&&&&&67H11101;&&&&&67H反011ggs300gg

16、9g0009g00甘0g90000011thenyunberis1?可以看出，r=5时，有12个m序列优选对，因为对角线上方和下方是代表一对优选对，在此就不重复写出:序歹U 1代表45E(q=1),序列2代表45E反(q=15),序歹13代表75G(q=3),序列4代表75G反(q=7),序歹U 5代表67H(q=5),序列6代表67H反(q=11)。表3-1r=5时所有序列互相关函数最大值序列1序列2J予列3序列4序列5序列6序列1119999日29999序列31199序列499序列511序列6那么，当r=5时，上表中数值不大于9时对应的两个m序列为一对优选对,一共有12对。分别为:J丁列1

17、J丁列3日2J丁列1序列4序列2J丁列1日5序列2J丁列1序列6序列2序列3序列4序列5序列6序列3序列5序列3序列6序列4序列5序列4序列6用每个节点表示一个序列，若两个序列是一对优选对,则用线连接起来,r=5时，所有优选对连接图如图所示:图3-1 r=5时优选对连接图由此看出，最多有三个节点能保证，彼此间是一对优选对，如：1-3-5,2-3-6等，故r=5时，最大连通集为3。(2)当r=4时，止匕时并没有定义值。根据定义，t(n)=9,所以未进行归一化时，如果两个m序列互相关是三值的分别为，-1,-9,7则该两个m序列为一对优选对。r=4时，只有两个m序列，其互相关值结果为:ini=101

18、1;1101;由此可以看出，除了上述三值外，还有-5,所以此两个m序列不是一对优选对。则不存在最大连通集。(3)当r=6时，可以得出依(工)1叼=17。将不同多项式产生的m序列依次求互相关函数，然后找到归一化后最大的数值乘以63,则是所需找的互相关值最大的数。r=6时，其互相关值最大值结果为：ini=100001;110000;110011;111001;110110;101101;0151?232317002317172300015IT2300002317aQQ0Q15thenujubecis£可以看出，r=6时，有9个m序列优选对，如表3-2所示，因为对角线上方和卜方是代表一对优

19、选对，在此就不重复写出表3-2r=6时所有序列互相关最大值序列1日2序列3序列4序列5日6序列11517232317日223171723J予列3151723序列42317序列515序列6那么，当r=6时，上表中数值不大于17时对应的两个m序列为一对优选对,一共有9对。分别为：J丁列1J丁列2J丁列2序列4J丁列3序列4J丁列1J丁列3J丁列2抒列5J丁列3抒列5J丁列1J丁列6J丁列4序列6日5序列6用每个节点表示一个序列，若两个序列是一对优选对，则用线连接起来，r=5时，所有优选对连接图如下图所示，：由此看出，最多有两个节点能保证，彼此间是一对优选对，如：1-3,2-4等,故r=6时，最大连

20、通集为2。通过以上（1）、（2）、（3）的分析，充分验证了图2-3中第三列的值，即:r=4时，没有最大连通集；r=5时，最大连通集为3;r=6时，最大连通集为2。3.2 m序列优选对的自相关及互相关函数m序列的自相关结果图形如下:图3-3m序列自相关函数图由此可见，m序列有良好的自相关特性。在0点处达到最大，即为1,其他地方几乎为0,这与理论值极为接近。r=5时，可以产生6组m序列，分别产生的自相关函数均是图3-3,验证了教材中的关于m序列的自相关特性的理论分析。按照节的分析，r=5时，有12对m序列优选对，因此应该有12组互相关特性曲线，并且每个曲线是应该有31个值。具体12组m序列优选对的

21、互相关函数图如下：图3-4序列1与序列3互相关函数图图3-5序列1与序列4互相关函数图图3-6序列1与序列5互相关函数图图3-7序列1与序列6互相关函数图图3-8序列2与序列3互相关函数图图3-9序列2与序列4互相关函数图图3-12序列3与序列5互相关函数图图3-13序列3与序列6互相关函数图图3-14序列4与序列5互相关函数图图3-15序列4与序列6互相关函数图由于优选对较多，只随机抽取三个优选对的互相关函数值，非归一化结果为:第2个优选对互相关函数值为：>>31*r(21,)axis=CoIluuls1七hxo>UR上l1>T =17 T T T7T-1T-1fCa

22、luums14thfaugh2dTT?-1CuIluiixl27thtuu&h31TT-917-1第7个优选对互相关函数值为:rr -9ColujiLiLbr 14-1 rcoIlibum zr7-1>>3I*irh±e1r_eairr(7*J3TLM=CoIijmtls1Through13-i?-1?-i-i?-1-rtHruuuiL267-f?-9-9-i-97-1-1-07through1111T-I第11个优选对互相关函数值为：>>31*.inter_corr(11;)Colvums14throughColumi-LS 27tlir ou

23、63;"h-1-1-9-1由以上可见，再次验证了m序列有选对的互相关值为三值，并且当r=5时,CdluiTirix1througli13此三值为：-9,-1,7。第4章Gold序列实验Gold序列的定义指出：给定移位寄存器级数r,总可以找到一对互相关函数值是最小的码序列,采用移位寄存器相加的方法构成新的码组，具互相关旁瓣都很小而且自相关函数和互相关函数均是有界的。这个新的码组被称为Gold码或Gold序列。Gold序列是m序列的复合码序列，由两个码长相等、码时钟速率相同的m序列优选对的模2和序列构成。每改变两个m序列相对位移就可得到一个新的Gold序列。加上两个m序列，共有2r+1个

24、Gold序歹1。并且产生的Gold码序列的周期都是2r-1。Gold序列的生成及特性4.2.1 Gold序列生成结果根据节，当r=5时，本实验可以产生33个周期为31的Gold序列。由于数目过多，本次实验只求出序列1和序列3（即45E和75G生成序列）这对优选对生成的所有Gold序列族。Gold序列具体组成下：ini1=10010;%45Eini3=11110;%74G图4-1第13组gold序列图4-2第46组gold序列图4-3第79组gold序列图4-4第1012组gold序列图4-5 第1315组gold序列图4-6第1618组gold序列图4-7第1921组gold序列图4-8 第2

25、224组gold序列图4-9第2527组gold序列图4-10第2830组gold序列图4-11第3133组gold序列根据图4-1到图4-11,r=5时，选取序列序列1和序列3（即45E和75G生成序列）这对优选对生成的所有Gold序列，生成结果总结如下表：表4-145E和75G生成所有的Gold序列统计表序列Gold序列11000010101110110001111100110100210000101101010001110111110010013000000001101111011010001111110140100011110100010010010011010000511100100

26、000111000000010110001106101101011100001100100011100110170001110100101100101100001001000811001001010110110111100100010109101000110110000010011101110101110000101100111110101101111101101111010011001111001110010110100001112011000011011010011101010000111113011101110001011101010100014011111000100011010001

27、011015111110011110111001100100116001110110011101000010011017010110100101000000101000118011010101110010100110101019111100101011111110111011120101111101001001011111100121100110001000010001011111022000010111000111100001110123110000100000101010100110024001001101100100001110010025110101001010100100011000

28、02600101101100110011010110101100012701010001000000011111011111101102811101111010011011101101010101012900110000011010111100110000001003011011111111110001100011101011003110101000001100010100001011110003210010011110101011000000000100103310001110001001111110000101001114.2.2 Gold码自相关及互相关特性由于本实验中，取序列1和序列3

29、生成的Gold序列有33个，为方便观察自相关函数及互相关函数，只需要从中选取几组Gold序列即可达到实验的目的，因此，从33组序列中随机选择了3组，分别是第5组，第14组和第29组(1)该三组的Gold序列自相关曲线分别如下图：图4-12抽取三组Gold序列自相关曲线图第5组gold序列的自相关函数值（归一化后）为:Colunui】through31.00000.2256F0323-0.0323-0.0323-0.0323-0.2303-0.0323oliimriM3,0.22asthrough-0.0323IfiD.22BS0.22E8Q.22E8-0.0323-0.0323Calunu-L

30、f17-0.0323thLDUh-0.290324-0.03230.22580.22580.2258F03230.2258Columns25-0.0323through-0.290331-tk0323-0.0323-0,0323-0.0323C,22E8第14组gold序列的自相关函数值（归一化后）为:Colujiuxs1through81.0000必D323-0.0323上口3230323-0,2003-0.03230.22&8Columns9through16-0.2Q030.225E-0.0323D.2258-0.D3230.225®一口.2003-0.2503Coil

31、umns1FthrougTi24一02903-Ol2900D,E2B8以0323.2258-0.03230.2258F29Q3ColUJftTLS26through310.2268-0.0323-0.2&03032303230323一0323第29组gold序列的自相关函数值(归一化后)为:Columns11,tbrouh81.ooooth2258-0.0323-040323-0.03230.2268口.2258Q.2258Columns9tlirouffh16-0.0323-0,2B013-0,2皿0323.225S5225S0323-0k2003CulmrirLa17thtoogl

32、L24-OU2903-00323(J,225S口.2258F口际0；2903一口，2的3-0.0321Culiuiixi25llii-uugli310.225B0.225EU.Z25E-LL1323-I.UJ323.2253可以看出这3组gold序列的自相关特性曲线都是在0点达到最大值，即为1,而在其他点，也就是旁瓣处呈现出三值特性。组与第29组，分别如下图所示:(2)该三组的Gold序列互相关相关曲线，即第5组与第14组，第5组与第29组，第14图413抽取三组Gold序列互相关曲线图第5组与第14组的互相关函数值（归一化结果乘以31）为：ColumnsIthrough13-1r?77-1-

33、9-0-1-1-977Cali±>ii314tRrougR26-1-1-1-1-9-9-9-1'（f71-1Colums27through319-1-1i-第5组与第29组的互相关函数值为（归一化结果乘以31）:Coluwns1through13-1-1-1-177-9-1-177-17CaluJtris14-througk25-1-17?r-177-177-1-9Cqlujins27throoigh31-9-17-1-1第14组与第29组的互相关函数值（归一化结果乘以31）为：Columnstthrough13Columns14 through 26-17一1一1-1

34、F77-1-9-1-1Coluitins27 through 31-I-1-0-1-1可以看出Gold序列互相关函数具有三值性，并且互相关函数的最大值仍然等于9。平衡Gold码相关实验4.3.1 平衡Gold码的定义按平衡性来分，Gold码序列可以分为平衡码序列和非平衡码序列。在一周期内，平衡码序列中1码元与0码元的个数之差为1，非平衡码中1码元与0码元的个数之差多于1。平衡码具有更好的频谱特性，解决了码不平衡时直接序列系统的载波泄露大的问题。r1在平衡gold码序列中，码序列1的数量为2r1个。而在该码族中，有这样1数量的序列数有2r11个，即2r11个平衡码4.3.2 平衡Gold码的判定

35、r1根据4.4.1所述理论，当r=5时，N=31,有21=17个平衡gold码序列，每个平衡gold码中1的数量为2r1个。将一个gold序列的每个数值加在一起，如果等于16，即表示1的个数正好比0的个数多一个，也就是平衡gold序列。在本实验中，只计算m序列1和m序列3（即45E和75G生成序列）生成的Gold码序列族中的平衡gold码序列。结果如下：ini1=10010;%45Eini3=11110;%75GColmansIthrough1616 IC 1212 Ifi 1216IS 201616 ifiCclumns14 through 26162012Ifi 202012 IS 121

36、21616Columns27 through 33162012 2C121216Columnst through 13Colimns14thr augh26Coluisns27throughThe number of balance is第一组输出结果是每个gold序列中的个数,可以看出值为16的码序列有17个第二组输出结果是1的个数与0的个数若是差1则输出返回值为1,否则为0,若为1则是平衡gold码序列，并且个数为17个第5章总结实验小结从第二章可以看出，已知r级移位寄存器及相应的一个本原多项式，就可以抽出所有的m序列，同时可以求出这些m序列彼此间的互反多项式。第三章给出了如何寻找m序列优

37、选对的方式，并且所求的优选对都具有良好的自相关特性，在0点处达到最大，即为1,其他地方几乎为0,但是m序列数目少，证明了m序列优选对的互相关函数值是三值的。第四章给出了Gold序列的生成方式及对应序列，从随机抽取的3组gold序列可以看出，其自相关特性曲线都是在0点达到最大值，即为1，而在其他点，即旁瓣处呈现出三值特性。而其互相关特性也是具有三值特性。同时，判定了平衡gold序列以及个数。实验心得最开始的时候，我仔细看了实验的内容，对照着书本看所需考虑的知识点，并将不太理解的地方与同学进行了交流。再对实验内容进行了基本了解之后，我也大概了有了编程的思路。但是，因为我对matlab编程并不是很熟

38、悉，而且编程技巧也比较缺乏，所以在写程序时，经常出现一些意想不到的错误或者是一些想法不能轻易的用编程来实现。在向同学请教了之后，又经过反复修改及添加适当内容，虽然用了很长时间才能写出来，但是达到了预期想要的结果，在此也要谢谢这些同学给我的帮助。在做完这次试验之后，对扩频序列尤其是对m序列及gold序列的产生抽取等基础知识有了更深的理解和认识，并且将自己的理解思路写入程序，不仅让我对matlab编程有了进一步的了解，也让我对理论基础知识的理解加深了。总之，我觉得经过这次的编程实验，自己收获非常大。(1)咙由取m序列，确定m列的个数else endclc;endclear;m=1;s=0;mm=1

39、;r=5;i=1;N=2Ar-1;%产生第一组 m序列，以q=1开头%可取值的个数for n=1:rfor i=1:N-1A(m,n)=i;b=gcd(i,N);i=mod(2*i,N);if b=1ends=s+1;for i=2:Nendb=gcd(i,N);enddisp (A)%以矩阵形式输出分组后的数%C=min(A',1);%求出矩阵每行最小值%disp(C)pp=0;qq=0;则 flag=0;k=(N-1)/2;%依次寻找本原多项式对应的反多项式%for m=1:s/rflag=0;for i=1:25for n=1:rq=mod(A(m,n)+i*N,k);if q=

40、0qq=(A(m,n)+i*N)/k;pp=A(m,1);%显示互反多项式fprintf('%d-%dn' ,pp,qq )flag=1;%只要找到反多项式，则跳出break ;endendifb=1%僻与NIE质，则进行分组c=0;formm=1:mforn=1:rifi=A(mm,n)c=1;%如果新取的数与之前已确定的数是重复的，跳出break;elsec=0;endendifc=1break;endendifc=0m=m+1;%新取数与之前已确定的数无重复，则作为%新一行第一位forn=1:rA(m,n)=i；i=mod(2*i,N);endendifflag=1bre

41、ak;endend(2) %生m5列的函数：produce_m_seqfunctionm_seq=produce_m_seq(ini)len=length(ini);%移位寄存器的长度，本实验r=5L=2Alen-1;%列周期，本实验为31位regi=1,zeros(1,len-1);%初始寄存器内容10000m_seq(1)=regi(1);fori=2:Lnewregi(1:len-1)=regi(2:len);temp1=ini.*regi;temp2=sum(temp1);newregi(len)=mod(temp2,2);regi=newregi;m_seq(i)=regi(1);e

42、ndendini3=11110;%75Gini4=10111;%75Gini5=11011;%67Hini6=11101;%67H%六组m列并显示出来m_seq1=produce_m_seq(ini1);fprintf('NO.%disn',1);fprintf('%d',m_seq1);fprintf('n');m_seq2=produce_m_seq(ini2);fprintf('NO.%disn',2);fprintf('%d',m_seq2);fprintf('n');m_seq3=prod

43、uce_m_seq(ini3);fprintf('NO.%disn',3);fprintf('%d',m_seq3);fprintf('n');m_seq4=produce_m_seq(ini4);fprintf('NO.%disn',4);fprintf('%d',m_seq4);fprintf('n');(3) %产生六组ml列程序produce_mclc;clear;ini1=10010;%45Eini2=10100;%45Em_seq5=produce_m_seq(ini5);fprintf

44、('NO.%disn',5);fprintf('%d',m_seq5);fprintf('n');m_seq6=produce_m_seq(ini6);fprintf( 'NO.%d is n ',6);fprintf('%d ' ,m_seq6);fprintf('n ' );%画出六组mf列figure(1);subplot(2,1,1);stem(m_seq1);axis(0 32 0 1);subplot(2,1,2);stem(m_seq2);axis(0 32 0 1);figure(2

45、);subplot(2,1,1);stem(m_seq3);axis(0 32 0 1);subplot(2,1,2);stem(m_seq4);axis(0 32 0 1);figure(3);subplot(2,1,1);stem(m_seq5);axis(0 32 0 1);subplot(2,1,2);stem(m_seq6);axis(0 32 0 1);(4)%求两个序列的相关性函数：seq_corr()functioncorr_values=seq_corr(m_seq1,m_seq2)m_seq1=2*m_seq1-1;%序列进行非零处理m_seq2=2*m_seq2-1;co

46、rr_values=zeros(1,length(m_seq1);fori=0:length(m_seq1)-1corr_values(i+1)=sum(m_seq2.*circshift(m_seq1,0i)/length(m_seq1);endend(5)%r=4寸，计算其互相关函数值clc;clear;ini=zeros(2,4);m_seq=zeros(4,15);ini=1011;1101;%以数组形式列出序列移位寄存器%产生mf?列fori=1:2m_seq(i,:)=produce_m_seq(ini(i,:);end%tm序列互相关函数归一化结果inter_corr(g,:)=

47、for j=(i+1):6inter_corr=seq_corr(m_seq(1,:),m_seq(2,:);four_corr=15*inter_corr;fprintf('%d',four_corr);fprintf('n');(6) %r=5时，找n列优选对、互相关函数clc;clear;ini=zeros(6,5);m_seq=zeros(5,31);ini=10010;10100;11110;10111;11011;11101;%以数组形式列出序列的移位寄存器%fori=1:6m_seq(i,:)=produce_m_seq(ini(i,:);end%

48、互相关函数%g=1;x=0;fori=1:5seq_corr(m_seq(i,:),m_seq(j,:);M(i,j)=max(31*abs(inter_corr(g,:);figure(g);ifM(i,j)<9|M(i,j)=9x=x+1;plot(inter_corr(g,:);axis(031);g=g+1;endendend%以矩阵形式输出每个互相关函数最大值disp(M);fprintf('thenumberis%dn',x)%自相关函数%self_corr=zeros(5,31);h=1;fori=1:6self_corr(h,:)=seq_corr(m_s

49、eq(i,:),m_seq(i,:);h=h+1;endforj=1:2figure(j+12);for k=1:3forj=(i+1):6for i=1:6subplot(3,1,k);plot(self_corr(k+(j-1)*3),:);axis(032-1);endend(7) %fr=6时，找mi列优选对、互相关函数clc;clear;ini=zeros(6,6);m_seq=zeros(6,63);ini=100001;110000;110011;111001;110110;101101;%,%用列fori=1:6m_seq(i,:)=produce_m_seq(ini(i,:);endx=0;g=1;%互相关函数%inter_corr(g,:)=seq_corr(m_seq(i,:),m_seq(j,:);M(i,j)=max(63*abs(inter_corr(g,:);ifM(i,j)<17|M(i,j)=17x=x+1;g=g+1;endendend%以矩阵形式输出每个互相关函数最大值disp(M);fprintf(&#