《數(shù)字信號處理》第三版課后實驗答案西安電子科技大學(xué)出版社

1、程序清單及波形顯示：clccloseallclearall%======內(nèi)容1：調(diào)用filter解差分方程，由系統(tǒng)對u(n)的響應(yīng)判斷穩(wěn)定性======A=[10.9]B=[0.050.05]%系統(tǒng)差分方程系數(shù)向量B和Ax1n=[11111111zeros(150)]%產(chǎn)生信號x1(n)=R8(n)x2n=ones(1128)%產(chǎn)生信號x2(n)=u(n)y1n=filter(BAx1n)%求系統(tǒng)對x1(n)的響應(yīng)y1(n)n=0:le

2、ngth(y1n)1subplot(221)stem(ny1n.)title((a)系統(tǒng)對R8(n)的響應(yīng)y1(n))xlabel(n)ylabel(y1(n))y2n=filter(BAx2n)%求系統(tǒng)對x2(n)的響應(yīng)y2(n)n=0:length(y2n)1subplot(222)stem(ny2n.)title((b)系統(tǒng)對u(n)的響應(yīng)y2(n))xlabel(n)ylabel(y2(n))hn=impz(BA58)%求系統(tǒng)單

3、位脈沖響應(yīng)h(n)n=0:length(hn)1subplot(223)y=hnstem(nhn.)title((c)系統(tǒng)單位脈沖響應(yīng)h(n))xlabel(n)ylabel(h(n))020406000.20.40.60.8(a)三三三R8(n)三三三y1(n)ny1(n)05010015000.51(b)三三三u(n)三三三y2(n)ny2(n)020406000.050.1(c)三三三三三三三三h(n)nh(n)%===內(nèi)容2：調(diào)

4、用conv函數(shù)計算卷積============================x1n=[11111111]%產(chǎn)生信號x1(n)=R8(n)h1n=ones(110)h2n=[12.52.51]y21n=conv(h1nx1n)y22n=conv(h2nx1n)figure(2)n=0:length(h1n)1subplot(221)stem(nh1n)title((d)系統(tǒng)單位脈沖響應(yīng)h1n)xlabel(n)ylabel(h1(n))

5、n=0:length(y21n)1subplot(222)stem(ny21n)title((e)h1(n)與R8(n)的卷積y21n)xlabel(n)ylabel(y21(n))n=0:length(h2n)1subplot(223)stem(nh2n)title((f)系統(tǒng)單位脈沖響應(yīng)h2n)xlabel(n)ylabel(h2(n))n=0:length(y22n)1subplot(224)stem(ny22n)title((g

6、)h2(n)與R8(n)的卷積y22n)xlabel(n)ylabel(y22(n))程序清單及波形顯示：%DTMF雙頻撥號信號產(chǎn)生6位電話號碼%clearallclctm=[1236545666789674203568]%DTMF信號代表的16個數(shù)N=205K=[1820222431343842]f1=[697770852941]%行頻率向量f2=[1209133614771633]%列頻率向量TN=input(鍵入6位電話號碼=)%

7、輸入6位數(shù)字TNr=0%接收端電話號碼初值為零fl=1:6d=fix(TN10^(6l))TN=TNd10^(6l)fp=1:4fq=1:4iftm(pq)==abs(d)breakend%檢測碼相符的列號qendiftm(pq)==abs(d)breakend%檢測碼相符的行號pendn=0:1023%為了發(fā)聲，加長序列x=sin(2pinf1(p)8000)sin(2pinf2(q)8000)%構(gòu)成雙頻信號sound(x8000)%

8、發(fā)出聲音pause(0.1)%接收檢測端的程序X=goertzel(x(1:205)K1)%用Goertzel算法計算八點DFT樣本val=abs(X)%列出八點DFT向量subplot(32l)stem(Kval.)gridxlabel(k)ylabel(|X(k)|)%畫出DFT(k)幅度axis([10500120])limit=80%fs=5:8ifval(s)limitbreakend%查找列號endfr=1:4ifval(r