【正文】 的輸出端，將輸出的數(shù)字信號轉(zhuǎn)換為模擬信號即為模擬系統(tǒng)的輸出。信號采樣包括兩個步驟：即將變量和值都連續(xù)的模擬信號先后轉(zhuǎn)換為在變量上離散的的離散信號和值上也離散的 數(shù)字信號 （ 量化 ）。自相關被定義為對信號本身在變化的時間和空間坐標上做 互相關 處理。 在頻域分析信號的最常見目的是分析信號屬性。 7. 在完成基本要求和專業(yè)知識上的擴展后，設計人機交流的窗口界面，實現(xiàn)輸出上的美觀。 clc：清屏。輸 入 A:39。 plot(n,x)：輸出 n— x的連續(xù) 信號圖形。 title(39。w39。)：標注 Y軸名稱。) % 設置信號的有關參數(shù) a=input(39。) T1= T2= T3= T0= x=A*exp(a*n*T0).*sin(w0*n*T0) y1=A*exp(a*n*T1).*sin(w0*n*T1) y2=A*exp(a*n*T2).*sin(w0*n*T2) y3=A*exp(a*n*T3).*sin(w0*n*T3) close all subplot(2,1,1) stem(n,x) % 繪制 x(n)的圖形 grid on title(39。) 桂林電子科技大學實訓說明書用紙 7 運行后，輸入 A=,a=,w0=，得到圖形如下： 0 5 10 15 20 25 30 35 40 45 50 5 0050100150離散時間信號0 5 10 15 20 25 30 35 40 45 50 5 0050100150連續(xù)時間信號 圖 接著輸入程序： figure(2) subplot(3,1,1) stem(n,y1) grid on title(39。) subplot(3,1,3) stem(n,y3) grid on title(39。w39。200Hz 理想采樣信號序列的幅度譜 39。幅角 39。*k) figure (4) subplot(2,1,1) plot(w,abs(Y2)) grid xlabel(39。) title(39。) ylabel(39。) 得到： 2 1 . 5 1 0 . 5 0 0 . 5 1 1 . 5 205001000w幅度5 0 0 H z 理想采樣信號序列的幅度譜 2 . 5 2 1 . 5 1 0 . 5 0 0 . 5 1 1 . 5 242024w幅角5 0 0 H z 理想采樣信號序列的相位譜 圖 輸入： Y3=y3*(exp(j*pi/)).^(n39。幅度 39。w39。1000Hz 理想采樣信號序列的相位譜 39。 桂林電子科技大學實訓說明書用紙 12 輸入： clc clear all close all n=1:50 hb=zeros(1,50) hb(1)=1 hb(2)=1 hb(3)=1 hb(4)=1 hb(5)=1 close all subplot(3,1,1) stem(hb) title(39。) y=conv(x,hb) subplot(3,1,3) stem(y) title(39。) subplot(3,1,2) plot(x) grid on title(39。) 運行得到： 桂林電子科技大學實訓說明書用紙 14 0 5 10 15 20 25 30 35 40 45 5000 . 51矩形序列的時域波形0 5 10 15 20 25 30 35 40 45 50101輸入信號 x [ n ] 的時域波形0 10 20 30 40 50 60 70 80 90 100101輸出信號 y [ n ] 的時域波形 圖 數(shù)字信號處理中經(jīng)常要進行卷積運算， conv 可以用來計算兩個有限長序列的卷積，該函數(shù)計算的兩個序列都是從 n=0 開始。) subplot(3,1,2) plot(abs(fft(x))) title(39。) 運行得到： 0 5 10 15 20 25 30 35 40 45 500246h ( n ) 的 FFT0 5 10 15 20 25 30 35 40 45 50012x ( n ) 的 FFT0 10 20 30 40 50 60 70 80 90 100012y ( n ) 的 FFT 圖 桂林電子科技大學實訓說明書用紙 16 matlab 中， 計算矢量 x 的 DFT 及其逆變換的函數(shù)分別為 fft 和 ifft，這兩個函數(shù)采用了混合算法，當 N為質(zhì)數(shù)時，采用的是原始的 DFT 算法。*W1) x=A*exp(a*n*T).*sin(W0*n*T) figure (1) subplot(4,1,1) plot(t*1000,xa) title(39。) ylabel(39。) axis([0 1000 0 1200]) subplot(4,1,3) stem(x) grid xlabel(39。) title(39。) ylabel(39。) grid xlabel(39。) axis([0 45 20 160]) 運行得到： 0 5 10 15 20 25 30 35 40 45 5 0050100150連續(xù)時間信號 x ( t )t: 毫秒x(t)0 100 200 300 400 500 600 700 800 900 100005001000連續(xù)時間信號頻譜 X a ( w 1 )0 5 10 15 20 25 30 35 40 45 50050100150t: 毫秒x(t)采樣序列 x ( n )0 5 10 15 20 25 30 35 40 45050100150由 x ( n ) 恢復 x 1 ( t )t: 毫秒x1(t) 圖 桂林電子科技大學實訓說明書用紙 18 繼續(xù)輸入： error=max(abs(x1xa)) k2=25:25 W2=(pi/)*k2 w2=W2/pi X=x*(exp(j*pi/)).^(n39。幅度 39。w239。輸入信號相位譜 39。 程序輸入： clear all。normalized39。visible39。 set(gcf,39。\fontname{隸書 }輸入時域信號的波形 39。 h_fig=get(H,39。unit39。position39。,h_fig,... 39。,39。,[0 50],39。,12)。text39。unit39。position39。left39。輸入 A39。 h_edit=uicontrol(h_fig,39。,... 39。,39。,10,. 39。,... 39。39。])。text39。unit39。position39。horizontal39。string39。a =39。,39。,39。,[,],? 39。,39。,[39。39。 h_text2=uicontrol(h_fig,39。,... 39。,39。,10,39。,39。,39。style39。unit39。position39。horizontal39。callback39。string39。,... 39。 h_push1=uicontrol(h_fig,39。,... 39。,39。,10,39。,39。)。push39。normalized39。fontsize39。離散 39。stem(n,x)39。我選擇的 題目是 基于 matlab 的時域信號采樣及頻譜分析 ， 通過程序的輸入運行實現(xiàn)信號的仿真 。 學好軟件，上機進行實際的操作是熟練軟件使用的有效方法 。本次實訓 最多兩 人做一個題目， 我的題目幾乎只是我自己一個人選 ，更 具有挑戰(zhàn)性 。確實，人的能力各不相同，但是努力程 度卻是我們所能控制的。所以我們不能完完全全復制人家的程序。這樣沒什么作用對我們。 實訓中也遇到了不少問題，比如 擴展創(chuàng)新時做了一個好看的輸入窗口。但是另一個問題出現(xiàn)了，就是英文水平不高，看不懂幫助的提示。我相信這也是不少同學遇到 問題 后能 進行下去的根本動力 。沒有老師的指導督促，本次實訓也不能按時按質(zhì)完成。 桂林電子科技大學實訓說明書用紙 23 謝 辭 首先，感謝應用科技學院為我 們電子信息工程專業(yè)的學生安排了這次實