版权声明：本文为CSDN博主「laomai」的原创文章

原文链接：https://blog.csdn.net/laomai/article/details/24342493

找了好多代码，这个步骤详细可以傻瓜操作，存着备用~

实验环境为matlab2013b

1、首先编写一个mseq.m文件,内容为:

function[mseq]=m_sequence(fbconnection)

n=length(fbconnection);

N=2^n-1;

register=[zeros(1,n-1) 1]; %移位寄存器的初始状态

mseq(1)=register(n); %m序列的第一个输出码元

for i=2:N

newregister(1)=mod(sum(fbconnection.*register),2);

for j=2:n

newregister(j)=register(j-1);

end;

register=newregister;

mseq(i)=register(n);

end

2、在matlab中打开这个文件

3、在matlab的命令行窗口中运行

primpoly(7,‘all’)

这个命令的目的是得到7阶M序列的所有本原多项式，如果想得到其他阶的m序列，自行修改数字7即可。

4、得到的结果为

Primitive polynomial(s) =

D7+D1+1

D7+D3+1

D7+D3+D2+D1+1

D7+D4+1

D7+D4+D3+D2+1

D7+D5+D2+D1+1

D7+D5+D3+D1+1

D7+D5+D4+D3+1

D7+D5+D4+D3+D2+D1+1

D7+D6+1

D7+D6+D3+D1+1

D7+D6+D4+D1+1

D7+D6+D4+D2+1

D7+D6+D5+D2+1

D7+D6+D5+D3+D2+D1+1

D7+D6+D5+D4+1

D7+D6+D5+D4+D2+D1+1

D7+D6+D5+D4+D3+D2+1

ans =

131

137

143

145

157

167

171

185

191

193

203

211

213

229

239

241

247

253

5、输入mseq([1 0 0 0 0 0 1])

6、便可得到一个127位的m序列，整理之后为

1000000111111101010

1001100111011101001

0110001101111011010

1101100100100011100

0010111110010101110

0110100010011110001

0100001100000

说明，

1 0 0 0 0 0 1这个输入参数的物理意义为各寄存器的初始状态，其数学意义是所选用的本原多项式D7+D1+1中D1到D7的各个系数。读者可以试试自行选用其他本原多项式，得到不同的m序列。

例如如果选用D7+D3+1，则只需调用

mseq([0 0 1 0 0 0 1])，得到的m序列为

1000000100100110100

1111011100001111111

0001110110001010010

1111101010100001011

0111100111001010110

0110000011011010111

0100011001000