通信系統(tǒng)課群綜合訓(xùn)練與設(shè)計(jì)

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

2、思考深入鉆研有關(guān)問(wèn)題，學(xué)會(huì)自己分析解決問(wèn)題的方法。設(shè)計(jì)內(nèi)容：模擬信源(模擬話音)數(shù)字化方式基帶碼信道碼調(diào)制方式信道解調(diào)信道解碼基帶解碼數(shù)模轉(zhuǎn)換自己構(gòu)造一時(shí)間函數(shù)PCMManchester碼漢明碼ASKAWGN與發(fā)送端對(duì)應(yīng)時(shí)間安排1根據(jù)設(shè)計(jì)任務(wù)，分析電路原理，確定實(shí)驗(yàn)方案2天2根據(jù)實(shí)驗(yàn)條件進(jìn)行電路的測(cè)試，并對(duì)結(jié)果進(jìn)行分析7天3撰寫課程設(shè)計(jì)報(bào)告1天指導(dǎo)教師簽名：年 月 日 系主任（或責(zé)任教師）簽名：年 月 日目錄TOC o 1-3 h u HYPERLINK l _Toc345883940 1 原理分析與方案論證 PAGEREF _Toc345883940 h 1 HYPERLINK l _Toc

3、345883941 1.1 通信系統(tǒng)架構(gòu) PAGEREF _Toc345883941 h 1 HYPERLINK l _Toc345883942 1.2 信源碼PCM碼 PAGEREF _Toc345883942 h 2 HYPERLINK l _Toc345883943 1.3 基帶碼Manchester碼 PAGEREF _Toc345883943 h 3 HYPERLINK l _Toc345883944 曼徹斯特碼簡(jiǎn)介及其編碼規(guī)則 PAGEREF _Toc345883944 h 3 HYPERLINK l _Toc345883945 1.3.2 曼徹斯特碼原理 PAGEREF _Toc

4、345883945 h 3 HYPERLINK l _Toc345883946 1.4 信道碼漢明碼 PAGEREF _Toc345883946 h 3 HYPERLINK l _Toc345883947 1.5 噪聲信道AWGN PAGEREF _Toc345883947 h 5 HYPERLINK l _Toc345883948 2 各模塊的MATLAB實(shí)現(xiàn) PAGEREF _Toc345883948 h 5 HYPERLINK l _Toc345883949 2.1 信號(hào)源的實(shí)現(xiàn) PAGEREF _Toc345883949 h 5 HYPERLINK l _Toc345883950 2.

5、2 信源編碼PCM編碼 PAGEREF _Toc345883950 h 6 HYPERLINK l _Toc345883951 2.2.1 PCM編碼原理 PAGEREF _Toc345883951 h 6 HYPERLINK l _Toc345883952 2.2.2 PCM編碼的實(shí)現(xiàn) PAGEREF _Toc345883952 h 8 HYPERLINK l _Toc345883953 2.3 基帶編碼Manchester編碼 PAGEREF _Toc345883953 h 9 HYPERLINK l _Toc345883954 2.4 信道編碼漢明碼編碼 PAGEREF _Toc3458

6、83954 h 10 HYPERLINK l _Toc345883955 2.4.1 漢明碼編碼分析 PAGEREF _Toc345883955 h 10 HYPERLINK l _Toc345883956 2.4.2 漢明碼的MATLAB實(shí)現(xiàn) PAGEREF _Toc345883956 h 11 HYPERLINK l _Toc345883957 2.5 ASK調(diào)制 PAGEREF _Toc345883957 h 12 HYPERLINK l _Toc345883958 2.5.1 ASK調(diào)制原理 PAGEREF _Toc345883958 h 12 HYPERLINK l _Toc3458

7、83959 2.5.2 ASK調(diào)制的MATLAB實(shí)現(xiàn) PAGEREF _Toc345883959 h 13 HYPERLINK l _Toc345883960 2.6 信道噪聲AWGN PAGEREF _Toc345883960 h 14 HYPERLINK l _Toc345883961 2.7 ASK解調(diào) PAGEREF _Toc345883961 h 14 HYPERLINK l _Toc345883962 2.8 漢明碼解碼 PAGEREF _Toc345883962 h 15 HYPERLINK l _Toc345883963 2.8.1 漢明碼校驗(yàn) PAGEREF _Toc3458

8、83963 h 15 HYPERLINK l _Toc345883964 2.8.2 漢明碼譯碼程序 PAGEREF _Toc345883964 h 16 HYPERLINK l _Toc345883965 2.9 manchester譯碼 PAGEREF _Toc345883965 h 17 HYPERLINK l _Toc345883966 2.10 PCM譯碼 PAGEREF _Toc345883966 h 17 HYPERLINK l _Toc345883967 2.11 誤碼率的計(jì)算 PAGEREF _Toc345883967 h 18 HYPERLINK l _Toc3458839

9、68 3 仿真結(jié)果分析 PAGEREF _Toc345883968 h 19 HYPERLINK l _Toc345883969 3.1 源信號(hào)與接收信號(hào)波形對(duì)比 PAGEREF _Toc345883969 h 19 HYPERLINK l _Toc345883970 3.2 誤碼率統(tǒng)計(jì) PAGEREF _Toc345883970 h 20 HYPERLINK l _Toc345883971 4 心得體會(huì) PAGEREF _Toc345883971 h 20 HYPERLINK l _Toc345883972 5 參考文獻(xiàn) PAGEREF _Toc345883972 h 21 HYPERLIN

10、K l _Toc345883973 6 附錄 PAGEREF _Toc345883973 h 22摘要 通信是通過(guò)某種媒體進(jìn)行的信息傳遞，目的是傳輸信息，通信系統(tǒng)是用以完成信息傳輸過(guò)程的技術(shù)系統(tǒng)的總稱，作用是將信息從信源發(fā)送到一個(gè)或多個(gè)目的地。調(diào)制與解調(diào)在信息的傳輸過(guò)程中占據(jù)著重要的地位，是不可或缺的，因此研究系統(tǒng)的調(diào)制和解調(diào)過(guò)程就極為重要。MATLAB是集數(shù)值計(jì)算、圖形繪制、圖像處理及系統(tǒng)仿真等強(qiáng)大功能于一體的科學(xué)計(jì)算語(yǔ)言，它強(qiáng)大的矩陣運(yùn)算和圖形可視化的功能以及豐富的工具箱，為通信系統(tǒng)的調(diào)制和解調(diào)過(guò)程的分析提供了極大的方便。 調(diào)制在通信系統(tǒng)中有十分重要的作用。通過(guò)調(diào)制，不僅可以進(jìn)行頻譜搬移，

11、把調(diào)制信號(hào)的頻譜搬移到所希望的位置上，從而將調(diào)制信號(hào)轉(zhuǎn)換成適合于傳播的已調(diào)信號(hào)，而且它對(duì)系統(tǒng)的傳輸有效性和傳輸?shù)目煽啃杂兄艽蟮挠绊懀{(diào)制方式往往決定了一個(gè)通信系統(tǒng)的性能。MATLAB軟件廣泛用于數(shù)字信號(hào)分析，系統(tǒng)識(shí)別，時(shí)序分析與建模，神經(jīng)網(wǎng)絡(luò)、動(dòng)態(tài)仿真等方面有著廣泛的應(yīng)用。本課題利用MATLAB軟件完成模擬信號(hào)進(jìn)行PCM編碼、Manchester編碼、漢明碼、ASK調(diào)制、經(jīng)過(guò)AWGN信道，再解調(diào)、譯碼的完整通信系統(tǒng)仿真，并通過(guò)統(tǒng)計(jì)誤碼率和對(duì)比前后波形，對(duì)這個(gè)通信系統(tǒng)進(jìn)行評(píng)估。關(guān)鍵字：MATLAB，編碼，調(diào)制，解調(diào)，譯碼，通信系統(tǒng)Abstract Communication is throug

12、h some of the information transmission media, the purpose is to transmit information, communication system is used to complete the process of information transmission technology system body, the function is the information from the source sent to one or more of the destination. Modulation and demodu

13、lation in information transmission process plays an important role, is indispensable, so the system of modulation and demodulation process is extremely important. MATLAB is set numerical calculation, graphics, image processing and system simulation and powerful function in the integration of scienti

14、fic computing language, its powerful matrix operation and graphics visual function and rich toolbox, as the communication system of the modulation and demodulation process analysis provides a great convenience. Modulation in a communication system has a very important role. Through the modulation, c

15、an not only move the modulation signal spectrum, and the spectrum of the move to want the position, and the modulation signal into suitable for transmission of the modulated signal, and it to the system transmission efficiency and the reliability of the transmission has a great influence on the modu

16、lation method, often determines a communication system performance. MATLAB software is widely used in digital signal analysis, system identification, time series analysis and modeling, neural networks, dynamic simulation for such applications. This topic using MATLAB software simulation signal PCM c

17、oding, Manchester coding, hamming code and ASK modulation, through AWGN channel, and demodulation, decoding the integrity of the communication system simulation, and through the statistical error rate and contrast before and after waveform, the communication system to evaluate.Key words：MATLAB，Codin

18、g, modulation and demodulation, decoding, communication system1 原理分析與方案論證1.1 通信系統(tǒng)架構(gòu)通信的目的是傳輸信息。通信系統(tǒng)的作用就是將信息從信息源發(fā)送到一個(gè)或多個(gè)目的地。對(duì)于任何一個(gè)通信系統(tǒng)，均可視為由發(fā)送端、信道和接收端三大部分組成（如圖1-1所示）。接收設(shè)備發(fā)送設(shè)備信息源編碼設(shè)備信 道信息源噪聲源發(fā)送端接收端信道調(diào)制設(shè)備譯碼碼設(shè)備解調(diào)設(shè)備圖1-1 通信系統(tǒng)一般模型1、信息源:把原始信息變換成原始電信號(hào)。2、信源編碼：實(shí)現(xiàn)模擬信號(hào)的數(shù)字化傳輸即完成A/D變化。提高信號(hào)傳輸?shù)挠行?。即在保證一定傳輸質(zhì)量的情況下，用盡可能

19、少的數(shù)字脈沖來(lái)表示信源產(chǎn)生的信息。信源編碼也稱作頻帶壓縮編碼或數(shù)據(jù)壓縮編碼。3、信道編碼：信源編碼的目的： 信道編碼主要解決數(shù)字通信的可靠性問(wèn)題。信道編碼的原理：對(duì)傳輸?shù)男畔⒋a元按一定的規(guī)則加入一些冗余碼（監(jiān)督碼），形成新的碼字，接收端按照約定好的規(guī)律進(jìn)行檢錯(cuò)甚至糾錯(cuò)。信道編碼又稱為差錯(cuò)控制編碼、抗干擾編碼、糾錯(cuò)編碼 。4、數(shù)字調(diào)制數(shù)字調(diào)制技術(shù)的概念：把數(shù)字基帶信號(hào)的頻譜搬移到高頻處，形成適合在信道中傳輸?shù)念l帶信號(hào)。數(shù)字調(diào)制的主要作用：提高信號(hào)在信道上傳輸?shù)男?，達(dá)到信號(hào)遠(yuǎn)距離傳輸?shù)哪康??；镜臄?shù)字調(diào)制方式：振幅鍵控ASK、頻移鍵控FSK、相移鍵控PSK。5、信道：信道是信號(hào)傳輸媒介的總稱，傳

20、輸信道的類型有無(wú)線信道（如電纜、光纖）和有線信道（如自由空間）兩種。6、噪聲源：通信系統(tǒng)中各種設(shè)備以及信道中所固有的，為了分析方便，把噪聲源視為各處噪聲的集中表現(xiàn)而抽象加入到信道。 信源碼PCM碼通常是把從模擬信號(hào)抽樣、量化，直到變換成為二進(jìn)制符號(hào)的基本過(guò)程，稱為脈沖編碼調(diào)制PCM，簡(jiǎn)稱脈碼調(diào)制。在編碼器中由沖激脈沖對(duì)模擬信號(hào)抽樣，得到在抽樣時(shí)刻上的信號(hào)抽樣值。這個(gè)抽樣值仍是模擬量。在量化之前，通常由保持電路將其作短暫保存，以便電路有時(shí)間對(duì)其量化。在實(shí)際電路中，常把抽樣和保持電路做在一起，稱為抽樣保持電路。圖中的量化器把模擬抽樣信號(hào)變成離散的數(shù)字量，然后在編碼器中進(jìn)行二進(jìn)制編碼。這樣每個(gè)二進(jìn)制

21、碼組就代表一個(gè)量化后的信號(hào)抽樣值。圖中的解碼器的原理和編碼過(guò)程相反。其中，量化與編碼的組合稱為模/數(shù)變換器（A/D變換器）；解碼與低通濾波的組合稱為數(shù)/模變換器（D/A變換器）。抽樣是對(duì)模擬信號(hào)進(jìn)行周期性的掃描，把時(shí)間上連續(xù)的信號(hào)變成時(shí)間上離散的信號(hào)。經(jīng)過(guò)抽樣的信號(hào)應(yīng)包含原信號(hào)的所有信息，即能無(wú)失真地恢復(fù)出原模擬信號(hào)。量化是把經(jīng)抽樣得到的瞬時(shí)值進(jìn)行幅度離散，即指定Q規(guī)定的電平，把抽樣值用最接近的電平表示。 編碼是用二進(jìn)制碼組表示有固定電平的量化值。實(shí)際上量化是在編碼過(guò)程中同時(shí)完成的。律和A律壓縮特性：律： （式1-1）A律： （式1-2）式中，x為歸一化輸入，y為歸一化輸出，A、為壓縮系數(shù)。

22、數(shù)字壓擴(kuò)技術(shù)：一種通過(guò)大量的數(shù)字電路形成若干段折線， 并用這些折線來(lái)近似A律或律壓擴(kuò)特性，從而達(dá)到壓擴(kuò)目的方法。即對(duì)數(shù)壓擴(kuò)特性的折線近似法。我國(guó)一般都采用A律。本文設(shè)計(jì)方案也采用A律-13折線法。1.3 基帶碼Manchester碼數(shù)字通信系統(tǒng)傳輸碼流中應(yīng)含有時(shí)鐘頻率分量，以便于接收端從碼流中提取同步信息，同步信息保證接收端按正確時(shí)序再生原始信息，減小誤碼率。因此信源信號(hào)進(jìn)行信源編碼調(diào)制后還需要進(jìn)行基帶編碼，使得碼流中“1”和“0”的統(tǒng)計(jì)概率應(yīng)各占1/2。曼徹斯特碼簡(jiǎn)介及其編碼規(guī)則 曼徹斯特碼 Manchester code (又稱裂相碼、雙向碼），一種用電平跳變來(lái)表示1或0的編碼，其變化規(guī)則

23、很簡(jiǎn)單，即每個(gè) HYPERLINK :/baike.baidu /view/601477.htm t _blank 碼元均用兩個(gè)不同相位的 HYPERLINK :/baike.baidu /view/1573140.htm t _blank 電平信號(hào)表示，也就是一個(gè)周期的方波，但0碼和1碼的相位正好相反。二進(jìn)制碼與曼徹斯特碼波形的對(duì)比關(guān)系如下。110 001 即把1變換為1/0序列中間下降沿代表1把0變換為0/1序列中間上升沿表示0。優(yōu)點(diǎn)在于其編碼后的時(shí)鐘信息都不受信源統(tǒng)計(jì)特性影響（連續(xù)0、連續(xù)1 時(shí)無(wú)法接受定時(shí)信號(hào)），有很強(qiáng)的定時(shí)分量。 1.3.2 曼徹斯特碼原理 用于數(shù)字基帶傳輸?shù)拇a型種類

24、Manchester碼是其中常用的一種。Manchester碼是一種用跳變沿而非電平來(lái)表示要傳輸?shù)亩M(jìn)制信息0或1一般規(guī)定在位元中間用下跳變表示“1”用上跳變表示“0”。曼徹斯特編碼被被認(rèn)為是一種自定時(shí)碼自定時(shí)意味著數(shù)據(jù)流的精確同步是可行的。每一個(gè)比特都準(zhǔn)確的在一預(yù)先定義時(shí)間時(shí)期的時(shí)間中被傳送。曼徹斯特編碼提供了一種簡(jiǎn)單的方法在長(zhǎng)時(shí)間段內(nèi)沒(méi)有電平跳變的情況下仍然能夠?qū)θ我獾亩M(jìn)制序列進(jìn)行編碼并且防止在這種情況下同步時(shí)鐘信號(hào)的丟失以及防止低通模擬電路中低頻直流飄移所引起的比特錯(cuò)誤。如果保證傳送的編碼交流信號(hào)的直流分量為零并且能夠防止中繼信號(hào)的基線漂移那么很容易實(shí)現(xiàn)信號(hào)的恢復(fù)和防止能量的浪費(fèi)。曼徹

25、斯特碼具有豐富的位定時(shí)信息。1.4 信道碼漢明碼漢明碼Hamming code：用于數(shù)據(jù)傳送，能檢測(cè)所有一位和雙位差錯(cuò)并糾正所有一位差錯(cuò)的二進(jìn)制代碼。與其他的錯(cuò)誤校驗(yàn)碼類似，漢明碼也利用了 HYPERLINK :/baike.baidu /view/492663.htm t _blank 奇偶校驗(yàn)位的概念，通過(guò)在 HYPERLINK :/baike.baidu /view/190231.htm t _blank 數(shù)據(jù)位后面增加一些比特，可以驗(yàn)證數(shù)據(jù)的有效性。利用一個(gè)以上的校驗(yàn)位，漢明碼不僅可以驗(yàn)證數(shù)據(jù)是否有效，還能在數(shù)據(jù)出錯(cuò)的情況下指明錯(cuò)誤位置。在接受端通過(guò)糾錯(cuò) HYPERLINK :/bai

26、ke.baidu /view/189742.htm t _blank 譯碼自動(dòng)糾正傳輸中的差錯(cuò)來(lái)實(shí)現(xiàn)碼糾錯(cuò)功能，稱為 HYPERLINK :/baike.baidu /view/1439048.htm t _blank 前向糾錯(cuò)FEC。在iew/1128788.htm t _blank 數(shù)據(jù)鏈路中存在大量噪音時(shí)，F(xiàn)EC可以增加數(shù)據(jù) HYPERLINK :/baike.baidu /view/4316.htm t _blank 吞吐量。通過(guò)在傳輸碼列中加入 HYPERLINK :/baike.baidu /view/7155908.htm t _blank 冗余位(也稱糾錯(cuò)位)可以實(shí)現(xiàn)43904

27、8.htm t _blank 前向糾錯(cuò)。但這種方法比簡(jiǎn)單重傳協(xié)議的成本要高。漢明碼利用奇偶?jí)K機(jī)制降低了 HYPERLINK :/baike.baidu /view/1439048.htm t _blank 前向糾錯(cuò)的成本?，F(xiàn)以數(shù)據(jù)碼1101為例講講漢明碼的編碼原理，此時(shí)D8=1、D4=1、D2=0、D1=1，在P1編碼時(shí)，先將D8、D4、D1的二進(jìn)制碼相加，結(jié)果為奇數(shù)3，漢明碼對(duì)奇數(shù)結(jié)果編碼為1，偶數(shù)結(jié)果為0（奇數(shù)位。若奇數(shù)結(jié)果編碼為0.偶數(shù)結(jié)果為1，則叫偶數(shù)位），因此P1值為1，D8+D2+D1=2，為偶數(shù)，那么P2值為0，D4+D2+D1=2，為偶數(shù)，P3值為0。這樣，參照上文的位置表，漢

28、明碼處理的結(jié)果就是1010101。在這個(gè)4位數(shù)據(jù)碼的例子中，我們可以發(fā)現(xiàn)每個(gè)漢明碼都是以三個(gè)數(shù)據(jù)碼為基準(zhǔn)進(jìn)行編碼的。從編碼形式上，我們可以發(fā)現(xiàn)漢明碼是一個(gè)校驗(yàn)很嚴(yán)謹(jǐn)?shù)木幋a方式。在這個(gè)例子中，通過(guò)對(duì)4個(gè) HYPERLINK :/baike.baidu /view/190231.htm t _blank 數(shù)據(jù)位的3個(gè)位的3次組合檢測(cè)來(lái)達(dá)到具體碼位的校驗(yàn)與修正目的（不過(guò)只允許一個(gè)位出錯(cuò)，兩個(gè)出錯(cuò)就無(wú)法檢查出來(lái)了，這從下面的糾錯(cuò)例子中就能體現(xiàn)出來(lái)）。在校驗(yàn)時(shí)則把每個(gè)漢明碼與各自對(duì)應(yīng)的 HYPERLINK :/baike.baidu /view/190231.htm t _blank 數(shù)據(jù)位值相加，如果

29、結(jié)果為偶數(shù)（糾錯(cuò)代碼為0）就是正確，如果為奇數(shù)（糾錯(cuò)代碼為1）則說(shuō)明當(dāng)前漢明碼所對(duì)應(yīng)的三個(gè)數(shù)據(jù)位中有錯(cuò)誤，此時(shí)再通過(guò)其他兩個(gè)漢明碼各自的運(yùn)算來(lái)確定具體是哪個(gè)位出了問(wèn)題。還是剛才的1101的例子，正確的編碼應(yīng)該是1010101，如果第三個(gè)數(shù)據(jù)位在傳輸途中因干擾而變成了1，就成了1010111。檢測(cè)時(shí)，P1+D8+D4+D1的結(jié)果是偶數(shù)4，第一位糾錯(cuò)代碼為0，正確。P2+D8+D2+D1的結(jié)果是奇數(shù)3，第二位糾錯(cuò)代碼為1，有錯(cuò)誤。P3+D4+D2+D1的結(jié)果是奇數(shù)3，第三位糾錯(cuò)代碼為1，有錯(cuò)誤。那么具體是哪個(gè)位有錯(cuò)誤呢？三個(gè)糾錯(cuò)代碼從高到低排列為 HYPERLINK :/baike.baidu /

30、view/802204.htm t _blank 二進(jìn)制編碼110，換算成十進(jìn)制就是6，也就是說(shuō)第6位數(shù)據(jù)錯(cuò)了，而數(shù)據(jù)第三位在漢明碼編碼后的位置正好是第6位。針對(duì)4位數(shù)據(jù)的漢明碼編碼示意圖漢明碼是一個(gè)在原有數(shù)據(jù)中插入若干校驗(yàn)碼來(lái)進(jìn)行錯(cuò)誤檢查和糾正的編碼技術(shù)。以典型的4位數(shù)據(jù)編碼為例，漢明碼將加入3個(gè)校驗(yàn)碼，從而使實(shí)際傳輸?shù)臄?shù)據(jù)位達(dá)到7個(gè)（位），它們的位置如果把上圖中的位置橫過(guò)來(lái)就是：數(shù)據(jù)位1234567代碼P1 P2 D8 P3 D4D2 D1說(shuō)明第1個(gè)漢明碼第2個(gè)漢明碼第1個(gè)數(shù)據(jù)碼第3個(gè)漢明碼第2個(gè)數(shù)據(jù)碼第3個(gè)數(shù)據(jù)碼第4個(gè)數(shù)據(jù)碼1.5 噪聲信道AWGN加性高斯白噪聲AWGN(Additive

31、 White Gaussian Noise) 是最基本的噪聲與干擾模型。加性噪聲：疊加在信號(hào)上的一種噪聲，通常記為n(t)，而且無(wú)論有無(wú)信號(hào)，噪聲n(t)都是始終存在的。因此通常稱它為加性噪聲或者加性干擾。白噪聲：噪聲的功率譜密度在所有的頻率上均為一常數(shù)，則稱這樣的噪聲為白噪聲。如果白噪聲取值的概率分布服從 HYPERLINK :/baike.baidu /view/573667.htm t _blank 高斯分布，則稱這樣的噪聲為高斯白噪聲。2 各模塊的MATLAB實(shí)現(xiàn)2.1 信號(hào)源的實(shí)現(xiàn)根據(jù)題目的要求，自己構(gòu)造一個(gè)時(shí)間函數(shù)，模擬話音信號(hào)，因?yàn)樵捯粜盘?hào)的頻率在300Hz3400Hz之間，取典

32、型值300Hz。構(gòu)造一個(gè)模擬正弦函數(shù)。設(shè)定采樣頻率為6000Hz，因而采樣點(diǎn)間步長(zhǎng)為1/6000s。程序?yàn)椋篺s=6000; %設(shè)定采樣頻率dt=1/fs; %設(shè)定步長(zhǎng)N=50; long=N; n=0:N-1;t=n*dt; %截止時(shí)間 fc=300; %時(shí)間函數(shù)頻率 y=sin(2*pi*fc*t); %構(gòu)造正弦函數(shù)figure(1); plot(t,y); %繪制時(shí)間函數(shù)圖象title(源信號(hào)時(shí)間函數(shù)波形)繪制出構(gòu)造的信號(hào)波形如圖2-1所示： 圖2-1 信號(hào)源時(shí)間函數(shù)波形2.2 信源編碼PCM編碼2. PCM編碼原理本設(shè)計(jì)采用A律折線法。實(shí)際中A律常用13折線法近似如圖2-2所示： 圖2

33、-2 A律13折線其具體分法如下：先將X軸的區(qū)間0，1一分為二，其中點(diǎn)為1/2,取區(qū)間1/2,1作為第八段;區(qū)間0,1/2再一分為二，其中點(diǎn)為1/4,取區(qū)間1/4,1/2作為第七段;區(qū)間0,1/4再一分為二，其中點(diǎn)為1/8,取區(qū)間1/8,1/4作為第六段;區(qū)間0,1/8一分為二，中點(diǎn)為1/16,取區(qū)間1/16,1/8作為第五段;區(qū)間0,1/16一分為二，中點(diǎn)為1/32,取區(qū)間1/32,1/16作為第四段; 區(qū)間0,1/32一分為二，中點(diǎn)為1/64,取區(qū)間1/64,1/32作為第三段;區(qū)間0,1/64一分為二，中點(diǎn)為1/128,區(qū)間1/128,1/64作為第二段;區(qū)間0,1/128作為第一段。然

34、后將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軸對(duì)應(yīng)。碼位的安排:目前國(guó)際上普遍采用8位非線性編碼。例如PCM 30/32路終端機(jī)中最大輸入信號(hào)幅度對(duì)應(yīng)4 096個(gè)量化單位(最小的量化間隔稱為一個(gè)量化單位)， 在4 096單位的輸入幅度范圍內(nèi)，被分成256個(gè)量化級(jí)，因此須用8位碼表示每一個(gè)量化級(jí)。用于13折線A律特性的8位非線性編碼的碼組結(jié)構(gòu)如表2-1所示：表2-1 8位非線性編碼的碼組結(jié)構(gòu)極性碼段落碼段內(nèi)碼M1M2M3M4M5M6M7M8其

35、中，第1位碼M1的數(shù)值“1”或“0”分別代表信號(hào)的正、負(fù)極性，稱為極性碼。從折疊二進(jìn)制碼的規(guī)律可知，對(duì)于兩個(gè)極性不同，但絕對(duì)值相同的樣值脈沖，用折疊碼表示時(shí)，除極性碼M1不同外，其余幾位碼是完全一樣的。因此在編碼過(guò)程中，只要將樣值脈沖的極性判出后，編碼器便是以樣值脈沖的絕對(duì)值進(jìn)行量化和輸出碼組的。這樣只要考慮13折線中對(duì)應(yīng)于正輸入信號(hào)的8段折線就行了。這8段折線共包含128個(gè)量化級(jí)，正好用剩下的7位碼(M2，, M8)就能表示出來(lái)。2. PCM編碼的實(shí)現(xiàn)設(shè)計(jì)處理段落碼子程序mdlGetDuanLuoMa，功能是對(duì)采樣得到的信號(hào)進(jìn)行比較，確定其PCM編碼對(duì)應(yīng)的段落碼。function duanL

36、uoMa= mdlGetDuanLuoMa( table, Is, i)switch i,case 2,duanLuoMa= ( Is= table( 4, 2) ) ;case 3,duanLuoMa= ( Is=table( 2, 2) ) |( Is= table( 6, 2) ) ;case 4,duanLuoMa= ( Is=table( 2, 2) & Is= table( 3, 2) ) | ( Is table ( 4, 2)& Is=table( 1, 2) | ( Is=table( 4, 2)&Is=table( 5, 2) | ( Is=table ( 4,2) &Is

37、 =table( 6, 2)& Is=table( 7,2);endend設(shè)計(jì)子程序pcmcoder，計(jì)算出相應(yīng)段落碼區(qū)間里的段內(nèi)號(hào)以及極性碼，完成PCM編碼。編碼結(jié)果如圖2-3所示： 圖2-3 PCM編碼結(jié)果2.3 基帶編碼Manchester編碼編寫實(shí)現(xiàn)Manchester編碼子程序，即對(duì)PCM編碼后的信號(hào)進(jìn)行處理，遇1則編為10，遇0則編為01即可。function mcst=mcstcoder(y,long)mcst=zeros(1,16*long);for n=1:8*long if(y(n)=1) mcst(2*n-1)=1; mcst(2*n)=0; else if(y(n)=0

38、) mcst(2*n-1)=0; mcst(2*n)=1; end endendendManchester編碼結(jié)果如圖2-4所示： 圖2-4 Manchester編碼2.4 信道編碼漢明碼編碼2. 漢明碼編碼分析現(xiàn)在以(7,4)分組碼為例來(lái)說(shuō)明漢明碼的特點(diǎn)。其主要參數(shù)如下：碼長(zhǎng)：信息位：校驗(yàn)位：，且最小距離：其生成矩陣G（前四位為信息位，后三位為冗余位）如下： （式2-1）系統(tǒng)碼可分為消息部分和冗余部分兩部分，根據(jù)生成矩陣，輸出碼字可按下 式計(jì)算： （式2-2）所以有信息位冗余位由以上關(guān)系可以得到(7,4)漢明碼的全部碼字如表2-2所示。表2-2 (7,4)漢明碼的全部碼字序號(hào)信息碼元冗余元序

39、號(hào)信息碼元冗余元000000008100011110001011910011002001010110101001030011110111011001401001101211000015010110113110101060110011141110100701110001511111112. 漢明碼的MATLAB實(shí)現(xiàn)function h= hanmincoder(s,L)if mod(L,4)=1; L=L+3; s(L)=0;s(L-1)=0;s(L-2)=0;elseif mod(L,4)=2; L=L+2; s(L)=0;s(L-1)=0; elseif mod(L,4)=3; L=L+1;

40、 s(L)=0; endh=zeros(1,L+L/4*3); N=L+L/4*3;j=1;for k=0:L/4-1 i=4*k+1; h(j)=s(i); h(j+1)=s(i+1); h(j+2)=s(i+2); h(j+3)=s(i+3); h(j+4)=xor(xor(s(i),s(i+1),s(i+2); h(j+5)=xor(xor(s(i),s(i+1),s(i+3); h(j+6)=xor(xor(s(i),s(i+2),s(i+3); j=j+7;endend由程序可見(jiàn)，程序先將輸入的信號(hào)序列補(bǔ)0至4的倍數(shù)，然后每4個(gè)相鄰數(shù)據(jù)作為一組，進(jìn)行處理，h序列即為編碼后的序列，h（

41、j）、h（j+1）、h（j+2）、h（j+3）為 信息位，h（j+4）、h（j+5）、h（j+6）為冗余位。編碼后的結(jié)果如圖2-5所示： 圖2-5 漢明碼編碼2.5 ASK調(diào)制2. ASK調(diào)制原理理論上數(shù)字調(diào)制與模擬調(diào)制在本質(zhì)上沒(méi)有什么不同，它們都屬于正弦波調(diào)制。但是，數(shù)字調(diào)制是源信號(hào)為離散型的正弦波調(diào)制，而模擬調(diào)制則是源信號(hào)為連續(xù)型的正弦波調(diào)制，因而，數(shù)字調(diào)制具有由數(shù)字信號(hào)帶來(lái)的一些特點(diǎn)。這些特點(diǎn)主要包括兩個(gè)方面：第一，數(shù)字調(diào)制信號(hào)的產(chǎn)生，除把數(shù)字的調(diào)制信號(hào)當(dāng)作模擬信號(hào)的特例而直接采用模擬調(diào)制方式產(chǎn)生數(shù)字調(diào)制信號(hào)外，還可以采用鍵控載波的方法。第二，對(duì)于數(shù)字調(diào)制信號(hào)的解調(diào)，為提高系統(tǒng)的抗噪聲性

42、能，通常采用與模擬調(diào)制系統(tǒng)中不同的解調(diào)方式。振幅鍵控是正弦載波的幅度隨數(shù)字基帶信號(hào)而變化的數(shù)字調(diào)制，即源信號(hào)為“1”時(shí)，發(fā)送載波，源信號(hào)為“0”時(shí)，發(fā)送0電平。 所以也稱這種調(diào)制為通、斷鍵控（OOK）。當(dāng)數(shù)字基帶信號(hào)為二進(jìn)制時(shí)，也稱為二進(jìn)制振幅鍵控（2ASK），2ASK信號(hào)的調(diào)制方法有模擬幅度調(diào)制方法和鍵控方法兩種。2ASK信號(hào)是數(shù)字調(diào)制方式中最早出現(xiàn)的，也是最簡(jiǎn)單的，但其抗噪聲性能較差，因此實(shí)際應(yīng)用并不廣泛，但經(jīng)常作為研究其它數(shù)字調(diào)制方式的基礎(chǔ)。 2ASK的時(shí)域特征 2ASK信號(hào)的時(shí)域表示式為： 其中， 為隨機(jī)的單極性矩形脈沖序列， 是經(jīng)過(guò)基帶成型處理之后的脈沖序列。2ASK信號(hào)的時(shí)域波形如

43、圖2-6所示： 圖2-6 2ASK時(shí)域波形2. ASK調(diào)制的MATLAB實(shí)現(xiàn)先將漢明碼序列進(jìn)行擴(kuò)展，每個(gè)數(shù)據(jù)都重復(fù)8次然后重組，與對(duì)應(yīng)的余弦載波進(jìn)行相乘，得到調(diào)制后的信號(hào)。調(diào)制后的部分波形圖如圖2-7所示 圖2-7 ASK調(diào)制信號(hào)2.6 信道噪聲AWGNMatlab本身自帶了于某信號(hào)中加入高斯白噪聲的函數(shù)：AWGNy = awgn(x,SNR) 在信號(hào)x中加入高斯白噪聲。SNR為信噪比，以dB為單位。x的強(qiáng)度假定為0dBW。設(shè)定SNR信噪比為13dB。繪制出收到噪聲干擾的傳輸信號(hào)波形圖如圖2-8所示： 圖2-8 加入高斯白噪聲的調(diào)制波2.7 ASK解調(diào)本設(shè)計(jì)采用相干解調(diào)，即采用與載波信號(hào)同步的

44、信號(hào)與已調(diào)信號(hào)相乘，再通過(guò)FIR低通濾波器，即可解調(diào)出調(diào)制信號(hào)。解調(diào)框圖如圖2-9所示：FIR濾波器乘法器raskz頻率為fc的相干波car 圖2-9 ASK解調(diào)框圖解調(diào)后的波形圖如圖2-10所示： 圖2-10 ASK相干解調(diào)后波形2.8 漢明碼解碼2.8.1 漢明碼校驗(yàn) 在發(fā)送端信息碼元M利用式（8-16），實(shí)現(xiàn)信道編碼，產(chǎn)生線性分組碼A；在傳輸過(guò)程中有可能出現(xiàn)誤碼，設(shè)接收到的碼組為B。則收發(fā)碼組之差為：（式2-3） 這里，表示i位有錯(cuò)，表示i位無(wú)錯(cuò)。基于這樣的原則接收端利用接收到的碼組B計(jì)算校正子：（式2-4） 因此，校正子僅與E有關(guān)，即錯(cuò)誤圖樣與校正子之間有確定的關(guān)系。對(duì)于上述（7，4）

45、碼，校正子S與錯(cuò)誤圖樣的對(duì)應(yīng)關(guān)系可由式（2-4）求得，其計(jì)算結(jié)果見(jiàn)表2-3所示。在接收端的譯碼器中有專門的校正子計(jì)算電路，從而實(shí)現(xiàn)檢錯(cuò)和糾錯(cuò)。 表2-3 （7，4）碼校正子與錯(cuò)誤圖樣的對(duì)應(yīng)關(guān)系序號(hào)錯(cuò)誤碼位ESe6e5e4e3e2e1e0S3S2S101234567/b0b1b2b3b4b5b6000000000000010000010000010000010000010000010000010000000000010101000111011101112.8.2 漢明碼譯碼程序 根據(jù)上節(jié)推算的校驗(yàn)子S，可以推算出出對(duì)應(yīng)的信息位，并能對(duì)發(fā)生一位錯(cuò)誤的信息進(jìn)行糾正。將接收到的漢明碼每7個(gè)為一組進(jìn)行解

46、碼，即可得到相應(yīng)的4位信息，重組后得到整個(gè)漢明碼組的譯碼結(jié)果。結(jié)果圖如圖2-11所示。 圖2-11 漢明碼解碼2.9 manchester譯碼Manchester譯碼即為其編碼的逆過(guò)程，將收到的序列每2個(gè)分為一組，遇“10”則譯為“1”，遇“01”則譯為“0”。function demcst=mcstdecoder(y,L)demcst=zeros(1,L/2);for n=1:L/2 if(y(2*n-1)=1 & y(2*n)=0 ) demcst(n)=1; else if(y(2*n-1)=0 & y(2*n)=1) demcst(n)=0; end endendend解碼后的結(jié)果如圖

47、2-12所示: 圖2-12 manchester解碼2.10 PCM譯碼PCM譯碼對(duì)接收到的二進(jìn)制PCM編碼序列進(jìn)行分組，每8個(gè)數(shù)據(jù)為一組，相應(yīng)的對(duì)應(yīng)著極性碼、段落碼、段內(nèi)碼。譯碼時(shí)，先由極性碼確定信號(hào)的符號(hào)，再由段落碼確定信號(hào)所在碼段，同時(shí)由段內(nèi)碼確定段內(nèi)量化單位數(shù)。由段落數(shù)確定段落起點(diǎn)電平和段內(nèi)量化值。最后譯碼得到的信號(hào)out= duanstartpoint+ duanneima.*duanneijianju.從而將二進(jìn)制序列譯碼出對(duì)應(yīng)的模擬信號(hào)采樣的值。PCM譯碼后的結(jié)果如圖2-13所示： 圖2-13 PCM譯碼后的波形2.11 誤碼率的計(jì)算 通過(guò)比較PCM編碼后與manchester譯

48、碼后的兩組數(shù)據(jù)，統(tǒng)計(jì)出不一致的數(shù)據(jù)的個(gè)數(shù)，除以總的PCM編碼個(gè)數(shù)，即可得到誤碼率。通過(guò)分析誤碼率的大小可以判斷此通信系統(tǒng)的好壞。len=length(demcst);errornum=0;for i=1:len if(demcst(i)=pcmy(i) errornum=errornum+1; endend disp(誤碼率為 );errornum/len3 仿真結(jié)果分析 3.1 源信號(hào)與接收信號(hào)波形對(duì)比對(duì)比輸入的源信號(hào)與經(jīng)過(guò)噪聲信道后接收到的信號(hào)波形，分別如圖3-1a和圖3-1b所示：圖3-1a 信源信號(hào) 圖3-1b 接收信號(hào)如圖可見(jiàn)，復(fù)原后的信號(hào)大致與源信號(hào)相同，但還是出現(xiàn)了變形，這是因?yàn)?/p>

49、ASK調(diào)制方式的抗噪聲性能并不是很好，當(dāng)高斯白噪聲同時(shí)干擾了一組碼中的二個(gè)以上的碼元時(shí)，便會(huì)出現(xiàn)誤碼。 3.2 誤碼率統(tǒng)計(jì) 通過(guò)誤碼率統(tǒng)計(jì)函數(shù)，可以精確的算出每次仿真的誤碼率，仿真的誤碼率結(jié)果如圖3-2所示： 圖3-2 誤碼率統(tǒng)計(jì)4 心得體會(huì)課程設(shè)計(jì)開始的時(shí)候是在準(zhǔn)備考研那段期之后，沒(méi)仔細(xì)認(rèn)真的看設(shè)計(jì)課題，當(dāng)時(shí)在做MATLAB綜合設(shè)計(jì)性實(shí)驗(yàn)的時(shí)候比較吃力，覺(jué)得使用MATLAB編寫程序好難?？偸且_好多Mfile文件進(jìn)行編寫函數(shù)，最麻煩的是寫完一個(gè)子函數(shù)在單獨(dú)對(duì)其進(jìn)行驗(yàn)證就存在一些函數(shù)輸入的問(wèn)題。當(dāng)時(shí)間過(guò)了2天之后，才拿起設(shè)計(jì)要求認(rèn)真看起來(lái)，想起老師在課堂上進(jìn)行MATLAB講解時(shí)還是很感興趣的，

50、于是按照老師的要求開始完成本次課程設(shè)計(jì)。對(duì)通信系統(tǒng)的理解：因?yàn)橹坝羞^(guò)課程設(shè)計(jì)的經(jīng)驗(yàn)，于是剛開始的時(shí)候沒(méi)有急于編寫代碼，而是仔細(xì)認(rèn)真研究書本上的知識(shí)，進(jìn)行消化和理解。定性了分析了整個(gè)通信系統(tǒng)；了解了通信系統(tǒng)研究的意義；通信系統(tǒng)課程設(shè)計(jì)是培養(yǎng)我們綜合運(yùn)用所學(xué)知識(shí),發(fā)現(xiàn)、提出、分析和解決實(shí)際問(wèn)題,鍛煉實(shí)踐能力的重要環(huán)節(jié),是對(duì)學(xué)生實(shí)際工作能力的具體訓(xùn)練和考察過(guò)程，通過(guò)課程設(shè)計(jì)我能夠比較系統(tǒng)的了解理論知識(shí)，把理論和實(shí)踐相結(jié)合，并且用到生活當(dāng)中。課程設(shè)計(jì)的收獲：本次課程設(shè)計(jì)結(jié)束了，受益良多。它帶來(lái)不僅是知識(shí)上的收獲，而且精神上的收獲也是無(wú)盡的。本次課程設(shè)計(jì)是我在大學(xué)寥寥可數(shù)的幾次設(shè)計(jì)性實(shí)驗(yàn)之一，是和同學(xué)

51、們一起奮斗的一場(chǎng)經(jīng)歷。看到同學(xué)們的努力，看到老師的辛勤，自己也都不敢懈怠了。實(shí)驗(yàn)開始的那一刻，告誡自己沒(méi)有失敗，老師既然出了設(shè)計(jì)項(xiàng)目，既然已經(jīng)選擇的這個(gè)設(shè)計(jì)項(xiàng)目，就不許中途退出，這不是游戲，我想那一刻要的是勇氣，是決心吧！實(shí)驗(yàn)成功的那一刻，收獲的是信心，是對(duì)自己本學(xué)期以來(lái)學(xué)習(xí)通信專業(yè)課的一個(gè)肯定，是對(duì)以后的學(xué)習(xí)的一種督促，也是對(duì)自己的未來(lái)大學(xué)校門外的生活的一個(gè)鼓勵(lì)。還有就是在設(shè)計(jì)肯定會(huì)有不甚理解的地方，這就需要我們自己動(dòng)手查找資料，增加的不僅是查找的資料上面的知識(shí)，更多的搜索的能力。編程體會(huì)：在編寫代碼的時(shí)候，不急于剛開始就對(duì)著電腦編寫程序，一定要理解和整理出來(lái)你所想要做什么；從前都沒(méi)有畫流程

52、圖的習(xí)慣，不過(guò)此次試驗(yàn)之前都有仔細(xì)分析，畫流程圖，發(fā)現(xiàn)在代碼編寫的過(guò)程中這樣不僅效率高，而且成功率還快；還有就是因?yàn)樵诔绦蛑羞^(guò)多的嵌套過(guò)多for循環(huán)語(yǔ)句和if選擇語(yǔ)句，所以很容易不小心就會(huì)導(dǎo)致代碼錯(cuò)誤，所以程序中應(yīng)盡量避免。還有就是一定要仔細(xì)，仔細(xì)，再仔細(xì)。5 參考文獻(xiàn)1樊昌信，曹麗娜. 通信原理（第六版）M,北京：國(guó)防工業(yè)出版社，20092程相君，陳生潭 信號(hào)與系統(tǒng) M,西安： 西安電子科技大學(xué)出版社,19903曹志剛，錢壓生 現(xiàn)代通信原理M,北京： 清華大學(xué)出版社，19924劉敏，魏玲 MATLAB通信仿真與應(yīng)用M，北京： 國(guó)防工業(yè)出版社，20015唐向宏，岳恒立，鄭學(xué)峰 MATLAB及在

53、電子信息類課程中的應(yīng)用M，北京： 電子工業(yè)出版社，2006 6陳杰 MATLAB寶典M，北京： 電子工業(yè)出版社7張德豐 MATLAB/Simulink建模與仿真M，北京 電子工業(yè)出版社6 附錄clc;clear all;fs=6000; %設(shè)定采樣頻率 dt=1/fs; %設(shè)定步長(zhǎng)N=50; %點(diǎn)數(shù)long=N; n=0:N-1;t=n*dt; %截止時(shí)間 fc=300; %時(shí)間函數(shù)的頻率 y=sin(2*pi*fc*t); %時(shí)間函數(shù) 余弦函數(shù) figure(1); plot(t,y); %繪制時(shí)間函數(shù)圖象title(源信號(hào)時(shí)間函數(shù)波形)L=length(y);pcmy=pcmcoder(y

54、,L,long); %調(diào)用PCM函數(shù)figure(2)stairs(pcmy); %繪制PCM編碼后的圖像title(PCM編碼)ylim(-0.2 1.2);%axis(0 300 -0.2 1.2);mcstz=mcstcoder(pcmy,long);figure(3)stairs(mcstz);title(曼徹斯特編碼)ylim(-0.2 1.2);axis(0 100 -0.2 1.2);L=length(mcstz);h=zeros(1,L+L/4*3); %產(chǎn)生編碼序列N=L+L/4*3;h=hanmincoder(mcstz,L); %調(diào)用漢明碼編碼程序figure(4);st

55、airs(h);title(漢明碼編碼)ylim(-0.2 1.2)axis(0 20 -0.2 1.2);q=askcoder(N,h);figure(5);plot(q);title(ASK調(diào)制后的波形)ylim(-2 2)xlim(0 500);r=awgn(q,14); %加入噪聲figure(6);plot(r); title(加入AWGN噪聲后的波形)xlim(0 500);rr=askdecoder(r);figure(7);plot(rr);title(ASK相干解調(diào)后的波形)xlim(0 1000);figure(8)cq=chouqu(rr);stairs(cq);ylim

56、(-0.2 1.2);xlim(0 20)title(抽取濾波后的波形)y=hanmingdecoder(cq,L) ;figure(9);stairs(y);title(漢明碼解碼)ylim(-0.2 1.2)axis(0 100 -0.2 1.2);L=length(y);demcst=mcstdecoder(y,L);figure(10)stairs(demcst);title(曼徹斯特解碼)ylim(-0.2 1.2)%axis(0 200 -0.2 1.2);pcmyout=pcmdecoder(demcst,long);figure(11)plot(pcmyout);title(P

57、CM譯碼后的波形);len=length(demcst);errornum=0;for i=1:len if(demcst(i)=pcmy(i) errornum=errornum+1; endenddisp(誤碼率為 );errornum/lenfunction pcmy= pcmcoder(y,L,long)pcmy=zeros(1,8*long);for n=1:L u=y(n); sys=mdlOutputs(u); pcmy(n-1)*8+1:8*n)=sys;endendfunction sys= mdlOutputs(u)%初始化變量y1= 0;y2= 0;y3= 0;y4= 0

58、;y5= 0;y6= 0;y7= 0;y8= 0;y1= ( u= 0) * 1; %求極性碼uabs=abs( u) ;table= 1, 1/ 128; 2 1/ 64; 3, 1/ 32; 4, 1/ 16; 5, 1/ 8; 6, 1/ 4;7, 1/ 2; 8, 1 ;uabs= uabs.*(uabs= 1).*0.9999;%限幅%求段落碼y2= mdlGetDuanLuoMa( table,uabs,2)* 1;y3= mdlGetDuanLuoMa( table,uabs,3) * 1;y4= mdlGetDuanLuoMa( table, uabs, 4) * 1;Duan

59、Hao= mdlGetDuanHao( y2, y3, y4) ;DuanNeiHao= mdlGetDuanNeiHao(DuanHao, uabs) ;y8= ( floor( DuanNeiHao./2) =DuanNeiHao./2) * 1;DuanNeiHao= floor(DuanNeiHao./ 2) ;y7= ( floor( DuanNeiHao./2)= DuanNeiHao./2) * 1;DuanNeiHao= floor(DuanNeiHao./2) ;y6= ( floor( DuanNeiHao./2)= DuanNeiHao./2) * 1;DuanNeiHa

60、o= floor(DuanNeiHao./2) ;y5= ( floor( DuanNeiHao./2)= DuanNeiHao./2) * 1;sys= y1 y2 y3 y4 y5 y6 y7 y8 ;endfunction duanHao= mdlGetDuanHao( y2, y3, y4)duanHao= y2*4+ y3*2+ y4;endfunction duanLuoMa= mdlGetDuanLuoMa( table, Is, i)switch i,case 2,duanLuoMa= ( Is= table( 4, 2) ) ;case 3,duanLuoMa= ( Is=t