基于卷積編碼的擴頻通信系統(tǒng)軟件平臺設(shè)計課程設(shè)計

1、<p><b>　　摘 要 </b></p><p>　　MATLAB是一種用于科學(xué)計算的高性能語言。它可將計算、可視化和程序設(shè)計集成在一個非常容易使用的環(huán)境中。MATLAB是一個交互系統(tǒng)，它的基本數(shù)據(jù)元素是數(shù)組，尤其適合解決用矩陣和向量組織數(shù)據(jù)的科學(xué)技術(shù)計算問題。此次課設(shè)就是利用MATLAB軟件編寫程序?qū)崿F(xiàn)卷積編碼和Viterbi譯碼的算法，并且實現(xiàn)擴頻通信系統(tǒng)的擴頻與解擴，從

2、而得出信噪比與誤碼率之間的關(guān)系。</p><p><b>　　關(guān)鍵詞</b></p><p>　　卷積 擴頻 信噪比 誤碼率</p><p><b>　　前 言</b></p><p>　　擴頻通信是現(xiàn)代通信系統(tǒng)中新的通信方式，它具有較強的抗干擾、抗衰落和抗多徑性能，頻譜利用率高。擴頻信號是

3、不可預(yù)測的、偽隨機的寬帶信號，其帶寬遠大于要傳輸?shù)臄?shù)據(jù)(信息)帶寬，同時接收機中必須有與寬帶載波同步的副本。在發(fā)端輸人的信息先調(diào)制形成數(shù)字信號，然后由擴頻碼發(fā)生器產(chǎn)生的擴頻碼序列去調(diào)制數(shù)字信號以展寬信號的頻譜，展寬后的信號再調(diào)制到射頻發(fā)送出去。在接收端收到的寬帶射頻信號，變頻至中頻，然后由本地產(chǎn)生的與發(fā)端相同的擴頻碼序列去相關(guān)解擴，再經(jīng)信息解調(diào)，恢復(fù)成原始信息輸出。通信系統(tǒng)都要進行3次調(diào)制和相應(yīng)的解調(diào)。一次調(diào)制為信息調(diào)制。與一般通信系統(tǒng)

4、比較，多了擴頻調(diào)制和解擴部分。擴頻通信具備如下特征：(1)數(shù)字傳輸方式；(2)傳輸信號的帶寬遠大于被傳信息帶寬；(3)帶寬的展寬，是利用與被傳信息無關(guān)的函數(shù)(擴頻函數(shù))對被傳信息的信元重新進行調(diào)制實現(xiàn)的；(4)接收端用相同的擴頻函數(shù)進行相關(guān)解調(diào)(解擴)，求解出被傳信息的數(shù)據(jù)。用擴頻函數(shù)(也稱偽隨機碼)調(diào)制和對信號相關(guān)處理是擴頻通信有別于其他通信的兩大特點。 </p><p>　　設(shè)計題目：基于卷積編碼的

5、擴頻通信系統(tǒng)軟件平臺設(shè)計</p><p><b>　　一 設(shè)計目的</b></p><p>　　鍛煉學(xué)生綜合運用所學(xué)知識設(shè)計實際通信系統(tǒng)的能力，為學(xué)生走向工作崗位處理本專業(yè)方向的各種實際問題打好基礎(chǔ)，做到學(xué)以致用。具體到本次設(shè)計：</p><p>　　1、學(xué)習(xí)通信中的糾錯編碼技術(shù)及其主要應(yīng)用；</p><p>　　2、掌

6、握擴頻通信原理技術(shù)及其特點，了解其在通信系統(tǒng)中特別是第三代移動通信系統(tǒng)所采用的碼分多址（CDMA）技術(shù)中的應(yīng)用；</p><p>　　3、學(xué)會用MATLAB軟件方法設(shè)計一般通信系統(tǒng)；</p><p>　　4、學(xué)習(xí)CDMA手機的信息處理過程，掌握其基帶電路主要組成部分和工作原理及功能；模擬CDMA手機的信息處理過程。完成基帶電路主要組成部分的軟件平臺設(shè)計并評估系統(tǒng)的可靠性。</p>

7、;<p><b>　　二 設(shè)計任務(wù)</b></p><p>　　1、建立CDMA手機處理信息的卷積編碼擴頻通信系統(tǒng)</p><p>　　2、軟件實現(xiàn)卷積編碼和Viterbi譯碼的算法</p><p>　　3、軟件實現(xiàn)擴頻通信系統(tǒng)的擴頻與解擴</p><p>　　4、評估卷積編碼擴頻通信系統(tǒng)的可靠性</

8、p><p><b>　　三 設(shè)計內(nèi)容</b></p><p>　　模擬CDMA手機的信息處理過程，完成基帶電路主要組成部分的軟件平臺設(shè)計并評估系統(tǒng)的可靠性。</p><p>　　1、系統(tǒng)設(shè)計的總體框圖</p><p>　　圖1 卷積編碼擴頻通信系統(tǒng)軟件仿真平臺總體設(shè)計框圖</p><p>　　2.、

9、總體設(shè)計框圖中各模塊的實現(xiàn)</p><p>　　(1) 數(shù)據(jù)源：隨機產(chǎn)生{1，-1}的二進制數(shù)據(jù)</p><p><b>　　(2) 卷積編碼：</b></p><p>　?、賹?shù)據(jù)源數(shù)據(jù)經(jīng)卷積編碼輸出。</p><p>　　②選定卷積編碼器，如碼率、生成矩陣、約束長度等。</p><p>　　(

10、3) m序列發(fā)生器：產(chǎn)生擴頻中用到的m序列</p><p>　　(4) 擴頻與解擴：</p><p>　　將編碼輸出與擴頻碼相乘，即完成擴頻；將擴頻后的信號（疊加有噪聲）與擴頻碼相乘，即為解擴。</p><p>　　(5) Viterbi譯碼：為卷積碼的譯碼算法</p><p>　　(6) 可靠度評估：用誤碼率（BER）指標(biāo)，計算公式

11、為：</p><p>　　比較兩個數(shù)據(jù) </p><p><b>　　3，設(shè)計原理</b></p><p>　　3.1 卷積碼概述</p><p>　　數(shù)字通信系統(tǒng)進行數(shù)據(jù)傳輸時，不可避免地會在接收端產(chǎn)生差錯，為了降低誤碼率，通常采用前向糾錯編碼(分組碼或卷積碼)

12、的方法來糾正在傳輸過程中產(chǎn)生的誤碼。卷積編碼是現(xiàn)代通信系統(tǒng)中最為常用的信道編碼方式，它可以大大提高在加性高斯白噪聲情況下的抗噪聲能力，被廣泛用于衛(wèi)星通信和空間無線通信上。</p><p>　　3.1.1　卷積編碼的原理</p><p>　　卷積碼(又稱連環(huán)碼)，是由伊萊亞斯(P.Elis)提出的一種非分組碼。它把k比特信息段編成n比特的碼組，該碼組不僅同當(dāng)前的k比特信息段有關(guān)，而且還同前面

13、的(N-1)個信息段有關(guān)聯(lián)(N為大于1的整數(shù))。通常，把卷積碼記作(n,k,N)，其中k為輸入碼元數(shù)，n為輸出碼元數(shù)，N為約束長度，表示編碼器的存儲器級數(shù)。卷積編碼屬于信道編碼，主要用來糾正碼元的隨機差錯，它是以犧牲效率來換取可靠性，利用增加監(jiān)督位，進行檢錯和糾錯。</p><p>　　卷積碼編碼器是一個由k個輸入端、n個輸出端，且具有(N-1)或m節(jié)移位寄存器構(gòu)成的有限狀態(tài)記憶系統(tǒng)，通常稱為時序網(wǎng)絡(luò)，其原理如圖

14、2所示。 </p><p>　　圖2　卷積碼編碼器原理圖 </p><p>　　通常卷積碼的編碼器由K級（每級K比特）的移位寄存器和n個線性代數(shù)函數(shù)發(fā)生器（這里是模2加法器）組成。需要編碼的二進制數(shù)據(jù)串行輸入移位寄存器，每次移入K比特數(shù)據(jù)。每個K比特的輸入序列對應(yīng)一個n比特的輸出序列。因此卷積碼的編碼效率定義為。參數(shù)K被稱作卷積碼的約束長度，它表示當(dāng)前的n比特輸出序列與多少個K比特輸入序列

15、有關(guān)系，同時也是一個決定編碼復(fù)雜程度的重要參數(shù)。</p><p>　　3.1.2、描述卷積碼的方法</p><p>　　描述卷積編碼的方法之一是給出其生成矩陣。一般說來卷積碼的生成矩陣是一個半無限矩陣，因為輸入序列是半無限長度的。</p><p>　　另一種描述卷積碼的方法是用一組n個矢量--稱為卷積碼的生成多項式來表示。對應(yīng)于n個代數(shù)函數(shù)發(fā)生器（這里是模2加法器）

16、與移位寄存器的連接方式，我們得到n個生成多項式，這里是n個長度為 的矢量。某個矢量的第1個元素為1，表示寄存器相應(yīng)的位置與該矢量對應(yīng)的模2加法器相連；反之，如果該元素為0，則表示寄存器相應(yīng)的位置與該矢量對應(yīng)的模2加法器不連接。</p><p>　　3.1.3　CDMA手機卷積碼編碼器 CDMA系統(tǒng)是在擴頻通信基礎(chǔ)上發(fā)展起來的一種新的無線通信系統(tǒng)，其反向信道(又稱上行信道)是由手機發(fā)射，基站接收。為了提高

17、反向信道信號的抗干擾能力，采用碼率為1/3的卷積編碼器，即輸入一個碼元,編碼器相應(yīng)輸出3個碼元。約束長度為9，即前后9個碼元有關(guān)聯(lián)。因此,它包含8級移位寄存器和3個模2加法器,電路組成如圖3所示。 </p><p>　　圖3　CDMA手機卷積編碼器</p><p>　　由圖3可知，該編碼器是一個(3,1,9)卷積編碼器,即k=1(1個輸入端)、n=3(3個輸出端)、N=9(8級移位寄存器)

18、。 若輸入信息序列為:U=(u0u1u2…),則對應(yīng)輸出為3個碼字序列：</p><p>　　C0=(c(0)0c(0)1c(0)2…),　　C1=(c(1)0c(1)1c(1)2…),　　C2=(c(2)0c(2)1c(2)2…). 其相應(yīng)編碼方程可寫為：　C0=U*G(0),　　C1=U*G(1),　　C2=U*G(2). </p><p>　　式中“*”表示卷積運算,G(

19、0),G(1)和G(2)表示編碼器的3個沖激響應(yīng)。編碼輸出可由輸入信息序列U和編碼器的3個脈沖沖激響應(yīng)的卷積得到,故稱卷積碼。由于編碼器有8級寄存器,所以沖激響應(yīng)至多可持續(xù)到9位,圖2所示卷積碼編碼器的3個沖激響應(yīng)可寫成：</p><p>　　G(0)=(101101111),　　　G(1)=(110110011),　　　G(2)=(111001001)1 若輸入信息序列為:U=(10001)，則</p&

20、gt;<p>　　C0=(10001)*(101101111)=(1011110011111), C1=(10001)*(110110011)=(1101010000011), C2=(10001)*(111001001)=(1110101011001).</p><p>　　經(jīng)過并串轉(zhuǎn)換,最后輸出的碼字為:C=(111011101110101110001000101101100110111)。&l

21、t;/p><p>　　3.1.4.維特比譯碼</p><p>　　(1)Viterbi譯碼原理簡述</p><p>　　卷積碼是1955年埃里亞斯(Elias)最早提出，1957年伍成克拉夫(Wozencraft)提出了一種有效譯碼方法，即序列譯碼。1963年梅西(Massey)提出了一種性能稍差，但比較實用的門限譯碼方法，由于這一實用性進展使卷積碼從理論走向?qū)嵱没６?/p>

22、后1967年維特比(Viterbi)提出了最大似然譯碼法。它對存儲器級數(shù)較小的卷積碼的譯碼很容易實現(xiàn)，人們后來稱它為維特比算法或維特比譯碼，并被廣泛地應(yīng)用于現(xiàn)代通信中。</p><p>　　卷積碼的 Viterbi 譯碼是根據(jù)接收碼字序列尋找編碼時通過網(wǎng)格圖最佳路徑的過程 ，找到最佳路徑即完成了譯碼過程 ，并可以糾正接收碼字中的錯誤比特。所謂“最佳”，是指最大后驗條件概率：</p><p>

23、;　　P( C/ R) = maxj[ P ( Cj/ R) ] ，</p><p>　　一般來說，信道模型并不使用后驗條件概率，因此利用Beyes公式，根據(jù)信道特性得出結(jié)論maxj[ P ( Cj/ R) ]與maxj[ P ( R/ Cj) ]等價。考慮到在系統(tǒng)實現(xiàn)中往往采用對數(shù)形式的運算，以求降低運算量，并且為求運算值為整數(shù)加入了修正因子 a1、 a2。經(jīng)過一系列轉(zhuǎn)化之后可以得出得出，尋找最佳路徑，就變成尋

24、找最大分支路徑量度。</p><p>　　(2)Viterbi譯碼流程</p><p>　?、俑鶕?jù)接收碼符號計算出相應(yīng)的分支量度值；</p><p>　?、谶M入某一狀態(tài)的分支量度與其前狀態(tài)路徑量度累加求和；</p><p>　?、郾容^到達當(dāng)前狀態(tài)的新的路徑量度的大小，選擇最大者作為新的狀態(tài)路徑量度存儲起來，并保存與此路徑對應(yīng)的碼字；</

25、p><p>　?、軐λ袪顟B(tài)都實施上述加、比、選運算；</p><p>　?、菰诿恳蛔g碼時刻，滿足延時就從 256 條留存路徑中，選擇路徑量度最大的一條路徑作為譯碼數(shù)輸出；</p><p>　　⑥進入下一譯碼時刻，重復(fù)以上步驟，直至譯碼結(jié)束。</p><p>　　(3) 對譯碼算法及實現(xiàn)的優(yōu)化</p><p>　　由于卷積

26、碼譯碼的復(fù)雜度隨著約束長度的增加以非線性方式迅速增加，在實際應(yīng)用中卷積碼的實際應(yīng)用性能往往受限于存儲器容量和系統(tǒng)運算速度，尤其是對約束長度比較大的卷積碼。為了在有限的硬件或軟件資源條件下保證系統(tǒng)較高的譯碼性能，需要對算法進行留存路徑更新算法優(yōu)化。</p><p><b>　　3.2 擴頻通信</b></p><p>　　3.2.1 擴頻通信的理論基礎(chǔ)</p&g

27、t;<p>　　擴頻通信的基本特點是傳輸信號所占用的頻帶寬度(W)遠大于原始信息本身實際所需的最小(有效)帶寬(DF)，其比值稱為處理增益Gp，其中 Gp = W/DF。</p><p>　　眾所周知，任何信息的有效傳輸都需要一定的頻率寬度，為了充分利用有限的頻率資源，增加通路數(shù)目，人們廣泛選擇不同調(diào)制方式，采用寬頻信道(同軸電纜、微波和光纖等)和壓縮頻帶等措施，同時力求使傳輸?shù)拿浇橹袀鬏數(shù)男盘栒加?/p>

28、盡量窄的帶寬。因現(xiàn)今使用的電話、廣播系統(tǒng)中，無論是采用調(diào)幅、調(diào)頻或脈沖編碼調(diào)制制式，Gp值一般都在十多倍范圍內(nèi)，統(tǒng)稱為“窄帶通信”，而擴頻通信的Gp值，高達數(shù)百、上千，稱為  “寬帶通信”。 而考慮到擴頻通信的可行性，是從信息論和抗干擾理論的基本公式中引伸而來的。 </p><p>　　信息論中關(guān)于信息容量的仙農(nóng)(Shannon)公式為： </p><p>　　C ＝ WLog2

29、(1十P/N) ...... (1) </p><p>　　式中：C指信道容量(用傳輸速率度量) ，W指信號頻帶寬度 ，P指信號功率 ，N指白噪聲功率。</p><p>　　公式(1)說明，在給定的傳輸速率C不變的條件下，頻帶寬度W和信噪比P／N是可以互換的。即可通過增加頻帶寬度的方法，在較低的信噪比P／N(S／N)情況下，傳輸信息.擴展頻譜換取信噪比要求的降低，正是擴頻通信的重

30、要特點，并由此為擴頻通信的應(yīng)用奠定了基礎(chǔ)。 </p><p>　　擴頻通信可行性的另一理論基礎(chǔ)，為柯捷爾尼可夫關(guān)于信息傳輸差錯概率的公式。公式說明了對于一定帶寬 DF的信息而言，用Gp值較大的寬帶信號來傳輸，可以提高通信抗干擾能力，保證強干擾條件下，通信的安全可靠。 </p><p>　　總之，我們用信息帶寬的100倍，甚至1000倍以上的寬帶信號來傳輸信息，就是為了提高通信的抗干擾能力，

31、即在強干擾條件下保證可靠安全地通信。這就是擴展頻譜通信的基本思想和理論依據(jù)。 </p><p><b>　　擴頻通信主要特點</b></p><p>　　①易于重復(fù)使用頻率，提高了無線頻譜利用率 </p><p>　　無線頻譜十分寶貴，雖然從長波到微波都得到了開發(fā)利用，仍然滿足不了社會的需求。在窄帶通信中，主要依靠波道劃分來防止信道之間發(fā)生干擾

32、。為此，世界各國都設(shè)立了頻率管理機構(gòu)，用戶只能使用申請獲準(zhǔn)的頻率。 </p><p>　　擴頻通信發(fā)送功率極低(1 --- 650mW)，采用了相關(guān)接收這一高技術(shù)，且可工作在信道噪聲和熱噪聲背景中，易于在同一地區(qū)重復(fù)使用同一頻率， 也可與現(xiàn)今各種窄道通信共享同一頻率資源。所以，在美國及世界絕大多數(shù)國家，擴頻通信不需申請頻率，任何個人與單位可以無執(zhí)照使用。 </p><p>　?、诳垢蓴_性強

33、，誤碼率低 </p><p>　　擴頻通信在空間傳輸時所占有的帶寬相對較寬，而收端又采用相關(guān)檢測的辦法來解擴，使有用寬帶信息信號恢復(fù)成窄帶信號，而把非所需信號擴展成寬帶信號，然后通過窄帶濾波技術(shù)提取有用的信號。這祥，對于各種干擾信號，因其在收端的非相關(guān)性，解擴后窄帶信號中只有很微弱的成份，信噪比很高，因此抗干擾性強。 </p><p>　　③隱蔽性好，對各種窄帶通信系統(tǒng)的干擾很小 <

34、/p><p>　　由于擴頻信號在相對較寬的頻帶上被擴展了，單位頻帶內(nèi)的功率很小，信號湮沒在噪聲里，一般不容易被發(fā)現(xiàn)，而想進一步檢測信號的參數(shù)(如偽隨機編碼序列)就更加困難，因此說其隱蔽性好。 再者，由于擴頻信號具有很低的功率譜密度，它對目前使用的各種窄帶通信系統(tǒng)的干擾很小。 </p><p><b>　?、芤子趯崿F(xiàn)碼分多址</b></p><p>

35、　　擴頻通信提高了抗干擾性能，但付出了占用頻帶寬的代價。 </p><p>　　如果讓許多用戶共用這一寬頻帶，則可大為提高頻帶的利用率。由于在擴頻通信中存在擴頻碼序列的擴頻調(diào)制，充分利用各種不同碼型的擴頻碼序列之間優(yōu)良的自相關(guān)特性和互相關(guān)特性，在接收端利用相關(guān)檢測技術(shù)進行解擴，則在分配給不同用戶碼型的情況下可以區(qū)分不同用戶的信號，提取出有用信號。這樣一來，在一寬頻帶上許多對用戶可以同時通話而互不干擾。 </

36、p><p>　　直擴通信占用寬帶頻譜資源通信，改善了抗干擾能力，是否浪費了頻段？其實正相反，擴頻通信提高了頻帶的利用率。正是由于直擴通信要用擴頻編碼進行擴頻調(diào)制發(fā)送，而信號接收需要用相同的擴頻編碼作相關(guān)解擴才能得到，這就給頻率復(fù)用和多址通信提供了基礎(chǔ)。充分利用不同碼型的擴頻編碼之間的相關(guān)特性，分配給不同用戶不同的擴頻編碼，就可以區(qū)別不同的用戶的信號，眾多用戶，只要配對使用自己的擴頻編碼，就可以互不干擾地同時使用同一頻

37、率通信，從而實現(xiàn)了頻率復(fù)用，使擁擠的頻譜得到充分利用。發(fā)送者可用不同的擴頻編碼，分別向不同的接收者發(fā)送數(shù)據(jù)；同樣，接收者用不同的擴頻編碼，就可以收到不同的發(fā)送者送來的數(shù)據(jù)，實現(xiàn)了多址通信。美國國家航天管理局（NASA）的技術(shù)報告指出：采用擴頻通信提高了頻譜利用率。另外，擴頻碼分多址還易于解決隨時增加新用戶的問題。 </p><p><b>　?、菘苟鄰礁蓴_ </b></p>&

38、lt;p>　　在無線通信的各個頻段，長期以來，多徑干擾始終是一個難以解決的問題之一。在以往的窄帶通信中，采用兩種方法來提高抗多徑干擾的能力： </p><p>　　一是把最強的有用信號分離出來，排除其他路徑的干擾信號，即采用分集／接收技術(shù)；二是設(shè)法把不同路徑來的不同延遲、不同相位的倍號在接收端從時域上對齊相加，合并成較強的有用信號，即采用梳狀濾波器的方法。 </p><p>　　這兩

39、種技術(shù)在擴頻通信中都易于實現(xiàn)。利用擴頻碼的自相關(guān)特性，在接收端可以用相關(guān)技術(shù) 從多徑信號中提取和分離出最強的有用信號，或把多個路徑來的同一碼序列的波形相加合成使之得到加強，這相當(dāng)于梳狀濾波器的作用。另外，采用頻率跳變擴頻調(diào)制方式的擴頻系統(tǒng)中，由于用多個頻率的信號傳送同一個信息，實際上起到了頻率分集的作用。 </p><p>　?、弈芫_地定時和測距 </p><p>　?、哌m合數(shù)字話音和數(shù)

40、據(jù)傳輸，以及開展多種通信業(yè)務(wù) </p><p>　　擴頻通信一般都采用數(shù)字通信、碼分多址技術(shù)，適用于計算機網(wǎng)絡(luò)，適合于數(shù)據(jù)和圖象傳輸。 </p><p>　　⑧安裝簡便，易于維護 </p><p>　　3.2.2 擴頻通信的實現(xiàn)方法</p><p>　　一個典型通信系統(tǒng)框圖如圖4所示：</p><p>　　圖4

41、典型通信系統(tǒng)框圖</p><p>　　圖5 擴頻通信系統(tǒng)一般原理框圖</p><p>　　擴頻通信系統(tǒng)一般原理框圖如圖5所示，在發(fā)端輸入的信息先經(jīng)信息調(diào)制形成數(shù)字信號，然后由擴頻碼發(fā)生器產(chǎn)生的擴頻碼序列去調(diào)制數(shù)字信號以展寬信號的頻譜。展寬后的信號再調(diào)制到射頻發(fā)送出去。 </p><p>　　一、擴頻解擴工作過程</p><p>　　在接收端

42、收到的寬帶射頻信號，變頻至中頻，然后由本地產(chǎn)生的與發(fā)端相同的擴頻碼序列去相關(guān)解擴。再經(jīng)信息解調(diào)、恢復(fù)成原始信息輸出。</p><p>　　由此可見，—般的擴頻通信系統(tǒng)都要進行三次調(diào)制和相應(yīng)的解調(diào)。一次調(diào)制為信息調(diào)制，二次調(diào)制為擴頻調(diào)制，三次調(diào)制為射頻調(diào)制，以及相應(yīng)的信息解調(diào)、解擴和射頻解調(diào)。 </p><p>　　與一般通信系統(tǒng)比較，擴頻通信就是多了擴頻調(diào)制和解擴部分。 </p>

43、;<p>　　二、擴頻通信的幾種工作方式 </p><p>　　按照擴展頻譜的方式不同，現(xiàn)有的擴頻通信系統(tǒng)可為</p><p>　?、僦苯有蛄袛U頻(Direct Sequence Spread Spectrum) </p><p>　　所謂直接序列(DS-Direct Scquency)擴頻，就是直接用具有高碼率的擴頻碼序列在發(fā)端去擴展信號的頻譜。而在

44、收端，用相同的擴頻碼序列去進行解擴，把展寬的擴頻信號還原成原始的信息。</p><p>　?、谔冾l率(Frequency Hopping)工作方式，簡稱跳頻(FH)方式 </p><p>　　另外一種擴展信號頻譜的方式稱為跳頻(FH－Frequency Hopping)。所謂跳頻，比較確切的意思是：用一定碼序列進行選擇的多頻率頻移鍵控。也就是說，用擴頻碼序列去進行頻移鍵控調(diào)制，使載波頻率

45、不斷地跳變，所以稱為跳頻。 </p><p>　?、厶儠r間(Time Hopping)工作方式，簡稱跳時(TH)方式 </p><p>　　與跳頻相似，跳時(TH－Time Hopping)是使發(fā)射信號在時間軸上跳變。首先把時間軸分成許多時片。在一幀內(nèi)哪個時片發(fā)射信號由擴頻碼序列去進行控制?？梢园烟鴷r理解為：用一定碼序列進行選擇的多時片的時移鍵控。 </p><p&g

46、t;　　由于采用了窄得很多的時片去發(fā)送信號，相對說來，信號的頻譜也就展寬了。在發(fā)端，輸入的數(shù)據(jù)先存儲起來，由擴頻碼發(fā)生器的擴頻碼序列去控制通－斷開關(guān)，經(jīng)二相或四相調(diào)制后再經(jīng)射頻調(diào)制后發(fā)射。在收端，由射頻接收機輸出的中頻信號經(jīng)本地產(chǎn)生的與發(fā)端相同的擴頻碼序列控制通－斷開關(guān)，再經(jīng)二相或四相解調(diào)器，送到數(shù)據(jù)存儲器和再定時后輸出數(shù)據(jù)。只要收發(fā)兩端在時間上嚴(yán)格同步進行，就能正確地恢復(fù)原始數(shù)據(jù)。 </p><p>　?、軐拵?/p>

47、線性調(diào)頻(Chirp Modulation)工作方式，簡稱Chirp方式 </p><p><b>　?、莞鞣N混合方式 </b></p><p>　　在上述幾種基本的擴頻方式的基礎(chǔ)上，可以組合起來，構(gòu)成各種混合方式。</p><p>　　本次設(shè)計采用直接序列擴頻。</p><p>　　所謂直接序列(DS)擴頻，就是直接用

48、具有高碼率的擴頻碼序列在發(fā)端去擴展信號的頻譜。而在收端，用相同的擴頻碼序列去進行解擴，把展寬的擴頻信號還原成原始的信息。</p><p>　　信號在射頻信道傳輸過程中必然受到各種外來信號的干擾。因此，在收端，進入接收機的除有用信號外還存在干擾信號。假定干擾為功率較強的窄帶信號，寬帶有用信號與干擾信號同時經(jīng)變頻至中心頻率為中頻fI輸出。不言而喻，對這一中頻寬帶信號必須進行解擴處理才能進行信息解調(diào)。 解擴實際上就是擴

49、頻的反變換，通常也是用與發(fā)端相同的調(diào)制器，并用與發(fā)端完全相同的偽隨機碼序列對收到的寬帶信號再一次進行二相相移鍵控。 </p><p>　　再一次的相移鍵控正好把擴頻信號恢復(fù)成相移鍵控前的原始信號。從頻譜上看則表現(xiàn)為寬帶信號被解擴壓縮還原成窄帶信號。這一窄帶信號經(jīng)中頻窄帶濾波器后至信息解調(diào)器再恢復(fù)成原始信息。但是對于進入接收機的變窄帶干擾信號，在收端調(diào)制器中同樣也受到偽隨機碼的雙相相移鍵控調(diào)制，它反而使窄帶干擾變成

50、寬度干擾信號。由于干擾信號頻譜的擴展，經(jīng)過中頻窄帶通濾波作用，只允許通帶內(nèi)的干擾通過，使干擾功率大為減少。由此可見，接收機輸入端的信號與噪聲經(jīng)過解擴處理，使信號功率集中起來通過濾波器，同時使干擾功率擴散后被濾波器大量濾除，結(jié)果便大大提高了輸出端的信號噪聲功率比。擴頻解擴過程可示意如圖7所示。</p><p>　　圖7擴頻解擴過程示意圖</p><p>　　這一過程說明了直擴系統(tǒng)的基本原理和

51、它是怎樣通過對信號進行擴頻與解擴處理從而獲得提高輸出信噪比的好處的。它體現(xiàn)了直擴系統(tǒng)的抗干擾能力。</p><p>　　綜上所述，直擴系統(tǒng)的特點是： </p><p>　　a．頻譜的擴展是直接由高碼率的擴頻碼序列進行調(diào)制而得到的。 </p><p>　　b.擴頻碼序列多采用偽隨機碼，也稱為偽噪聲(PN)碼序列。 </p><p>　　c.擴頻

52、調(diào)制方式多采用BPSK或QPSK等幅調(diào)制。擴頻和解擴的調(diào)制解調(diào)器多采用平衡調(diào)制器，制作簡單又能抑制載被。</p><p>　　d.模擬信息調(diào)制多采用頻率調(diào)制(FM)，而數(shù)字信息調(diào)制多采用脈沖編碼調(diào)制(PCM)或增量調(diào)制(DM)。 </p><p>　　e.接收端多采用產(chǎn)生本地偽隨機碼序列對接收信號進行相關(guān)解擴，或采用匹配濾波器來解擴信號。 </p><p>　　f.

53、擴頻和解擴的偽隨機碼序列應(yīng)有嚴(yán)格的同步，碼的搜捕和跟蹤多采用匹配濾波器或利用偽隨機碼的優(yōu)良的相關(guān)特性在延遲鎖定環(huán)中實現(xiàn)。 </p><p>　　g.一般需要用窄帶通濾波器來排除干擾，以實現(xiàn)其抗干擾能力的提高。 </p><p><b>　　3.2.3.解擴 </b></p><p><b>　　(1)概述</b></

54、p><p>　　正如在一般的窄帶通信中，已調(diào)信號在接收端都要進行解調(diào)來恢復(fù)所傳的信息。在擴頻通信中接收端則用與發(fā)送端相同的擴頻碼序列與收到的擴頻信號進行相關(guān)解調(diào)，恢復(fù)所傳的信息。換句話說，這種相關(guān)解調(diào)起到解擴的作用。即把擴展以后的信號又恢復(fù)成原來所傳的信息。這種在發(fā)端把窄帶信息擴展成寬帶信號，而在收端又將其解擴成窄帶信息的處理過程，會帶來一系列好處。弄清楚擴頻和解擴處理過程的機制，是理解擴頻通信本質(zhì)的關(guān)鍵所在。<

55、;/p><p><b>　　(2)原理分析</b></p><p>　　在現(xiàn)實應(yīng)用中，一般采用相關(guān)檢測或匹配濾波的方法來解擴。 </p><p>　　所謂相關(guān)檢測，是指當(dāng)你想檢測出所需要的有用信號，有效的方法是在本地產(chǎn)生一個相同的信號，然后用它與接收到的信號對比，求其相似性。換句話說，就是用本地產(chǎn)生的相同的信號與接收到的信號進行相關(guān)運算，其中相關(guān)函

56、數(shù)最大的就最可能是所要的有用信號。 </p><p>　　基本的擴頻過程就是在收端產(chǎn)生與發(fā)端完全相同的PN碼，對收到的擴頻信號，在平衡調(diào)制器中再一次進行二相相移鍵控調(diào)制?？梢钥闯霭l(fā)端相移鍵控調(diào)制后的信號在收端又被恢復(fù)成原來的載波信號。當(dāng)然一個必要的條件是本地的PN碼信號的相位必須和收到的相移后的信號在相移點對準(zhǔn)，才能正確地將相移后的信號再翻轉(zhuǎn)過來。由此可見，收發(fā)兩端信號的同步十分重要。 </p>&

57、lt;p><b>　　(3)解擴過程</b></p><p>　　平衡調(diào)制器把收到的展寬的信號解擴成信息調(diào)制的載波。最后經(jīng)帶通濾波器輸出。以上所述就是所謂的相關(guān)解擴過程。通常為了處理方便，大多在中頻進行。也就是接收到的擴頻信號，先在變頻器中先變換到中頻，再進入到平衡調(diào)制器中解擴。其后接中頻帶通濾波器輸出。有時為了避免強干擾信號從平衡調(diào)制器的輸入端繞過它而泄漏到輸出端去，可以來用外差相關(guān)

58、解擴。</p><p>　　本地產(chǎn)生的PN碼先與本地振蕩器產(chǎn)生的與接收信號差一個中頻信號的本地振蕩信號在下面一個平衡調(diào)制器進行調(diào)制，產(chǎn)生本地參考信號。它是一個展寬了的信號。然后，此本地參考信號與接收的信號在上面一個平衡調(diào)制器調(diào)制成中頻輸出信號。這時平衡調(diào)制器實際上起的是混頻器的作用。由于它的輸入信號與輸出信號不同，也就不會發(fā)生強干擾信號直接繞過去的泄漏了。并且后面還有一個中頻帶通濾波器，可以起到濾除干擾的作用。&

59、lt;/p><p>　　(4)效果分析及實現(xiàn)方法</p><p>　　相關(guān)解擴過程對擴頻通信至關(guān)重要。正是這一解擴過程大大提高了系統(tǒng)的抗擾能力。 </p><p>　　相關(guān)解擴在性能上固然很好，但總是需要在接收端產(chǎn)生本地PN碼。這一點有時帶來許多不方便。例如，解決本地信號與接收信號的同步問題就很麻煩，還不能做到實時把有用信號檢測出來。因為匹配濾波和相關(guān)檢測的作用在本質(zhì)上

60、是一樣的，我們可以用匹配濾波器來解擴直擴信號。 </p><p>　　所謂匹配濾波器，就是與信號相匹配的濾波器，它能在多種信號或干擾中把與之匹配的信號檢測出來。對于視頻矩形脈沖序列來說，無源匹配濾波器就是抽頭延遲線上加上加法累加器。有時稱為橫向濾波器。</p><p><b>　　4，心得體會</b></p><p>　　不知不覺已經(jīng)到了大四期

61、末，按照教學(xué)安排，我們進行了最后一次課程設(shè)計，此次課程設(shè)計內(nèi)容是用MATLAB軟件編寫程序，并實現(xiàn)卷積編碼和Viterbi譯碼的算法，并且實現(xiàn)擴頻通信系統(tǒng)的擴頻與解擴，從而得出信噪比與誤碼率之間的關(guān)系。在大一時候接觸過MATLAB，但是當(dāng)時只是知道該軟件是一款相當(dāng)強大的數(shù)學(xué)語言，可以實現(xiàn)復(fù)雜的數(shù)學(xué)運算，經(jīng)過為期三周的課設(shè)讓我們對MATLAB有了進一步的了解和認識。</p><p>　　在第一周，*老師先講解了課設(shè)

62、的相關(guān)背景知識。在接下來的時間里，我們就開始自己操作軟件，編寫程序。說實話專業(yè)基礎(chǔ)知識不是很扎實，而且之前沒有接觸過MATLAB語言，要用一門新語言編寫程序并實現(xiàn)仿真很困難。于是在老師的推薦下去圖書館借閱了參考書，請教老師，編寫模塊…….最終課設(shè)有了一些進展。</p><p>　　在編寫過程中，出現(xiàn)了一系列錯誤，請教同學(xué)，咨詢*老師，一個一個解決，功夫不負有心人，在程序成功運行出結(jié)果看到對比曲線心中十分高興。在運

63、行程序---提示錯誤---解決錯誤—總結(jié)原因的過程中，再次意識到做工程，各個方面都要考慮到，尤其是那種細節(jié)的錯誤，更要得到重視。最后幾天我們更是弄懂了每一個程序，將每個程序進行仿真，分析誤碼率與信噪比之間的關(guān)系，同時將每個仿真結(jié)果進行對比，分析性能，得出結(jié)論，最后順利的通過了驗收，完成了大學(xué)的最后一此課程設(shè)計。</p><p><b>　　5，參考文獻</b></p><

64、p>　　【1】興雯等.擴展頻譜通信及其多址技術(shù)[M].西安:西安電子科技大學(xué)出版社, 2004.</p><p>　　【2】蓉,李署堅. 擴頻通信技術(shù)[J]. 無線電工程, 2001,(S1) . </p><p>　　【3】淑君,柳鐸. 淺談擴頻通信技術(shù)及其應(yīng)用[J]. 山東電子, 2004,(01) . </p><p>　　【4】孟令彪. 無線

65、擴頻通信技術(shù)述略[J]. 商場現(xiàn)代化, 2004,(15) . 1998,(06) . </p><p>　　【6】徐明元，邵玉斌.matlab仿真在通信與電子工程中的應(yīng)用[M].西安電子科技大學(xué)出版社，2001</p><p>　　【7】郭海燕，畢紅軍.MATLAB在偽隨機碼的生成及仿真中的應(yīng)用[J].計算機仿真，21(3),2004.3</p><p><

66、b>　　6，附錄</b></p><p>　　程序一：本程序未加卷積碼，未加擴頻（文件名DS0）</p><p>　　EN = input(' 請輸入信噪比（用分貝表示）: ');</p><p>　　L= input(' 請輸入數(shù)據(jù)幀長度(= 信息比特數(shù) + 尾比特數(shù), 默認值: 400) : ');</

67、p><p>　　if isempty(L)</p><p><b>　　L = 400;</b></p><p><b>　　end</b></p><p>　　ferrlim=input('請輸入程序終止的幀錯誤數(shù),默認值: 15 : ');</p><p>

68、　　if isempty(ferrlim)</p><p>　　ferrlim = 15;</p><p><b>　　end </b></p><p>　　c=[-1 -1 -1 1 1 1 1 -1 1 -1 1 1 -1 -1 1];</p><p>　　N=length(c);</p><p&

69、gt;　　g=[1 0 1;1 1 1];k0=1;</p><p>　　for j=1:length(EN);</p><p>　　errs(j) = 0;</p><p>　　nferr(j) = 0;</p><p>　　nframe = 0;</p><p>　　while nferr(j)<ferrli

70、m</p><p>　　nframe = nframe + 1;</p><p>　　tem=rand(1,L);</p><p><b>　　for i=1:L</b></p><p>　　if tem(i)<=0.5</p><p><b>　　d(i)=0;</b>

71、;</p><p><b>　　else</b></p><p><b>　　d(i)=1;</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p>　　n=(1/10^(EN

72、(j)/10))*randn(1,L);</p><p><b>　　b=2*d-1;</b></p><p><b>　　for i=1:L</b></p><p>　　y(i)=b(i)+n(i);</p><p>　　if y(i)>=0</p><p><

73、b>　　Db(i)=1;</b></p><p><b>　　else</b></p><p><b>　　Db(i)=0;</b></p><p><b>　　end</b></p><p><b>　　end</b></p>

74、;<p>　　err=length(find(Db~=d));</p><p><b>　　if err>0</b></p><p>　　nferr(j)=nferr(j)+1;</p><p><b>　　end</b></p><p>　　errs(j) = errs(j)

75、+ err;</p><p>　　if rem(nframe,3)==0 | nferr(j)==ferrlim </p><p>　　% Bit error rate</p><p>　　ber(j) = errs(j)/nframe/(L);</p><p>　　% Frame error rate</p><p&g

76、t;　　fer(j) = nferr(j)/nframe;</p><p>　　% Display intermediate results in process </p><p>　　fprintf('************** Eb/N0 = % 5.2f db **************\n', EN(j));</p><p>　　fpr

77、intf('Frame size = %d, rate 1/%d. \n', L, 2);</p><p>　　fprintf('%d frames transmitted, %d frames in error.\n', nframe, nferr(j));</p><p>　　fprintf('%8.4e ', ber(j));

78、 </p><p>　　fprintf('\n'); </p><p>　　fprintf('%8.4e ', fer(j));</p><p>　　fprintf('\n');</p><p>　　fprintf('***************************

79、********************\n\n');</p><p><b>　　end</b></p><p>　　% Save intermediate results </p><p>　　end%while</p><p>　　end %nEN</p><p>　　sem

80、ilogy(EN,ber,'.');</p><p><b>　　end</b></p><p>　　程序二：本程序無卷積，擴頻周期L=15（文件名為DS1）</p><p>　　EN = input(' 請輸入信噪比（用分貝表示）: ');</p><p>　　L=input(

81、9;請輸入數(shù)據(jù)幀長度(=信息比特數(shù)+尾比特數(shù),默認值: 400):');</p><p>　　ferrlim=input(' 請輸入程序終止的幀錯誤數(shù), 默認值: 15 : ');</p><p>　　if isempty(ferrlim)</p><p>　　ferrlim = 15;</p><p><b&

82、gt;　　end </b></p><p>　　if isempty(L)</p><p><b>　　L = 400;</b></p><p><b>　　end</b></p><p>　　c=[-1 -1 -1 1 1 1 1 -1 1 -1 1 1 -1 -1 1];</

83、p><p>　　N=length(c);</p><p>　　g=[1 0 1;1 1 1];k0=1;</p><p>　　for j=1:length(EN);</p><p>　　errs(j) = 0;</p><p>　　nferr(j) = 0;</p><p>　　nframe = 0

84、;</p><p>　　while nferr(j)<ferrlim</p><p>　　nframe = nframe + 1;</p><p>　　tem=rand(1,L);</p><p><b>　　for i=1:L</b></p><p>　　if tem(i)<=0.5

85、</p><p><b>　　d(i)=0;</b></p><p><b>　　else</b></p><p><b>　　d(i)=1;</b></p><p><b>　　end</b></p><p><b>　

86、　end</b></p><p>　　n=(1/10^(EN(j)/10))*randn(1,L*N);</p><p><b>　　b=2*d-1;</b></p><p><b>　　for i=1:L</b></p><p>　　y(i)=N*b(i)+n((i-1)*N+1:i*

87、N)*c';</p><p>　　if y(i)>=0</p><p><b>　　Db(i)=1;</b></p><p><b>　　else</b></p><p><b>　　Db(i)=0;</b></p><p><b&g

88、t;　　end</b></p><p><b>　　end</b></p><p>　　err=length(find(Db~=d));</p><p><b>　　if err>0</b></p><p>　　nferr(j)=nferr(j)+1;</p><

89、;p><b>　　end</b></p><p>　　errs(j) = errs(j) + err;</p><p>　　if rem(nframe,3)==0 | nferr(j)==ferrlim</p><p>　　% Bit error rate</p><p>　　ber(j) = errs(j)/nf

90、rame/(L);</p><p>　　% Frame error rate</p><p>　　fer(j) = nferr(j)/nframe;</p><p>　　% Display intermediate results in process </p><p>　　fprintf('**************Eb/N0=

91、%5.2f db **************\n', EN(j));</p><p>　　fprintf('Frame size = %d, rate 1/%d. \n', L, 2);</p><p>　　fprintf('%d frames transmitted, %d frames in error.\n', nframe, nferr(j

92、));</p><p>　　fprintf('%8.4e ', ber(j)); </p><p>　　fprintf('\n'); </p><p>　　fprintf('%8.4e ', fer(j));</p><p>　　fprintf('\n'

93、);</p><p>　　fprintf('***********************************************\n\n');</p><p><b>　　end</b></p><p>　　% Save intermediate results </p><p>　　end

94、%while</p><p>　　end %nEN</p><p>　　semilogy(EN,ber,'.');</p><p><b>　　end</b></p><p>　　程序三：本程序無卷積，擴頻周期L=31（文件名為DS2）</p><p>　　EN = input

95、(' 請輸入信噪比（用分貝表示）: ');</p><p>　　L=input('請輸入數(shù)據(jù)幀長度(=信息比特數(shù)+尾比特數(shù),默認值: 400):');</p><p>　　ferrlim=input(' 請輸入程序終止的幀錯誤數(shù), 默認值: 31: ');</p><p>　　if isempty(ferrlim

96、)</p><p>　　ferrlim = 31;</p><p><b>　　end </b></p><p>　　if isempty(L)</p><p><b>　　L = 400;</b></p><p><b>　　end</b></

97、p><p>　　c=[-1 1 1 1 -1 -1 -1 -1 1 -1 -1 -1 -1 1 1 -1 -1 1 -1 -1 1 -1 1 1 1 1 -1 -1 -1 -1 -1 ];</p><p>　　N=length(c);</p><p>　　g=[1 0 1;1 1 1];k0=1;</p><p>　　for j=1:lengt

98、h(EN);</p><p>　　errs(j) = 0;</p><p>　　nferr(j) = 0;</p><p>　　nframe = 0;</p><p>　　while nferr(j)<ferrlim</p><p>　　nframe = nframe + 1;</p><p&

99、gt;　　tem=rand(1,L);</p><p><b>　　for i=1:L</b></p><p>　　if tem(i)<=0.5</p><p><b>　　d(i)=0;</b></p><p><b>　　else</b></p><

100、;p><b>　　d(i)=1;</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p>　　n=(1/10^(EN(j)/10))*randn(1,L*N);</p><p><b>　　b=2*d-1;

101、</b></p><p><b>　　for i=1:L</b></p><p>　　y(i)=N*b(i)+n((i-1)*N+1:i*N)*c';</p><p>　　if y(i)>=0</p><p><b>　　Db(i)=1;</b></p>&l

102、t;p><b>　　else</b></p><p><b>　　Db(i)=0;</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p>　　err=length(find(Db~=d));&

103、lt;/p><p><b>　　if err>0</b></p><p>　　nferr(j)=nferr(j)+1;</p><p><b>　　end</b></p><p>　　errs(j) = errs(j) + err;</p><p>　　if rem(nfr

104、ame,3)==0 | nferr(j)==ferrlim</p><p>　　% Bit error rate</p><p>　　ber(j) = errs(j)/nframe/(L);</p><p>　　% Frame error rate</p><p>　　fer(j) = nferr(j)/nframe;</p>

105、<p>　　% Display intermediate results in process </p><p>　　fprintf('**************Eb/N0=%5.2f db **************\n', EN(j));</p><p>　　fprintf('Frame size = %d, rate 1/%d. \n'

106、, L, 2);</p><p>　　fprintf('%d frames transmitted, %d frames in error.\n', nframe, nferr(j));</p><p>　　fprintf('%8.4e ', ber(j)); </p><p>　　fprintf('\n&

107、#39;); </p><p>　　fprintf('%8.4e ', fer(j));</p><p>　　fprintf('\n');</p><p>　　fprintf('***********************************************\n\n');</p>&l

108、t;p><b>　　end</b></p><p>　　% Save intermediate results </p><p>　　end%while</p><p>　　end %nEN</p><p>　　semilogy(EN,ber,'.');</p><p&g

109、t;<b>　　end</b></p><p>　　程序四 ：本程序包含卷積，擴頻周期L=15</p><p>　　EN = input(' 請輸入信噪比（用分貝表示）: ');</p><p>　　L= input(' 請輸入數(shù)據(jù)幀長度(= 信息比特數(shù) + 尾比特數(shù), 默認值: 400) : ');<

110、/p><p>　　if isempty(L)</p><p><b>　　L = 400;</b></p><p><b>　　end</b></p><p>　　ferrlim=input(' 請輸入程序終止的幀錯誤數(shù), 默認值: 15 : ');</p><p

111、>　　if isempty(ferrlim)</p><p>　　ferrlim = 15;</p><p><b>　　end </b></p><p>　　c=[-1 -1 -1 1 1 1 1 -1 1 -1 1 1 -1 -1 1]; % 產(chǎn)生m序列</p><p>　　N=length(c);

112、</p><p>　　g=[1 0 1;1 1 1];k0=1; % 卷積生成矩陣</p><p>　　for j=1:length(EN);</p><p>　　errs(j) = 0; % 賦初值 </p><p>　　nferr(j) = 0;&

113、lt;/p><p>　　nframe = 0;</p><p>　　while nferr(j)<ferrlim</p><p>　　nframe = nframe + 1;</p><p>　　tem=rand(1,L);</p><p>　　for i=1:L %判決生成單極性

114、序列d(t)</p><p>　　if tem(i)<=0.5</p><p><b>　　d(i)=0;</b></p><p><b>　　else</b></p><p><b>　　d(i)=1;</b></p><p><b>

115、　　end</b></p><p><b>　　end</b></p><p>　　b1=cnv_encd(g,k0,d); %卷積生成矩陣b1</p><p>　　L1=length(b1);</p><p>　　n=(1/10^(EN(j)/10))*randn(1,L1*N)

116、; %生成噪聲</p><p>　　b=2*b1-1; %通過公式計算生成雙極性的序列b(t)</p><p>　　for i=1:L1 %判決</p><p>　　y(i)=N*b(i)+n((i-1)*N+1:i*N)*c'; %擴頻輸出y(i)</p>

117、<p>　　if y(i)>=0 %判決生成擴頻后的序列Db(t)</p><p><b>　　Db(i)=1;</b></p><p><b>　　else</b></p><p><b>　　Db(i)=0;</b></p><p>

118、<b>　　end</b></p><p><b>　　end</b></p><p>　　[Dd,survivor_state,cumulated_metric]=viterbi(g,k0,Db); </p><p>　　%viterbi譯碼</p><p>　　err=length(find(D

119、d~=d)); %比較Db與d</p><p><b>　　if err>0</b></p><p>　　nferr(j)=nferr(j)+1;</p><p><b>　　end</b></p><p>　　errs(j) = errs(j) + err;</p

120、><p>　　if rem(nframe,3)==0 | nferr(j)==ferrlim</p><p>　　% Bit error rate</p><p>　　ber(j) = errs(j)/nframe/(L);</p><p>　　% Frame error rate</p><p>　　fer(j) = n

121、ferr(j)/nframe;</p><p>　　% Display intermediate results in process </p><p>　　fprintf('**************Eb/N0=%5.2fdb **************\n', EN(j));</p><p>　　fprintf('Frame siz

122、e = %d, rate 1/%d. \n', L, 2);</p><p>　　fprintf('%d frames transmitted, %d frames in error.\n', nframe, nferr(j));</p><p>　　fprintf('%8.4e ', ber(j)); </p>&