5阶移位寄存器的m序列伪随机数的生成

 5 阶移位寄存器的本原多项式为 f(x) = x^5 x^2 1，各系数是 C5 = C2 = C0 = 1， C1 = C3 = C4 = 0。
 a5 ⊕ a3 ⊕ a0 = 0, 即 a5 = a3 ⊕ a0，
其中 ⊕ 表示“模二加”运算，或“异或”运算。

from matplotlib import pyplot as plt

a = [0 0 0 0 1]#5级寄存器初始化
m = []#m序列周期为31（0～30），但为了看重复性，多看1位。

print（“a4a3a2a1a0“）

for i in range（032）:
    a.reverse（）
    print（a）#倒序打印，即a4a3a2a1a0的顺序
    a.reverse（）
    a.append（a[3]^a[0]）#a5=a3^a0
    temp = a.pop（0）#把a0输出
    m.append（temp）#a0输出填入m序列

print（“m序列为：（多加一位）“）
print（m）

m_str = [str（x） for x in m]#首先每个数字变成字符
m_bin = ‘‘.join（m_str）#然后列表字符串化
m_bin = “0b“ + m_bin
print（“\n m序列的二进制数为“ + m_bin）

m_int = int（m_bin 2）
print（“\n m序列的十进制数为“）
print（m_int）

s0 = 0b10011001100110011001100110011001#信源码举例
print（“信源码的二进制表达“）
print（bin（s0））
print（“信源码