在matlab中生成m序列

               

实验环境为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) = 
 
D^7+D^1+1
D^7+D^3+1
D^7+D^3+D^2+D^1+1
D^7+D^4+1
D^7+D^4+D^3+D^2+1
D^7+D^5+D^2+D^1