Berlekamp Massey(BM)算法
下面是BM算法的步骤
主要是实现过程:
这里的函数都是关于x的多项式,我们很自然的想到用数组去表达x的多项式,例如:a[]={1,0,1}表达的就是y=1*(x^0)+0*(x^1)+1*(x^2);
因此在这里我们把T(x),f(x),B(x)分别用整数数组 t,f,b 去承接。
代码解释:
next是求下一步离差的,按照图片中的写就可以;
主函数中主要是执行step2.2的
注意!!!
N是一开始输入的,在代码中是用i去表示n不断增加的这个过程。代码中所有的小n就是图片中原理中的N,偷了个懒。
还需要注意的是,数组交接时如何取范围,因为我们初始时给出的是N长的二元序列,而且求的是最短级数的LFSR,所以数组交接时循环的次数不必太大,为N即可。
#include<stdio.h>