版權(quán)聲明：本文為CSDN博主「laomai」的原創(chuàng)文章

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

找了好多代碼，這個步驟詳細可以傻瓜操作，存著備用~

實驗環(huán)境為matlab2013b

1、首先編寫一個mseq.m文件,內(nèi)容為:

function[mseq]=m_sequence(fbconnection)

n=length(fbconnection);

N=2^n-1;

register=[zeros(1,n-1) 1]; %移位寄存器的初始狀態(tài)

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,‘a(chǎn)ll’)

這個命令的目的是得到7階M序列的所有本原多項式，如果想得到其他階的m序列，自行修改數(shù)字7即可。

4、得到的結(jié)果為

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這個輸入?yún)?shù)的物理意義為各寄存器的初始狀態(tài)，其數(shù)學(xué)意義是所選用的本原多項式D7+D1+1中D1到D7的各個系數(shù)。讀者可以試試自行選用其他本原多項式，得到不同的m序列。

例如如果選用D7+D3+1，則只需調(diào)用

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

1000000100100110100

1111011100001111111

0001110110001010010

1111101010100001011

0111100111001010110

0110000011011010111

0100011001000

總結(jié)

以上是生活随笔為你收集整理的mseq matlab,Matlab生成M序列的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。