MATLAB 下的數(shù)字信號處理實現(xiàn)示例
本部分內(nèi)容是本講義中數(shù)據(jù)信號處理實驗部分實驗項目在 MatLab 下實現(xiàn)代碼。之所以
提供這些代碼，是希望通過研究以下代碼，能夠更快、更好地掌握用MatLab 進(jìn)行數(shù)據(jù)信號處理實驗的方法；提高實驗質(zhì)量。希望同學(xué)們在閱讀代碼的時候，注意學(xué)習(xí)方法，在最短的時間內(nèi)熟悉MatLab，提高應(yīng)用能力。示例中有些部分是實驗項目中的內(nèi)容實現(xiàn)，有些是一些典型例題的實現(xiàn)。研究示例代碼，倡導(dǎo)個性化編程是我們的目標(biāo)，希望同學(xué)們能在在進(jìn)行實驗項目的過程中提高M(jìn)atLab 的應(yīng)用能力；在學(xué)習(xí)MatLab 編程的同時加強(qiáng)對數(shù)字信號處理有關(guān)實驗項目的理解。
以下代碼段均在 MatLab5.3 下調(diào)試通過，但是由于排版或其他一些原因，可能有部分代
碼段不能得到正常結(jié)果。您可以在“http://202.38.75.33/dsp/matlab/”得到本講義的修訂內(nèi)容，同時可以在這個網(wǎng)址獲取所有代碼。
附錄 C1 信號、系統(tǒng)和系統(tǒng)響應(yīng)
1、理想采樣信號序列
（1）首先產(chǎn)生信號x(n),0<=n<=50
n=0:50; %定義序列的長度是50
A=444.128; %設(shè)置信號有關(guān)的參數(shù)
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 %清除已經(jīng)繪制的x(n)圖形
subplot(3,1,1);stem(x); %繪制x(n)的圖形
title(‘理想采樣信號序列’); %設(shè)置結(jié)果圖形的標(biāo)題
（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）改變參數(shù)為： 1, 0.4, 2.0734, 1 0 A = α = Ω = T =
n=0:50; %定義序列的長度是50
A=1; a=0.4; w0=2.0734; T=1; %設(shè)置信號有關(guān)的參數(shù)和采樣率T
x=A*exp(-a*n*T).*sin(w0*n*T); %pi 是MATLAB 定義的π，信號乘可采用“.*”