【正文】 桂林電子科技大學實訓說明書用紙 24 參考文獻 [1] 薛年喜 .MATLAB在數(shù)字信號處理中的應(yīng)用 [M].北京 : 清華大學 ,2020 [2] 吳瑛 . 數(shù)字信號處理 [M].西安 ： 電子科技大學出版社 ， 2020 [3] 管愛紅 . MATLAB基礎(chǔ)及其應(yīng)用教程 [M].北京：電子工業(yè)出版社， 2020 。本次實訓中 要 長時間占用計算機機房 ，有機會進行實際操作訓練的機會十分難得，感謝學院領(lǐng)導對我們本次實訓的支持和關(guān)心， 使我們有了這次難得的機會加強理論知識 ，增強動手能力 。通過這樣一點點的積累，我相信我們一桂林電子科技大學實訓說明書用紙 22 定能為祖國的電子事業(yè)做出一點貢獻。雖然有時候老師的要求苛刻了一點，但是用意我們都理解。 最后 ，在這里要感謝老師的指導。 在這里，也順便說一下， 遇到問題有時候是好事，它可以讓我們學得更多。 當圖形出來的那一瞬間，前面一次又一次的調(diào)試我覺得都是值得的，這種調(diào)試成功的興奮就是不斷推動我們把實訓進行下去的動力。這時只有找書本或上網(wǎng)查找了。剛開始，遇到不懂的函數(shù)找?guī)椭藛?。有關(guān)圖形的編程涉及的函數(shù)很多，而且對于學習數(shù)字信號處理的我們也不要求掌握。 但是后面的擴展創(chuàng)新就比較 難了 。這次程序的編寫 主要要實現(xiàn) 時域信號采樣及頻譜分析 ，相對比較簡單，所以程序中也是 復(fù)雜的指令 語句 不多 。模仿學習并不可恥我覺得，可恥的是用而不想。改人家的變成自己的這是我們大部分人采用的方法。但是一般跟自己的要求總差一點點。 這次實訓的難度可能就在于程序的編寫了。就像實訓動員時老師說的那樣，只要認真去做，無論做得多少。雖然題目的難易程度不同，但是我覺得只要做了都會學到很多東西。把以前在 課堂上的理論知識變成 計算機的仿真 的過程是艱辛而漫長的，但是做出作品的那一瞬間也是幸福的。以前上 MATLAB 理論課 的時 候，我們都是 只是了解一些操作的理論知識，現(xiàn)在要直接運用軟件解決實際問題 ，這對我們的鍛煉相當大。要想更好地發(fā)揮硬件的功能，軟件的 仿真作用顯而易見。硬件和軟件都是我們電子信息工程專業(yè)必須要掌握的知識。 本次實訓 時間兩周，主要是熟悉 MATLAB 和運用 MATLAB 進行數(shù)字信號處理的仿真。)。,39。,39。,39。,1039。,[,],? 39。,39。,39。,... 39。,39。 h_push2=uicontrol(h_fig,39。plot(n,x)39。callback39。連續(xù) 39。string39。fontsize39。position39。normalized39。unit39。push39。style39。])。 x=A*exp(a*n*T0).*sin(w0*n*T0)。39。39。39。,[39。,... 39。,39。,10,39。,[,],? 39。,39。,39。,... 39。,39。 桂林電子科技大學實訓說明書用紙 20 h_edit2=uicontrol(h_fig,39。w0 =39。輸入 w039。string39。left39。horizontal39。fontsize39。position39。normalized39。unit39。text39。style39。,])。))。string39。a=str2num(get(gcbo,39。callback39。left39。horizontal39。fontsize39。position39。normalized39。unit39。edit39。style39。})。,39。,{39。,39。,39。,10,39。,[,],? 39。,39。,39。,8,... 39。,39。,39。 h_text1=uicontrol(h_fig,39。))。string39。A=str2num(get(gcbo,39。callback39。left39。horizontal39。fontsize39。position39。normalized39。unit39。edit39。style39。})。,39。,{39。,39。,39。,[,],... 39。,39。,39。,10,... 39。,39。,39。 h_text=uicontrol(h_fig,39。fontsize39。ylim39。xlim39。position39。normalized39。unit39。parent39。,[,])。,39。,39。 set(h_fig,39。parent39。,13)。 text(,str,39。 str=39。currentaxes39。)。,39。,[0,0,1,1],39。,39。,39。 H=axes(39。 3 發(fā)揮擴展 為了使輸出更美觀，設(shè)計 人機 交流的窗口界面。) axis([2 2 5 5]) 2 1 . 5 1 0 . 5 0 0 . 5 1 1 . 5 205001000w2幅度輸入信號幅度譜2 1 . 5 1 0 . 5 0 0 . 5 1 1 . 5 2505w2幅角輸入信號相位譜 圖 桂林電子科技大學實訓說明書用紙 19 恢復(fù)曲線與原信號曲線相同，說明恢復(fù)誤差很小，如果采樣頻率減小，誤差增大，采樣頻率增大，則恢復(fù)誤差更小。) title(39。) ylabel(39。) axis([2 2 0 1000]) subplot(2,1,2) plot(w2,angle(X)) grid xlabel(39。) title(39。) ylabel(39。*k2) % 序列的付里葉變換函數(shù) figure (2) subplot(2,1,1) plot(w2,abs(X)) grid xlabel(39。x1(t)39。t:毫秒 39。由 x(n)恢復(fù) x1(t)39。x(t)39。t:毫秒 39。采樣序列 x(n)39。x(n)39。n39。連續(xù)時間信號頻譜 Xa(w1)39。x(t)39。t:毫秒 39。連續(xù)時間信號 x(t)39。 程序輸入： n=0:49 T=1/fs t0=10/a Dt=1/(5*a) t=0:Dt:t0 xa=A*exp(a*t).*sin(W0*t) K1=50 k1=0:1:K1 W1max=2*pi*500 W1=W1max*k1/K1 w1=W1/pi Xa=xa*exp(j*t39。如果 x為一個矩陣時，則調(diào)用后計算出每列的 N點 FFT。y(n)的 FFT39。x(n)的 FFT39。h(n)的 FFT39。 FFT 對信號 X(n),Y(n),h(n)進行譜分析，觀察與 4中結(jié)果有無差別。輸出信號 y[n]的時域波形 39。輸入信號 x[n]的時域波形 39。矩形序列的時域波形 39。輸出信號 y[n]39。輸入序列 x[n]39。系統(tǒng) hb[n]39。 2. 設(shè)系統(tǒng)單位抽樣響應(yīng)為 h(n)，求解當輸入為 X(n)時的系統(tǒng)響應(yīng) Y(n)，畫出X(n),Y(n),h(n)的時域波形及幅頻特性曲線，并利用結(jié)果驗證卷積定理的正確性 .（此內(nèi)容將參數(shù)設(shè)置為 A=， a=， w0=， f=1000）。) 得到： 2 1 . 5 1 0 . 5 0 0 . 5 1 1 . 5 205001000w幅度1 0 0 0 H z 理想采樣信號序列的幅度譜 2 . 5 2 1 . 5 1 0 . 5 0 0 . 5 1 1 . 5 242024w幅角1 0 0 0 H z 理想采樣信號序列的相位譜 圖 觀察以上圖形，采樣頻率為 1000Hz 時沒有失真， 500Hz 時有橫線，產(chǎn)生失真， 200Hz時橫線加長，失真加大。) title(39。) ylabel(39。) 桂林電子科技大學實訓說明書用紙 11 axis([2 2 0 1000]) subplot(2,1,2) plot(w,angle(Y3)) grid xlabel(39。) title(39。) ylabel(39。*k) figure (5) subplot(2,1,1) plot(w,abs(Y3)) grid xlabel(39。500Hz 理想采樣信號序列的相位譜 39。幅角 39。w39。500Hz 理想采樣信號序列的幅度譜 39。幅度 39。w39。) 得到： 2 1 . 5 1 0 . 5 0 0 . 5 1 1 . 5 205001000w幅度2 0 0 H z 理想采樣信號序列的幅度譜 2 . 5 2 1 . 5 1 0 . 5 0 0 . 5 1 1 . 5 242024w幅角2 0 0 H z 理想采樣信號序列的相位譜 圖 輸