mseq matlab,Matlab生成M序列
版權聲明:本文為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
總結
以上是生活随笔為你收集整理的mseq matlab,Matlab生成M序列的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: mcdlc.exe - mcdlc是什么
- 下一篇: cpu影响matlab仿真速度吗,Pro