2023年全國碩士研究生考試考研英語一試題真題(含答案詳解+作文范文)_第1頁
已閱讀1頁,還剩17頁未讀, 繼續(xù)免費閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認(rèn)領(lǐng)

文檔簡介

1、<p>  通信系統(tǒng)課程設(shè)計報告</p><p><b>  單邊帶電臺仿真</b></p><p><b>  院 (系):</b></p><p><b>  專業(yè)年級(班):</b></p><p>  學(xué) 生:</p><

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

3、調(diào)。首先,利用相移法從語音信號的雙邊帶信號得到單邊帶信號,再編寫MATLAB程序,使單邊帶信號得到調(diào)制和調(diào)解。分析調(dià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è)計任務(wù)與要求3</p><p><b>  二.開發(fā)環(huán)境3</b></p><p>  三.電臺設(shè)計原理及方案3</p><p>  3.1發(fā)送模塊的設(shè)計方案

7、3</p><p>  3.2信道模塊的設(shè)計方案4</p><p>  3.3接收模塊的設(shè)計方案4</p><p>  四.電臺詳細(xì)設(shè)計4</p><p>  4.1對音頻輸入與調(diào)制的仿真4</p><p>  4.2對指定信噪比信道的仿真5</p><p>  4.3話音信號的解調(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 原信號和解調(diào)后的信號對比7</p><p>  4.4.1 原信號和解調(diào)后信號波形對比7</p>

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

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

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

12、求加入高斯噪聲。</p><p>  3.接收機采用相干解調(diào)方式。為了模擬載波頻率誤差對解調(diào)話音音質(zhì)的影響,相干解調(diào),設(shè)本地載波頻率分別為為 9.8kHz,9.9kHz,10kHz。解調(diào)濾波器設(shè)計為 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>  三.電臺設(shè)計原理及方案</p><p>  3.1發(fā)送模塊的設(shè)計方案</p><p>  有限語音信號可以采用直接用程序讀取,當(dāng)讀取一次之后,可以把音頻文件保存起來,

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

33、.系統(tǒng)仿真學(xué)報,2007,19(8):1886-1889.DOI:10.3969/j.issn.1004-731X.2007.08.058. </p><p>  [2]劉素心,王汝芳,張廣森等.基于Matlab的通信系統(tǒng)的實驗仿真設(shè)計[J].實驗室科學(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è)計[J].南通職業(yè)大學(xué)學(xué)報,2013,27(1):96-100.DOI:10.3969/j.issn.1008-5327.2013.01.026. </p><p>  [4]陳蕾,姚遠程,秦明偉等.自適應(yīng)抗干擾通信系統(tǒng)中頻譜感知技術(shù)研究[J].電視技術(shù),2014,38(5):101-104. </p><p

35、>  [5]梅志紅,鄧文華.基于Matlab與VB集成技術(shù)的通信系統(tǒng)仿真軟件的開發(fā)[J].計算機仿真,2004,21(6):195-199.DOI:10.3969/j.issn.1006-9348.2004.06.056.</p><p>  [6]范偉,翟傳潤,戰(zhàn)興群等.基于MATLAB的擴頻通信系統(tǒng)仿真研究[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].計算機技術(shù)與發(fā)展,2006,16(9):166-168.DOI:10.3969/j.issn.1673-629X.2006.09.057. </p><p>  [8]解偉俊,胡修林,張?zhí)N玉等.基于MATLAB的仿真通信語音庫實現(xiàn)[J].華中理工大學(xué)學(xué)報,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>  %功能,采樣點數(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>  %計算聲音時間長度</b></p><p>  t_end = 1/fs * length(wav);</p><p>  %仿真系統(tǒng)采樣時間點</p><p>  Fs = 50000;</p><

40、;p>  t = 1/Fs:1/Fs:t_end;</p><p>  %設(shè)計300~3400hz的帶通濾波器</p><p>  [fenzi, fenmu] = butter(3, [300 3400]/(fs/2));</p><p>  %對音頻信號進行濾波</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>  %音頻信號的希爾伯特變換</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('基帶信號波形');xlabel('時間(樣值數(shù))');</p><p>  subplot(2, 2, 2); psd(wav, 10000, Fs); axis([0 2500 -20 10]);</p

45、><p>  title('基帶信號功率譜');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信號波形');xlabel('時間(樣值數(shù))');</p><p>  subplot(2, 2, 4); psd(SSB_OUT, 10000, Fs);axis([0 25000 -20 10]);</p><p>  title('SSB信號功率譜');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>  %功能;計算實際信噪比</p><p><b>

48、;  %clear;</b></p><p>  [in, Fs] = wavread('SSB_OUT.wav');</p><p>  %設(shè)定信噪比為20時的新到輸出</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 輸入信號序列</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>  % 測量得出信噪比</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>  % 讀入信道輸出信號數(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è)計 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)信號波形');xlabel('時間(樣值數(shù))');</p><p>  subplot(1, 2, 2); psd(5*demod_out1, 10000, Fs);

61、axis([0 2500 -20 10]);</p><p>  title('解調(dià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>  % 讀入信道輸出信號數(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è)計 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)信號波形');xlabel('時間(樣值數(shù))');</p><p>  subplot(1, 2, 2); psd(5*demod_out2, 10000, Fs); axis([0 2500 -20 10]);</p><p>  title('解調(dià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、信號數(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è)計 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)信號波形');xlabel('時間(樣值數(shù))');</p><p

74、>  subplot(1, 2, 2); psd(5*demod_out3, 10000, Fs); axis([0 2500 -20 10]);</p><p>  title('解調(dià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、>  %將原來信號的頻譜和三次解調(diào)出來的頻譜結(jié)果作對比,并分析結(jié)果</p><p>  figure(5);</p><p>  subplot(2, 2, 1); psd(wav, 10000, Fs); axis([0 2500 -20 10]);</p><p>  title('基帶信號功率譜');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)信號功率譜');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)信號功率譜');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)信號功率譜');xlabel('頻率/Hz');ylabel('功率譜/(dB)');</p><p>  %將原來信號的波形和三次

81、解調(diào)出來的信號波形結(jié)果作對比,并分析結(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('基帶信號波形');xlabel('時間(樣值數(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)信號波形');xlabel('時間(樣值數(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)信號波形');xlabel('時間(樣值數(shù))');</p><p>  subplot(2, 2, 4); plot(5*demod_out3(53550:53750)); axis([0 200 -0.3 0.3

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 眾賞文庫僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論