matlab在伪随机码,MATLAB在伪随机码的生成及仿真中的应用

收稿日期:2003 - 01 - 09 　　第 21 卷 　第 3 期 计 　算 　机 　仿 　真 2004 年 3 月 　　 文章编号:1006 - 9348(2004)03 - 0149 - 04 MATLAB在伪随机码的生成及仿真中的应用 郭海燕 ,毕红军 (北方交通大学电子信息工程学院 ,北京 100044) 摘要:主要介绍 m 序列、Gold 序列、Kasami 序列和 JPL 序列的生成原理及其 MATLAB 的实现 ,并介绍如何用 MATLAB 中Simulink 对其进行硬件仿真 ,给出了相应自相关和互相关结果。关键词:序列;自相关;互相关 中图分类号:TN914. 42 ;TN911. 22 　　文献标识码:A 1 　引言 在 CDMA 数字移动通信和其他扩频通信中 ,地址码序列几乎都还有扩展频谱的作用 ,并且要求其有良好的伪随机特性和相关特性。本文中涉及到的 m 序列、Gold 序列、Kasami 序列和JPL 序列理论已经很成熟 ,而且 m序列和 Gold 序列作为地址码序列 ,已广泛应用于扩频 CDMA 系统。Kasami 序列由于自、互相关特性均较好的特点 ,且序列数量也很可观 ,正逐步得到应用。JPL 测距码 ,其构造特征非常适合于远距测距 ,码周期很长但同步很快。本文在简单叙述了伪随机码基本原理上 ,从硬件实现和仿真的角度 ,利用 MATLAB 的Simulink 工具箱对伪随机码进行软件实现及硬件仿真 ,得出了相应结果。 2 　m 序列、Gold 序列、Kasami 序列的基本原理与线性移位寄存器实现 2. 1 　m序列[1]发生器使用了 n 级移位寄存器作为延迟线 , 移位寄存器的输入是移位寄存器状态的函数 ,可以用特征多项式来表示 ,其一般形式为 : F(X) = C0 + C1X1 + C2X2 + ⋯+ Cn - 1Xn - 1 　　得到 m序列的本原多项式后[1] ,即可计算移位寄存器的初始化状态(C0 ,C1 ,C2 , ⋯,Cn) ,并得出其硬件发生电路。 2.2 　Gold 序列[1]是由两个 m 序列优选对[2]逐位模 2 加得到 ,当改变其中一个 m 序列相位(向后或向前移位) 时 ,可得到一新的 Gold 序列。 2. 3 　Kasami 序列[3]分为小集合 Kasami 序列和大集合 Kasami 序列 ,都是在 m 序列的基础上生成的。Kasami 序列的实现 , 关键在于如何得到取样后产生的序列。 2. 3. 1 　小集合 Kasami 序列的产生原理 :选定一周期为 2n - 1 (n 为偶数)的 m 序列 a ,因为 :2n - 1 = (2n/ 2 + 1) (2n/ 2 - 1) ,对序列 a 每隔 2n/ 2 + 1 个数进行抽样 ,得到一长度为 2n/ 2 - 1 的序列 ,可证明此序列仍为 m 序列 ,将此序列重复 2n/ 2 + 1 遍 , 得到一与 a 序列同长的序列 a′。把 a 与 a′序列逐位模 2 加就可以得到一个小集合 Kasami 序列 ,当改变其中一个序列相位(向后或向前移动)时 ,可得到一新的小集合 Kasami 序列。 Ks = ( a , b) = { a , ab, aTb , agT2 b, ⋯, agT2 n/ 2- 2 b}　　其中 T ib表示对b进行i位循环移位后得到的序列 ,g表示模二加。 2. 3. 2 　大集合 Kasami 序列的产生原理 :同小集合一样 ,要先得到 a 及 a’序列 ,再对序列 a 每隔 2(n + 2) /