【正文】 分別打印其幅頻特性曲線。在MATLAB中提供了roots子函數(shù)，用于求多項(xiàng)式的根。調(diào)用格式：zplane(z，p)；繪制由列向量z確定的零點(diǎn)、列向量p確定的極點(diǎn)構(gòu)成的零極點(diǎn)分布圖。觀察分236。x(n+iN)]RN(n)i=165。165。四 思考題(1)比較線性卷積和循環(huán)卷積在序列長(zhǎng)度不同時(shí)兩者的聯(lián)系？(2)試著寫出循環(huán)卷積的源代碼？實(shí)驗(yàn)三 時(shí)域采樣理論與頻域采樣定理驗(yàn)證一、實(shí)驗(yàn)?zāi)康?時(shí)域采樣理論與頻域采樣理論是數(shù)字信號(hào)處理中的重要理論。循環(huán)卷積運(yùn)用到離散傅立葉變換的循環(huán)移位性質(zhì)，即時(shí)域循環(huán)移位定理。)。%顯示原連續(xù)信號(hào) title(39。復(fù)指數(shù)序列的虛部39。時(shí)間（n）39。(2)掌握MATLAB產(chǎn)生常用時(shí)域離散信號(hào)的方法。title(39。n=m*。+Qk。Pk162。)。Pxx12=abs(fft(w.*xn(353:384),Ns).^2)/norm(w)^2。%抽樣頻率,單位為kHZ n=0:N1。subplot(2,2,4)。功率(dB)39。plot(f1,Pxx1)。xn1=a*sin(2*pi*f*n1./fs)+wn1。 由于觀測(cè)數(shù)據(jù)有限，所以周期圖法估計(jì)分辨率低，估計(jì)誤差大。178。178。)[numa,dena]=impinvar(b,a,Fs)。經(jīng)濾波處理的心電圖信號(hào)的頻譜39。xlabel(39。)。[b,a] = butter(N,Wn,39。x(n)=[4,2,0,4,6,4,2,4,6,6,4,4,6,6,2,6,12,8,0,16,38,60,84,90,66,32,4,2,4,8,12,12,10,6,6,6,4,0,0,0,0,0,2,4,0,0,0,2,2,0,0,2,2,2,2,0] 對(duì)序列x(n)用FFT做頻譜分析，生成x(n)的頻譜圖。k39。k39。%subplot(2,2,3)。理論計(jì)算 X(k)=DFT[x(n)]2N。、難點(diǎn)、特點(diǎn)快速傅里葉變換（FFT）的原理及應(yīng)用。%subplot(2,2,1)。)。)。實(shí)驗(yàn)儀器及材料計(jì)算機(jī)，MATLAB軟件實(shí)驗(yàn)內(nèi)容及要求1．設(shè)計(jì)巴特沃斯低通數(shù)字濾波器對(duì)人體心電信號(hào)進(jìn)行濾波(1)人體心電圖信號(hào)在測(cè)量過(guò)程中會(huì)受到工業(yè)高頻干擾，所以必須經(jīng)過(guò)低通濾波處理，才能作為判斷心臟功能的有用信息。[N,Wn] = buttord(Wp,Ws,Rp,Rs,39。)。title(39。)。e=sqrt(10^(Rp/10)1)。a2=[2 3 1]。 [BZ,AZ] = impinvar(B,A,Fs)：沖激響應(yīng)不變法，返回值為數(shù)字濾波器的系數(shù)。n=0N120163。n4=0:N41。Pxx4=10*log10(abs(fft(xn4).^2)/N4)。頻率39。title(39。%數(shù)據(jù)長(zhǎng)度 Ns=32。Pxx9=abs(fft(w.*xn(257:288),Ns).^2)/norm(w)^2。)。H=P(P+1)239。p(1)=20。%產(chǎn)生和H（k）維數(shù)相同的單位矩陣； p(k)=(IH(k)*Ck)*p1(k)。,39。本講義在第三版的基礎(chǔ)上編寫了五個(gè)實(shí)驗(yàn)，前2個(gè)實(shí)驗(yàn)屬基礎(chǔ)性的驗(yàn)證性實(shí)驗(yàn)，第5個(gè)實(shí)驗(yàn)屬基本應(yīng)用綜合性實(shí)驗(yàn)。單位脈沖序列39。)。離散信號(hào)39。 subplot(2,1,2),plot(tn,f)。h(m)x(nm)N1MATLAB信號(hào)處理工具箱提供了conv函數(shù)，該函數(shù)用于計(jì)算兩個(gè)有限序列的卷積。因此取其主值序列就滿足=。[xa(t)229。jwnjwX(e)，即 上式的右邊就是序列的傅立葉變換?(jW)=X(ejw)Xaw=WT上式說(shuō)明理想采樣信號(hào)的傅立葉變換可用相應(yīng)的采樣序列的傅立葉變換得到，只要將自變量ω用WT代替即可。X(k)=FFT[x(n)]，k=0,1,2,3,,M1式中k代表的頻率為wk=2pkM。點(diǎn)），再對(duì)其進(jìn)行16點(diǎn)DFT(FFT),得到的就是四 思考題如果序列x(n)的長(zhǎng)度為M，希望得到其頻譜()jXeω在]2,0[π上的N點(diǎn)等間隔采樣，當(dāng)N實(shí)驗(yàn)四 離散系統(tǒng)的零極點(diǎn)分析一、實(shí)驗(yàn)?zāi)康?． 熟悉MATLAB的仿真及應(yīng)用環(huán)境2． 在MATLAB的環(huán)境下研究控制系統(tǒng)穩(wěn)定性二、實(shí)驗(yàn)內(nèi)容和要求。系統(tǒng)極點(diǎn)的位置對(duì)系統(tǒng)響應(yīng)有著非常明顯的影響。如果是模擬周期信號(hào)，也應(yīng)該選取整數(shù)倍周期的長(zhǎng)度，經(jīng)過(guò)采樣后形成周期序列，按照周期序列的譜分析進(jìn)行。x3(n)=cospn4選擇FFT的變換區(qū)間N為8和16 兩種情況分別對(duì)以上序列進(jìn)行頻譜分析。四 思考題1結(jié)合本次實(shí)驗(yàn)與書本上相關(guān)原理，對(duì)書本后面的習(xí)題進(jìn)行相關(guān)的matlab軟件仿真？2因果穩(wěn)定的離散系統(tǒng)必須滿足的充分必要條件是什么？MATLAB提供了哪些進(jìn)行零極點(diǎn)求解的子函數(shù)？如何使用？實(shí)驗(yàn)五 基于MATLAB的快速傅里葉變換一 實(shí)驗(yàn)?zāi)康膶W(xué)習(xí)用FFT對(duì)連續(xù)信號(hào)和時(shí)域離散信號(hào)進(jìn)行譜分析的方法，了解可能出現(xiàn)的誤差及其原因，以便正確應(yīng)用FFT。 功能：求多項(xiàng)式的根。13239。對(duì)比上面敘述的時(shí)域采樣原理和頻域采樣原理，得到一個(gè)有用的結(jié)論，這兩個(gè)采樣理論具有對(duì)偶性：“時(shí)域采樣頻譜周期延拓，頻域采樣時(shí)域信號(hào)周期延拓”。n=165。公式1165。n和m的變化區(qū)間均是[0,L1]，直接計(jì)算該式比較麻煩。離散信號(hào)39。在上述語(yǔ)句輸入完成之后，敲擊回車鍵，彈出圖形窗口，顯示出如下圖形，即已經(jīng)滿足題干所述條件，顯示了原連續(xù)信號(hào)和采樣獲得的離散信號(hào)波形(4)在matlab命令窗口中逐行輸入下列語(yǔ)句 f=200。Um=1。幅度x（n）39。(3)一個(gè)連續(xù)的周期性方波信號(hào)頻率為200Hz,信號(hào)幅度在1~+1V之間，在窗口上顯示2個(gè)周期信號(hào)波形，用Fs=4kHz的頻率對(duì)連續(xù)信號(hào)進(jìn)行采樣，試顯示原連續(xù)信號(hào)和采樣獲得的離散信號(hào)波形。fontsize39。,n,x,39。%設(shè)信號(hào)初值為0；H(1)=p1(1)*Ck39。(112a)(112b)(112c)(112d)假設(shè)初始條件Ak,Ck,Qk,Rk,yk,xk1,Pk1已知，其中x0=E[x0],P0=var[x0]。、自適應(yīng)濾波器和各種譜估計(jì)方法。Pxx15=abs(fft(w.*xn(449:480),Ns).^2)/norm(w)^2。xn=a*sin(2*pi*f1*n./fs)+wn。頻率39。功率譜Pxx，N=25639。)。wn2=randn(1,N2)。四、實(shí)驗(yàn)程序與實(shí)驗(yàn)結(jié)果（1）用周期圖法進(jìn)行譜估計(jì) A、實(shí)驗(yàn)程序： %用周期法進(jìn)行譜估計(jì) clear all。我認(rèn)識(shí)到了自己理論知識(shí)的不足，也認(rèn)識(shí)到了我們學(xué)習(xí)的基礎(chǔ)知識(shí)究竟能運(yùn)用于什么領(lǐng)域，如何運(yùn)用。 [n,Wn] = buttord(Wp,Ws,Rp,Rs,39。plot(w*Fs/pi,20*log10(abs(h)/norm))title(39。編程設(shè)計(jì)教材習(xí)題62，求模擬濾波器Ha(s)的表達(dá)式。subplot(2,1,2)stem(l,xx)。xlabel(39。[numa,dena]=impinvar(b,a,1)。jw|H(e)|。2*NFFT39。2*NFFT39。%subplot(2,2,2)。n2=0:1:2*N1。 編程環(huán)境。xk2=1/2*(wkconj(wk(mod(kk,N)+1)))/i。subplot(2,1,1)。subplot(2,1,2)。八、心得體會(huì)第三篇：數(shù)字信號(hào)處理實(shí)驗(yàn)二北京信息科技大學(xué)數(shù)字信號(hào)處理實(shí)驗(yàn)報(bào)告題 目：數(shù)字信號(hào)處理課程設(shè)計(jì)實(shí)驗(yàn)學(xué) 院： 信息與通信工程學(xué)院 專 業(yè)： 通信工程專業(yè) 姓 名： 班 級(jí)：學(xué) 號(hào)： 指導(dǎo)老師：實(shí)驗(yàn)?zāi)康氖煜IR數(shù)字濾波器的設(shè)計(jì)原理與方法。Ws=*pi。grid。xx=filter(b,a,x)。figure(5)subplot(2,1,1)plot(l,abs(y))。[N,Wn]=cheb1ord(Wp,Ws,Rp,Rs,39。源程序close all clear all clc b1=[0 0 1]。178。n=0x*(n)x(n+m)(11)229。%抽樣頻率,單位為kHZ n1=0:N11。f3=((0:length(Pxx3)1))/length(Pxx3)。subplot(2,2,2)。ylabel(39。)。Pxx6=abs(fft(w.*xn(161:192),Ns).^2)/norm(w)^2。)。x?k=+Hk()239。Ak=exp()。*inv(Ck*p1(k)*Ck39。,16)。四．實(shí)驗(yàn)結(jié)果五．實(shí)驗(yàn)小結(jié)通過(guò)卡