MATLAB 下的數字信號處理實現示例
本部分內容是本講義中數據信號處理實驗部分實驗項目在 MatLab 下實現代碼。之所以
提供這些代碼，是希望通過研究以下代碼，能夠更快、更好地掌握用MatLab 進行數據信號處理實驗的方法；提高實驗質量。希望同學們在閱讀代碼的時候，注意學習方法，在最短的時間內熟悉MatLab，提高應用能力。示例中有些部分是實驗項目中的內容實現，有些是一些典型例題的實現。研究示例代碼，倡導個性化編程是我們的目標，希望同學們能在在進行實驗項目的過程中提高MatLab 的應用能力；在學習MatLab 編程的同時加強對數字信號處理有關實驗項目的理解。
以下代碼段均在 MatLab5.3 下調試通過，但是由于排版或其他一些原因，可能有部分代
碼段不能得到正常結果。您可以在“http://202.38.75.33/dsp/matlab/”得到本講義的修訂內容，同時可以在這個網址獲取所有代碼。
附錄 C1 信號、系統和系統響應
1、理想采樣信號序列
（1）首先產生信號x(n),0<=n<=50
n=0:50; %定義序列的長度是50
A=444.128; %設置信號有關的參數
a=50*sqrt(2.0)*pi;
T=0.001; %采樣率
w0=50*sqrt(2.0)*pi; %ω符號在MatLab 中不能輸入，用w 代替
x=A*exp(-a*n*T).*sin(w0*n*T); %pi 是MATLAB 定義的π，信號乘可采用“.*”
close all %清除已經繪制的x(n)圖形
subplot(3,1,1);stem(x); %繪制x(n)的圖形
title(‘理想采樣信號序列’); %設置結果圖形的標題
（2）繪制信號x(n)的幅度譜和相位譜
k=-25:25;
W=(pi/12.5)*k;
X=x*(exp(-j*pi/12.5)).^(n’*k);
magX=abs(X); %繪制x(n)的幅度譜
subplot(3,1,2);stem(magX);title(‘理想采樣信號序列的幅度譜’);
angX=angle(X); %繪制x(n)的相位譜
subplot(3,1,3);stem(angX) ; title (‘理想采樣信號序列的相位譜’)
（3）改變參數為： 1, 0.4, 2.0734, 1 0 A = α = Ω = T =
n=0:50; %定義序列的長度是50
A=1; a=0.4; w0=2.0734; T=1; %設置信號有關的參數和采樣率T
x=A*exp(-a*n*T).*sin(w0*n*T); %pi 是MATLAB 定義的π，信號乘可采用“.*”