通信系統(tǒng)學(xué)科基礎(chǔ)課群綜合訓(xùn)練

1、精選優(yōu)質(zhì)文檔-傾情為你奉上課程設(shè)計任務(wù)書學(xué)生姓名： 專業(yè)班級： 指導(dǎo)教師： 工作單位： 題 目: 通信系統(tǒng)學(xué)科基礎(chǔ)課群綜合訓(xùn)練初始條件：PC機(jī)一臺、Matlab軟件要求完成的主要任務(wù): 掌握以下各種電路與通信技術(shù)的基本原理，掌握實驗的設(shè)計、電路調(diào)試與測量的方法。1.培養(yǎng)學(xué)生根據(jù)需要選學(xué)參考書，查閱手冊，圖表和文獻(xiàn)資料的自學(xué)能力，通過獨(dú)立思考深入鉆研有關(guān)問題，學(xué)會自己分析解決問題的方法。2.通過對實驗電路的分析計算，了解簡單實用電路的分析方法和工程設(shè)計方法。3.掌握示波器，頻譜儀，失真度儀的正確使用方法，學(xué)會簡單電路的實驗調(diào)試和整機(jī)指標(biāo)測試方法，提高動手能力。設(shè)計內(nèi)容：模擬信源(模擬話音)數(shù)字化

2、方式基帶碼信道碼調(diào)制方式信道解調(diào)信道解碼基帶解碼數(shù)模轉(zhuǎn)換自己構(gòu)造一時間函數(shù)PCMMiller碼漢明碼ASKAWGN與發(fā)送端對應(yīng)時間安排：1根據(jù)設(shè)計任務(wù)，分析電路原理，確定實驗方案2天2根據(jù)實驗條件進(jìn)行電路的測試，并對結(jié)果進(jìn)行分析7天3撰寫課程設(shè)計報告1天指導(dǎo)教師簽名： 年 月 日 系主任（或責(zé)任教師）簽名： 年 月 日專心-專注-專業(yè)目錄摘 要通信的目的是傳輸信息。通信系統(tǒng)的作用就是將信息從信源發(fā)送到一個或多個目的地。對于電通信來說，首先要把消息轉(zhuǎn)變成電信號，然后經(jīng)過發(fā)送設(shè)備，將信號送入信道，在接收端利用接收設(shè)備對接收信號作相應(yīng)的處理后，送給信宿再轉(zhuǎn)換為原來的消息。數(shù)字通信系統(tǒng)是利用數(shù)字信號來

3、傳遞信息的通信系統(tǒng)，數(shù)字通信涉及的技術(shù)問題很多，其中主要有信源編碼與譯碼、信道編碼與譯碼、數(shù)字調(diào)制與解調(diào)、同步以及加密與解密等。本次課程設(shè)計主要是利用Matlab仿真軟件設(shè)計完成一個典型的通信系統(tǒng)。按照此次課程設(shè)計的要求，將一模擬話音信號經(jīng)過數(shù)字化，信源編碼，信道編碼，數(shù)字調(diào)制后再經(jīng)過相應(yīng)的解碼調(diào)制后，得到原始信號。其中數(shù)字化方式為PCM編碼，基帶碼為Miller碼，信道碼為漢明碼，數(shù)字調(diào)制方式為ASK調(diào)制，信道為AWGN信道。按照這一流程利用Matlab的腳本語句編寫出各系統(tǒng)模塊，然后按照逆過程實現(xiàn)信號的解調(diào)與解碼，最后對系統(tǒng)進(jìn)行分析。關(guān)鍵字：通信系統(tǒng)，ASK仿真， PCM仿真，Miller

4、碼仿真，漢明碼仿真AbstractThe purpose of communication is the transmission of information. The role of communication systems that send information from the source to one or more destinations. For electrical communication, the message should first converted into electrical signals, and then through the transm

5、ission device, the signal is fed into the channel, the receiver uses the receiving apparatus corresponding to the received signal processing, and then sent to the sink into the original message .The course is primarily designed using Matlab simulation software designed to complete a typical communic

6、ation system. After the course in accordance with the design requirements, an analog voice signal will be digitized, source coding, channel coding, digital modulation and then through the corresponding decoding modulation, the original signal. Way in which the digital PCM encoded baseband code Mille

7、r codes, Hamming code channel, digital modulation of ASK modulation, the channel is AWGN. In accordance with this process to use Matlab script statement to write each system module, and then follow the reverse process to achieve demodulation and decoding signals the end of the system for analysis.Ke

8、y words：Communication systems, ASK simulation, PCM simulation, Miller code simulation，Hamming code simulation1 課程設(shè)計要求要求：掌握以上各種電路與通信技術(shù)的基本原理，掌握實驗的設(shè)計、電路調(diào)試與測量的方法。1.培養(yǎng)學(xué)生根據(jù)需要選學(xué)參考書，查閱手冊，圖表和文獻(xiàn)資料的自學(xué)能力，通過獨(dú)立思考深入鉆研有關(guān)問題，學(xué)會自己分析解決問題的方法。2.通過對實驗電路的分析計算，了解簡單實用電路的分析方法和工程設(shè)計方法。3.掌握示波器，頻譜儀，失真度儀的正確使用方法，學(xué)會簡單電路的實驗調(diào)試和整機(jī)指標(biāo)測試方

9、法，提高動手能力。2 數(shù)字通信系統(tǒng)模型數(shù)字通信系統(tǒng)是利用數(shù)字信號來傳遞信息的通信系統(tǒng)，數(shù)字通信涉及的技術(shù)問題很多，其中主要有信源編碼與譯碼、信道編碼與譯碼、數(shù)字調(diào)制與解調(diào)、同步以及加密與解密等。1）信源編碼與譯碼信源編碼有兩個基本功能：一是提高信息傳輸?shù)挠行裕赐ㄟ^某種數(shù)據(jù)壓縮技術(shù)設(shè)法減少碼元數(shù)目和降低碼元速率。碼元速率決定傳輸所占的帶寬，而傳輸帶寬反映了通信的有效性。二是完成模/數(shù)轉(zhuǎn)換，即當(dāng)信息源給出的是模擬信號時，信源編碼器將其轉(zhuǎn)換成數(shù)字信號，以實現(xiàn)模擬信號的數(shù)字話傳輸。信源譯碼是信源編碼的逆過程。2）信道編碼與譯碼信道編碼的目的是增強(qiáng)數(shù)字信號的抗干擾能力。數(shù)字信號在信道傳輸時受到噪聲等

10、影響后將會引起差錯。為了減小差錯，信道編碼器對傳輸?shù)男畔⒋a元按一定的規(guī)則加入保護(hù)成分（監(jiān)督元），組成所謂“抗干擾碼”。接收端的信道譯碼器按相應(yīng)的逆規(guī)則進(jìn)行解碼，從中發(fā)現(xiàn)錯誤或糾正錯誤，提高通信系統(tǒng)的可靠性。3）加密與解密在需要實現(xiàn)保密通信的場合，為了保證所傳信息的安全，人為地將被傳輸?shù)臄?shù)字序列擾亂，即加上密碼，這種處理過程叫加密。在接收端利用與發(fā)送端相同的密碼復(fù)制品對收到的數(shù)字序列進(jìn)行解密，恢復(fù)原來信息。4）數(shù)字調(diào)制與解調(diào)數(shù)字調(diào)制就是把數(shù)字基帶信號的頻譜搬移到高頻處，形成適合在信道中傳輸?shù)膸ㄐ盘?。基本的?shù)字調(diào)制方式有振幅鍵控、頻移鍵控、絕對相移鍵控、相對相移鍵控。在接收端可以采用相干解調(diào)或者

11、非相干解調(diào)還原數(shù)字基帶信號。對高斯噪聲下的信號檢測，一般采用相關(guān)器或匹配濾波器來實現(xiàn)。5）同步同步是使收發(fā)兩端的信號在時間上保持一致，是保證數(shù)字通信系統(tǒng)有序、準(zhǔn)確、可靠工作的前提條件。按照同步的功用不同，分為載波同步、位同步、群同步和網(wǎng)同步。數(shù)字通信系統(tǒng)的特點有：（1）抗干擾能力強(qiáng)，且噪聲不積累（2）傳輸差錯可控（3）便于用現(xiàn)代數(shù)字信號處理技術(shù)對數(shù)字信息進(jìn)行處理、變換、存儲。（4）利于集成，是通信設(shè)備微型化，重量輕。（5）易于加密處理，且保密性好。3 通信系統(tǒng)原理3.1 PCM編譯碼原理3.1.1 PCM概述脈沖編碼調(diào)制（PCM）簡稱脈碼調(diào)制，它是一種用二進(jìn)制數(shù)字代碼來代替連續(xù)信號的抽樣值，從

12、而實現(xiàn)通信的方式。由于這種通信方式抗干擾能力強(qiáng)，因此在光釬通信、數(shù)字微波通信、衛(wèi)星通信中均獲得了極為廣泛的運(yùn)用。PCM信號的形成是模擬信號經(jīng)過“抽樣、量化、編碼”三個步驟實現(xiàn)的。分別完成時間上離散、幅度上離散、及量化信號的二進(jìn)制表示。根據(jù)CCITT的建議，為改善小信號量化性能，采用壓擴(kuò)非均勻量化，有兩種建議方式，分別為A律和律方式，我國采用了A律方式，由于A律壓縮實現(xiàn)復(fù)雜，常使用 13 折線法編碼。3.1.2 抽樣所謂抽樣，就是對模擬信號進(jìn)行周期性掃描，把時間上連續(xù)的信號變成時間上離散的信號。該模擬信號經(jīng)過抽樣后還應(yīng)當(dāng)包含原信號中所有信息，也就是說能無失真的恢復(fù)原模擬信號。它的抽樣速率的下限是

13、由抽樣定理確定的。在一個頻帶限制在內(nèi)的時間連續(xù)信號，如果以的時間間隔對它進(jìn)行抽樣，那么根據(jù)這些抽樣值就能完全恢復(fù)原信號?；蛘哒f，如果一個連續(xù)信號的頻譜中最高頻率不超過，當(dāng)抽樣頻率時，抽樣后的信號就包含原連續(xù)的全部信息。這就是抽樣定理。3.1.3 量化方式從數(shù)學(xué)上來看，量化就是把一個連續(xù)幅度值的無限數(shù)集合映射成一個離散幅度值的有限數(shù)集合。量化器Q輸出L個量化值，k=1，2，3，L。常稱為重建電平或量化電平。當(dāng)量化器輸入信號幅度落在與之間時，量化器輸出電平為。這個量化過程可以表達(dá)為： （式3-1）這里稱為分層電平或判決閾值。通常稱為量化間隔。量化后的抽樣信號于量化前的抽樣信號相比較，當(dāng)然有所失真，

14、且不再是模擬信號。這種失真在接收端還原模擬信號是變現(xiàn)為噪聲，并稱為量化噪聲。量化噪聲的大小取決于把樣值分級“取整”的方式，分的級數(shù)越多，即量化極差或間隔越小，量化噪聲也越小。模擬信號的量化分為均勻量化和非均勻量化。由于均勻量化存在的主要缺點是：無論抽樣值大小如何，量化噪聲的均方根值都固定不變。因此，當(dāng)信號較小時，則信號量化噪聲功率比也就很小，這樣的話化信噪比就難以達(dá)到給定的要求。通常，把滿足信噪比要求的輸入信號取值范圍定義為動態(tài)范圍，可見，對于弱信號時，均勻量化時的信號動態(tài)范圍將受到較大的限制。為了克服這個缺點，實際中，往往采用非均勻量化。非均勻量化是根據(jù)信號的不同區(qū)間來確定量化間隔的。對于信

15、號取值小的區(qū)間，其量化間隔也??；反之，量化間隔就大。它與均勻量化相比，有兩個突出的優(yōu)點。首先，當(dāng)輸入量化器的信號具有非均勻分布的概率密度（實際中常常是這樣）時，非均勻量化器的輸出端可以得到較高的平均信號量化噪聲功率比；其次，非均勻量化時，量化噪聲功率的均方根值基本上與信號抽樣值成比例。因此量化噪聲對大、小信號的影響大致相同，即改善了小信號時的量化信噪比。實際中，非均勻量化的實際方法通常是將抽樣值通過壓縮再進(jìn)行均勻量化。通常使用的壓縮器中，大多采用對數(shù)式壓縮。廣泛采用的兩種對數(shù)壓縮律是 壓縮律和A壓縮律。美國采用壓縮律，我國和歐洲各國均采用A壓縮律，因此，PCM編碼方式采用的也是A壓縮律。所謂A

16、壓縮律也就是壓縮器具有如下特性的壓縮律： （式3-2）式中，x為歸一化輸入，y為歸一化輸出，A、為壓縮系數(shù)。A律壓擴(kuò)特性是連續(xù)曲線，A值不同壓擴(kuò)特性亦不同，在電路上實現(xiàn)這樣的函數(shù)規(guī)律是相當(dāng)復(fù)雜的。實際中，往往都采用近似于A律函數(shù)規(guī)律的13折線（A=87.6）的壓擴(kuò)特性。這樣，它基本上保持了連續(xù)壓擴(kuò)特性曲線的優(yōu)點，又便于用電路實現(xiàn)，本設(shè)計中所用到的PCM編碼正是采用這種壓擴(kuò)特性來進(jìn)行編碼的。實際中A律常采用13折線近似如圖2-2所示圖3-1區(qū)間0,1/16一分為二，中點為1/32,取區(qū)間1/32,1/16作為第四段; 區(qū)間0,1/32一分為二，中點為1/64,取區(qū)間1/64,1/32作為第三段;

17、區(qū)間0,1/64一分為二，中點為1/128,區(qū)間1/128,1/64作為第二段;區(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軸對應(yīng)。3.1.4 PCM編碼所謂編碼就是把量化后的信號變換成代碼，其相反的過程稱為譯碼。當(dāng)然，這里的編碼和譯碼與差錯控制編碼和譯碼是完全不同的，前者是屬于信源編碼的范疇。量化后的抽樣信號在一定的取值范圍內(nèi)僅有有限個可取的樣值，且信號正、負(fù)幅度分布的對稱性使正、負(fù)樣值的個數(shù)相等，正、負(fù)向

18、的量化級對稱分布。若將有限個量化樣值的絕對值從小到大依次排列，并對應(yīng)的依次賦予一個十進(jìn)制數(shù)字代碼，在碼前以“+”、“”號為前綴，來區(qū)分樣值的正負(fù)，則量化后的抽樣信號就轉(zhuǎn)化為按抽樣時序排列的一串十進(jìn)制數(shù)字碼流，即十進(jìn)制數(shù)字信號。把量化的抽樣信號變換成給定字長的二進(jìn)制碼流的過程為編碼。目前國際上普遍采用8位非線性編碼。例如PCM 30/32路終端機(jī)中最大輸入信號幅度對應(yīng)4 096個量化單位(最小的量化間隔稱為一個量化單位)在4 096單位的輸入幅度范圍內(nèi)，被分成256個量化級，因此須用8位碼表示每一個量化級。用于13折線A律特性的8位非線性編碼的碼組結(jié)構(gòu)如表2-1所示：表3-1 8位非線性編碼的碼

19、組結(jié)構(gòu)極性碼段落碼段內(nèi)碼M1M2M3M4M5M6M7M8其中，第1位碼M1的數(shù)值“1”或“0”分別代表信號的正、負(fù)極性，稱為極性碼。從折疊二進(jìn)制碼的規(guī)律可知，對于兩個極性不同，但絕對值相同的樣值脈沖，用折疊碼表示時，除極性碼M1不同外，其余幾位碼是完全一樣的。因此在編碼過程中，只要將樣值脈沖的極性判出后，編碼器便是以樣值脈沖的絕對值進(jìn)行量化和輸出碼組的。這樣只要考慮13折線中對應(yīng)于正輸入信號的8段折線就行了。用M2至M4位表示段落碼，它的8種可能狀態(tài)來分別代表8個段落的起點電平。其它四位表示段內(nèi)碼，它的16種可能狀態(tài)來分別代表每一段落的16個均勻劃分的量化級。這樣處理的結(jié)果，8個段落被劃分成1

20、28個量化級。3.2 Miller編譯碼原理Miller碼也稱延遲調(diào)制碼，是一種變形雙向碼。其編碼規(guī)則為：對原始符號“1”碼元起始不躍變，中心點出現(xiàn)躍變來表示，即用10或01表示。對原始符號“0”則分成單個“0”還是連續(xù)“0”予以不同處理；單個“0”時，保持0前的電平不變，即在碼元邊界處電平不躍變，在碼元中間點電平也不躍變；對于連續(xù)“0”，則使連續(xù)兩個“0”的邊界處發(fā)生電平躍變。Miller碼具有這樣的性質(zhì)：若兩個“1”碼中間有一個“0”碼時，Miller碼流中會出現(xiàn)最大寬度為兩倍碼元周期的波形，根據(jù)這一性質(zhì)可以用來進(jìn)行宏觀檢錯。3.3 漢明碼3.3.1 漢明碼編譯碼原理一般來說，若漢明碼長為

21、n，信息位數(shù)為k，則監(jiān)督位數(shù)r=n-k。若希望用r個監(jiān)督位構(gòu)造出r個監(jiān)督關(guān)系式來指示一位錯碼的n種可能位置，則要求或 （式3-3）下面以（7，4）漢明碼為例說明原理： 設(shè)漢明碼（n,k）中k=4，為了糾正一位錯碼，由式（3-1）可知，要求監(jiān)督位數(shù)r3。若取r=3,則n=k+r=7。我們用來表示這7個碼元，用的值表示3個監(jiān)督關(guān)系式中的校正子，則的值與錯誤碼元位置的對應(yīng)關(guān)系可以規(guī)定如表3-1所列。表3-2 校正子和錯碼位置的關(guān)系錯碼位置錯碼位置001101010110 100111011000無錯碼則由表1可得監(jiān)督關(guān)系式： （式3-4）邊緣檢測的實質(zhì)是采用某種算法來提取出圖像中對象與背景間的交界線

22、。我們將邊緣定義為圖像中灰度發(fā)生急劇變化的區(qū)域邊界。圖像灰度的變化情況可以在發(fā)送端編碼時，信息位的值決定于輸入信號，因此它們是隨機(jī)的。監(jiān)督位、應(yīng)根據(jù)信息位的取值按監(jiān)督關(guān)系來確定，即監(jiān)督位應(yīng)使（式3-4）中、的值為0（表示編成的碼組中應(yīng)無錯碼） （式3-5）式（4-3）經(jīng)過移項運(yùn)算，解出監(jiān)督位 （式3-6）式（4-4）其等價形式為： （式3-7）式（3-7）還可以簡記為 或 （式3-8）其中 （式3-9） （式3-10） （式3-11） （式3-12） （式3-13）所以有 （式3-14）式（3-14）等價于 （式3-15）其中Q為P的轉(zhuǎn)置，即 （式3-16）式（3-9）表示,信息位給定后，用信

23、息位的行矩陣乘矩陣Q就產(chǎn)生出監(jiān)督位。我們將Q的左邊加上一個k×k階單位方陣，就構(gòu)成一個矩陣G （式3-17）G稱為生成矩陣，因為由它可以產(chǎn)生整個碼組，即有 （式3-18）或者 （式3-19）式(3-19)即漢明碼的編碼原理。3.3.2 解碼和糾錯原理當(dāng)數(shù)字信號編碼成漢明碼形式（本文中即A）后在信道中傳輸，由于信道中噪聲的干擾，可能由于干擾引入差錯，使得接收端收到錯碼，因此在接收端進(jìn)行漢明碼糾錯，以提高通信系統(tǒng)的抗干擾能力及可靠性。 一般來說接收碼組與A不一定相同。若設(shè)接收碼組為一n列的行矩陣B，即 （式3-20）則發(fā)送碼組和接收碼組之差為 （式3-21）E就是傳輸中產(chǎn)生的錯碼行矩陣

24、（式3-22）若ei=0，表示接收碼元無錯誤，若ei=1，則表示該接收碼元有錯。式（3-21）可改寫成 （式3-23）若E=0，即接收碼組無錯,則 ，將它代人式（3-19），該是仍成立，即有 （式3-24）當(dāng)接收碼組有錯時，E0，將B帶入式（3-8）后，該式不一定成立。在未超過檢錯能力時，式（3-24）不成立。假設(shè)此時式（3-24）的右端為S,即 （式3-25）將代入式（3-25），可得 （式3-26）由式（3-6）可知，所以 （式3-25）此處S與前面的有著一一對應(yīng)關(guān)系，則S能代表錯碼位置。因此，糾錯原理即，接收端收到碼組后按式（3-25）計算出S,再根據(jù)表1判斷錯碼情況，進(jìn)行差錯糾正。3.

25、4 二進(jìn)制幅移鍵控（2ASK）調(diào)制與解調(diào)原理幅移鍵控（ASK）相當(dāng)于模擬信號中的調(diào)幅，只不過與載頻信號相乘的是二進(jìn)數(shù)碼而已。幅移就是把頻率、相位作為常量，而把振幅作為變量，信息比特是通過載波的幅度來傳遞的。由于調(diào)制信號只有0或1兩個電平，相乘的結(jié)果相當(dāng)于將載頻或者關(guān)斷，或者接通，它的實際意義是當(dāng)調(diào)制的數(shù)字信號"1時，傳輸載波；當(dāng)調(diào)制的數(shù)字信號為"0"時，不傳輸載波。 圖3-3 2ASK調(diào)制波形圖由圖可以看出2ASK信號的時間波形e2ASK（t）隨二進(jìn)制基帶信號s（t）通斷變化。所以又被稱為通斷鍵控信號圖3-4 2ASK解調(diào)波形圖3.5 AWGN信道加性高斯白噪聲A

26、WGN(Additive White Gaussian Noise)是最基本的噪聲與干擾模型。它的幅度分布服從高斯分布，而功率譜密度是均勻分布的。加性高斯白噪聲 AWGN(Additive White Gaussian Noise) 是最基本的噪聲與干擾模型。加性噪聲：疊加在信號上的一種噪聲，通常記為n(t)，而且無論有無信號，噪聲n(t)都是始終存在的。因此通常稱它為加性噪聲或者加性干擾。白噪聲：噪聲的功率譜密度在所有的頻率上均為一常數(shù)，則稱這樣的噪聲為白噪聲。如果白噪聲取值的概率分布服從高斯分布，則稱這樣的噪聲為高斯白噪聲。AWGN，在通信上指的是一種通道模型（channel model）

27、，此通道模型唯一的信號減損是來自于寬帶（Wideband）的線性加成或是穩(wěn)定譜密度（以每赫茲瓦特的帶寬表示）與高斯分布振幅的白噪聲。白噪聲是指功率譜密度在整個頻域內(nèi)均勻分布的噪聲，即其功率譜密度為常數(shù)。AWGN從統(tǒng)計上而言是隨機(jī)無線噪聲，其特點是其通信信道上的信號分布在很寬的頻帶范圍內(nèi)。高斯白噪聲的概念：“白”指功率譜恒定；高斯指幅度取各種值時的概率p(x)是高斯函數(shù)。功率譜密度恒定的話，自相關(guān)系數(shù)則是功率譜密度的反變換，高斯白噪聲的自相關(guān)系數(shù)為無延時的沖擊函數(shù)，則在時間差不等于零的時候，自相關(guān)等于0，也就是不同時間的高斯白噪聲的幅度是不相關(guān)的。4 系統(tǒng)仿真分析原信號為單一頻率400Hz信號，

28、是一般定義的話音信號頻率的最低值。圖4-1 原始信號圖4-2 PCM編碼信號圖4-3 Miller編碼信號圖4-4 Hamming編碼信號圖4-5 ASK調(diào)制信號圖4-6 經(jīng)過AWGN信道后的信號圖4-7 ASK解調(diào)信號圖4-8 Hamming解碼信號圖4-9 Miller解碼信號圖4-10 PCM解碼信號（原始信號）圖4-11 漢明碼信號與恢復(fù)出的漢明碼信號對比圖4-12 Miller信號與恢復(fù)出的Miller信號對比圖4-13 PCM信號與恢復(fù)出的PCM信號對比圖4-14 原始信號與恢復(fù)出的原始信號對比由輸出的圖像可以看出，編碼前和解碼后的數(shù)據(jù)是一致的，可以判斷PCM編解碼，Miller編

29、解碼，漢明編解碼程序是正確的。由圖4-11，并比較PCM編碼信號、Miller碼編碼信號、漢明碼編碼信號、ASK調(diào)制信號與2ASK解調(diào)信號可看出，2ASK調(diào)制與解調(diào)代碼正確，可實現(xiàn)2ASK調(diào)制與解調(diào)。比較原始正弦信號與輸出信號波形可看出，最終在接收端能夠近似地再現(xiàn)信源，實現(xiàn)整個通信系統(tǒng)的仿真。5 總結(jié)從上大學(xué)以來，我接觸Matlab軟件已經(jīng)很多次，當(dāng)時也學(xué)了部分這個軟件基本操作并且能很熟練的使用這個軟件。通過這次本次課程設(shè)計，我更加系統(tǒng)的學(xué)習(xí)了Matlab語言，發(fā)現(xiàn)其要比C語言相對簡單多。但是Matlab語言中的函數(shù)還是比較多的，在編程的時候，一定要借助Matlab幫助。另外當(dāng)一個函數(shù)看不懂時

30、，要把它代到程序中去，單獨(dú)實現(xiàn)結(jié)果，通過結(jié)果來分析函數(shù)，象像很多函數(shù)是單憑看理論是看不懂的。另外，當(dāng)看懂一個函數(shù)時，一定要在程序中注釋，否則過一段時間就忘記了其功能。在這次課程設(shè)計中，我主要學(xué)習(xí)了如何通過MATLAB程序進(jìn)行通信系統(tǒng)的仿真，對任意的信號進(jìn)行信源編碼，信道編碼，通過信道，再經(jīng)過相應(yīng)的解碼來恢復(fù)原信號。通過仿真，更加深刻的理解了通信原理中學(xué)到的信源編碼、信道編碼、調(diào)制等的相關(guān)理論知識，對這些知識點有了進(jìn)一步的掌握，對通信系統(tǒng)的流程更有了宏觀的把握。參考文獻(xiàn)1 樊昌信等.通信原理（第6版）.國防工業(yè)出版社.20082 趙鴻圖等.通信原理MATLAB仿真教程.人民郵電出版社.20103

31、 劉保柱等.MATLAB7.0從入門到精通.人民郵電出版社.20104 陳懷琛等.MATLAB及在電子信息課程中的應(yīng)用.電子工業(yè)出版社.2007 5 郭文彬等. 通信原理基于Matlab的計算機(jī)仿真. 北京郵電大學(xué)出版 社.20066 張德豐等.MATLAB通信工程仿真.機(jī)械工業(yè)出版社.20107 西瑞克斯.無線通信的MATLAB和FPGA實現(xiàn).人民郵電出版社. 2009附錄: 程序clcclose allclear all%*%原信號產(chǎn)生 % ( y )%*fs=8000;dt=1/fs;N=21;n=0:N-1;t=n*dt;fc=400;y=sin(2*pi*fc*t);figure (

32、1)stem(t,y)axis(0 N*dt -1.2 1.2);title('原始輸入信號')xlabel('Time(S)');ylabel('Amplitude');%*%PCM編碼 % （PCMcode1）%*PCMcode1=PCMENC(y,1);PCMcode2=PCMcode1'figure (2)stairs(0:length(PCMcode2(1:N*8)-1,PCMcode2(1:N*8);axis(0 length(PCMcode2(1:N*8) -0.2 1.2);title('PCM編碼信號')

33、%*%Miller編碼 % (millercode)%*millercode=MillerENC(PCMcode1,N);figure(3)stem(millercode)title('Miller編碼信號')axis(0 length(millercode) -0.2 1.2);%*%Hamming碼編碼(7,4)% (hammingcode)%*n=7;k=4;hammingcode = encode(millercode,n,k,'hamming/binary');figure(4)stem(hammingcode)title('Hamming編碼

34、信號')axis(0 length(hammingcode) -0.2 1.2);%*%ASK調(diào)制% (askcode)%*fc=1000;Tc=1/fc;baudnum=length(hammingcode);re=8; %原信號每比特重復(fù)re次bit=; c=1:baudnum*re; for i=1:baudnum if hammingcode(i)=0 bit1=zeros(1,re); else bit1=ones(1,re); end bit=bit,bit1; endt1=linspace(1.25*10-4,baudnum*Tc,baudnum*re);c=sin(2*

35、pi*t1*fc); askcode=c.*bit; figure(5);plot(askcode);axis(0 baudnum*re -1.2 1.2);title('ASK調(diào)制信號');%*%AWGN% (awgn)%*noise=4;awgn=AWGN(askcode,noise);figure(6)stairs(awgn)title('加AWGN后信號')%*%ASK Demodulation% (askdecode)%*yn1=askcode.*c;for i=1:length(yn1) if(yn1(i)>=0.2) %抽樣判決 yn1(i)

36、=1; else yn1(i)=0; endendaskdecode=zeros(1,length(hammingcode);i1=1;for i2=1:8:length(yn1) %每隔8個點 askdecode(i1)=yn1(i2); i1=i1+1;endfigure(7);stem(askdecode)title('ASK解調(diào)后信號')axis(0 length(askdecode) -0.2 1.2);%*%Hamming Decode% (hammingdecode)%*n=7;k=4;hammingdecode = decode(askdecode,n,k,&#

37、39;hamming/binary');figure(8)stem(hammingdecode)title('Hamming解碼信號')axis(0 length(hammingdecode) -0.2 1.2);%*%Miller Decode% (millerdecode)%*millerdecode=MillerDEC(hammingdecode);figure(9)stem(millerdecode)title('Miller解碼信號')axis(0 length(millerdecode) -0.2 1.2);%*%PCM Decode% (P

38、CMdecode)%*figure(10) PCMdecode=PCMDEC(millerdecode,N);stem(t,PCMdecode)axis(0 N*dt -1.2 1.2);title('PCM解碼信號')xlabel('Time(S)');ylabel('Amplitude');%*%對比%*figure(11)subplot(2,1,1);stem(t,y)axis(0 N*dt -1.2 1.2);title('原始輸入信號')xlabel('Time(S)');ylabel('Ampl

39、itude');subplot(2,1,2);stem(t,PCMdecode)axis(0 N*dt -1.2 1.2);title('PCM解碼信號(恢復(fù)的原始信號)')xlabel('Time(S)');ylabel('Amplitude');figure(12)subplot(2,1,1);stairs(0:length(PCMcode2(1:N*8)-1,PCMcode2(1:N*8);axis(0 length(PCMcode2(1:N*8) -0.2 1.2);title('PCM編碼信號')subplot(

40、2,1,2);stairs(0:length(millerdecode)-1,millerdecode);title('Miller解碼信號(PCM碼)')axis(0 length(millerdecode) -0.2 1.2);figure(13)subplot(2,1,1);stem(millercode)title('Miller編碼信號')axis(0 length(millercode) -0.2 1.2);subplot(2,1,2);stem(hammingdecode)title('Hamming解碼信號(Miller碼)')a

41、xis(0 length(hammingdecode) -0.2 1.2);figure(14)subplot(2,1,1);stem(hammingcode)title('Hamming編碼信號')axis(0 length(hammingcode) -0.2 1.2);subplot(2,1,2);stem(askdecode)title('ASK解調(diào)后信號（Hamming碼）')axis(0 length(askdecode) -0.2 1.2);%PCM編碼程序function code=PCMENC(S,range) %S為輸入信號z=sign(S);

42、 %判斷S的正負(fù)MaxS=range; %求S的最大值 S=abs(S/MaxS); %歸一化Q=2048*S; %量化code=zeros(length(S),8); %代碼存儲矩陣（全零）%段落碼判斷程序for i=1:length(S) if (Q(i)>=128)&(Q(i)<=2048) code(i,2)=1; %在第五段與第八段之間，段位碼第一位都為“1” end if (Q(i)>32)&(Q(i)<128)|(Q(i)>=512)&(Q(i)<=2048) code(i,3)=1; %在第三四七八段內(nèi)，段位碼第二位為

43、“1” endif (Q(i)>=16)&(Q(i)<32)|(Q(i)>=64)&(Q(i)<128)|.(Q(i)>=256)&(Q(i)<512)|(Q(i)>=1024)&(Q(i)<=2048) code(i,4)=1; %在二四六八段內(nèi)，段位碼第三位為“1” endend %段內(nèi)碼判斷程序N=zeros(length(S);for i=1:length(S)N(i)=bin2dec(num2str(code(i,2:4)+1; %找到code位于第幾段enda=0,16,32,64,128,256,512,1024; %量化間隔b=1,1,2,4,8,16,32,64; %除以16，得到每段