移动通信中m序列的Matlab仿真

Gilgamesh Emiya

已于 2022-11-05 17:58:03 修改

阅读量5.4k

点赞数 11

分类专栏：仿真实验文章标签： matlab 开发语言

于 2022-08-06 12:26:09 首次发布

本文链接：https://blog.csdn.net/qq_45068436/article/details/126191966

版权

仿真实验专栏收录该内容

2 篇文章 2 订阅

订阅专栏

本篇文章简要介绍了移动通信中的m序列，并通过Matlab进行了简单的仿真实验。

m序列的介绍：二进制m序列是一种重要的伪随机序列，有优良的自相关特性，有时称为伪噪声(PN)序列。“伪”的意思是说这种码是周期性的序列，易于产生和复制，但其随机性接近于噪声或随机序列。m序列在扩展频谱及码分多址技术中有着广泛的应用，并且在m序列基础上还能构成其它的码序列。

m序列的含义：m序列是最长线性移位寄存器序列的简称。顾名思义，m序列是由多级移位寄存器或其延迟元件通过线性反馈产生的最长的码序列。在二进制移位寄存器中，若n为移位寄存器的级数，n级移位寄存器共有2n个状态，除去全0状态外还剩下2n-1中状态，因此它能产生的最大长度的码序列为2n-1位，也就是说，一个n级线性反馈移位寄存器产生的最长周期等于 2n-1 。

m序列的产生：如下图1所示是由n级移位寄存器构成的码序列发生器。寄存器的状态决定于时钟控制下输入的信息（“0”或“1”），例如第i级移位寄存器状态决定于前一时钟脉冲后的第i－1级移位寄存器的状态。

图1：n级循环序列发生器的模型

图中C0，C1，…，Cn均为反馈线，其中C0＝C1＝1，表示反馈连接。因为m序列是由循环序列发生器产生的，因此C0和Cn肯定为1，即参与反馈。而反馈系数C1，C2，…，Cn－1，若为1，参与反馈；若为0，则表示断开反馈线，即开路，无反馈连线。

一个线性反馈移动寄存器能否产生m序列，决定于它的反馈系数Ci (i=1,2,…,n)，下表中列出了部分m序列的反馈系数，按照如下表一中的系数来构造移位寄存器，就能产生相应的m序列。

表1：部分m序列反馈系数表

根据表1中的八进制的反馈系数，可以确定m序列发生器的结构。以7级m序列反馈系数Ci=(211)8为例，首先将八进制的系数转化为二进制的系数即Ci=(010001001)2，由此我们可以得到各级反馈系数分别为：C0=1，C1=0，C2=0，C3=0，C4=1，C5=0，C6=0，C7=1，由此就很容易地构造出相应的m序列发生器。根据反馈系数，其他级数的m序列的构造原理与上述方法相同。

m序列的性质的分析

m序列具有以下性质：

m序列的随机性：在m序列中，码元为“1”的数目和码元为“0”的数目只相差1个。比如m序列1010011，其中码元为“1”的有4个，为“0”的有三个，即“1”和“0”相差一个。

m序列的游程分布：m序列中取值相同的那些相继的元素合称为一个“游程”。游程中元素的个数称为游程长度。n级的m序列中，总共有2^(n-1)个游程，其中长度为1的游程占总游程数的1/2，长度为2的游程占总游程数的1/4，长度为k的游程占总游程数的2^k。且长度为k的游程中，连0与连1的游程数各占一半。

m序列的自相关函数：自相关函数只用于表征一个信号与延迟 $\tau$ 后自身信号的相似性。采用二进制的码序列，长度(周期为)P的码序列x的自相关函数为

$R_{x}\left ( \tau \right )= \sum_{i=1}^{P}x_{i}\cdot x_{i+\tau }$

，由该公式可知，在二进制序列情况下，只需要比较原序列与移位后序列对应位的码元即可，自相关函数为 $R_{x}\left ( \tau \right )=A-D$ ，A为对应位码元相同数目，D为对应位码元不同的数目，自相关系数为 $\rho \left ( \tau \right )=\frac{A-D}{P}=\frac{A-D}{A+D}$ ，对于m序列，其码长为P=2^n-1，在这里P也等于码序列中码元数，即“0”和“1”的综合。m序列的自相关系数为

$\rho \left ( \tau \right )=\left\{\begin{matrix} 1&,\tau =0 \\ -\frac{1}{P}&,\tau \neq 0,\tau =1,2,...,P-1 & \end{matrix}\right.$

m序列的互相关函数：两个码序列的互相关函数是两个不同码序列一致程度(相似度)的度量，它也是位移量的函数。研究表明，两个长度周期相同，由不同反馈系数产生的m序列，其互相关函数(或互相关系数)与自相关函数相比，没有尖锐的二值特性，是多值的。在二进制情况下，假设码序列的周期为P的两个m序列，其互相关函数为

$R_{xy}\left ( \tau \right )=A-D$

Matlab仿真主要流程：

（1）设计产生两个阶数n>5的m序列(反馈系数自选)

根据实验要求，我设计了两个阶数为6的m序列，6级寄存器的初始状态定义为数组an=[1 1 1 0 1 1]，两个移位寄存器的反馈系数分别为(103)8和(147)8，转化为二进制数分别为(1000011)2和(1100111)2，（此处括号后的数字为进制数），在matlab程序中我们用数组cn和cn2存放反馈系数，通过length函数算出寄存器的级数，定义变量L根据初始状态算出寄存器的级数，由m(1)=an(len)和n(1)=an(len)存放两个序列第一个输出码元，如图2所示。