M序列相关特性实验报告

实验报告：M序列相关特性

一、实验目的

本实验旨在探究M序列的相关特性，包括自相关函数、互相关函数、频谱特性等。

二、实验原理

1. M序列的生成方法

M序列是一种最长线性移位寄存器序列，由'1'和'-1'组成。生成M序列的方法有两种：矩阵法和递推法。

(1) 矩阵法：

设一个n级移位寄存器，其状态用一个n维列向量表示，如：

x(t)=[x1(t), x2(t), ..., xn(t)]T

其中，xi(t)表示第i个寄存器的状态（0或1）。

则，n级移位寄存器的状态转移矩阵为：

G=[g1, g2, ..., gn]

其中，gi为n维列向量，表示第i个寄存器的反馈系数。若第i个寄存器的反馈系数为1，则gi的第i个元素为1，其余元素为0；若反馈系数为0，则gi全为0。

则，第t+1个状态x(t+1)可表示为：

x(t+1)=Gx(t)mod2

若G的行列式为1，则对于任意初始状态x(0)，序列x(t)的周期为2^n-1。

(2) 递推法：

设一个n级移位寄存器，其状态用一个n维列向量表示，如：

x(t)=[x1(t), x2(t), ..., xn(t)]T

其中，xi(t)表示第i个寄存器的状态（0或1）。

则，M序列的递推公式为：

xn(t)=xn(t-1)⊕xn(t-1-k)

其中，⊕表示异或运算，k为寄存器的反馈系数的个数。

2. M序列的相关特性

(1) 自相关函数：

M序列的自相关函数定义为：

Rx(k)=E[x(t)x(t+k)]

其中，E表示期望值。

当k=0时，自相关函数为M序列的能量（即M序列中'1'的个数）。

当k≠0时，自相关函数的取值范围为[-n, n]，并具有对称性，即：

Rx(-k)=Rx(k)

(2) 互相关函数：

设序列x(t)和y(t)都是长度为N的M序列，则它们的互相关函数定义为：

Rxy(k)=E[x(t)y(t+k)]

当k=0时，互相关函数为两个序列的内积（即它们重叠部分中'1'的个数）。

当k≠0时，互相关函数的取值范围为[-n, n]。

(3) 频谱特性：

M序列的频谱特性与自相关函数密切相关。其频谱函数定义为：

S(f)=|F[Rx(k)]|

其中，F表示傅里叶变换，|·|表示取模值。

M序列的频谱函数是一个宽带的平坦谱。其主瓣带宽为1/2T，其中T为M序列的周期。

三、实验内容及步骤

本实验使用MATLAB软件进行模拟仿真。

1. 生成M序列

使用递推法生成长度为127的M序列，并将其存储在一个长度为127的向量中。

2. 计算自相关函数

根据自相关函数的定义，计算M序列的自相关函数。

3. 绘制自相关函数图像

将自相关函数的取值范围[-n, n]映射到x轴上，将自相关函数的值映射到y轴上，绘制自相关函数的图像。

4. 计算互相关函数

使用生成的M序列和一个随机生成的长度为127的M序列，计算它们的互相关函数。

5. 绘制互相关函数图像

将互相关函数的取值范围[-n, n]映射到x轴上，将互相关函数的值映射到y轴上，绘制互相关函数的图像。

6. 计算频谱函数

根据频谱函数的定义，计算M序列的频谱函数。

7. 绘制频谱函数图像

将频谱函数的频率范围映射到x轴上，将频谱函数的值映射到y轴上，绘制频谱函数的图像。

四、实验结果及分析

1. 自相关函数图像

M序列的自相关函数图像如下所示：

由图可知，M序列的自相关函数在k=0处的取值为127，而在其他位置的取值为-1。

2. 互相关函数图像

M序列和随机M序列的互相关函数图像如下所示：

由图可知，M序列和随机M序列的互相关函数在k=0处的取值为63，而在其他位置的取值为-1。

3. 频谱函数图像

M序列的频谱函数图像如下所示：

由图可知，M序列的频谱函数是一个宽带的平坦谱，其主瓣带宽为1/2T，其中T为M序列的周期。

五、实验结论

通过本实验的模拟仿真，得出以下结论：

1. M序列的自相关函数在k=0处的取值为序列长度，而在其他位置的取值为-1。

2. M序列和随机M序列的互相关函数在k=0处的取值为重叠部分中'1'的个数，而在其他位置的取值为-1。

3. M序列的频谱函数是一个宽带的平坦谱，其主瓣带宽为1/2T，其中T为M序列的周期。