基本信息

---

下载

   友情提示：（无需注册或充值，赞助后即可获取资源下载链接）

---

     嘿，亲！知识可是无价之宝呢，但咱这精心整理的资料也耗费了不少心血呀。小小地破费一下，绝对物超所值哦！如有下载和支付问题，请联系我们QQ(微信同号)：813200300

本次赞助数额为： 2 元　

支付宝支付

微信支付

   源码介绍

---

 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("a4,a3,a2,a1,a0")

for i in range(0,32):
    a.reverse()
    print(a)#倒序打印，即a4,a3,a2,a1,a0的顺序
    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)

#获得绘图的(x,y)数据
mx, my = [], []
for j in range(0,32):
    for k in range(1,101):
        my.append(m[j])
        mx.append(j k*0.01)

#根据数据绘制图形
fig = plt.figure(dpi=128, figsize=(10,6))
plt.plot(mx, my, c='red')

#设置图形的格式
plt.title('m sequence (31 1)', fontsize=24)
plt.xlabel('', fontsize=16)
#fig.autofmt_xdate()
plt.ylabel('', fontsize=16)
plt.tick_params(axis='both', which='major', labelsize=16)

plt.show()