在matlab中生成m序列

2021-09-20 06:45:41 字數 3751 閱讀 9912

實驗環境為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+1

d^7+d^5+d^3+d^1+1

d^7+d^5+d^4+d^3+1

d^7+d^5+d^4+d^3+d^2+d^1+1

d^7+d^6+1

d^7+d^6+d^3+d^1+1

d^7+d^6+d^4+d^1+1

d^7+d^6+d^4+d^2+1

d^7+d^6+d^5+d^2+1

d^7+d^6+d^5+d^3+d^2+d^1+1

d^7+d^6+d^5+d^4+1

d^7+d^6+d^5+d^4+d^2+d^1+1

d^7+d^6+d^5+d^4+d^3+d^2+1

ans =

131137

143145

157167

171185

191193

203211

213229

239241

247253

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這個輸入引數的物理意義為各暫存器的初始狀態,其數學意義是所選用的本原多項式d^7+d^1+1中d^1到d^7的各個係數。讀者可以試試自行選用其他本原多項式,得到不同的m序列。

例如如果選用d^7+d^3+1,則只需呼叫

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

1000000100100110100

1111011100001111111

0001110110001010010

1111101010100001011

0111100111001010110

0110000011011010111

0100011001000

實驗環境為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+1

d^7+d^5+d^3+d^1+1

d^7+d^5+d^4+d^3+1

d^7+d^5+d^4+d^3+d^2+d^1+1

d^7+d^6+1

d^7+d^6+d^3+d^1+1

d^7+d^6+d^4+d^1+1

d^7+d^6+d^4+d^2+1

d^7+d^6+d^5+d^2+1

d^7+d^6+d^5+d^3+d^2+d^1+1

d^7+d^6+d^5+d^4+1

d^7+d^6+d^5+d^4+d^2+d^1+1

d^7+d^6+d^5+d^4+d^3+d^2+1

ans =

131137

143145

157167

171185

191193

203211

213229

239241

247253

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這個輸入引數的物理意義為各暫存器的初始狀態,其數學意義是所選用的本原多項式d^7+d^1+1中d^1到d^7的各個係數。讀者可以試試自行選用其他本原多項式,得到不同的m序列。

例如如果選用d^7+d^3+1,則只需呼叫

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

1000000100100110100

1111011100001111111

0001110110001010010

1111101010100001011

0111100111001010110

0110000011011010111

0100011001000

在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 reg...

matlab生成m序列的方法

2.2利用idi nput idinput idinpu t函式引言m序列屬於偽隨機序列的一種。在通訊領域應用較為廣泛。由於其具有很好的偽雜訊特性,因此在保密通訊中,可以將其應用在高可靠性的保密通道中。另外,雷達領域也有廣泛應用,偽隨機相位編碼脈衝 連續波雷達中用到的調製訊號正是m序列。m序列的週期...

在delphi中生成GUID

什麼是 guid guid 的格式為8 4 4 4 12 xx x x x 其中每個 x 是 0 9 或 a f 範圍內的乙個十六進製制的數字。例如 6f9619ff 8b86 d011 b42d 00c04fc964ff 即為有效的 guid 值。為什麼要使用 guid 世界上的任何兩台計算機都不...