利用matlab實現(xiàn)連續(xù)信號的采樣與重構(gòu)仿真課程設計

1、<p><b>　　課程設計</b></p><p>　　課程名稱：連續(xù)信號的采樣與重構(gòu)</p><p><b>　　專業(yè)班級：通信工程</b></p><p><b>　　目 錄</b></p><p><b>　　1、摘要1</b><

2、;/p><p><b>　　2、正文2</b></p><p>　　2.1、設計目的2</p><p>　　2.2、設計原理 </p><p>　　(1)、連續(xù)時間信號2</p><p>　　(2)、采樣定理3</p><p>　　(3)、信號重構(gòu)5</p

3、><p>　　2.3、信號采樣與恢復的程序5</p><p>　　（1）設計連續(xù)信號6</p><p>　　(2)設計連續(xù)信號的頻譜7</p><p>　?。?）設計采樣信號8</p><p>　　（4）設計采樣信號的頻譜圖9</p><p>　?。?）設計低通濾波器10</p&g

4、t;<p>　　（6）恢復原信號12</p><p>　　3、總結(jié)與致謝13</p><p><b>　　4、參考文獻14</b></p><p><b>　　1.摘要</b></p><p>　　本次課程設計應用MATLAB實現(xiàn)連續(xù)信號的采樣與重構(gòu)仿真，了解MATLAB軟件，學

5、習應用MATLAB軟件的仿真技術(shù)。它主要側(cè)重于某些理論知識的靈活運用，以及一些關鍵命令的掌握，理解，分析等。初步掌握線性系統(tǒng)的設計方法，培養(yǎng)獨立工作能力。</p><p>　　加深理解采樣與重構(gòu)的概念，掌握利用MATLAB分析系統(tǒng)頻率響應的方法和掌握利用MATLAB實現(xiàn)連續(xù)信號采用與重構(gòu)的方法。計算在臨界采樣、過采樣、欠采樣三種不同條件下重構(gòu)信號的誤差，并由此總結(jié)采樣頻率對信號重構(gòu)誤差的影響。</p>

6、<p>　　要做到以下基本要求：</p><p>　　1. 掌握利用MATLAB分析系統(tǒng)頻率響應的方法，增加對仿真軟件MATLAB的感性認識，學會該軟件的操作和使用方法。</p><p>　　2. 掌握利用MATLAB實現(xiàn)連續(xù)信號采用與重構(gòu)的方法，加深理解采樣與重構(gòu)的概念。</p><p>　　3 . 初步掌握線性系統(tǒng)的設計方法，培養(yǎng)獨立工作能力。&

7、lt;/p><p>　　4. 學習MATLAB中信號表示的基本方法及繪圖函數(shù)的調(diào)用，實現(xiàn)對常用連續(xù)時間信號的可視化表示，加深對各種電信號的理解。</p><p>　　5. 加深理解采樣對信號的時域和頻域特性的影響；驗證信號與系統(tǒng)的基本概念、基本理論，掌握信號與系統(tǒng)的分析方法。</p><p>　　6. 加深對采樣定理的理解和掌握，以及對信號恢復的必要性；掌握對連續(xù)信號在

8、時域的采樣與重構(gòu)的方法。</p><p><b>　　2.正文</b></p><p>　　2.1 設計目的與要求</p><p>　　對連續(xù)信號進行采樣，在滿足采樣定理和不滿足采用定理兩種情況下對連續(xù)信號和采樣信號進行FFT頻譜分析。</p><p><b>　　2.2 設計原理</b></

9、p><p>　?。?） 連續(xù)時間信號</p><p>　　連續(xù)信號是指自變量的取值范圍是連續(xù)的，且對于一切自變量的取值，除了有若干個不連續(xù)點以外，信號都有確定的值與之對應。嚴格來說，MATLAB并不能處理連續(xù)信號，而是用等時間間隔點的樣值來近似表示連續(xù)信號。當取樣時間間隔足夠小時，這些離散的樣值就能較好地近似連續(xù)信號。</p><p>　　在一定條件下，一個連續(xù)時間信號

10、完全可以用該信號在等時間間隔上的瞬時值來表示，并且可以用這些樣本值把信號完全恢復過來。這樣，抽樣定理為連續(xù)時間信號與離散時間信號的相互轉(zhuǎn)換提供了理論依據(jù)。通過觀察采樣信號的頻譜，發(fā)現(xiàn)它只是原信號頻譜的線性重復搬移，只要給它乘以一個門函數(shù)，就可以在頻域恢復原信號的頻譜，在時域是否也能恢復原信號時，利用頻域時域的對稱關系，得到了信號。</p><p><b>　　(2)采樣定理</b></

11、p><p>　　模擬信號經(jīng)過 (A/D) 變換轉(zhuǎn)換為數(shù)字信號的過程稱為采樣，信號采樣后其頻譜產(chǎn)生了周期延拓，每隔一個采樣頻率 fs，重復出現(xiàn)一次。為保證采樣后信號的頻譜形狀不失真，采樣頻率必須大于信號中最高頻率成分的兩倍，這稱之為采樣定理。</p><p>　　時域采樣定理從采樣信號恢復原信號必需滿足兩個條件：</p><p>　　a、必須是帶限信號，其頻譜函數(shù)在 ＞

12、各處為零；（對信號的要求，即只有帶限信號才能適用采樣定理。）</p><p>　　b、 取樣頻率不能過低，必須 ＞2 （或 ＞2）。（對取樣頻率的要求，即取樣頻率要足夠大，采得的樣值要足夠多，才能恢復原信號。）</p><p>　　如圖1所示，給出了信號采樣原理圖</p><p>　　圖1 信號采樣原理圖</p><p>　　由圖1可見，，其

13、中，沖激采樣信號的表達式為： </p><p>　　其傅立葉變換為，其中。設，分別為，的傅立葉變換，由傅立葉變換的頻域卷積定理，可得</p><p>　　（1） </p><p>　　若設是帶限信號，帶寬為， 經(jīng)過采樣后的頻譜就是將在頻率軸上搬移至處（幅度為原頻譜的倍）。因此，當時

14、，頻譜不發(fā)生混疊；而當時，頻譜發(fā)生混疊。</p><p>　　一個理想采樣器可以看成是一個載波為理想單位脈沖序列的幅值調(diào)制器，即理想采樣器的輸出信號，是連續(xù)輸入信號調(diào)制在載波上的結(jié)果，如圖2所示。</p><p><b>　　圖2 信號的采樣</b></p><p>　　用數(shù)學表達式描述上述調(diào)制過程，則有 </p><

15、;p>　　理想單位脈沖序列可以表示為 </p><p>　　其中是出現(xiàn)在時刻，強度為1的單位脈沖。由于</p><p>　　的數(shù)值僅在采樣瞬時才有意義，同時，假設</p><p>　　所以又可表示為 </p><p><b>　?。?） 信號重構(gòu)</b&g

16、t;</p><p>　　設信號被采樣后形成的采樣信號為，信號的重構(gòu)是指由經(jīng)過內(nèi)插處理后，恢復出原來信號的過程,又稱為信號恢復。</p><p>　　若設是帶限信號，帶寬為，經(jīng)采樣后的頻譜為。設采樣頻率，則由式（1）知是以為周期的譜線?，F(xiàn)選取一個頻率特性（其中截止頻率滿足）的理想低通濾波器與相乘，得到的頻譜即為原信號的頻譜。</p><p>　　2.3 信號采樣與恢

17、復的程序</p><p>　　現(xiàn)在以正弦函數(shù)為例，進行MATLAB仿真實驗。</p><p><b>　　設計連續(xù)信號。</b></p><p>　　先制作一個程序，使之產(chǎn)生一個正弦連續(xù)信號。所用程序如下所示：</p><p><b>　　f1=50;</b></p><p>

18、;　　t=(1:50)/2000; %時間軸步距</p><p>　　x=sin(2*pi*t*f1);</p><p>　　figure(1);</p><p>　　plot(x); %繪制x(t)的圖形圖片號加底框</p><p>　　xlabel('t');ylabe

19、l('x(t)');</p><p>　　title('連續(xù)時間信號波形'); %圖片命名 </p><p><b>　　grid;</b></p><p><b>　　產(chǎn)生的圖形如下：</b></p><

20、;p>　?。?）設計連續(xù)信號的頻譜</p><p>　　設計一頻譜程序，使其產(chǎn)生頻譜波形圖。程序如下：</p><p>　　n=0:511; %長度</p><p>　　N=512; %設采樣點的N值</p><p>　　Xk=abs(fft(x,N)); </p>

21、<p>　　figure(2); %頻域波形</p><p>　　plot(n,Xk); </p><p>　　axis([0 N 1.1*min(Xk) 1.1*max(Xk)]); </p><p>　　%可用axis函數(shù)來調(diào)整圖軸的范圍</p><p>　　xla

22、bel('時域頻譜波形圖'); </p><p>　　ylabel('|Xk|');</p><p><b>　　波形如下：</b></p><p><b>　?。?）設計采樣信號</b></p><p>　　設計一采樣程序，使之輸出采樣波形。程序如下：</p&

23、gt;<p>　　X=fft(x,512);w=(0:255)/256*500; </p><p>　　T=4*t; </p><p>　　x=sin(2*pi*T*f1);</p><p>　　figure(3); </p><p>　　stem(x

24、) ; %圖形x（n）的繪制</p><p>　　xlabel('n');ylabel('x(n)');</p><p>　　title('采樣信號波形圖'); %圖形命名</p><p><b>　　grid;</b&

25、gt;</p><p><b>　　波形如下：</b></p><p>　?。?）設計采樣信號的頻譜圖</p><p>　　設計出該采樣信號的頻譜程序，程序如下：figure(4);plot(w,abs([X(1:256)])); %頻譜圖的繪制</p><p>　　xlabel('Hz');yla

26、bel('頻率響應幅度');</p><p>　　title('采樣頻譜波形圖'); %命名</p><p><b>　　grid;</b></p><p><b>　　波形如下：</b></p><p>　?。?）設計低通濾波器</

27、p><p>　　設計一低通濾波器，使之具有濾波作用。程序如下：</p><p>　　[B,A]=butter(8,350/500); %巴特沃斯低通濾波器的設計</p><p>　　[H,w]=freqz(B,A,512,2000);</p><p>　　figure(5);subplot(2,1,1);</p><p>

28、　　plot(w*2000/(2*pi),abs(H)); %低通頻譜圖的繪制</p><p>　　xlabel('Hz');ylabel('頻率響應幅度');</p><p>　　title('低通濾波器波形圖'); %命名</p><p><b>　　grid;</b>&

29、lt;/p><p>　　當采樣頻率f=350Hz時，波形圖如下：</p><p><b>　?。?）恢復原信號。</b></p><p>　　設計程序，對采樣信號頻譜進行濾波，并輸出該信號所恢復頻譜信號與連續(xù)信號，程序如下：</p><p>　　y=filter(B,A,x);figure(6);subplot(2,1,1)

30、;</p><p><b>　　plot(y);</b></p><p>　　xlabel('t');ylabel('x(t)');</p><p>　　title('連續(xù)信號波形');</p><p><b>　　grid;</b></p>

31、;<p>　　Y=fft(y,512);w=(0:255)/256*500;</p><p>　　subplot(2,1,2);</p><p>　　plot(w,abs([Y(1:256)])); %頻譜圖的繪制</p><p>　　xlabel('Hz');ylabel('頻率響應幅度');</p>

32、;<p>　　title('恢復后的頻譜波形圖');</p><p><b>　　grid;</b></p><p><b>　　波形如下：</b></p><p>　　下圖為采樣f=150Hz時的圖形</p><p><b>　　4.總結(jié)與致謝</b&

33、gt;</p><p>　　經(jīng)過此次MATLAB課程設計我學到了一些軟件知識和學習方法。我現(xiàn)有的知識還不足以完成這次課程設計，所以為了這次的課程設計，我查閱了一些資料，并上網(wǎng)搜索了與此有關的知識。</p><p>　　在這次設計中，同樣也學到了對信號的采樣定理的應用，以及信號的重構(gòu)，并通過觀察MATLAB所生成的頻譜圖，進一步了解了有關信號的采樣與重構(gòu)。</p><p&

34、gt;　　同時，感謝本組成員的熱心幫助下，使我能夠順利的完成課程設計。</p><p><b>　　參考資料</b></p><p>　　1．董長虹. Matlab信號處理與應用[M].北京：國防工業(yè)出版社，2005.01。</p><p>　　2．甘俊英. 基于MATLAB的信號與系統(tǒng)實驗指導[M].北京：清華大學出版社，2007.8。<

35、;/p><p>　　3.吳大正. 信號與線性系統(tǒng)分析[M].北京：高等教育出版社，2005.08。</p><p>　　4．樓順天，劉小東，李博菡.基于MATLAB7.X的系統(tǒng)分析與設計——信號處理[M].西安：西安電子科技大學出版社，2005.05。</p><p>　　5.丁志中 葉中付.頻譜無混疊采樣和信號完全可重構(gòu)采樣[J].數(shù)據(jù)采集與處理,2005,20(3)