通信原理課程設(shè)計(jì)報(bào)告--單邊帶電臺(tái)仿真

1、<p>　　通信系統(tǒng)課程設(shè)計(jì)報(bào)告</p><p><b>　　單邊帶電臺(tái)仿真</b></p><p><b>　　院 （系）：</b></p><p><b>　　專業(yè)年級(jí)（班）：</b></p><p>　　學(xué) 生：</p><

2、p>　　學(xué) 號(hào)：</p><p><b>　　指 導(dǎo) 教 師：</b></p><p><b>　　完 成 時(shí) 間：</b></p><p><b>　　摘 要</b></p><p>　　這篇文章介紹了使用MATLAB軟件編制程序，以實(shí)現(xiàn)單邊帶信號(hào)的調(diào)制和解

3、調(diào)。首先，利用相移法從語音信號(hào)的雙邊帶信號(hào)得到單邊帶信號(hào)，再編寫MATLAB程序，使單邊帶信號(hào)得到調(diào)制和調(diào)解。分析調(diào)制前后的時(shí)域和頻域波形圖，以更加深入理解單邊帶信號(hào)的調(diào)制和解調(diào)的原理。</p><p>　　關(guān)鍵詞：單邊帶；調(diào)制和解調(diào)；MATLAB</p><p>　　Single sideband radio simulation</p><p><b>

4、;　　Abstract：</b></p><p>　　This article introduces the way to modulate and demodulate the single side band with the soft program from MATLAB . First, get the single side band signal from the speech sign

5、al`s couple sides band , then wright the MATLAB program to modulate and demodulate the signal. Analyze the pictures ,and understand the theory furtherly.</p><p>　　Keywords：single side band , modulation and d

6、emodulation , MATLAB</p><p><b>　　目錄</b></p><p>　　一．設(shè)計(jì)任務(wù)與要求3</p><p><b>　　二．開發(fā)環(huán)境3</b></p><p>　　三．電臺(tái)設(shè)計(jì)原理及方案3</p><p>　　3.1發(fā)送模塊的設(shè)計(jì)方案

7、3</p><p>　　3.2信道模塊的設(shè)計(jì)方案4</p><p>　　3.3接收模塊的設(shè)計(jì)方案4</p><p>　　四．電臺(tái)詳細(xì)設(shè)計(jì)4</p><p>　　4.1對(duì)音頻輸入與調(diào)制的仿真4</p><p>　　4.2對(duì)指定信噪比信道的仿真5</p><p>　　4.3話音信號(hào)的解調(diào)

8、5</p><p>　　4.3.1 載頻9.8KHz解調(diào)5</p><p>　　4.3.2 載頻9.9KHz解調(diào)6</p><p>　　4.3.3 載頻10KHz解調(diào)7</p><p>　　4.4 原信號(hào)和解調(diào)后的信號(hào)對(duì)比7</p><p>　　4.4.1 原信號(hào)和解調(diào)后信號(hào)波形對(duì)比7</p>

9、<p>　　4.4.2 原信號(hào)和解調(diào)后信號(hào)頻率對(duì)比8</p><p><b>　　五．設(shè)計(jì)總結(jié)9</b></p><p><b>　　參考文獻(xiàn)：9</b></p><p><b>　　附錄：10</b></p><p><b>　　一．設(shè)計(jì)任務(wù)與要求

10、</b></p><p>　　這個(gè)仿真以真實(shí)的音頻信號(hào)作為輸入，設(shè)計(jì)一個(gè)單邊帶發(fā)信機(jī)。將基帶信號(hào)調(diào)制為SSB信號(hào)后送入帶通型高斯噪聲信道，加入給定功率的噪聲之后，再送入單邊帶接收機(jī)。單邊帶接收機(jī)將型號(hào)解調(diào)下來，通過計(jì)算機(jī)聲卡將解調(diào)信號(hào)播放出來試聽效果，從而對(duì)信道信噪比與解調(diào)音質(zhì)之間的關(guān)系進(jìn)行主觀測(cè)試?，F(xiàn)設(shè)計(jì)一個(gè)單邊帶發(fā)信機(jī)、帶通信道和相應(yīng)的接收機(jī)，參數(shù)定位如下值。</p><p&g

11、t;　　1.輸入信號(hào)為一個(gè)話音信號(hào)，采樣率為 8000Hz。話音輸入后首先進(jìn)行預(yù)濾波，濾波器是一個(gè)頻率范圍在[300, 3400]Hz的帶通濾波器，其目的是將話音頻頻譜限制在3400Hz以下。單邊帶調(diào)制的載波頻率設(shè)計(jì)為 10kHz，調(diào)制輸出上邊帶。要求觀測(cè)單邊帶調(diào)制前后的信號(hào)功率譜。</p><p>　　2.信道是一個(gè)帶限高斯噪聲信道，其通帶頻率范圍是[10000, 13500]Hz。能夠根據(jù)信噪比 SNR 的要

12、求加入高斯噪聲。</p><p>　　3.接收機(jī)采用相干解調(diào)方式。為了模擬載波頻率誤差對(duì)解調(diào)話音音質(zhì)的影響，相干解調(diào)，設(shè)本地載波頻率分別為為 9.8kHz，9.9kHz,10kHz。解調(diào)濾波器設(shè)計(jì)為 300～3400Hz的帶通濾波器。</p><p><b>　　二．開發(fā)環(huán)境</b></p><p>　　操作系統(tǒng)：Microsoft Windo

13、ws 7 Ultimate（64位）</p><p>　　交互工具：鍵盤/鼠標(biāo)</p><p>　　開發(fā)工具：MATLAB2010a</p><p>　　三．電臺(tái)設(shè)計(jì)原理及方案</p><p>　　3.1發(fā)送模塊的設(shè)計(jì)方案</p><p>　　有限語音信號(hào)可以采用直接用程序讀取，當(dāng)讀取一次之后，可以把音頻文件保存起來，

14、當(dāng)用到?jīng)]有麥克風(fēng)的電腦時(shí)可以備用，發(fā)送的語音信號(hào)首先進(jìn)行希爾比特變化，然后用函數(shù)可以進(jìn)行單邊帶調(diào)制，保存單邊帶語音信號(hào)。</p><p>　　3.2信道模塊的設(shè)計(jì)方案</p><p>　　讀取發(fā)送模塊保存的語音信號(hào)，然后加入高斯白噪聲，然后改變信噪比，通常情況下加入的信噪比為20，然后濾波，模擬現(xiàn)實(shí)當(dāng)中的信道傳輸，并保存得到的語音信號(hào)。</p><p>　　3.3接

15、收模塊的設(shè)計(jì)方案</p><p>　　首先讀取信道傳輸過來的語音信號(hào)，加入不同頻率的載波，然后解調(diào)，并觀察分析不同載波解調(diào)出來的語音信號(hào)和原語音信號(hào)的區(qū)別，分析失真度。</p><p><b>　　四．電臺(tái)詳細(xì)設(shè)計(jì)</b></p><p>　　4.1對(duì)音頻輸入與調(diào)制的仿真</p><p>　　根據(jù)設(shè)置的參數(shù)，系統(tǒng)中信號(hào)最

16、高頻率約為 14kHz。為了較好地顯示調(diào)制波形，系統(tǒng)仿真采樣率設(shè)為 50kHz，滿足取樣定理。由于話音信號(hào)的采樣率為 8000Hz，與系統(tǒng)仿真采樣率不等，因此，在進(jìn)行信號(hào)處理之前，必須將話音的采樣率提高到 50kHz，用插值函數(shù)來做這一任務(wù)。</p><p>　　先編寫程序?qū)⒒鶐б纛l信號(hào)讀入，進(jìn)行[300，3400]Hz的帶通濾波，并將信號(hào)采樣率提高到 50kHz，進(jìn)行單邊帶調(diào)制之后，將調(diào)制輸出結(jié)果保存為 wav

17、 文件，文件名為 SSB_OUT.wav。</p><p>　　圖1 基帶信號(hào)和SSB的波形、功率譜</p><p>　　分析調(diào)制前后的波形頻率，可以得到SSB調(diào)制只是把原來波形搬移到10KHz位置，并沒有改變?cè)瓉聿ㄐ蔚男螤詈皖l率范圍，說明調(diào)制的效果很好。</p><p>　　4.2對(duì)指定信噪比信道的仿真</p><p>　　信道是一個(gè)帶限高

18、斯噪聲信道，其帶通頻率范圍是[10000,13500]Hz.，目前輸入的信噪比SNR為20。仿真指定信噪比信道計(jì)算信噪比為20dB時(shí)的信道輸出，將結(jié)果保存為 Chanel_out.wav 文件。方便以后的實(shí)用。</p><p>　　運(yùn)行結(jié)束：計(jì)算出的信噪比結(jié)果SNR_dB=19.9126</p><p>　　4.3話音信號(hào)的解調(diào)</p><p>　　4.3.1 載頻

19、9.8KHz解調(diào) </p><p>　　接收機(jī)采用想干解調(diào)方式。為了模擬載波頻率誤差對(duì)解調(diào)話音音質(zhì)的影響，社本地載波的頻率為9.8KHz，與發(fā)信機(jī)載波頻率相差200Hz。解調(diào)濾波器設(shè)計(jì)為300Hz到3400Hz。</p><p>　　圖2 9.8KHz解調(diào)后波形和功率譜</p><p>　　解調(diào)輸出信號(hào)被保存為音頻文件 SSBDemo_OUT9.8KHz.wav，并

20、由 sound 函數(shù)播放。聆聽播放解調(diào)輸出信號(hào)的聲音可知，在20dB信道信噪比條件下，即使解調(diào)本地載波頻率誤差達(dá)到 200Hz，聲音仍然是清晰可懂的。</p><p>　　圖3 9.8KHz解調(diào)后波形和功率譜</p><p>　　解調(diào)輸出信號(hào)被保存為音頻文件 SSBDemo_OUT9.8KHz.wav，并由 sound 函數(shù)播放。聆聽播放解調(diào)輸出信號(hào)的聲音可知，在10dB信道信噪比條件下，

21、即使解調(diào)本地載波頻率誤差達(dá)到 200Hz，聲音仍然清晰地，但是也略有失真，效果明顯比信噪比為20dB差。</p><p>　　4.3.2 載頻9.9KHz解調(diào) </p><p>　　接收機(jī)采用想干解調(diào)方式。為了模擬載波頻率誤差對(duì)解調(diào)話音音質(zhì)的影響，社本地載波的頻率為9.9KHz，與發(fā)信機(jī)載波頻率相差100Hz。解調(diào)濾波器設(shè)計(jì)為300Hz到3400Hz。將得到的波形和波形的功率譜輸出。<

22、;/p><p>　　圖4 9.9KHz解調(diào)后波形和功率譜</p><p>　　解調(diào)輸出信號(hào)被保存為音頻文件 SSBDemo_OUT9.9KHz.wav，并由 sound 函數(shù)播放。聆聽播放解調(diào)輸出信號(hào)的聲音可知，在20dB信道信噪比條件下，即使解調(diào)本地載波頻率誤差達(dá)到 100Hz，聲音仍然是清晰可懂的。</p><p>　　圖5 9.9KHz解調(diào)后波形和功率譜</

23、p><p>　　解調(diào)輸出信號(hào)被保存為音頻文件 SSBDemo_OUT9.9KHz.wav，并由 sound 函數(shù)播放。聆聽播放解調(diào)輸出信號(hào)的聲音可知，在10dB信道信噪比條件下，即使解調(diào)本地載波頻率誤差達(dá)到 100Hz，聲音仍然是清晰可懂的，但是效果比信噪比為20dB效果要差一點(diǎn)，不過基本分不出來。</p><p>　　4.3.3 載頻10KHz解調(diào) </p><p>

24、　　接收機(jī)采用想干解調(diào)方式。為了與上述不同載波解調(diào)出的結(jié)果相對(duì)比，這次的在載波和原來調(diào)制的載波相同。解調(diào)濾波器設(shè)計(jì)為300Hz到3400Hz。將得到的波形和波形的功率譜輸出</p><p>　　圖6 10KHz解調(diào)后波形和功率譜</p><p>　　解調(diào)輸出信號(hào)被保存為音頻文件 SSBDemo_OUT10KHz.wav，并由 sound 函數(shù)播放。聆聽播放解調(diào)輸出信號(hào)的聲音可知，在解調(diào)波和

25、載波頻率相同，20dB信道信噪比條件下，聲音是很清晰的。</p><p>　　圖7 10KHz解調(diào)后波形和功率譜</p><p>　　解調(diào)輸出信號(hào)被保存為音頻文件 SSBDemo_OUT10KHz.wav，并由 sound 函數(shù)播放。聆聽播放解調(diào)輸出信號(hào)的聲音可知，在解調(diào)波和載波頻率相同，10dB信道信噪比條件下，聲音是很清晰的，并且和信噪比為20dB基本沒有差別，耳朵幾乎分辨不出。<

26、;/p><p>　　4.4 原信號(hào)和解調(diào)后的信號(hào)對(duì)比</p><p>　　4.4.1 原信號(hào)和解調(diào)后信號(hào)波形對(duì)比</p><p>　　圖8 原信號(hào)波形與解調(diào)后信號(hào)波形對(duì)比</p><p>　　從波形上來看，三個(gè)頻率解調(diào)出來的音頻信號(hào)波形整體差別不大，要是要仔細(xì)看，在低頻部分最后一個(gè)頻率解調(diào)出來的波形效果最好。</p><p&g

27、t;　　4.4.2 原信號(hào)和解調(diào)后信號(hào)頻率對(duì)比</p><p>　　圖9 原信號(hào)波形與解調(diào)后信號(hào)頻譜對(duì)比</p><p>　　從頻率上來看，三個(gè)頻率解調(diào)出來的音頻信號(hào)波形整體差別不大，在細(xì)微處還是有差別的，尤其是在高頻和低頻處，有部分的失真，但是對(duì)音頻的整體效果影響不大，最后一個(gè)借條出來的頻譜和原信號(hào)的頻譜最接近，因此也是最好的解調(diào)頻率。</p><p><

28、b>　　五．設(shè)計(jì)總結(jié)</b></p><p>　　通過這門實(shí)驗(yàn)使我學(xué)習(xí)掌握了許多知識(shí)。首先是對(duì)matlab有了一個(gè)全新的認(rèn)識(shí)，其次是對(duì)matlab的更多操作和命令的使用有了更高的掌握，最重要的事對(duì)matlab的處理能力有了一個(gè)更高的飛躍尤其是對(duì)相關(guān)函數(shù)的使用及相關(guān)問題的處理。 </p><p>　　就對(duì)matlab相關(guān)的命令操作而言，通過這次實(shí)驗(yàn)的親身操作和實(shí)踐

29、，學(xué)習(xí)掌握了許多原本不知道的或者不太熟悉的命令。比如說相關(guān)m文件的建立，畫圖用到的標(biāo)注，配色，坐標(biāo)控制，同一張圖里畫幾幅不同的圖像，相關(guān)參數(shù)的設(shè)置以及相關(guān)函數(shù)的調(diào)用格式等等。 就拿建立一個(gè)數(shù)學(xué)方程而言，通過設(shè)置不同的參數(shù)達(dá)到所需要的要求和結(jié)果，而且還可以在不同的窗口建立不同的函數(shù)而達(dá)到相同的效果。而自己對(duì)于矩陣及閉環(huán)傳遞函數(shù)的建立原本所掌握的知識(shí)幾乎為零，而通過這次實(shí)驗(yàn)使我徹底的掌握了相關(guān)的命令操作和處理的方法，在這里我們不僅

30、可以通過建立函數(shù)和參數(shù)來達(dá)到目標(biāo)效果，而且還可以通過可視化的編程達(dá)到更快更方便，更簡潔的效果。就拿可視化編程而言原本根本就只是聽說而已罷了，從來就沒有親身去嘗試過，然而現(xiàn)在自己卻可以和容易的通過搭建不同功能木塊來實(shí)現(xiàn)相關(guān)的函數(shù)及功能。這些在原本根本就不敢相信，然而通過學(xué)習(xí)和實(shí)驗(yàn)親身操作這些原本看似不可能的操作在此就變的輕而易舉的事了。</p><p>　　通過對(duì)同一個(gè)模塊分析其對(duì)應(yīng)的不同的參數(shù)分析圖的建立去分析和

31、解釋其對(duì)應(yīng)的相關(guān)功能和技術(shù)指標(biāo)和性能分析是非常重要的，我們不可能只需要建立相關(guān)的模塊和功能就說自己掌握了所有的相關(guān)知識(shí)和技術(shù)，真正的技術(shù)和知識(shí)是怎么去分析和解釋相關(guān)的技術(shù)指標(biāo)和功能參數(shù)才是重中之重。就此而言，我坦誠的說自己所掌握的還是十分的有限的，但是老師給我們介紹的相關(guān)方法和技巧還是十分有效果的，如果自己真的想在這方面有什么建樹對(duì)自己以后的要求還是需要更改的要求的，萬不可以就此止步不前，自命不凡，我們還需掌握和了解還有許多許多，我們真

32、正所掌握的只是皮毛，要想取得更大的成績就得不斷的去努力學(xué)習(xí)和汲取相關(guān)的知識(shí)和技巧。萬不可自以為傲，裹足不前，matlab真的是個(gè)非常強(qiáng)大和有用的工具我們真正的能把它學(xué)懂學(xué)透的話還是需要下非常大的功夫和努力的。</p><p><b>　　參考文獻(xiàn)：</b></p><p>　　[1]張渭濱,陳方,蘇武潯等.Chen-Mobius多路數(shù)字通信系統(tǒng)的Matlab仿真[J]

33、.系統(tǒng)仿真學(xué)報(bào),2007,19(8):1886-1889.DOI:10.3969/j.issn.1004-731X.2007.08.058. </p><p>　　[2]劉素心,王汝芳,張廣森等.基于Matlab的通信系統(tǒng)的實(shí)驗(yàn)仿真設(shè)計(jì)[J].實(shí)驗(yàn)室科學(xué),2008,(3):101-103,105.DOI:10.3969/j.issn.1672-4305.2008.03.043. </p><p

34、>　　[3]王峰,丁金林.基于MATLAB/GUI的數(shù)字通信系統(tǒng)虛擬實(shí)驗(yàn)平臺(tái)的設(shè)計(jì)[J].南通職業(yè)大學(xué)學(xué)報(bào),2013,27(1):96-100.DOI:10.3969/j.issn.1008-5327.2013.01.026. </p><p>　　[4]陳蕾,姚遠(yuǎn)程,秦明偉等.自適應(yīng)抗干擾通信系統(tǒng)中頻譜感知技術(shù)研究[J].電視技術(shù),2014,38(5):101-104. </p><p

35、>　　[5]梅志紅,鄧文華.基于Matlab與VB集成技術(shù)的通信系統(tǒng)仿真軟件的開發(fā)[J].計(jì)算機(jī)仿真,2004,21(6):195-199.DOI:10.3969/j.issn.1006-9348.2004.06.056.</p><p>　　[6]范偉,翟傳潤,戰(zhàn)興群等.基于MATLAB的擴(kuò)頻通信系統(tǒng)仿真研究[J].微計(jì)算機(jī)信息,2006,22(19):242-244.DOI:10.3969/j.issn

36、.1008-0570.2006.19.086. </p><p>　　[7]周曉蘭,張杰.MATLAB在通信系統(tǒng)仿真中的應(yīng)用[J].計(jì)算機(jī)技術(shù)與發(fā)展,2006,16(9):166-168.DOI:10.3969/j.issn.1673-629X.2006.09.057. </p><p>　　[8]解偉俊,胡修林,張?zhí)N玉等.基于MATLAB的仿真通信語音庫實(shí)現(xiàn)[J].華中理工大學(xué)學(xué)報(bào),20

37、00,28(2):77-79.DOI:10.3321/j.issn:1671-4512.2000.02.027.</p><p><b>　　附錄：</b></p><p>　　%FileName:ForSSB.m</p><p><b>　　clc;</b></p><p>　　clear all

38、;</p><p>　　%功能，采樣點(diǎn)數(shù)40000，采樣率為8000</p><p>　　%jilu = wavrecord(5*8000, 8000, 'double');</p><p>　　%wavwrite(jilu, 'GDGvoice8000.wav');</p><p>　　[wav, fs] =

39、 wavread('GDGvoice8000.wav');</p><p><b>　　%計(jì)算聲音時(shí)間長度</b></p><p>　　t_end = 1/fs * length(wav);</p><p>　　%仿真系統(tǒng)采樣時(shí)間點(diǎn)</p><p>　　Fs = 50000;</p><

40、;p>　　t = 1/Fs:1/Fs:t_end;</p><p>　　%設(shè)計(jì)300~3400hz的帶通濾波器</p><p>　　[fenzi, fenmu] = butter(3, [300 3400]/(fs/2));</p><p>　　%對(duì)音頻信號(hào)進(jìn)行濾波</p><p>　　wav = filter(fenzi, fenmu

41、, wav);</p><p><b>　　%輸出濾波后的聲音</b></p><p>　　wavwrite(wav, 'LVBO_OUT.wav');</p><p>　　%利用插值函數(shù)將音頻的采樣率提升為50khz</p><p>　　wav = interp1([1/fs:1/fs:t_end],

42、wav, t, 'spline');</p><p>　　%音頻信號(hào)的希爾伯特變換</p><p>　　wav_hilbert = imag(hilbert(wav));</p><p><b>　　%載波頻率</b></p><p>　　fc = 10000;</p><p>&

43、lt;b>　　%單邊帶調(diào)制</b></p><p>　　SSB_OUT = wav.*cos(2*pi*fc*t) - wav_hilbert.*sin(2*pi*fc*t);</p><p><b>　　%輸出波形</b></p><p>　　figure(1);</p><p>　　subplot(

44、2, 2, 1); plot(wav(53550:53750)); axis([0 200 -0.5 0.5]);</p><p>　　title('基帶信號(hào)波形');xlabel('時(shí)間（樣值數(shù)）');</p><p>　　subplot(2, 2, 2); psd(wav, 10000, Fs); axis([0 2500 -20 10]);</p

45、><p>　　title('基帶信號(hào)功率譜');xlabel('頻率/Hz');ylabel('功率譜/(dB)');</p><p>　　subplot(2, 2, 3); plot(SSB_OUT(53550:53750)); axis([0 200 -0.5 0.5]);</p><p>　　title('S

46、SB信號(hào)波形');xlabel('時(shí)間（樣值數(shù)）');</p><p>　　subplot(2, 2, 4); psd(SSB_OUT, 10000, Fs);axis([0 25000 -20 10]);</p><p>　　title('SSB信號(hào)功率譜');xlabel('頻率/Hz');ylabel('功率譜/(dB)

47、');</p><p>　　%將 SSB 調(diào)制輸出存盤備用</p><p>　　wavwrite(0.5*SSB_OUT,Fs, 'SSB_OUT.wav');</p><p>　　%FileName:ForSNB.m</p><p>　　%功能;計(jì)算實(shí)際信噪比</p><p><b>

48、;　　%clear;</b></p><p>　　[in, Fs] = wavread('SSB_OUT.wav');</p><p>　　%設(shè)定信噪比為20時(shí)的新到輸出</p><p>　　SNR_db = 20;</p><p>　　out = ChanelSimulink(in, SNR_db);</p

49、><p>　　wavwrite(out, Fs, 'Chanel_OUT.wav');</p><p>　　%FileName:ChanelSimulink.m</p><p>　　function out = ChanelSimulink(in, SNR_db)</p><p>　　% SNR_db 設(shè)定信噪比</p>

50、;<p>　　% in 輸入信號(hào)序列</p><p>　　% out 信道輸出序列</p><p><b>　　% 系統(tǒng)采樣率</b></p><p>　　Fs = 50000;</p><p>　　Power_of_in = var(in);</p><p>　　Power_of_

51、noise = Power_of_in/(10.^(SNR_db/10));</p><p><b>　　% 信道帶寬</b></p><p>　　bandwidth = 13500 - 10000;</p><p>　　% 噪聲功率譜密度值 W/Hz</p><p>　　NO = Power_of_noise/band

52、width;</p><p>　　Gause_noise = sqrt(NO*Fs/2) .* randn(size(in));</p><p>　　% 噪聲通道 10～13.5kHz</p><p>　　[num, den] = butter(4, [10000 13500]/(Fs/2));</p><p>　　signal_of_fil

53、ter_out = filter(num, den, in);</p><p>　　noise_of_filter_out = filter(num, den, Gause_noise);</p><p>　　SNR_dB = 10*log10(var(signal_of_filter_out)/var(noise_of_filter_out));</p><p>

54、<b>　　% 測(cè)量得出信噪比</b></p><p><b>　　% 信道輸出</b></p><p>　　out = signal_of_filter_out + noise_of_filter_out;</p><p><b>　　SNR_dB</b></p><p>　

55、　% FileName:Forjietiao 9.8K.m</p><p><b>　　% clear;</b></p><p><b>　　Fs=50000;</b></p><p>　　% 讀入信道輸出信號(hào)數(shù)據(jù)</p><p>　　[recvsignal, Fs] = wavread('

56、Chanel_OUT.wav');</p><p>　　t = (1/Fs:1/Fs:length(recvsignal)/Fs)';</p><p><b>　　% 本地載波頻率</b></p><p>　　fc_local1 = 10000 -200;</p><p>　　% fc_local = 1

57、0000 -100;</p><p>　　% fc_local = 10000;</p><p><b>　　% 本地載波</b></p><p>　　local_carrier1 = cos(2*pi*fc_local1.*t);</p><p><b>　　% 相干解調(diào)</b></p>

58、;<p>　　xianggan_out1 = recvsignal.*local_carrier1;</p><p>　　% 設(shè)計(jì) 300～3400Hz 的帶通濾波器</p><p>　　[fenzi, fenmu] = butter(3, [300 3400]/(Fs/2));</p><p>　　demod_out1 = filter(fenzi,

59、 fenmu, xianggan_out1);</p><p>　　%sound(demod_out1/max(demod_out1), Fs);</p><p>　　wavwrite(demod_out1, Fs, 'SSBDemod_OUT_9.8k.wav');</p><p>　　figure(2);</p><p>

60、　　subplot(1, 2, 1); plot(5*demod_out1(53550:53750)); axis([0 200 -0.3 0.3]);</p><p>　　title('解調(diào)信號(hào)波形');xlabel('時(shí)間（樣值數(shù)）');</p><p>　　subplot(1, 2, 2); psd(5*demod_out1, 10000, Fs);

61、axis([0 2500 -20 10]);</p><p>　　title('解調(diào)信號(hào)功率譜');xlabel('頻率/Hz');ylabel('功率譜/(dB)');</p><p>　　%FileName:Forjietiao 9.9K.m</p><p><b>　　%clear;</b>

62、</p><p><b>　　Fs=50000;</b></p><p>　　% 讀入信道輸出信號(hào)數(shù)據(jù)</p><p>　　[recvsignal, Fs] = wavread('Chanel_OUT.wav');</p><p>　　t = (1/Fs:1/Fs:length(recvsignal)/F

63、s)';</p><p><b>　　% 本地載波頻率</b></p><p>　　%fc_local = 10000 -200;</p><p>　　fc_local2 = 10000 -100;</p><p>　　%fc_local = 10000;</p><p><b>

64、;　　% 本地載波</b></p><p>　　local_carrier2 = cos(2*pi*fc_local2.*t);</p><p><b>　　% 相干解調(diào)</b></p><p>　　xianggan_out2 = recvsignal.*local_carrier2;</p><p>　　%

65、 設(shè)計(jì) 300～3400Hz 的帶通濾波器</p><p>　　[fenzi, fenmu] = butter(3, [300 3400]/(Fs/2));</p><p>　　demod_out2 = filter(fenzi, fenmu, xianggan_out2);</p><p>　　%sound(demod_out2/max(demod_out2),

66、Fs);</p><p>　　wavwrite(demod_out2, Fs, 'SSBDemod_OUT_9.9k.wav');</p><p>　　figure(3);</p><p>　　subplot(1, 2, 1); plot(5*demod_out2(53550:53750)); axis([0 200 -0.3 0.3]);</

67、p><p>　　title('解調(diào)信號(hào)波形');xlabel('時(shí)間（樣值數(shù)）');</p><p>　　subplot(1, 2, 2); psd(5*demod_out2, 10000, Fs); axis([0 2500 -20 10]);</p><p>　　title('解調(diào)信號(hào)功率譜');xlabel('

68、;頻率/Hz');ylabel('功率譜/(dB)');</p><p>　　%FileName:Forjietiao 10K.m</p><p><b>　　%clear;</b></p><p><b>　　Fs=50000;</b></p><p>　　% 讀入信道輸出

69、信號(hào)數(shù)據(jù)</p><p>　　[recvsignal, Fs] = wavread('Chanel_OUT.wav');</p><p>　　t = (1/Fs:1/Fs:length(recvsignal)/Fs)';</p><p><b>　　% 本地載波頻率</b></p><p>　　%

70、fc_local = 10000 -200;</p><p>　　%fc_local = 10000 -100;</p><p>　　fc_local3 = 10000;</p><p><b>　　% 本地載波</b></p><p>　　local_carrier3 = cos(2*pi*fc_local3.*t);

71、</p><p><b>　　% 相干解調(diào)</b></p><p>　　xianggan_out3 = recvsignal.*local_carrier3;</p><p>　　% 設(shè)計(jì) 300～3400Hz 的帶通濾波器</p><p>　　[fenzi, fenmu] = butter(3, [300 3400]/

72、(Fs/2));</p><p>　　demod_out3 = filter(fenzi, fenmu, xianggan_out3);</p><p>　　%sound(demod_out/max(demod_out), Fs);</p><p>　　wavwrite(demod_out3, Fs, 'SSBDemod_OUT_10k.wav');

73、</p><p>　　figure(4);</p><p>　　subplot(1, 2, 1); plot(5*demod_out3(53550:53750)); axis([0 200 -0.3 0.3]);</p><p>　　title('解調(diào)信號(hào)波形');xlabel('時(shí)間（樣值數(shù)）');</p><p

74、>　　subplot(1, 2, 2); psd(5*demod_out3, 10000, Fs); axis([0 2500 -20 10]);</p><p>　　title('解調(diào)信號(hào)功率譜');xlabel('頻率/Hz');ylabel('功率譜/(dB)');</p><p>　　%插值函數(shù)interp1的用法</p&

75、gt;<p>　　hours=1:12;</p><p>　　temps=[5 8 9 15 25 29 31 30 22 25 27 24];</p><p>　　h=1:0.1:12;</p><p>　　t=interp1(hours,temps,h,'spline');%(直接輸出將是很多的，spline三次樣條插值)</

76、p><p>　　plot (hours,temps,'+',h,t,hours,temps,'r');%作圖</p><p>　　xlabel('hour'),ylabel('degrees celsius');</p><p>　　%filename:Forduibi.m</p><p

77、>　　%將原來信號(hào)的頻譜和三次解調(diào)出來的頻譜結(jié)果作對(duì)比，并分析結(jié)果</p><p>　　figure(5);</p><p>　　subplot(2, 2, 1); psd(wav, 10000, Fs); axis([0 2500 -20 10]);</p><p>　　title('基帶信號(hào)功率譜');xlabel('頻率/Hz&#

78、39;);ylabel('功率譜/(dB)');</p><p>　　subplot(2, 2, 2); psd(5*demod_out1, 10000, Fs); axis([0 2500 -20 10]);</p><p>　　title('9.8KHz解調(diào)信號(hào)功率譜');xlabel('頻率/Hz');ylabel('功率譜/(

79、dB)');</p><p>　　subplot(2, 2, 3); psd(5*demod_out2, 10000, Fs); axis([0 2500 -20 10]);</p><p>　　title('9.9KHz解調(diào)信號(hào)功率譜');xlabel('頻率/Hz');ylabel('功率譜/(dB)');</p>

80、<p>　　subplot(2, 2, 4); psd(5*demod_out3, 10000, Fs); axis([0 2500 -20 10]);</p><p>　　title('10KHz解調(diào)信號(hào)功率譜');xlabel('頻率/Hz');ylabel('功率譜/(dB)');</p><p>　　%將原來信號(hào)的波形和三次

81、解調(diào)出來的信號(hào)波形結(jié)果作對(duì)比，并分析結(jié)果</p><p>　　figure(6);</p><p>　　subplot(2, 2, 1); plot(wav(53550:53750)); axis([0 200 -0.3 0.3]);</p><p>　　title('基帶信號(hào)波形');xlabel('時(shí)間（樣值數(shù)）');</p

82、><p>　　subplot(2, 2, 2); plot(5*demod_out1(53550:53750)); axis([0 200 -0.3 0.3]);</p><p>　　title('9.8KHz解調(diào)信號(hào)波形');xlabel('時(shí)間（樣值數(shù)）');</p><p>　　subplot(2, 2, 3); plot(5*de

83、mod_out2(53550:53750)); axis([0 200 -0.3 0.3]);</p><p>　　title('9.9KHz解調(diào)信號(hào)波形');xlabel('時(shí)間（樣值數(shù)）');</p><p>　　subplot(2, 2, 4); plot(5*demod_out3(53550:53750)); axis([0 200 -0.3 0.3