【正文】 expresscase test1 {mands} case {test2}{mands} ……otherwise {mands} End m序列仿真實(shí)現(xiàn)第三章中介紹了特征多項(xiàng)式，根據(jù)特征多項(xiàng)式就可以知道移位寄存器的結(jié)構(gòu)，從而產(chǎn)生m序列。本論文中m序列生成選擇的本原多項(xiàng)式為 （）則可以確定移位寄存器的結(jié)構(gòu)如下： 利用Matlab編寫(xiě)該移位寄存器的源程序，其流程圖如下： 寄存器初始化寄存器逐個(gè)移位采樣生成m序列相關(guān)特性仿真 軟件流程圖m序列仿真的源程序如下：% clear。 K=3。 %周期a=[1 0 1 0 0]。 % 初始狀態(tài)m=length(a)。 for k=1:K for i=1:2.^m1 %寄存移位 s((k1)*(2.^m1)+i)=a(5)。 temp=xor(a(3),a(5))。 a(5)=a(4)。 a(4)=a(3)。 a(3)=a(2)。 a(2)=a(1)。 a(1)=temp。 end end n=length(s)。 for i=1:n %采樣 for j=1:16 ss((i1)*16+j)=s(i)。 end end T=length(ss)。 %長(zhǎng)度t=[0:1:T1]。 plot(t,ss)。title(39。偽隨機(jī)序列39。)。axis([0 T1 ])。%畫(huà)圖 s1=s。 s1=s1*21。 L=length(s1)。 for delta=1:n s2=[zeros(1,delta) s1(1:Ldelta)]。 s2(1:delta)=s1(Ldelta+1:L)。 corr1(delta)=sum(s1.*s2)。 end figure plot(corr1)。title(39。m序列的自相關(guān)函數(shù)39。)K=3。 a=[1 0 1 0 0]。 m=length(a)。 for k=1:K for i=1:2.^m1 l((k1)*(2.^m1)+i)=a(5)。 temp1=xor(a(5),a(4))。temp2=xor(a(2),temp1)temp=xor(temp2,a(1))。 a(5)=a(4)。 a(4)=a(3)。 a(3)=a(2)。 a(2)=a(1)。 a(1)=temp。 end end n=length(s)。 for i=1:n for j=1:16 ll((i1)*16+j)=l(i)。 end end T=length(ll)。 t=[0:1:T1]。 figure(3)plot(t,ss) axis([0 T1 ]) l1=l。l1=2*l11。L=length(s1)。 for delta=1:n s2=[zeros(1,delta) s1(1:Ldelta)]。 s2(1:delta)=s1(Ldelta+1:L)。 corr2(delta)=sum(l1.*s2)。 end figure(3)。plot(1:length(corr2),corr2,39。*39。)。title(39。兩個(gè)隨機(jī)序列的互相關(guān)函數(shù)39。)。figure(4)。plot(1:length(corr1),corr1,39。r*39。,1:length(corr2),corr2,39。b39。)。title(39。m序列的性質(zhì)39。) 仿真生成m序列m序列生成一個(gè)周期如下：…0010110011111000110111010100001…。第三章介紹了m序列游程分布的隨機(jī)性，一般說(shuō)來(lái)，在q元n級(jí)m序列的一個(gè)周期p中，游程總數(shù)為。其中長(zhǎng)度為n的游程數(shù)目q1個(gè)，長(zhǎng)度為k的游程數(shù)目占總游程數(shù)目的（）。上面生成的m序列為一個(gè)q=2的5級(jí)移位寄存器生成周期為31的m序列。一個(gè)周期為31個(gè)碼元中，0的游程為8個(gè)，1的游程為8個(gè)，共有個(gè)游程，其中長(zhǎng)度為5的游程有一個(gè)，即“11111” ；長(zhǎng)度為4的游程有一個(gè)，即“0000” ；長(zhǎng)度為3的游程有兩個(gè)，即“000” 和“111” ；長(zhǎng)度為2的游程有四個(gè)，即“00”和“11”各兩個(gè)，長(zhǎng)度為1的游程有8個(gè)。同時(shí)在上面的序列中1的個(gè)數(shù)有16個(gè)比0的個(gè)數(shù)多一個(gè)，這也印證了m序列的均衡性。而該序列的移位序列與原序列模二加又是一個(gè)新的m序列，必定是該m序列的一個(gè)狀態(tài)，同時(shí)也滿(mǎn)足了m序列的移位相加性。由于m序列的自相關(guān)函數(shù)表示為： （）該m序列的自相關(guān)函數(shù)仿真如下，由于m序列具有周期性，其自相關(guān)函數(shù)也相應(yīng)地具有周期性，下面給出的是3個(gè)周期的碼長(zhǎng)為31的m序列的自相關(guān)函數(shù)曲線。 m序列的自相關(guān)函數(shù)不同的m序列之間也存在著不同的關(guān)系，把這種相互之間的關(guān)系一般用互相關(guān)來(lái)表示。由數(shù)字“1”和“0”構(gòu)成的二進(jìn)制序列的互相關(guān)函數(shù)就可以改寫(xiě)為 （）式中，A表示兩個(gè)序列在一個(gè)周期中對(duì)應(yīng)元素相同的個(gè)數(shù)，或兩序列模二加后0的個(gè)數(shù)；D表示兩序列在一個(gè)周期中對(duì)應(yīng)元素不同的個(gè)數(shù)，或兩序列模二加后為1的個(gè)數(shù)；p表示一個(gè)序列中元素的總數(shù)，即序列長(zhǎng)度，p=A+D。 計(jì)算時(shí)可先求兩序列的模二和序列，然后將和序列中0的個(gè)數(shù)減去1的個(gè)數(shù)得到差值，將該差值再除以序列長(zhǎng)度p就可以得到互相關(guān)函數(shù)。，其值具有三值性。 兩個(gè)m序列的互相關(guān)函數(shù) m序列的自相關(guān)和互相關(guān)函數(shù)生成m序列后可進(jìn)行直序擴(kuò)頻的仿真，將基帶信息序列加在m序列上即得到直序擴(kuò)頻信號(hào)。以下是直擴(kuò)序列實(shí)現(xiàn)的源程序和圖形。clear。 K=3。 a=[1 0 1 0 0]。 m=length(a)。 for k=1:K for i=1:2.^m1 s((k1)*(2.^m1)+i)=a(5)。 temp=xor(a(3),a(5))。 a(5)=a(4)。 a(4)=a(3)。 a(3)=a(2)。 a(2)=a(1)。 a(1)=temp。 end end n=length(s)。 for i=1:n for j=1:16 ss((i1)*16+j)=s(i)。 end end T=length(ss)。 t=[0:1:T1]。 plot(t,ss) axis([0 T1 ]) d=round(rand(1,16))。 %產(chǎn)生基帶信息序列signal=[]。 for o=1:16 if d(1,o)==0 sig=ones(1,T)。 else sig=ones(1,T)。 end signal=[signal sig]。 end figure(2)。plot(signal)。 axis([0 14000 ])。 title(39?；鶐盘?hào)39。)。ss=[ss ss ss ss ss ss ss ss ss ss]。 q=signal(1:length(ss)).*ss。 figure(3)plot(q)。axis([0 length(ss) ])。title(’直接序列擴(kuò)頻’) 基帶信號(hào)基帶信號(hào)即要傳送的信息序列，下面給出基帶信號(hào)疊加在m序列上的直擴(kuò)信號(hào)仿真圖形。第五章 結(jié)論 本論文主要針對(duì)直擴(kuò)信號(hào)發(fā)生器的設(shè)計(jì)，介紹了擴(kuò)頻通信的基本原理，在Matlab平臺(tái)上實(shí)現(xiàn)了m序列生成和直序擴(kuò)頻調(diào)制，并對(duì)m序列的自相關(guān)和互相關(guān)特性進(jìn)行了仿真分析。第一章對(duì)擴(kuò)頻技術(shù)的發(fā)展歷程做了一個(gè)概要的闡述，簡(jiǎn)單介紹了擴(kuò)頻通信的基本理論、特點(diǎn)和應(yīng)用。第二章主要介紹了擴(kuò)頻通信系統(tǒng)的良好性能。第三章主要介紹了與本次畢業(yè)設(shè)計(jì)的核心內(nèi)容有關(guān)的原理，如移位寄存器、m序列的生成方法、特性及其相關(guān)函數(shù)。第四章詳細(xì)說(shuō)明了在Matlab環(huán)境中對(duì)m序列和直序擴(kuò)頻信號(hào)的仿真及對(duì)m序列相關(guān)特性進(jìn)行的分析，給出了源程序和仿真結(jié)果。參考文獻(xiàn)1 [M].西安: 西安電子科技大學(xué)出版社，2007.2 李建新等．現(xiàn)代通信系統(tǒng)分析與仿真一MATLAB通信工具箱[M].西安：西安電子科技大學(xué)出版杜，2000.3 [M].合肥:中國(guó)科學(xué)技術(shù)大學(xué)出版社,1995.4 張國(guó)云，[J].現(xiàn)代電子技術(shù)，2007，12.5 [J].裝備制造技術(shù)，2008.6 ：北京郵電大學(xué)出版社，2002，10.7 張平等．第三代蜂窩移動(dòng)通信系統(tǒng)WCDMA．北京：北京郵電大學(xué)出版社，2000.8 維特比A J，李世鶴等譯．CDMA擴(kuò)頻通信原理．北京：人民郵電出版社，1997.9 桑林 郝建軍 劉丹譜. ：北京郵電大學(xué)出版社，2002.10 曾一凡 李暉. ：機(jī)械工業(yè)出版社，2005. 11 ：人民郵電出版社，2001.致 謝 三個(gè)月的畢業(yè)設(shè)計(jì)告一段落，由于直擴(kuò)信號(hào)和Matla b對(duì)我來(lái)說(shuō)是一個(gè)全新的概念，要在這三個(gè)月時(shí)間內(nèi)單靠我一個(gè)人完成畢業(yè)設(shè)計(jì)有點(diǎn)困難，在設(shè)計(jì)期間周?chē)耐瑢W(xué)老師給與了大大的支持和鼓勵(lì)，在此向他們表示感謝。首先，我要由衷的感謝鄒瓊導(dǎo)師，在畢業(yè)設(shè)計(jì)的整個(gè)過(guò)程中，鄒老師給我提供和介紹了很多很有價(jià)值的參考書(shū)籍和文獻(xiàn)，使我很快掌握了理論知識(shí)，在畢業(yè)設(shè)計(jì)過(guò)程中，對(duì)我的整個(gè)畢設(shè)做了宏觀的計(jì)劃和精心的指導(dǎo)，對(duì)我在課題設(shè)計(jì)中涉及到的一些沒(méi)弄明白的問(wèn)題耐心地講解，同時(shí)給我提出了很多寶貴的建議。在此我對(duì)鄒導(dǎo)師給予的支持和幫助表示衷心的感謝。其次，同是做Matlab仿真的廖峰同學(xué)在設(shè)計(jì)的過(guò)程當(dāng)中經(jīng)常和我交流討論，也幫助我解決了一些問(wèn)題，在此一并表示感謝。附 錄程序清單：m序列生成及其特性：clear。 K=3。 %周期a=[1 0 1 0 0]。 % 初始狀態(tài)m=length(a)。 for k=1:K for i=1:2.^m1 %寄存移位 s((k1)*(2.^m1)+i)=a(5)。 temp=xor(a(3),a(5))。 a(5)=a(4)。 a(4)=a(3)。 a(3)=a(2)。 a(2)=a(1)。 a(1)=temp。 end end n=length(s)。 for i=1:n %采樣 for j=1:16 ss((i1)*16+j)=s(i)。 end end T=length(ss)。 %長(zhǎng)度t=[0:1:T1]。 plot(t,ss)。title(39。偽隨機(jī)序列39。)。axis([0 T1 ])。%畫(huà)圖 s1=s。 s1=s1*21。 L=length(s1)。 for delta=1:n s2=[zeros(1,delta) s1(1:Ldelta)]。 s2(1:delta)=s1(Ldelta+1:L)。 corr1(delta)=sum(s1.*s2)。 end figure plot(corr1)。title(39。m序列的自相關(guān)函數(shù)39。)K=3。 a=[1 0 1 0 0]。 m=length(a)。 for k=1:K for i=1:2.^m1 l((k1)*(2.^m1)+i)=a(5)。 temp1=xor(a(5),a(4))。 temp2=xor(a(2),temp1) temp=xor(temp2,a(1))。 a(5)=a(4)。 a(4)=a(3)。 a(3)=a(2)。 a(2)=a(1)。 a(1)=temp。 end end n=length(l)。 for i=1:n for j=1:16 ll((i1)*16+j)=l(i)。 e