【文章內容簡介】 的采樣序列進行頻譜分析，繪制其幅頻和相頻曲線，對各頻率下采樣序列 X(t)和 X(n)的幅頻曲線 進行觀察 ； 4. 設系統(tǒng)單位抽樣響應為 h(n)，求解當輸入為 X(n)時的系統(tǒng)響應 Y(n)，畫出X(n),Y(n),h(n)的時域波形及幅頻特性曲線，并利用結果驗證卷積定理的正確性 。 5. 用 FFT 對信號 X(n),Y(n),h(n)進行頻譜分析，觀察與 4 中結果有無差別； 6. 由采樣序列 X(n)恢復出連續(xù)時間信號 X1(t)，畫出其時域波形，對比 X1(t)與原來的連續(xù)時間信號 X(t)的時域波形，計算并記錄兩者最大誤差。 7. 在完成基本要求和專業(yè)知識上的擴展后，設計人機交流的窗口界面，實現輸出上的美觀。 設計方法 涉及的指令 程序中用到的指令不多，多為單純的輸入函數，輸出圖形 。主要用到的指令有以下一些。在完成第七項的人機界面時涉及的指令和調用的函數較多，比較麻煩。 clc：清屏。 clear all： 清除了所有的變量，包括全局變量 global。 close all：關閉所有窗口。 input(39。輸 入 A:39。)：輸入指令，在屏幕上顯示括號內的文字，同時要操作者輸入一個 A的值。 subplot(a,b,c)： 中 a代表所畫圖形的行數 ， b代表所畫圖形的列數 ， c 代表所畫圖形的序號 。 figure（ 1） ： 與 figure 的作用一樣，建立一個新的圖形 。 plot(n,x)：輸出 n— x的連續(xù) 信號圖形。 桂林電子科技大學實訓說明書用紙 6 stem(n,x)：輸出 n— x的離散信號圖形。 grid on：顯示格子標尺。 grid of：取消格子標尺。 title(39。連續(xù)時間信號 39。)：輸出“連續(xù)時間信號”標題。 xlabel(39。w39。)：標注 X 軸名稱。 ylabel(39。幅度 39。)：標注 Y軸名稱。 conv(x,hb)：計算 x 和 hb 的卷積。 運行程序 X(t)及其 200HZ、 500HZ 和 1000HZ 頻率抽樣信號函數 X(n)程序： clc clear all close all n=0:50 % 定義序列的長度是 50 A=input(39。請入 A 的值 A:39。) % 設置信號的有關參數 a=input(39。請入 a 的值 a:39。) w0=input(39。請入 w0的值 w0: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。離散時間信號 39。) subplot(2,1,2) plot(n,x) grid on title(39。連續(xù)時間信號 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。200Hz 理想采樣信號序列 39。) subplot(3,1,2) stem(n,y2) grid on title(39。500Hz 連續(xù)時間信號 39。) subplot(3,1,3) stem(n,y3) grid on title(39。1000Hz 連續(xù)時間信號 39。) 運行得到 桂林電子科技大學實訓說明書用紙 8 0 5 10 15 20 25 30 35 40 45 50 1 0 001002002 0 0 H z 理想采樣信號序列0 5 10 15 20 25 30 35 40 45 50 2 0 002005 0 0 H z 連續(xù)時間信號0 5 10 15 20 25 30 35 40 45 50 2 0 002001 0 0 0 H z 連續(xù)時間信號 圖 再輸入： k=25:25 W=(pi/)*k w=W/pi Y1=y1*exp(j*pi/).^(n39。*k) figure (3) subplot(2,1,1) plot(w,abs(Y1)) grid xlabel(39。w39。) ylabel(39。幅度 39。) title(39。200Hz 理想采樣信號序列的幅度譜 39。) axis([2 2 0 1000]) subplot(2,1,2) plot(w,angle(Y1)) grid xlabel(39。w39。) ylabel(39。幅角 39。) 桂林電子科技大學實訓說明書用紙 9 title(39。200Hz 理想采樣信號序列的相位譜 39。) 得到： 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 理想采樣信號序列的相位譜 圖 輸入： Y2=y2*(exp(j*pi/)).^(n39。*k) figure (4) subplot(2,1,1) plot(w,abs(Y2)) grid xlabel(39。w39。) ylabel(39。幅度 39。) title(39。500Hz 理想采樣信號序列的幅度譜 39。) axis([2 2 0 1000]) subplot(2,1,2) plot(w,angle(Y2)) grid 桂林電子科技大學實訓說明書用紙 10 xlabel(39。w39。) ylabel(39。幅角 39。) title(39。500Hz 理想采樣信號序列的相位譜 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。*k) figure (5) subplot(2,1,1) plot(w,abs(Y3)) grid xlabel(39。w39。) ylabel(39。幅度 39。) title(39。1000Hz 理想采樣信號序列的幅度譜 39。) 桂林電子科技大學實訓說明書用紙 11 axis([2 2 0 1000]) subplot(2,1,2) plot(w,angle(Y3)) grid xlabel(39。w39。) ylabel(39。幅角 39。) title(39。1000Hz 理想采樣信號序列的相位譜 39。) 得到： 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 時有橫線，產生失真， 200Hz時橫線加長，失真加大。說明采樣頻率越大，失真越小。 2. 設系統(tǒng)單位抽樣響應為 h(n)，求解當輸入為 X(n)時的系統(tǒng)響應 Y(n)，畫出X(n),Y(n),h(n)的時域波形及幅頻特性曲線，并利用結果驗證卷積定理的正確性 .（此內容將參數設置為 A=， a=，