資源描述：
《基于matlab的偽隨機(jī)序列實現(xiàn).ppt》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫。
1、基于matlab的偽隨機(jī)序列實現(xiàn)如果一個序列，一方面它是可以預(yù)先確定的，并且是可以重復(fù)地生產(chǎn)和復(fù)制的；一方面它又具有某種隨機(jī)序列的隨機(jī)特性（即統(tǒng)計特性），我們便稱這種序列為偽隨機(jī)序列。偽隨機(jī)序列具有良好的隨機(jī)性和接近于白噪聲的相關(guān)函數(shù)，并且有預(yù)先的可確定性和可重復(fù)性。這些特性使得偽隨機(jī)序列在通信加密、雷達(dá)信號設(shè)計、通信系統(tǒng)中得到了廣泛的應(yīng)用。一、偽隨機(jī)序列簡介二、偽隨機(jī)序列產(chǎn)生的原理m序列是最長線性反饋移位寄存器序列的簡稱，它是由帶線性反饋的移位寄存器產(chǎn)生的周期最長的一種序列。m序列的生成可用移
2、位寄存器序列發(fā)生器的本原多項式來確定，一個本原多項式對應(yīng)一個最大長度序列，即對應(yīng)一個m序列。本原多項式可在matlab中生成，打開Matlab程序，輸入primpoly(n,’all’)能得到n階移位寄存器所對應(yīng)的所有的本原多項式。輸入primpoly(4,’all’)能得到4階移位寄存器所對應(yīng)的所有的本原多項式。輸出結(jié)果為：Primitivepolynomial(s)=D^4+D^1+1D^4+D^3+1以多項式為例：（下圖為對應(yīng)的移位寄存器反饋連線圖）y2y3y1y4由上圖可以看出，寄存器狀
3、態(tài)滿足：y1(i)=y1(i-1)y4(i-1);y2(i)=y1(i-1);y3(i)=y2(i-1);y4(i)=y3(i-1);給寄存器組賦予一個非全零的初始狀態(tài)，y4的輸出即可作為碼周期為15的一個m序列。輸入mseq=mseries([1001])可得mseq=Columns1through100001111010Columns11through1511001偽隨機(jī)碼生成函數(shù)mseries.mfunction[mseq]=mseries(fbconnection)n=length(fb
4、connection);%所需要的移位寄存器的長度N=2^n-1;%偽隨機(jī)碼的周期register=[zeros(1,n-1)1];%初始化寄存器的狀態(tài)fori=1:Nnewregister(1)=mod(sum(fbconnection.*register),2）forj=2:nnewregister(j)=register(j-1);endregister=newregister;mseq(i)=register(n);end三、偽隨機(jī)序列的自相關(guān)性衡量一個m序列與它的j次移位序列之間的相關(guān)
5、程度，并把它叫做m序列()的自相關(guān)函數(shù)。記作，當(dāng)采用二進(jìn)制數(shù)字0和1代表碼元的可能取值時，（其中A表示，D表示）由m序列的均衡性和移位相加特性可知：自相關(guān)程序：a=[000111101011001];b=a;m=length(b);N=100;x=1:N;fork=1:Nc=xor(a,b);D=sum(c);A=m-D;R(k)=(A-D)/(A+D);plot(k,R(k),'ro')b=[b(m),b(1:m-1)];holdonendplot(x,R,'b')gridon從圖中可以看出在
6、15、30、45、60、75、90處得到了最大的相關(guān)特性，其他值處自相關(guān)值接近于零，此仿真驗證了m序列具有良好的自相關(guān)特性。謝謝
