通信系統(tǒng)綜合課程設計

1、課程設計任務書學生姓名： 專業(yè)班級：指導教師： 工作單位： 題 目：通信系統(tǒng)課群綜合訓練與設計1課程設計的目的通過課程設計，使學生加強對電子電路的理解，學會對電路分析計算以及 設計。進一步提高分析解決實際問題的能力，通過完成綜合設計型和創(chuàng)新性實驗 及訓練，創(chuàng)造一個動腦動手、獨立開展電路實驗的機會， 鍛煉分析、解決電子電 路問題的實際本領，實現(xiàn)由課本知識向?qū)嶋H能力的轉(zhuǎn)化；加深對通信原理的理解, 提高學生對現(xiàn)代通信系統(tǒng)的全面認識，增強學生的實踐能力。2課程設計要求要求：掌握以上各種電路與通信技術(shù)的基本原理，掌握實驗的設計、電路 調(diào)試與測量的方法。1. 培養(yǎng)學生根據(jù)需要選學參考書，查閱手冊，圖表和文

2、獻資料的自學能力， 通過獨立思考、深入鉆研有關(guān)問題，學會自己分析解決問題的方法。2. 通過對實驗電路的分析計算，了解簡單實用電路的分析方法和工程設計 方法。3. 掌握示波器，頻譜儀，失真度儀的正確使用方法，學會簡單電路的實驗 調(diào)試和整機指標測試方法，提高動手能力。3課程設計進度安排序號設計 內(nèi) 容所用時間1根據(jù)設計任務，分析電路原理，確定實驗方案2天2根據(jù)實驗條件進行電路的測試，并對結(jié)果進行分析7天3 1撰寫課程設計報告1天合計2周指導教師簽名：系主任（或責任教師）簽名:目 錄摘要IAbstractII1設計任務11.1設計目的11.2課程設計要求11.3課程設計任務12實驗原理分析22.1

3、實驗原理框圖 22.2 PCM原理介紹32.2.1 扌由樣（Sampling） 32.2.2 量化（quantizing ） 32.2.3 A壓縮律：42.2.4 編碼（Coding ） 92.2基帶傳輸PST碼92.3循環(huán)碼10循環(huán)碼介紹10循環(huán)碼編碼原理12循環(huán)碼的糾錯原理 132.4二進制頻移鍵控（2PSK調(diào)制與解調(diào)原理 152.4.1 數(shù)字調(diào)制技術(shù)15二進制相移鍵控（2PSK）基本原理 16實現(xiàn)2PSK調(diào)制與解調(diào) 182.5衰落信道183 MATLAB軟件仿真及結(jié)果 193.1 MATLAB 簡介193.2 發(fā)送端仿真結(jié)果與分析 19仿真結(jié)果圖193.2.2 實驗結(jié)果分析 244.小結(jié)

4、25參考文獻26附錄：27摘要通信系統(tǒng)的作用就是將信息從信源發(fā)送到一個或多個目的地。 一個數(shù)字通信 系統(tǒng)是利用數(shù)字信號來傳遞信息的通信系統(tǒng)， 主要有信源編碼與譯碼、 信道編碼 與譯碼、數(shù)字調(diào)制與解調(diào)。 具體過程是把信息源的模擬信號經(jīng)過模數(shù)轉(zhuǎn)換變成數(shù) 子信號然后，再進行基帶編碼、信道編碼、載波調(diào)制然后發(fā)送到信道中，然后造 接收端進行對應的解調(diào)、譯碼、數(shù)模轉(zhuǎn)換得到最終的接收信號。本設計結(jié)合PCM勺抽樣、量化、編碼原理把模擬信號轉(zhuǎn)換數(shù)字信號、 漢明碼 編譯碼原理實現(xiàn)信道編譯碼、PST編譯碼原理用于實現(xiàn)信源編譯碼、PSK調(diào)制解 調(diào)原理用于實現(xiàn)數(shù)字調(diào)制與解調(diào)、 衰落信道原理用于信道和噪聲源勺近似， 利用

5、 MATLAB件編程，完成了整個通信系統(tǒng)編碼仿真分析。關(guān)鍵詞：MATLAB件編程、數(shù)字通信系統(tǒng)、PCM PST碼、循環(huán)碼、PSKAbstractCommunication system is the role of the information from the source sent to one or more of the desti nati on. A digital com muni cati on system is to use digital sig nal to tran sfer in formatio n com muni cati on system, basica

6、lly have source codi ng and decoding, channel coding and decoding, digital modulation and demodulation. The specific process is the source of the an alog sig nal after a/d convert into nu mber son signal and then, again carries on the baseband coding, channel coding, carrier modulation and then sent

7、 to the channel, then made on the receiving end corresponding demodulation, decoding, d/a conversion to get the final receiving sig nal.This design with PCM sampling and quantization and coding principle to an alog sig nal into the digital sig nal, ham ming code kn itt ing decod ing theory to realiz

8、e channel knitting decoding, PST knitting decoding principle used to implement the source knitting decoding, PSK modulation demodulation principle used to implement the digital modulati on and demodulati on, Fadi ng cha nnel and the prin ciple for no ise approximation, using MATLAB software programm

9、ing, complete the whole com muni cati on system cod ing simulatio n an alysis.Keywords: MATLAB software program ming, digital com muni catio n system, PCM, PST, Cyclic code, PSK1設計任務1.1設計目的通過課程設計，使學生加強對電子電路的理解，學會對電路分析計算以及設計。進一 步提高分析解決實際問題的能力，通過完成綜合設計型和創(chuàng)新性實驗及訓練，創(chuàng)造一個動 腦動手、獨立開展電路實驗的機會，鍛煉分析、解決電子電路問題的實際本

10、領，實現(xiàn)由課 本知識向?qū)嶋H能力的轉(zhuǎn)化；加深對通信原理的理解，提高學生對現(xiàn)代通信系統(tǒng)的全面認識， 增強學生的實踐能力。1.2課程設計要求要求：掌握以上各種電路與通信技術(shù)的基本原理，掌握實驗的設計、電路調(diào)試與測量 的方法。1. 培養(yǎng)學生根據(jù)需要選學參考書，查閱手冊，圖表和文獻資料的自學能力，通過獨立 思考、深入鉆研有關(guān)問題，學會自己分析解決問題的方法。2. 通過對實驗電路的分析計算，了解簡單實用電路的分析方法和工程設計方法。3. 掌握示波器，頻譜儀，失真度儀的正確使用方法，學會簡單電路的實驗調(diào)試和整機 指標測試方法，提高動手能力。1.3課程設計任務利用仿真軟件（如Matlab或SystemView

11、），或硬件實驗系統(tǒng)平臺上設計完成一個典型 的通信系統(tǒng)。學生要完成整個系統(tǒng)各環(huán)節(jié)以及整個系統(tǒng)的仿真，最終在接收端或者精確或者近似地 再現(xiàn)輸入（信源），計算失真度，并且分析原因。（信源為模擬信源，數(shù)字化方式為增量調(diào) 制，基帶碼為AMI碼，信道碼為循環(huán)碼，調(diào)制方式為 ASK調(diào)制，信道類型為衰落信道）。完成整個系統(tǒng)各環(huán)節(jié)以及整個系統(tǒng)的仿真，最終在接收端或者精確或者近似地再現(xiàn)輸入（信源），計算失真度，并且分析原因。信源：自己構(gòu)造一時間函數(shù)，數(shù)字化方式 ：PCM， 基帶碼:PST，信道碼:循環(huán)碼，調(diào)制方式:PSK，信道類型:衰落信道；解調(diào)，信道解碼，基 帶解碼，數(shù)模轉(zhuǎn)換的過程與輸入端對應。2實驗原理分析2

12、.1實驗原理框圖圖2.1通信系統(tǒng)的原理框圖模數(shù)轉(zhuǎn)換中有三個基本過程：抽樣、量化、編碼。根據(jù)編碼的不同有幾種調(diào)制方法， 如脈沖振幅調(diào)制（PAM）,脈沖編碼調(diào)制（PCM），差分脈沖編碼調(diào)制（DPCM）,增量調(diào) 制。基帶傳輸?shù)某Ｓ么a型有 AMI碼，HDB3碼，PST碼，曼徹斯特碼，密勒碼，CMI碼。 信道編碼的常用編碼有漢明碼，卷積碼，循環(huán)碼，BCH碼。常見的數(shù)字調(diào)制方式有振幅鍵控（ASK），頻移鍵控（FSK）,相移鍵控（PSK）。 常見的信道有加性高斯白噪聲信道，多徑衰落信道。根據(jù)題目要求，數(shù)字化方式為 PCM調(diào)制，基帶碼為PST碼，信道碼為循環(huán)碼，數(shù)字 調(diào)制方式為PSK調(diào)制，信道為衰落信道，則有

13、以下原理框圖圖2.2實驗原理框圖2.2 PCM原理介紹模擬信號數(shù)字化必須經(jīng)過三個過程，即抽樣、量化和編碼，以實現(xiàn)話音數(shù)字化的脈沖編碼調(diào)制(PCM，Pulse Coding Modulation)技術(shù)。2.2.1 抽樣(Sampli ng)離散時間信號通常是有連續(xù)時間信號經(jīng)周期采樣得到的。完成采樣功能的器件稱為采樣器，下圖所示為采樣器的示意圖。圖中Xa(t)表示模擬信號，Xa(nt)表示采樣信號，T為采樣周期，n=0，1，2，。一般可以把采樣器視為一個每隔 T秒閉合一次的電子開關(guān)So 在理想情況下，開關(guān)閉合時間 t滿足t v<T。實際采樣過程可視為脈沖調(diào)幅過程，Xa(t)為調(diào)制信號，被調(diào)脈

14、沖載波p(t)是周期為T、脈寬為t的周期脈沖串。當T 0時的理想 采樣情況是實際采樣的一種科學的、本質(zhì)的抽象，同時可使數(shù)學推導得到簡化。下面主要 討論理想采樣。2.2.2 量化(quantizing )抽樣信號雖然是時間軸上離散的信號，但仍然是模擬信號，其樣值在一定的取值范圍 內(nèi)，可有無限多個值。顯然，對無限個樣值給出數(shù)字碼組來對應是不可能的。為了實現(xiàn)以數(shù)字碼表示樣值，必須采用“四舍五入”的方法把樣值分級“取整”，分的級數(shù)越多,即量化級差或間隔越小，量化噪聲也越小。223 A壓縮律：所謂A壓縮律也就是壓縮器具有如下特性的壓縮律:1 十 In / y - <1十血川穴.1 + ln A其中

15、，A為壓縮系數(shù)；y為歸一化的壓縮器輸出電壓；x為歸一化的壓縮器輸入電壓。圖畫出了 A為某一取值的歸一化壓縮特性。A律壓縮特性是以原點奇對稱的，為了簡便,圖中只給出了正半軸部分。圖2.4 A壓縮律特性上圖中，x和y都在-1和+1之間，取量化級數(shù)為N（在y方向上從-1到+1被均勻劃分 為N個量化級），則量化間隔為Aj = 當N很大時，在每一量化級中壓縮特性曲線可看作是直線，因此有2 dx式中，xi為第i個量化級間隔的中間值 因此(2.1)為了使量化信噪比不隨信號 而變小，即應使各量化級間隔與x變化，也就是說在小信號時的量化信噪比不因 x成線性關(guān)系，即Ax, oc Xjx的減小則式（2.1）可寫成其

16、中k為比例常數(shù)。當量化級數(shù)很大時，可以將它看成連續(xù)曲線，因而式(2.2)成為線性微分方程解此微分方程其中C為常數(shù)。為了滿足歸一化要求，當x=1時,k H-c 0y=1,代入式(2.3)可得故所得結(jié)果為In工二炒一氏如果壓縮特性滿足上式，就可獲得理想的壓縮效果,(2.2)(2.3)(2.4)其量化信噪比和信號幅度無關(guān)。滿圖2.5理想壓縮特性曲線足上式的曲線如下圖所示,定的修改。A律壓縮特性就是對式(2.4)修改后的函數(shù)。在上圖中，通過原點作理想壓縮特性曲線 的切線0C,將oc、cd作為實際的壓縮特性。修改以后，必須用兩個不同的方程來描述這 段曲線，以切點c為分界點，線段0C的方程:設切點c的坐標

17、為(x1，y1)斜率為ck則由式(2.4)可得(2.5)所以線段oc的方程為1所以當x=x1時，y1=1/k時，有1 , L=1 + 一n x, k k因此有芯=嚴宀所以，切點坐標為(exp-(k-1),1/k)，令可二丄二嚴1丄則上三1+ In將它代入式(2.5),就可得到以切點c為邊界的段的方程為y =, 0 c x < (2.6)(2.7)l+lii JA因cd段的方程，滿足式(3.4),所以由該式可得1.1 + In 貝 xV = 1 +In 芹匚1+1hj4 Uln A由以上分析可見，經(jīng)過修改以后的理想壓縮特性與圖5中所示的曲線近似，而式(2.6)式(2.7)和式(2.4)完全

18、一樣13折線：實際中，A壓縮律通常采用13折線來近似，13折線法如圖2.4所示，圖中 先把丫軸的0,1區(qū)間分為8個不均勻段。7/86/84/fl2/8/81/41/2*二 u圖2.6折線示意圖1/64 i/32其具體分法如下：a將區(qū)間0，1一分為二，其中點為1/2,取區(qū)間1/2,1 作為第八段;b. 將剩下的區(qū)間0,1/2再一分為二，其中點為1/4,取區(qū)間1/4,1/2作為第七段；c. 將剩下的區(qū)間0,1/4再一分為二，其中點為1/8,取區(qū)間1/8,1/4作為第六段；d. 將剩下的區(qū)間0,1/8再一分為二，其中點為1/16,取區(qū)間1/16,1/8作為第五段;e. 將剩下的區(qū)間0,1/16再一分

19、為二，其中點為1/32,取區(qū)間1/32,1/16作為第四段;f. 將剩下的區(qū)間0,1/32再一分為二，其中點為1/64,取區(qū)間1/64,1/32作為第三段；g. 將剩下的區(qū)間0,1/64再一分為二，其中點為1/128取區(qū)間1/128,1/64作為第二段;h.最后 剩下的區(qū)間0,1/128作為第一段。然后將y軸的0,1區(qū)間均勻地分成八段，從第一段到第八段分別為 0,1/8,(1/8,2/8,(2/8,3/8,(3/8,4/8,(4/8,5/8,(5/8,6/8,(6/8,7/8,(7/8,1。分別與 x 軸的八段 對應。采用上述的方法就可以作出由八段直線構(gòu)成的一條折線，該折線和A壓縮律近似，圖2

20、.4中的八段線段的斜率分別為：表1各段落的斜率段落12345678斜率161684211/21/4從上表中可以看出，除一、二段外，其他各段折線的斜率都不相同。圖2.4中只畫出了第一象限的壓縮特性，第三象限的壓縮特性的形狀與第一象限的壓縮特性的形狀相同，且它們以原點為奇對稱，所以負方向也有八段直線，總共有16個線段。但由于正向一、二兩段和負向一、二兩段的斜率相同，所以這四段實際上為一條直線，因此，正、負雙向 的折線總共由13條直線段構(gòu)成，這就是13折線的由來。從A律壓縮特性中可以看出，取 A=87.6主要基于下述兩個原因：1使壓縮特性曲線在原點附近的斜率為16;2當用13折線逼近時，的八段量化分

21、界點近似為1/2八n(n=0,1,2,7)。從表1可以看出，當要求滿足x=1/2An時，相應有y=1-n/8代入式中，有1 _ 1 _ 1 _ 1f =冏亠=嚴=用f因此有2 =將上式代入式(2.7),就可以得到對應A=94.4時的壓縮特性(2.8)1 + ln (j1)In(256a)l + lii Z山宓) In 256此壓縮特性如果用13折線逼近，除了第一段落起始點外，其余各段落的分界點的 x、 y都應滿足式(2.8)。在13折線中，第一段落起始點要求的 x、y都應該為零，而若按照式 (2.8)計算時，當x=0時，Ls;而當y=0, 乂=1/2八8。因此，需要對式(2.8)的壓縮特性曲線

22、 作適當?shù)男拚?，我們可以在原點和點(1/2A7,1/8)之間用一段直線代替原來的曲線，這段直線 的斜率是1/8- 1/2八7=16。為了找到一個能夠表示修正后的整個壓縮特性曲線的方程，將式(2.8)變成_ ln(1 + 255z)-(2.9)從上式中可以看出，它滿足x=0時，y=0;x=1時，y=1。雖然式(2.9)在其他點上會有誤 差，但x在區(qū)間(1/128,1 內(nèi)，1+255X都能和原來的256x比較接近。所以，在絕大部分范 圍內(nèi)的壓縮特性仍和A律壓縮特性非常接近，只有在x-0的小信號部分和A律壓縮特性 有些差別若在式(2.9)中，令 卩=255，則式(2.9)可寫成（2.10）ha 1

23、+ 占）血1+山式(2.10)的壓縮特性與 卩律壓縮特性完全一致。(2)按照量化的維數(shù)分，量化分為標量量化和矢量量化。標量量化是一維的量化， 一個幅度對應一個量化結(jié)果。而矢量量化是二維甚至多維的量化，兩個或兩個以上的幅度 決定一個量化結(jié)果。以二維情況為例，兩個幅度決定了平面上的一點。而這個平面事先按照概率已經(jīng)劃分 為N個小區(qū)域，每個區(qū)域?qū)粋€輸出結(jié)果(碼數(shù)，codebook)。由輸入確定的那一點落在了哪個區(qū)域內(nèi)，矢量量化器就會輸出那個區(qū)域?qū)拇a字(codeword)。矢量量化的好處是引入了多個決定輸出的因素，并且使用了概率的方法，一般會比標量量化效率更高。編碼（Coding）量化后的抽樣

24、信號在一定的取值范圍內(nèi)僅有有限個可取的樣值，且信號正、負幅度分布的對稱性使正、負樣值的個數(shù)相等，正、負向的量化級對稱分布。若將有限個量化樣值的絕對值從小到大依次排列，并對應地依次賦予一個十進制數(shù)字代碼(例如，賦予樣值 0的十進制數(shù)字代碼為0),在碼前以“ +”、“一”號為前綴，來 區(qū)分樣值的正、負，則 量化后的抽樣信號就轉(zhuǎn)化為按抽樣時序排列的一串十進制數(shù)字碼流，即十進制數(shù)字信號。簡單高效的數(shù)據(jù)系統(tǒng)是二進制碼系統(tǒng)，因此，應將十進制數(shù)字代碼變換成二進制編碼。根據(jù)十進制數(shù)字代碼的總個數(shù)，可以確定所需二進制編碼的位數(shù)，即字長。這種把量化的 抽樣信號變換成給定字長的二進制碼流的 過程稱為編碼。2.2基帶

25、傳輸 PST碼PST:Pair Selected Ternary成 對選擇三進碼。在實際的基帶傳輸系統(tǒng)中，并不是所有代碼的電波形都能在信道中傳輸。含有直流分 量和較豐富低頻分量的單極性基帶波形就不適宜在低頻傳輸特性差的信道中傳輸。PST碼是成對選擇三進碼。其編碼過程是：先將二進制代碼兩兩分組，然后再把每一碼組編碼成 兩個三進制數(shù)字(+、-、0)。因為兩位三進制數(shù)字共有9種狀態(tài)，故可靈活地選擇其中的 4種狀態(tài)。為防止PST碼的直流漂移，當在一個碼組中僅發(fā)送單個脈沖時，兩個模式應交 替變換。例如：表1 PST碼二進制代碼+模式-模式0 0-+-+0 10 +0 -1 0+ 0-01 1+ -+ -

26、PST碼能提供足夠的定時分量，且無直流成分，編碼過程也較簡單。但這種碼在識別時需要提供 分組”信息，即需要建立幀同步。 每位二進制信碼都被變換成1位三電平取值 （+1、0、-1）的碼，因而有時把這類碼稱為IB/ 1T碼。2.3循環(huán)碼231循環(huán)碼介紹循環(huán)碼是線性分組碼的一種，所以它具有線性分組碼的一般特性，此外還具有循環(huán)性。 循環(huán)碼的編碼和解碼設備都不太復雜，且檢 （糾）錯能力強。它不但可以檢測隨機的錯誤， 還可以檢錯突發(fā)的錯誤。（n,k）循環(huán)碼可以檢測長為n-k或更短的任何突發(fā)錯誤，包括首 尾相接突發(fā)錯誤。循環(huán)碼是一種無權(quán)碼，循環(huán)碼編排的特點是相鄰兩個數(shù)碼之間符合卡諾圖中的鄰接條 件，即相鄰兩

27、個數(shù)碼之間只有一位碼元不同，碼元就是組成數(shù)碼的單元。符合這個特點的 有多種方案，但循環(huán)碼只能是表中的那種。循環(huán)碼的優(yōu)點是沒有瞬時錯誤，因為在數(shù)碼變 換過程中，在速度上會有快有慢，中間經(jīng)過其它一些數(shù)碼形式，稱它們?yōu)樗矔r錯誤。這在 某些數(shù)字系統(tǒng)中是不允許的，為此希望相鄰兩個數(shù)碼之間僅有一位碼元不同，即滿足鄰接 條件，這樣就不會產(chǎn)生瞬時錯誤。循環(huán)碼就是這樣一種編碼，它可以在卡諾圖中依次循環(huán) 得到。循環(huán)碼最大的特點就是碼字的循環(huán)特性，所謂循環(huán)特性是指：循環(huán)碼中任一許用碼組 經(jīng)過循環(huán)移位后，所得到的碼組仍然是許用碼組。若（an 1 an 2aj a。）為一循環(huán)碼組，則（an 2 an 3a。a. 1）、

28、（an 3 an 4an 1 an 2） 還是許用碼組。也就是說，不論是左移還是右移，也不 論移多少位，仍然是許用的循環(huán)碼組(1) 循環(huán)碼的多項式表示設碼長為n的循環(huán)碼表示為(an1務2a：印a。)，其中ai為二進制數(shù)，通常把碼組中各碼元當做二進制的系數(shù)，即把上式中長為n的各個分量看做多項式： Tx anixan2xaiXaix a。(2.11)的各項系數(shù)，則碼字與碼多項式一一對應，這種多項式中，x僅表示碼元位置的標記，因此我們并不關(guān)心x的取值，這種多項式稱為碼多項式。(2) (n,k)循環(huán)碼的生成多項式(n,k)循環(huán)碼的生成多項式寫為g(x),它是(n,k)循環(huán)碼碼集中唯一的，幕次為n-k的

29、碼多項式,則xkg(x)是一個幕次為n的碼多項式。按模(xn 1)運算，此時:Q ( x) xkg(x)R (x)nR(x)(2.12)(2.13)且因xk g(x)也是n階幕，故Q(x)=1。由于它是循環(huán)碼，故xkg(x)按模(xn 1)運算后的 余 式”也是循環(huán)碼的一個碼字，它必能被g(x)整除，即：3 F(x)(214)由以上兩式可以得到：g(x)xkg(x)Q(x)(xn1) R(x) (xn 1) f(x)g(x)( 2.15)和xn1 xk f (x) g(x) h(x)g(x)(2.16)從上式中可以看出，生成多項式g(x)應該是xn 1的一個因式，即循環(huán)碼多項式應該是xn 1的

30、一個n-k次因式。(3)循環(huán)碼的生成矩陣和一致校驗矩陣對所有的i=0,1,2, -1,用生成多項式g(x)除xn ki，有：xn k ：a： (x)g(x) bi(x)(2.17)式中bi(x)是余式，表示為：bi (x)b,n k 1Xn k 1b,1Xbi,0(2.18)因此，xn k i b(x)是g(x)的倍式，即xn k 1bi(x)是碼多項式，由此得到系統(tǒng)形式的生成矩陣為：1oobk 1,n k 1bk 1,1bk 1,oGo1obk 2,n k 1bk 2,1bkx,ooooobo, n k 1bo,1bo,o(2.19)它是一個 k n 階的矩陣。同樣，由 G HT=0 可以得

31、到系統(tǒng)形式的一致校驗矩陣為：bk 1, n k 1bk 2,n k 1bo,n k 11ooH( 2.2o)bk 1,1bk 2,1bo,1o1obk 1,obk 2,obo,ooo1g(x) x3 x 1 ，42h(x) xxx 1 。寫得其生成矩陣和校驗矩陣分別為1o11oooo1o11ooGoo1o11oooo1o112.21)如已知( 7， 4)循環(huán)碼的生成多項式和校驗多項式分別為：111o1ooHo111o1ooo111o12.22)2.3.2 循環(huán)碼編碼原理有信息碼構(gòu)成信息多項式 m(x) mk 1xk 1mo，其中最高幕次為k-1 ；用xn k乘以信息多項式m(x)，得到的xn

32、km(x)，最高幕次為n-1 ，該過程相當于把信息碼mk 1 ，mk 2,m!, mo)移位到了碼字德前k個信息位，其后是r個全為零的監(jiān)督位；用g(x)除xn km(x)得到余式r(x),其次數(shù)必小于g(x)的次數(shù)，即小于(n-k)，將此r(x)加于信息位后做監(jiān)督位，即將r(x)于xn km(x)相加，得到的多項式必為一碼多項式。(7, 4)循環(huán)碼編碼的程序框圖如圖 2.5 所示：圖2.7編碼程序框圖循環(huán)碼的糾錯原理糾錯碼的譯碼是該編碼能否得到實際應用的關(guān)鍵所在。譯碼器往往比編碼較難 實現(xiàn)，對于糾錯能力強的糾錯碼更復雜。根據(jù)不同的糾錯或檢錯目的，循環(huán)碼譯碼器 可分為用于糾錯目的和用于檢錯目的的

33、循環(huán)碼譯碼器。通常，將接收到的循環(huán)碼組進行除法運算，如果除盡，則說明正確傳輸；如果未 除盡，則在寄存器中的內(nèi)容就是錯誤圖樣，根據(jù)錯誤圖樣可以確定一種邏輯，來確定 差錯的位置，從而達到糾錯的目的。用于糾錯目的的循環(huán)碼的譯碼算法比較復雜，感 興趣的話可以參考一些參考書。而用于檢錯目的循環(huán)碼，一般使用ARQ通信方式。檢測過程也是將接受到的碼組進行除法運算，如果除盡，則說明傳輸無誤；如果未除盡，則表明傳輸出現(xiàn)差錯， 要求發(fā)送端重發(fā)。用于這種目的的循環(huán)碼經(jīng)常被成為循環(huán)冗余 校驗碼，即CRC校驗碼。CRC校驗碼由于編碼電路、檢錯電路簡單且易于實現(xiàn)，因此 得到廣泛的應用。在通過MODEMS輸文件的協(xié)議如 Z

34、MODEMXMODE協(xié)議中均用到了CRC校驗技術(shù)。在磁盤、光盤介質(zhì)存儲技術(shù)中也使用該方法。當碼字c通過噪聲信道傳送時，會受到干擾而產(chǎn)生錯誤。如果信道產(chǎn)生的錯誤圖樣是e,譯碼器收到的n重接受矢量是y,則表示為:yc e(2.23)上式也可以寫成多項式形式：y(x)c(x)e(x)(2.24)譯碼器的任務就是從y(x)中得到Ae(x)，然后求的估值碼字Ac(x)y(x)Ae(x)(2.25)并從中得到信息組m(x)循環(huán)碼譯碼可按以下三個步驟進行：(1) 有接收到的y(x)計算伴隨式s(x);(2) 根據(jù)伴隨式s(x)找出對應的估值錯誤圖樣e(x)；AAA(3) 計算c(x) y(x) e(x)，得

35、到估計碼字c(x)。若c(x) c(x)，則譯碼正確，否則，A若c(x)c(x)，則譯碼錯誤。由于g(x)的次數(shù)為n - k次，g(x)除E(x)后得余式(即伴隨式)的最高次數(shù)為 n-k-1次， 故S(x)共有2n-k個可能的表達式，每一個表達式對應一個錯誤格式?？梢灾?(7,4)循環(huán) 碼的S(x)共有2(7-4)= 8個可能的表達式，可根據(jù)錯誤圖樣表來糾正(7,4)循環(huán)碼中的一位 錯誤，其伴隨式如表2所示。表2 BCH( 7，4)循環(huán)碼錯誤圖樣表錯誤圖樣錯誤圖樣碼字伴隨式S(x)伴隨式E6(X)=X 61000000x2100E5(X)=X 501000002 X +x110E4(X)=X

36、 400100002x2+x+1111E3(X)=X 30001000x+1011E2(X)=X 20000100x2+1101E1 (x)=x 10000010X010Eo(x)=x 000000011001E(x)=000000000000上式指出了系統(tǒng)循環(huán)碼的譯碼方法：將收到的碼字R(x)用g(x)去除，如果除盡則無錯；否則有錯。如果有錯，可由余式 S(x) 一一找出對應圖樣，然后將錯誤圖樣E(x)與R(x)模2和，即為所求碼字C(x)，從而實現(xiàn)糾錯目的。根據(jù)前面的討論，可得(7, 4)循環(huán)碼譯碼的程序框圖如圖2.6所示：圖2.8譯碼程序框圖2.4二進制頻移鍵控(2PSK調(diào)制與解調(diào)原理數(shù)

37、字調(diào)制技術(shù)為了使數(shù)字信號在帶通信道中傳輸，必須用數(shù)字基帶信號對載波進行調(diào)制以使信號與 信道的特性相匹配。這種用數(shù)字信號控制載波，把數(shù)字基帶信號變換為數(shù)字帶通信號(已 調(diào)信號)的過程成為數(shù)字調(diào)制。在接收端通過解調(diào)器把帶通信號還原成數(shù)字基帶信號的過 程成為數(shù)字解調(diào)。通常把包括調(diào)制和解調(diào)過程的數(shù)字傳輸系統(tǒng)叫做數(shù)字帶通傳輸系統(tǒng)。一般來說，數(shù)字調(diào)制與模擬調(diào)制的基本原理相同，但是數(shù)字調(diào)制有離散取值的特點。 因此數(shù)字調(diào)制技術(shù)有兩種方法： 利用模擬調(diào)制的方法去實現(xiàn)數(shù)字式調(diào)制； 通過開關(guān)鍵控載波，通常稱為鍵控法。基本鍵控方式有振幅鍵控、頻移鍵控、相移 鍵控振幅鍵控頻移鍵控相移鍵控圖2.9正弦載波的三種鍵控波形二

38、進制相移鍵控(2PSK基本原理2PSK數(shù)字調(diào)制可分為二進制調(diào)制和多進制調(diào)制相移鍵控是利用載波的相位變化來傳遞數(shù)字信息，而振幅和頻率保持不變。在中，通常用初始相位0和n分別表示二進制“T和“ 0”。因此，2PSK信號的時域表達式為:e2PSK（t）Acos（ctn）式中，n表示第n個符號的絕對相位：0，發(fā)送,發(fā)送“因此，上式可以改寫為0”1”A cos ct, A cos ct, 由于表示信號的兩種碼元的波形相同，極性相反,e2PSK（t）概率為概率為1故2PSK信號可以表述為一個雙極性全占空矩形脈沖序列與一個正弦載波的相乘:e2psK(t) st cos ct式中s(t)ang(t nTs)n

39、這里，g(t)是脈寬為Ts的單個矩形脈沖，而an的統(tǒng)計特性為1, 概率為P1,概率為1 P即發(fā)送二進制符號“ 0”時，e2psk(t)取0相位；發(fā)送二進制符號“ 1”時，Qpsk(t)取n相位。 這種以載波的不同相位直接去表示相應二進制數(shù)字信號的調(diào)制方式，稱為二進制絕對相移 方式。雙極性(a)模擬調(diào)制的方法開關(guān)電路(b)鍵控法圖2.10 2PSK信號調(diào)制器原理方框圖2PSK信號的解調(diào)通常采用相干解調(diào)法，解調(diào)器原理框圖如下圖2.9所示。假設數(shù)字信息為“ 1001T,相干解調(diào)各點波形如下圖2.9所示。圖2.11 2PSK信號的解調(diào)原理框圖abcde圖2.12 2PSK信號相干解調(diào)時各點時間波形24

40、3 MATLAB實現(xiàn)2PSK調(diào)制與解調(diào)1、2PSK調(diào)制流程圖圖2.13 2psk調(diào)制2、2PSK解調(diào)流程圖圖2.14 2psk解調(diào)2.5衰落信道瑞利衰落信道(Rayleigh fading channe)是一種無線電信號傳播環(huán)境的統(tǒng)計模型。這種模型假設信號通過無線信道之后，其信號幅度是隨機的，表現(xiàn)為衰落”特性，并且多徑衰落的信號包絡服從瑞利分布。由此，這種多徑衰落也稱為瑞利衰落。這一信道模型能夠描述由電離層和對流層反射的短波信道，以及建筑物密集的城市環(huán)境。瑞利衰落只適用 于從發(fā)射機到接收機不存在直射信號的情況，否則應使用萊斯衰落信道作為信道模型。假設經(jīng)反射(或散射)到達接收天線的信號為N個幅值

41、和相位均隨機的且統(tǒng)計獨立的信號之和。信號振幅為r,相位為，則其包絡概率密度函數(shù)為(2.26)(2.27)rtP(r)= 2 e 2 (r 0)相位概率密度函數(shù)為：P( )=1/2( 02 )3 MATLAB軟件仿真及結(jié)果3.1 MATLAB 簡介MATLAB是一個高級的矩陣/陣列語言，它包含控制語句、函數(shù)、數(shù)據(jù)結(jié)構(gòu)、輸入和輸 出和面向?qū)ο缶幊烫攸c。用戶可以在命令窗口中將輸入語句與執(zhí)行命令同步，也可以先編 寫好一個較大的復雜的應用程序（M文件）后再一起運行，可移植性好、可拓展性極強。MATLAB是一個包含大量計算算法的集合。 其擁有600多個工程中要用到的數(shù)學運算 函數(shù)，可以方便的實現(xiàn)用戶所需的

42、各種計算功能。MATLAB的這些函數(shù)所能解決的問題大 致包括矩陣運算和線性方程組的求解、微分方程及偏微分方程的組的求解、符號運算、傅 立葉變換和數(shù)據(jù)的統(tǒng)計分析、工程中的優(yōu)化問題、稀疏矩陣運算、復數(shù)的各種運算、三角 函數(shù)和其他初等數(shù)學運算、多維數(shù)組操作以及建模動態(tài)仿真等。3.2發(fā)送端仿真結(jié)果與分析3.2.1 仿真結(jié)果圖圖3. 1原始信號與PCk編碼輸出圖3.2 PST編碼圖3.3二進制轉(zhuǎn)三進制后循環(huán)碼信號輸出P圖3.4 PSK調(diào)制信號帶迫遊液后國埶2F2D01002003(10400500600 7DC別 0900n&叵1” 圖3.6衰落信道信號圖3.9 PST解碼后的信號圖解嗎信號圖3

43、.10 PCM解調(diào)后的信號322實驗結(jié)果分析實驗最后的結(jié)果與預想的還是差不多的， 調(diào)制，編碼和解調(diào)，解碼的過程都是正確的 有幾點可能對本設計存在較大影響：1.是取樣的點數(shù)我取了 50個不是很多會帶來誤差。 2.再就是信道噪聲會造成誤差。3.PSK解調(diào)的過程中通過低通濾波器和抽樣判決的過程也 都是有誤差的。4.最后用A律曲線逼近過程逼近時以一小段一小段的曲線去近似逼近。4.小結(jié)在本次課程設計中我學習了整個通信系統(tǒng)的工作原理，了解了PCM編解碼、循環(huán)碼編解碼、PST編解碼、PSK調(diào)制解調(diào)的原理和 MATLAB®程，也學習很多 MATLA的函數(shù)。雖然之前有接觸MATLA這個軟件，但是也只是

44、會調(diào)用一些簡單的函數(shù)來做計算， 或用 一些簡單的畫圖命令來作圖，一直覺得這個仿真軟件的功能很強大，但沒有怎么深入地去 運用。這次的課設讓我知道，我之前對它的了解簡直只是冰山一角。而要用它來實現(xiàn)自己 的專業(yè)目的，更是增加了難度。將在課堂上學到的理論知識用它來仿真實現(xiàn)，這就要求我 們對理論知識有全面深入透徹的理解， 并且很熟悉仿真軟件，但我在這兩個方面都有欠缺， 所以實際操作起來也沒有剛開始想象的那樣容易。這次的課程設計讓我收獲頗豐，同時也讓我發(fā)現(xiàn)了自身的不足。在此間發(fā)現(xiàn)的不足， 我將努力改善，通過學習、實踐等方式不斷提高，克服那些不應成為學習、獲得知識的障 礙。這里要感謝幫助我解決編程過程中遇到

45、問題的同學和老師。參考文獻1 樊昌信 曹麗娜，通信原理第六版，國防工業(yè)出版社，20072 周開利，鄧春暉主編 MATLAB基礎及其應用教程，北京大學出版社，20073 董振海，精通MATLAB 7編程與數(shù)據(jù)庫應用，電子工業(yè)出版社,20074 陳懷琛，MATLAB在電子信息課程中的應用（第二版），電子工業(yè)出版社，2006 張德豐主編，MATLAB通信工程仿真，機械工業(yè)出版社，2010附錄：% 主函數(shù) %clc close all clear all% 生成時間函數(shù) %t=1:0.1:5; s=sin(2*3.14*1000*t);%pcm 調(diào)制 %pcmcode,sMax=pcm(s);figu

46、re(1)subplot(2,1,1);plot(s); title(' 原始正弦信號 '); subplot(2,1,2);stairs(pcmcode);zoom xon;axis(0 750 -1 2);%pst 編碼 %pstcode=pst_c(pcmcode); figure(2);stairs(pstcode);zoom xon;axis(0 750 -1 2);% 三進制轉(zhuǎn)為二進制 %code=ZH_C(pstcode);figure(3);stairs(code);zoom xon;axis(0 750 -1 2);% 循環(huán)碼編碼 %xh=XH_C(code)

47、;figure(4);stairs(xh);zoom xon;axis(0 750 -1 2);%psk 調(diào)制 %psk,bit,mod=psk_encode(xh);figure(5); plot(psk); zoom xon;axis(0 750 -1 2);% figure(12);plot(psk);axis(0 1000 -1.2 1.2);title('PSK Modulation'); % 衰落信道 %chan el=rayleighcha n(1/64000,100);%生產(chǎn)瑞麗衰落信道 noisig0=filter(chanel,psk);% 把處理過的信號送

48、入信道 noisig=abs(noisig0);figure(6); plot(noisig);zoom xon;title(' 通過衰落信道后的波形 ');% figure;plot(noisig);axis(0 1000 -1.2 1.2)%psk 解調(diào) %pskdec=psk_decode(psk,mod);figure(10);stairs(pskdec);zoom xon;axis(0 750 -1 2);% 循環(huán)碼解碼 %xm=XH_D(pskdec);figure(11);stairs(xm);zoom xon;axis(0 750 -1 2);% 二進制轉(zhuǎn)為三進

49、制 %dcode=ZH_D(xm);figure(12);stairs(dcode);zoom xon;axis(0 750 -1 2);%pst 解碼 %pstdec=repst(dcode);figure(13);stairs(pstdec);zoom xon;axis(0 750 -1 2);%pcm 解調(diào) %pcmdec=pcm_de(pstdec,sMax); plot(pcmdec); title(' 還原的正弦信號 '); %PCM 調(diào)制 %判斷 signal 的正負%求 signal 的最大值%歸一化量化%存儲矩陣(全零)function code_lin,sM

50、ax=pcm(signal) sSig=sign(signal); sMax=max(abs(signal);Q=abs(signal/sMax)*2048; code=zeros(length(signal),8);%段落碼判斷 %for m=1:length(signal)if Q(m)>128 && Q(m)<=2048 %在第五段與第八段之間，段位碼第一位都為“ 1”code(m,2)=1;endif (Q(m)>32 && Q(m)<128) | (Q(m)>512 && Q(m)<=2048)cod

51、e(m,3)=1;%在第三四七八段內(nèi)，段位碼第二位為“ 1endif (Q(m)>16&&Q(m)<32)|(Q(m)>64&&Q(m)<128)|(Q(m)>256&&Q(m)<512)|(Q(m)>1024&& Q(m)<=2048)code(m,4)=1;%在二四六八段內(nèi)，段位碼第三位為“ 1”end%符號位的判斷if sSig(m)>0code(m,1)=1;elseif sSig(m)<0code(m,1)=0;endend%段內(nèi)碼判斷程序N=zeros(len

52、gth(signal),1);for m=1:length(signal)N(m)=code(m,2)*4+code(m,3)*2+code(m,4)+1;% 找到 Y 位于第幾段enda=0,16,32,64,128,256,512,1024;%量化間隔b=1,1,2,4,8,16,32,64;%除以 16，得到每段的最小量化間隔for m=1:length(signal)q=ceil(Q(m)-a(N(m)/b(N(m);if q=0code(m,(5:8)=0,0,0,0;else k=num2str(dec2bin(q-1,4);code(m,5)=str2num(k(1);code(m,6)=str2num(k(2);code(m,7)=str2num(k(3);code(m,8)=str2num(k(4);end%求出在段內(nèi)的位置%如果輸入為零則輸出“ 0”%編碼段內(nèi)碼為二進制end%將 N 行 8 列矩陣轉(zhuǎn)換為 1行 8*N 列的矩陣 code_lin=zeros(size(code,1)*8,1);c=1;for i=1:size(code,1)for j=1:8code_lin(c)=code(i,j); c=c+1;endendcode_lin=code_lin' %code_lin End%PST 編碼 % functiony=pst_c(