通信原理(陈启兴版)第10章课后习题答案.doc

第10章 正交编码与伪随机序列
10.1 学习指导
10.1.1 要点
正交编码与伪随机序列的要点主要包括正交编码的概念、常见的正交编码和伪随机序列。
1. 正交编码的概念
对于二进制信号，用一个数字序列表示一个码组。这里，我们只讨论二进制且码长相同的编码。两个码组的正交性可用它们的互相关系数来表述。
设码长为n的编码中码元只取值+1和-1。如果x和y是其中的两个码组：x = (x1, x2, …, xn)，y = (y1, y2, …, yn)，其中，xi, yi ∈ {+1, -1}，i = 1, 2, …, n，则码组x和y的互相关系数被定义为
2.
如果码组x和y正交，则r(x, y) = 0。两两正交的编码称为正交编码。
类似地，我们还可以定义一个码组的自相关系数。一个长为n的码组x的自相关系数被定义为

其中，x的下标按模n运算，即xn＋k º xk。
在二进制编码理论中，常采用二进制数字“0”和“1”表示码元的可能取值。若规定用二进制数字“0”代替上述码组中的“-1”，用二进制数字“1”代替“+1”，则码组x和y的互相关系数被定义为

其中，a表示码组 x和y中对应码元相同的个数，b表示码组x和y中对应码元不同的个数。例如，对于4个码组：x1 = (1,1, 1, 1)，x2 = (1, 1, 0,0)，x3 = (1, 0, 0, 1)，x4 = (1, 0, 1, 0)，它们任意两者之间的相关系数都为0。
对于采用二进制数字“0”和“1”表示的码元，若用x的j次循环移位代替y，就得到x的自相关系数rx (j)。比如，如果一个长为n的码组x = (x1, x2, …, xn)，则y = (x1 + j, x2 + j, …, xn, x1, x2, …, xj)。根据上式计算出码组x和y的互相关系数就是码组x的自相关系数。
显然，无论是采用二进制数字“0”和“1”表示的码元，还是采用二进制数字“+1”和“-1”表示的码元，互相关系数和自相关系数都是在1与-1之间取值。若两个码组间的互相关系数r < 0，则称这两个码组互相超正交。如果一种编码中任意两码组之间均超正交，则称这种编码为超正交码。例如，对于3个码组：x1 = (+1, +1, +1)，x2 = (+1, -1, -1)，x3 = (-1, -1, +1)，由它们构成的编码是超正交码。
由正交编码和其反码构成的编码就是双正交编码。例如，4个码组：x1 = (+1, +1, +1, +1)，x2 = (+1, +1, -1, -1)，x3 = (+1, -1, -1, +1)，x4 = (+1, -1, +1, -1)，其反码为：y1 = (-1, -1, -1, -1)，x2 = (-1, -1, +1, +1)，x3 = (-1, +1, +1 -1)，x4 = (-1, +1, -1, +1)。这8个码组构成的编码就是双正交编码，任意两个码组之间的互相关系数r为0或-1。
2．常见的正交编码
常见的正交编码有Hadamard码矩阵、Walsh矩阵和伪随机序列等。
Hadamard码矩阵是法国数学家M. J. Hadamard于1893年首先构造出来的一种方阵，仅由元素+1和-1构成，而且其任意两行(列)之间是互相正交的，简记为H矩阵。H矩阵的最低阶数为2，即

为了简便起见，把上式中的＋1和－1简写为＋和－，上式就表示为

阶数为2k的高阶H矩阵从下列递推关系得出

其中，k为正整数，是直积运算。上式的直积运算是指将矩阵Hk / 2中的每一个元素用矩阵H2代替，比如，

H2矩阵、H4矩阵和H8矩阵都是对称矩阵，而且第一行和第一列的元素全为“＋”，我们把这样的H矩阵称为Hadamard码矩阵的正规形式，或称为正规Hadamard码矩阵。
在H矩阵中，交换任意两行或两列，或改变任一行或列中每个元素的符号，都不会影响矩阵的正交性质。因此，正规H矩阵经过上述各种交换或改变后仍为H矩阵，但不一定是正规的了。
按照递推关系式可以构造出所有2k阶的H矩阵。可以证明，高于2阶的H矩阵的阶数一定是4的倍数。不过，以4的倍数作为阶数是否一定存在H矩阵，这一问题并未解决。
H矩阵是正交方阵。如果把其中每一行看作是一个码组，则这些码组也是互相正交的，而整个H矩阵就是一种码长为n的正交编码，它包含n个码组。因为长度为n的编码共有2n个不同码组，如果只将这n个码组作为许用码组，其余(2n - n)个为禁用码组，则可以将其多余度用来纠错。这种编码在纠错编码理论中称为Reed-Muller码。
Walsh函数的定义常用三角函数法、拉德马赫函数乘积表示法、Hadamard矩阵表示法和递推公式法等。这里介绍Walsh函数的递推公式形式，其定义为

其中，j = 0,1,2, …, q = 0或1，[j/2]表示j/2的整数部分。
为了便于理解，做以下几点说明：
(1) 当把Wal(j, t)改成Wal(j, 2t)时，表示保持波形相对形状不变，只是将时基从-1/2 ≤ t ≤ 1/2压缩到-1/4 ≤ t ≤ 1/4；
(2) 当把Wal(j, 2t)改成Wal[j, 2(t ± 1/4)]时，表示保持波形相对形状不变，只是将波形向左(对应“+”号)或向右(对应“-”号)平移 1/4。
例如，Wal(5, t)应该根据Wal(2, t)递推出来，此时，k = 5, j = 2, q = 1, [j/2] = 1。

其中，

前八个Walsh函数中的任意两个函数都是正交的。将前N个Walsh函数在等距的N个点抽样，再将抽样值写成矩阵形式，即得N × N矩阵。例如，N = 8时，可以得到8 × 8矩阵：

如果把Walsh矩阵的每一行作为一个码组，就得到Walsh编码。
3. 伪随机序列
一方面，由于随机噪声的存在，通信系统的性能会变坏；另一方面，有时为了达到特殊的目的，通信系统中又需要噪声，比如保密通信中，让有用信号隐藏在随机噪声之中，以达到防止有用信号被截获。相同的随机噪声难以重复产生，这给通信接收方造成了困难，因为难以从隐藏在随机噪声之中提取出有用信号。后来，人们发明了伪随机序列，才使得这个难题得以解决。
伪随机序列具有类似于随机噪声的某些统计特性，同时又能够重复产生。目前广泛应用的伪随机噪声都是由周期性数字序列经过滤波等处理后得出的。我们将这种周期性数字序列称为伪随机序列，它有时又被称为伪随机信号和伪随机码。
常用的伪随机序列有m序列、M序列、二次剩余序列和双素数序列。
10.1.2 难点
正交编码与伪随机序列的难点主要是m序列的产生。
m序列是最长线性反馈移位寄存器序列的简称，它是由带线性反馈移存器产生的周期最长的一种序列。
一般来说，希望用尽可能少的级数产生尽可能长的序列。一个n级线性反馈移存器可能产生的最长周期等于(2n - 1)。这种最长的序列称为最长线性反馈移存器序列就是m序列。
反馈电路连接方法不同，移存器产生的序列的周期长度也会不同。
一般的线性反馈移存器原理方框图如图10-1所示，其中，各级移存器的状态用ai表示，ai = 0或1，i为非负整数。反馈线的连接状态用ci表示，ci＝1表示此线接通，ci＝0表示此线断开。反馈线的连接状态不同，输出序列的周期就可能不同。

根据上图，可以得到寄存器an-1的新状态为

事实上，上式是一个递推公式。显然，其它寄存器的新状态等于前一级寄存器的旧状态。
ci的取值决定了移存器的反馈连接和序列的周期长度，故ci是一个很重要的参量。现在将它们与一个n阶方程一一对应，让它们在为n阶方程的系数，即

这个n阶方程被称为特征方程或特征多项式。例如，特征方程

对应的连接系数取值为：c0＝c1＝c4＝1，c2＝c3＝0按照此特征方程构成的反馈移存器就图10-3所示。
任何一个寄存器的输出都可以作为一个伪随机序列。如果我们把寄存器an-1的输出序列{an, n = 0, 1, 2, …}的每个元素与一个代数方程建立一一对应的关系，即

这个代数方程被称为母函数。
关于递推方程、特征方程和母函数之间的关系，由几个定理来阐述。
定理10-1 如果多项式u(x)的阶数低于特征方程f(x)的阶数，该特征方程f(x)对应的母函数为G(x)，则

定理10-2 一个n级线性反馈移存器的状态具有周期性，且周期p £ 2n－1。
定理10-3 如果序列A = { ak, k = 0, 1, 2, …}具有最长周期p = 2n - 1，则其特征多项式f(x)应为阶数位n的既约多项式 (所谓既约多项式是指不能分解因子的多项式)。
定理10-4 如果n级线性反馈移存器的特征多项式f (x)是既约的，则由其产生的序列A = { ak, k = 0, 1, 2, …}的周期等于使(xp + 1)被f (x)整除的最小正整数p。
10.2 习题详解
10-1 已知一个4级线性反馈移位寄存器的特征方程为f(x) = x4 + x3 + 1，假设4个移位寄存器的初始状态为(a3, a2, a1, a0) = (1, 0, 1, 0)，试画出其组成方框图，并列出4个移位寄存器状态更新表。
解 组成方框图和状态更新表如图答10-1所示。

10-2 某3级线性反馈移位寄存器的特征方程为f(x) = x3 + x + 1，试证明该特征方程是本原多项式。
证明f(x)为3阶多项式，如果它能分解因子，则其因子只有x, (x + 1)，(x2 + 1)，和(x2 + x + 1)四种可能。不难验证，f(x)不能被上述四种因子整除，所以，f(x)是既约的。
3级线性反馈移位寄存器产生的序列的最长周期为p = 23 - 1 = 7。由于

即(2p+ 1)能被f(x)整除。
不难验证，(26 + 1)、(25 + 1)、(24 + 1)和(23 + 1)都不能被f(x)整除。
综上所述，f(x)是本原多项式。
10-3 某4级线性反馈移位寄存器的特征方程为f(x) = x4 + x3 + x2 + x + 1，试证明该4级线性反馈移位寄存器产生的序列不是m序列。
证明 4级线性反馈移位寄存器产生的序列的最长周期为p = 24 - 1 = 15。由于

即(25 + 1)能被f(x)整除，而5 < p，因此，f(x)不是本原多项式，也就是说，该4级线性反馈移位寄存器产生的序列不是m序列。
10-4 由一个8级线性反馈移位寄存器产生的m序列中，在一个周期内有可能产生哪些长度的游程？每种长度游程的数目为多少？
解 8级线性反馈移位寄存器产生的m序列的周期为

该m序列的游程总数为

游程长度k = 1的数目为

游程长度k = 2的数目为

游程长度k = 3的数目为

游程长度k = 4的数目为

游程长度k = 5的数目为

游程长度k = 6的数目为

游程长度k = 7的数目为

游程长度k = 8的数目为1。