【文章內(nèi)容簡(jiǎn)介】 (210) （5）m序列發(fā)生器下面給出產(chǎn)生m序列的條件: ① r級(jí)移位寄存器產(chǎn)生的碼, 周期n＝2r1, 其特征多項(xiàng)式必然是不可約的， 即不能再因式分解而產(chǎn)生最長(zhǎng)序列。 因此, 反饋抽頭不能隨便決定, 否則將會(huì)產(chǎn)生短碼。② 所有的次數(shù)r＞1的不可約多項(xiàng)式f(x)必然能除盡1＋Xn, 因?yàn)閍N(x)=(1+xn)／f(x)。 ③ 如果2r1是一個(gè)素?cái)?shù), 則所有r次不可約多項(xiàng)式產(chǎn)生的線性移位寄存器序列, 一定是m序列, 產(chǎn)生這個(gè)m序列的不可約多項(xiàng)式稱為本原多項(xiàng)式。 ④ 除了第r階以外, 如果還有偶數(shù)個(gè)抽頭的反饋結(jié)構(gòu), 則產(chǎn)生的序列就不是最長(zhǎng)線性移位寄存器序列。2．2 m序列的反饋系數(shù) 一個(gè)線性反饋移位寄存器能否產(chǎn)生m序列, 決定于它的電路反饋系數(shù)ci, 也就是它的遞歸關(guān)系式。 不同的反饋系數(shù), 產(chǎn)生不同的移位寄存器序列。 表21列出了不同級(jí)數(shù)的最長(zhǎng)線性移位寄存器序列的反饋系數(shù)。 r≥9時(shí), 由于m序列的條數(shù)很多, 不可能在此一一列出, 故只列出了一部分, 反饋系數(shù)Ci是以八進(jìn)制有示的。使用該表時(shí)，首先將每位八進(jìn)制數(shù)寫成二進(jìn)制形式。最左邊的1就是C0（C0恒為1），從此向右，依次用二進(jìn)制數(shù)表示C1，C2，……，C2，……Cn值后，就可構(gòu)成M序列發(fā)生器。例如，表中N=5，反饋系數(shù)Ci=（45）8，將它們化成二進(jìn)制數(shù)為100101，即相應(yīng)的反饋系數(shù)依次為C0=1，C1=0，C2=0，C3=1，C4=0，C5=1。表21 部分m序列反饋系統(tǒng)數(shù)表級(jí)數(shù)n周期P反饋系數(shù)（八進(jìn)制）37134152353145，67，75，663103，147，1557127203，211，217，235，277，313，325，345，3678255435，453，537，543，545，551，703，74795111021，1055，1131，1157，1167，11751010232011，2033，2157，2443，2745，34711120474005，4445，5023，5263，6211，736312409610123，11417，12515，13505，14127，1505313819120033，23261，24633，30741，32535，37505141638342103，51761，55753，60153，71147，674011532797100003，110013，120265，133663，1423051665535210013，233303，307572，311405，34743317131071400011，411335，444257，527427，646775表中的m序列的反饋系數(shù)只列出了一部分。 通過(guò)這些反饋系數(shù), 還可以求出對(duì)應(yīng)的鏡像序列的反饋抽頭和特征多項(xiàng)式。 所謂的鏡像序列是與原序列相反的序列。 如r＝3的序列為1110100, 鏡像序列為0010111。 可以通過(guò)下式, 由原序列的特征多項(xiàng)式f(x)求鏡像序列的特征多項(xiàng)式f(R)(x)， 即 （211）2．3 m序列發(fā)生器結(jié)構(gòu) m序列發(fā)生器的結(jié)構(gòu)一般有兩種形式, 簡(jiǎn)單型(SSRG)和模件抽頭型(MSRG)。 SSRG的結(jié)構(gòu)如圖23所示。 這種結(jié)構(gòu)的反饋邏輯由特征多項(xiàng)式確定, 這種結(jié)構(gòu)的缺點(diǎn)在于反饋支路中的器件時(shí)延是疊加的, 即等于反饋支路中所有模2加法器時(shí)延的總和。 因此限制了偽隨機(jī)序列的工作速度。 提高SSRG工作速率的辦法之一是選用抽頭數(shù)目少的m序列, 這樣, 還可簡(jiǎn)化序列產(chǎn)生器的結(jié)構(gòu)。 圖23 SSRG結(jié)構(gòu)發(fā)生器原理圖2．4 m序列的基本性質(zhì)（1）均衡性在m序列的一個(gè)周期內(nèi), “1”和“0”的數(shù)目基本相等。 準(zhǔn)確地說(shuō), “1”的個(gè)數(shù)比“0”的個(gè)數(shù)多一個(gè)。（2）游程分布把一個(gè)序列中取值相同的那些相繼元素合稱一個(gè)游程。 在一個(gè)游程中, 元素的個(gè)數(shù)稱為游程長(zhǎng)度。 （3）移位相加性 一個(gè)序列{an}與其經(jīng)m次遲延移位產(chǎn)生的另一不同序列{an+m}模2加, 得到的仍然是{an}的某次遲延移位序列{an+k}, 即 {an}+{an+m}={an+k} (213)（4）周期性 m序列的周期為N＝2r1, r為反饋移位寄存器的級(jí)數(shù)。 （5）偽隨機(jī)性 如果對(duì)一正態(tài)分布白噪聲取樣, 若取樣值為正, 記為“＋”。 若取樣值為負(fù), 記為“－”,則將每次取樣所得極性排成序列, 可以寫成 …+++++++++…這是一個(gè)隨機(jī)序列, 具有如下基本性質(zhì)： (1) 序列中“＋”和“－”的出現(xiàn)概率相等。 (2) 序列中長(zhǎng)度為1的游程約占1／2, 長(zhǎng)度為2的游程約占1／4, 長(zhǎng)度為3的游程約占1/8。 (3) 由于白噪聲的功率譜為常數(shù), 自相關(guān)函數(shù)為一沖激函數(shù)δ(τ)。 2．5 m序列的相關(guān)性信號(hào)的自相關(guān)函數(shù)和功率譜之間形成一傅里葉變換對(duì), 即 （214） 由于m序列的自相關(guān)函數(shù)是周期性的, 則對(duì)應(yīng)的頻譜是離散的。 自相關(guān)函數(shù)的波形是三角波, 對(duì)應(yīng)的離散譜的包絡(luò)為Sa2(x)。 由此可得m序列的功率譜G(ω)為 (215)圖25給出G(ω)的頻譜圖, Tc為偽碼chip的持續(xù)時(shí)間。 圖25 m序列的頻譜圖由此可得: (1) m序列的功率譜為離散譜, 譜線間隔ω1=2π/(NTc)。 (2) 功率譜的包絡(luò)為Sa2(Tcω/2N), 每個(gè)分量的功率與周期N成反比。 (3) 直流分量與N2成反比, N越大, 直流分量越小, 載漏越小。 (4) 帶寬由碼元寬度Tc決定, Tc越小, 即碼元速率越高, 帶寬越寬。 (5) 第一個(gè)零點(diǎn)出現(xiàn)在2π／Tc。 (6) 增加m序列的長(zhǎng)度N, 減小碼元寬度Tc, 將使譜線加密, 譜密度降低, 更接近于理想噪聲特性。第三章 詳細(xì)設(shè)計(jì)3．1 十階m序列的設(shè)計(jì)框圖十階m序列即n=10，本設(shè)計(jì)選取反饋系數(shù)為（2011）8 轉(zhuǎn)換為二進(jìn)制為10000001001。對(duì)應(yīng)多項(xiàng)式為 (31)則其反饋電路如（31）所示：C10=1C3=1C0=1D2D3D4D5D6D7D8D9D10輸出D1時(shí)鐘 圖31 十階m序列發(fā)生器原理圖圖中C0，C1，……C10均為反饋線（其中無(wú)反饋線，系數(shù)C為0，圖中省略），其中C0=Cn=1，表示反饋連接。因?yàn)閙序列是由循環(huán)序列發(fā)生器產(chǎn)生的，因此C0和C10肯定為1，即參于反饋。而反饋系數(shù)C3若為1，參于反饋；為0的，則表示斷開(kāi)反饋線，即開(kāi)路，無(wú)反饋連線。一個(gè)線性反饋移位寄存器能否產(chǎn)生m序列，決定于它的反饋系數(shù)Ci(C0,C1,……C10的總稱)。3．2工作流程圖為方便說(shuō)明，其工作原理如圖32流程所示：否是初始化輸入初始值Di設(shè)置P=0輸出值=D10Di=Di1周期P+1P是否等于1023結(jié)束D0=D3⊕D10 圖33 m序列生成器工作流程圖根據(jù)33所示流程圖，用MATLAB做出仿真程序?yàn)椋篺unction [seq]=gen_m_seq(connections)。 % connections:反饋系數(shù)% 寄存器初始化為[0......0 1]connections=randint(1,10)m=length(connections)。%碼長(zhǎng)L=2^m1。%寄存器初始化[0......0 1]registers=[zeros(1,m1) 1]。%序列初始化seq=zeros(1,L)。seq(1)=registers(m)。%找connections中不為零的系數(shù)index=find(connections==1)。for i=2:LM=0。for k=1:length(index)。M=xor(registers(index(k)),M)。end。%移位registers=[M,registers(1:m1)]。%輸出seq(i)=registers(m)。end自相關(guān)性與互相關(guān)性代碼如下：function [max_coor] = plot_pcf(a, b)% 計(jì)算等長(zhǎng)序列a,b的循環(huán)自相關(guān)和互相關(guān)并繪圖a=randint(1,10)b=randint(1,10)L = length(a)。K = floor(L/2)。a_span = [a(LK+1:L),a,a(1:K)]。b_span = [b(LK+1:L),b,b(1:K)]。for k = K : K % 周期自相關(guān)函數(shù)PACFPACFa(k+K+1) = 1/L * a(1:L)*a_span(k+K+1:k+K+L)39。PACFb(k+K+1) = 1/L * b(1:L)*b_span(k+K+1:k+K+L)39。endfigure。subplot(3,1,1)。plot([K : K], abs(PACFa), 39。k39。)。ylabel(39。a的自相關(guān)39。)。axis([K K ])。title(39。序列a,b的循環(huán)自、互相關(guān)39。)。subplot(3,1,2)。plot(