MATLAB仿真m序列和Gold序列自相關(guān)與互相關(guān)

1、MATLAB仿真m序列和Gold序列自相關(guān)與互相關(guān)MATLAB 是美國MathWorks公司出品的商業(yè)數(shù)學(xué)軟件，用于算法開發(fā)、數(shù)據(jù)可視化、數(shù)據(jù)分析以及數(shù)值計(jì)算的高級技術(shù)計(jì)算語言和交互式環(huán)境，主要包括MATLAB和Simulink兩大部分。MATLAB可以進(jìn)行矩陣運(yùn)算、繪制函數(shù)和數(shù)據(jù)、實(shí)現(xiàn)算法、創(chuàng)建用戶界面、連接其他編程語言的程序等，主要應(yīng)用于工程計(jì)算、控制設(shè)計(jì)、信號處理與通訊、圖像處理、信號檢測、金融建模設(shè)計(jì)與分析等領(lǐng)域。本章正是運(yùn)用MATLAB來仿真m序列和Gold序列的相關(guān)特性，以及OCDMA系統(tǒng)的誤碼率同用戶數(shù)的關(guān)系曲線。1 仿真過程在理論分析的基礎(chǔ)上,下面使用附錄上兩段程序,通過MAT

2、LAB仿真得出m序列和Gold序列的自己相關(guān)性。這段m序列產(chǎn)生程序采用了8個(gè)移位寄存器，將最后兩個(gè)移位寄存器的值進(jìn)行異或處理反饋給第一個(gè)移位寄存器，然后向前移位，輸出最后一個(gè)移位寄存器的值，Gold序列的產(chǎn)生只是將兩個(gè)m序列中的一個(gè)進(jìn)行延時(shí)移位，再進(jìn)行異或，產(chǎn)生的主要原理和m序列并無較大差異。在這兩段程序個(gè)前半部分m序列和Gold序列生成的基礎(chǔ)上，只要將y1=xcorr（x1）改為y1=xcorr（x1，x2）即可求出它們的互相關(guān)仿真。2 仿真結(jié)果在系統(tǒng)中采用上述序列仿真得到自相關(guān)和互相關(guān)特性曲線如圖1 , 2及圖3 所示。圖1 m 序列的自相關(guān)曲線圖4-2 m 序列與Gold序列的自相關(guān)曲線

3、圖4-3 m 序列與Gold 序列的互相關(guān)曲線附錄程序1X1=1;X2=0;X3=1;X4=0; %移位寄存器輸入Xi初T態(tài)（0101）， Yi為移位寄存器各級輸出m=120; %置M序列總長度for i=1:m %1# Y8=X8; Y7=X7; Y6=X6; Y5=X5; Y4=X4; Y3=X3; Y2=X2; Y1=X1; X8=Y7; X7=Y6; X6=Y5; X5=Y4; X4=Y3; X3=Y2; X2=Y1; X1=xor(Y7,Y8); %異或運(yùn)算 if Y8=0 U(i)=-1; else U(i)=Y8;end endM=U%繪圖i1=it=1:1:i1;x1=(2*M

4、)-1;%將運(yùn)行結(jié)果m序列M從單極性序列變?yōu)殡p極性序列y1=xcorr（x1）；%求自相關(guān)性t=1：1：i1；plot（t,y1(1:i1)）;axis(1,120,-12,288)%繪出信號的相關(guān)圖gridxlabel('t')ylabel('相關(guān)性')title('移位寄存器產(chǎn)生的M序列的相關(guān)性')程序2function c=gold()n=7;a=1 1 1 1 1 1 1 1;co=;for v=1:2n-1co=co,a(1);a(8)=mod(a(5)+a(1),2);a(1)=a(2);a(2)=a(3);a(3)=a(4);a(4

5、)=a(5);a(5)=a(6);a(6)=a(7);a(7)=a(8);endm1=co;b=1 0 1 0 0 0 0 1;co=;for v=1:2n-1co=co,b(1);m=mod(b(5)+b(1),2);p=mod(b(6)+m,2);b(8)=mod(b(5)+b(1),2);b(1)=b(2);b(2)=b(3);b(3)=b(4);b(4)=b(5);b(5)=b(6);b(6)=b(7);b(7)=b(8);endm2=co;c=xor(m1,m2);x2=(2*c)-1;%將運(yùn)行結(jié)果Gold序列c從單極性序列變?yōu)殡p極性序列y1=xcorr（x2）；%求自相關(guān)性t=1:1:120；plot（t,y1(1:120)）;axis(1,120,-12,288)%繪出信號的相