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