數(shù)字信號處理第三版課后實驗程序(高西全)

1、%實驗實驗1：系統(tǒng)響應(yīng)及系統(tǒng)穩(wěn)定性：系統(tǒng)響應(yīng)及系統(tǒng)穩(wěn)定性closeallclearall%======內(nèi)容內(nèi)容1：調(diào)用：調(diào)用filter解差分方程，由系統(tǒng)對解差分方程，由系統(tǒng)對u(n)的響應(yīng)判斷穩(wěn)定性的響應(yīng)判斷穩(wěn)定性======A=[10.9]B=[0.050.05]%系統(tǒng)差分方程系數(shù)向量系統(tǒng)差分方程系數(shù)向量B和Ax1n=[11111111zeros(150)]%產(chǎn)生信號產(chǎn)生信號x1(n)=R8(n)x2n=ones(1128)%產(chǎn)生信

2、號產(chǎn)生信號x2(n)=u(n)hn=impz(BA58)%求系統(tǒng)單位脈沖響應(yīng)求系統(tǒng)單位脈沖響應(yīng)h(n)subplot(221)y=h(n)tstem(hny)%調(diào)用函數(shù)調(diào)用函數(shù)tstem繪圖繪圖title((a)系統(tǒng)單位脈沖響應(yīng)系統(tǒng)單位脈沖響應(yīng)h(n))boxony1n=filter(BAx1n)%求系統(tǒng)對求系統(tǒng)對x1(n)的響應(yīng)的響應(yīng)y1(n)subplot(222)y=y1(n)tstem(y1ny)title((b)系統(tǒng)對系統(tǒng)對R

3、8(n)的響應(yīng)的響應(yīng)y1(n))boxony2n=filter(BAx2n)%求系統(tǒng)對求系統(tǒng)對x2(n)的響應(yīng)的響應(yīng)y2(n)subplot(224)y=y2(n)tstem(y2ny)title((c)系統(tǒng)對系統(tǒng)對u(n)的響應(yīng)的響應(yīng)y2(n))boxon%===內(nèi)容內(nèi)容2：調(diào)用：調(diào)用conv函數(shù)計算卷積函數(shù)計算卷積============================x1n=[11111111]%產(chǎn)生信號產(chǎn)生信號x1(n)=R8(

4、n)h1n=[ones(110)zeros(110)]h2n=[12.52.51zeros(110)]y21n=conv(h1nx1n)y22n=conv(h2nx1n)figure(2)subplot(221)y=h1(n)tstem(h1ny)%調(diào)用函數(shù)調(diào)用函數(shù)tstem繪圖繪圖title((d)系統(tǒng)單位脈沖響應(yīng)系統(tǒng)單位脈沖響應(yīng)h1(n))boxonsubplot(222)y=y21(n)tstem(y21ny)title((e)h

5、1(n)與R8(n)的卷積的卷積y21(n))boxonsubplot(223)y=h2(n)tstem(h2ny)%調(diào)用函數(shù)調(diào)用函數(shù)tstem繪圖繪圖title((f)系統(tǒng)單位脈沖響應(yīng)系統(tǒng)單位脈沖響應(yīng)h2(n))boxonsubplot(224)y=y22(n)tstem(y22ny)title((g)h2(n)與R8(n)的卷積的卷積y22(n))boxon%=========內(nèi)容內(nèi)容3：諧振器分析：諧振器分析===========

6、=============un=ones(1256)%產(chǎn)生信號產(chǎn)生信號u(n)n=0:255xsin=sin(0.014n)sin(0.4n)%產(chǎn)生正弦信號產(chǎn)生正弦信號A=[11.82370.9801]B=[1100.4901100.49]%系統(tǒng)差分方程系數(shù)向量系統(tǒng)差分方程系數(shù)向量B和Ay31n=filter(BAun)%諧振器對諧振器對u(n)的響應(yīng)的響應(yīng)y31(n)y32n=filter(BAxsin)%諧振器對諧振器對u(n)的響

7、應(yīng)的響應(yīng)y31(n)figure(3)subplot(211)y=y31(n)tstem(y31ny)title((h)諧振器對諧振器對u(n)的響應(yīng)的響應(yīng)y31(n))boxonsubplot(212)y=y32(n)tstem(y32ny)title((i)諧振器對正弦信號的響應(yīng)諧振器對正弦信號的響應(yīng)y32(n))boxontitle((e)32點頻域采樣)xlabel(k)ylabel(|X_3_2(k)|)axis([01602

8、00])n1=0:N1subplot(326)stem(n1x32n.)boxontitle((f)32點IDFT[X_3_2(k)])xlabel(n)ylabel(x_3_2(n))axis([032020])10.3.2實驗程序清單實驗程序清單%第10章實驗3程序exp3.m%用FFT對信號作頻譜分析clearallcloseall%實驗內(nèi)容(1)==========================================

9、=========x1n=[ones(14)]%產(chǎn)生序列向量x1(n)=R4(n)M=8xa=1:(M2)xb=(M2):1:1x2n=[xaxb]%產(chǎn)生長度為8的三角波序列x2(n)x3n=[xbxa]X1k8=fft(x1n8)%計算x1n的8點DFTX1k16=fft(x1n16)%計算x1n的16點DFTX2k8=fft(x2n8)%計算x1n的8點DFTX2k16=fft(x2n16)%計算x1n的16點DFTX3k8=fft

10、(x3n8)%計算x1n的8點DFTX3k16=fft(x3n16)%計算x1n的16點DFT%以下繪制幅頻特性曲線subplot(221)mstem(X1k8)%繪制8點DFT的幅頻特性圖title((1a)8點DFT[x_1(n)])xlabel(ωπ)ylabel(幅度)axis([0201.2max(abs(X1k8))])subplot(223)mstem(X1k16)%繪制16點DFT的幅頻特性圖title((1b)16點D

11、FT[x_1(n)])xlabel(ωπ)ylabel(幅度)axis([0201.2max(abs(X1k16))])figure(2)subplot(221)mstem(X2k8)%繪制8點DFT的幅頻特性圖title((2a)8點DFT[x_2(n)])xlabel(ωπ)ylabel(幅度)axis([0201.2max(abs(X2k8))])subplot(222)mstem(X2k16)%繪制16點DFT的幅頻特性圖tit

12、le((2b)16點DFT[x_2(n)])xlabel(ωπ)ylabel(幅度)axis([0201.2max(abs(X2k16))])subplot(223)mstem(X3k8)%繪制8點DFT的幅頻特性圖title((3a)8點DFT[x_3(n)])xlabel(ωπ)ylabel(幅度)axis([0201.2max(abs(X3k8))])subplot(224)mstem(X3k16)%繪制16點DFT的幅頻特性圖t