关于matlab生成M序列的代码

#关于matlab生成M序列的代码
参考了相关书籍和网上matlab生产m序列的代码，我发现和书上的定义有些不一样（系统辨识原理及其运用-李彦俊）于是我就自己学了一个m文件来供大家参考。使用起来简单方便。以下就是函数
function [seq]=msep(coef)
%==========================================================================
% 此函数用来生成最大长度m序列
% coef为特征多项式向量
%==========================================================================
m=length(coef);
B=coef;
len=2^m-1; % 得到序列的长度
seq=[];
p=0;
for i=1:len
p=double(xor(B(3),B(4)));
B=[B(4),B(1:3)];
B(1)=p;
seq(i)=B(1);
end
//希望大家印用的时候可以表明出处，学习，有什么不足请多指点