【正文】 入了解，有了更深刻的認(rèn)識。title(39。ylabel(39。plot(w/pi,magH)。ylabel(39。line([0,10],[0,0])。axis([0,8,])。*w)。i=i+(i==0)*eps。alpha=(N1)/2。w163。grid on。freqz(b,1,512)。w1=。)。b=fir1(N,ws/pi,Window)。Window=hamming(N+1)。窗函數(shù)設(shè)計FIR濾波器步驟如下：（1）給定理想頻率響應(yīng)的幅頻特性和相頻特性；（2）求理想單位脈沖響應(yīng)，在實際計算中，可對理想頻率響應(yīng)采樣。但是將h（n）作為濾波器的脈沖響應(yīng)有兩個問題：一是它是無限長的；二是它是非因果的。)。頻率39。)subplot(2,1,2)。xlabel(39。[H,W]=freqz(bz,az)。[Z,P,K]=buttap(N)。ws2=2*fs*tan(ws1/2)。ts=1/fs。ws=*pi。幅度39。xlabel(39。[P,Q]=freqs(H,W)。)。s39。246。233。217。%阻帶截止頻率 rp=1。三、實驗內(nèi)容及分析利用MATLAB編程方法或利用MATLAB中fdatool工具設(shè)計不同功能的IIR數(shù)字濾波器。所以從模擬傳遞函數(shù)可直接通過代數(shù)置換得到數(shù)字濾波器的傳遞函數(shù)。（2）雙線性變化法。本實驗通過使用MATLAB函數(shù)對數(shù)字濾波器進(jìn)行設(shè)計和和實現(xiàn)，要求掌握IIR數(shù)字巴特沃斯濾波器、數(shù)字切比雪夫濾波器的設(shè)計原理、設(shè)計方法和設(shè)計步驟；能根據(jù)給定的濾波器指標(biāo)進(jìn)行濾波器設(shè)計；同時也加深學(xué)生對數(shù)字濾波器的常用指標(biāo)和設(shè)計過程的理解。相位39。title(39。stem(w/pi,angx)。ylabel(39。stem(w/pi,magx)。title(39。magx=abs(x)。x1=+*exp(1*j*w)+*exp(2*j*w)+*exp(3*j*w)+*exp(4*j*w)+*exp(5*j*w)。三、實驗內(nèi)容與結(jié)果分析已知LTI離散時間系統(tǒng)，要求由鍵盤實現(xiàn)系統(tǒng)參數(shù)輸入，并繪出幅頻和相頻響應(yīng)曲線和零極點分布圖，進(jìn)而分析系統(tǒng)的濾波特性和穩(wěn)定性。利用系統(tǒng)的頻率響應(yīng)，可以分析系統(tǒng)對各種頻率成分的響應(yīng)特性，并推出系統(tǒng)的特性（高通、低通、帶通、帶阻等）。它的作用是在z平面上畫出單位圓、零點與極點。二、實驗原理離散時間系統(tǒng)的系統(tǒng)函數(shù)定義為系統(tǒng)零狀態(tài)響應(yīng)的Z變化與激勵的Z變化之比。64點39。stem(n,xn)。xn=cos(5*pi*n/16)。32點39。stem(n,xn)。xn=cos(5*pi*n/16)。)。imaginary39。xlabel(39。axis([0 11 1])。)。frequency in pi unit39。subplot(222)。title(39。)。plot(w/pi,angX)。幅度 part39。ylabel(39。grid。realX=real(H)。y=[1,]。三、實驗內(nèi)容及分析（１）編程計算并畫出下面DTFT的實部、虛部、幅度和相位譜。通過這樣的表示，可以將時域的離散序列映射到頻域以便于進(jìn)一步的處理。)。%學(xué)號01 x=a.^(n)。grid on。stem(n,x)。(1)N=10。)。)。w=1。grid on。stem(n,x)。w=2*pi。grid on。stem(n,x)。w=2*pi/1。)。x=A*cos(w*n+phi)。A=1。及232。2230。title(39。x=[zeros(1,N1)ones(1,N)]。grid on。stem(n,x)。(2)N=6。title(39。程序: 4(1)N=5。1n179。)。y=[zeros(1,NM+1)1 zeros(1,NM1)]。grid on。stem(n,x)。形。n=0：N1 x=exp（（a.+j*w）*n）stem(n,real(x))或stem(n,imag(x))三、實驗內(nèi)容及分析236。n=0：N1 x=A*cos（2*pi*f*n/Fs+phase）（5）單邊實指數(shù)序列x（n）長度為N的實指數(shù)序列x（n）可以通過下面的MATLAB命令實現(xiàn)。（4）正弦序列x（n）這里，正弦序列的參數(shù)都是實數(shù)。n=（N1）：N1 x=[zeros(1,N1)ones(1,N)]； stem(n,x)延遲的單位階躍序列可以使用類似于單位沖激序列的方法獲得。二、實驗原理離散時間信號是指在離散時刻才有定義的信號，簡稱離散信號，或者序列。八．參考文獻(xiàn)[1]數(shù)字信號處理教程(第三版)程佩青 清華大學(xué)出版社 [2]MATLAB信號處理 劉波 文忠 電子工業(yè)出版社 [3] 王宏 清華大學(xué)出版社[4]MATLAB基礎(chǔ)與編程入門 張威 西安電子科技大學(xué)出版社[5] 數(shù)字信號處理及其MATLAB實驗 趙紅怡 張常 化學(xué)工業(yè)出版社[6]MATLAB信號處理詳解 陳亞勇等 人民郵電出版社 [7] 數(shù)字信號處理錢同惠 機(jī)械工業(yè)出版社 29第三篇：數(shù)字信號處理實驗報告JIANGSUUNIVERSITY OF TECHNOLOGY數(shù)字信號處理實驗報告學(xué)院名稱： 電氣信息工程學(xué)院專 業(yè)：班 級： 姓 名： 學(xué) 號： 指導(dǎo)老師： 張維璽（教授）2013年12月20日實驗一 離散時間信號的產(chǎn)生一、實驗?zāi)康臄?shù)字信號處理系統(tǒng)中的信號都是以離散時間形態(tài)存在的，所以對離散時間信號的研究是數(shù)字信號的基本所在。整個過程中學(xué)到了很多新的知識，特別是對Matlab的使用終于有些了解。第2部分 課程設(shè)計總結(jié)通過本次課程設(shè)計，使我們對數(shù)字信號處理相關(guān)知識有了更深刻的理解，尤其是對各種濾波器的設(shè)計。：時域圖的振幅大小與性別無關(guān)，只與說話人音量大小有關(guān)，音量越大，振幅越大。[100 10000])，窗的長度可根據(jù)需要定義。圖形如下：噪聲信號的時域波形1100噪聲信號的頻譜圖幅值/ 10疊加信號的時域波形10101時間/t2200幅值/db34幅值/，常會出現(xiàn)實際發(fā)出聲音落后錄制動作半拍，可在[x,fs,bits]=wavread(39。)。頻率w/pi39。%繪制模值 title(39。)。低通濾波后信號的時域波形39。%將濾波后的語音信號進(jìn)行快速傅立葉變換 figure(5)。%buttord函數(shù)（n為階數(shù)，Wn為截至頻率）[num,den]=butter(N,Wn)。%通帶截至頻率 Ws=2*2000/fs。)。頻率w/pi39。%疊加信號的頻譜圖 title(39。grid。)。疊加信號的時域波形39。plot(yy,39。%采樣間隔 n1=1:N/2。%對噪聲信號進(jìn)行采樣yy=x+x1。)。)。幅值39。y0=fft(x1)。grid。plot(x1)。%采樣間隔n1=1:N/2。)。C:39。)。錄音結(jié)束39。x=wavrecord(time*Fs,Fs,39。%文字提示 pause。%抽樣頻率 time=3。幅值 39。%加標(biāo)題 xlabel(39。%在四行兩列的第三個窗口顯示圖形 plot(f,abs(y0(n1)))。)。原始男生信號的時域波形39。f=[(n11)*(2*pi/N)]/pi。% 對語音信號進(jìn)行采樣wavplay(y,fs)。%顯示縱坐標(biāo)表示的意義 grid。)。原始女生信號的頻譜圖39。%加網(wǎng)格y0=fft(x)。%加標(biāo)題 ylabel(39。plot(x)。%采樣間隔 n1=1:N/2。)。C:39。)。錄音結(jié)束39。x=wavrecord(time*Fs,Fs,39。%文字提示 pause。%抽樣頻率 time=3。%在Matlab中，函數(shù)sound可以對聲音進(jìn)行回放，其調(diào)用格式: sound(xhigh, Fs,bits)。%加標(biāo)題 grid。)。頻率w/pi39。%對濾波后的信號進(jìn)行快速傅立葉變換 subplot(4,2,8)。幅值/A39。%加標(biāo)題 xlabel(39。%在四行兩列的第六個窗口顯示圖形 plot(ys1)。[h1,w]=freqz(h)。n=ceil(*pi)。width=ws1wc1。wc1=2*pi*1000/fs。%標(biāo)題 grid。)。頻率w/pi39。%將濾波后的語音信號進(jìn)行快速傅立葉變換 subplot(4,2,7)。幅值/A39。%顯示標(biāo)表示的意義 title(39。)。%將信號送入高通濾波器濾波 subplot(4,2,5)。high39。width=wswc。%播放濾波后信號3．2高通濾波器濾波 fs=8000。低通濾波器濾波后信號的頻譜圖39。ylabel(39。%在四行兩列的第四個窗口顯示圖形 plot(f, abs(fftwave(n1)))。幅值/A39。%加標(biāo)題 xlabel(39。%在四行兩列的第二個窗口顯示圖形 plot(ys)。ys=fftfilt(h,x)。h=fir1(n,wn/pi,39。ws=2*pi*1200/fs。%加標(biāo)題 grid。)。頻率w/pi39。%原始信號的頻譜圖 title(39。%顯示縱坐標(biāo)的表示意義 grid。)。subplot(4,2,1)。%返回采樣點數(shù) df=fs/N。C:39。耗損（dB）39。xlabel(39。plot(w,abs(h))。)。s39。Fs1=1。rs=100。%帶通濾波fp=[1200,3000]。耗損（dB）39。xlabel(39。plot(w,abs(h))。)。[B,A]=butter(N,wc,39。was=2*tan(ws/2)。ws=2*pi*fs/Fs。Fs=22050。gridon。)。巴特沃斯低通濾波器39。figure(1)。[B,A]=butter(N,wc,39。was=2*tan(ws/2)。wp=2*pi*fp/Fs。fs=1200。)。頻率39。%畫出截取后的時域圖形 title(39。x=wavread(39。)。幅值 39。%加標(biāo)題 xlabel(39。plot(f,abs(y0(n1)))。%顯示縱坐標(biāo)的表示意義 grid。)。subplot(2,1,1)。%返回采樣點數(shù) df=fs/N。%按任意鍵播放語音信號wavwrite(x,Fs,39。按Enter鍵回放錄音39。%錄制語音信號 fprintf(39。)。,time)。程序如下：Fs=8000。Ws/2即連續(xù)信號是帶限的，且信號最高頻率不超過抽樣頻率的二分之一，則可不失真恢復(fù)。Xa(jW)=237。TTk=165。2p1165。?a(t)Xa(jW)171。165。xa(t)d(tnT)=229。=fft(x)對采集的信號做快速傅立葉變換，并用[h1,w]=freqz(h)進(jìn)行DTFT變換。實驗課主要是以MATLAB為平臺，充分利用MATLAB的數(shù)字信號處理各種功能讓學(xué)生親自動手將課堂所學(xué)進(jìn)行仿真實現(xiàn)。數(shù)字信號處理課程的特點是課程本身理論性強(qiáng)、公式推導(dǎo)較多、概念比較抽象，使我們感到有枯燥難學(xué)之感。三、課程評價和建議我們的數(shù)字信號處理課是羅老師教的，羅老師有豐富工作的經(jīng)驗，對于這門課的實際用途很了解，另外羅老師本身就很幽默，對于這門課采用多種教學(xué)方法，豐富教學(xué)內(nèi)容，偶爾給我們講些生活上的問題，吸引學(xué)生對課程的關(guān)注。對于電氣信息類專業(yè)的學(xué)生來說，這些內(nèi)容是學(xué)習(xí)后續(xù)專業(yè)課程的重要基礎(chǔ)，也是實際工作中必不可少的專業(yè)基礎(chǔ)知識?？梢詮氖掠布矫?，比如說PCB，別小看這門技術(shù)，平時我們在試驗時制作的簡單，這一技術(shù)難點就在于板的層數(shù)越多，要做的越穩(wěn)定就越難，這可是非常有難度的，如果學(xué)好了學(xué)精了，也是非常好找工作的。通過該專業(yè)的學(xué)習(xí)，可以掌握電話網(wǎng)、廣播電視網(wǎng)、互聯(lián)網(wǎng)等各種通信系統(tǒng)的原理，研究提高信息傳送速度的技術(shù)，根據(jù)實際需要設(shè)計新的通信系統(tǒng)，開發(fā)可迅速準(zhǔn)確地傳送各種信息的通信工具等。第四單元的課程我們重點理解基2 FFT算法——時域抽取法﹑頻域抽取法，F(xiàn)FT的編程方法，分裂基FFT算法。信息科學(xué)是研究信息的獲取、傳輸、處理和利用的一門科學(xué)，信息要用一定形式的信號來表示，才能被傳輸、處理、存儲、顯示和利用，可以說，信號是信息的表現(xiàn)形式。數(shù)字信號處理技術(shù)正飛速發(fā)展，它不但自成一門學(xué)科，更是以不同形式影響和滲透到其他學(xué)科：它與國民經(jīng)濟(jì)息息相關(guān)，與國防建設(shè)緊密相連；它影響或改變著我們的生產(chǎn)、生活方式，因此受到人們普遍的關(guān)注。第三單元的課程我們學(xué)習(xí)了離散傅立葉變換定義和性質(zhì)，離散傅立葉變換應(yīng)用——快速卷積，頻譜分析。第七單元的課程我們學(xué)習(xí)了線性相位有限脈沖響應(yīng)（FIR）數(shù)字濾波器，窗函數(shù)法設(shè)計有限脈沖響應(yīng)（FIR）數(shù)字濾波器，頻率采樣法設(shè)計有限脈沖響應(yīng)（FIR）數(shù)字濾波器二、專業(yè)認(rèn)識和未來規(guī)劃通信工程是一門工程學(xué)科，主要是在掌握通信基本理論的基礎(chǔ)上，運用各種工程方法對通信中的一些實際問題進(jìn)行處理。當(dāng)然,就業(yè)形勢每年都會變化，所以關(guān)鍵還是要看自己。數(shù)字信號處理課程屬于專業(yè)基礎(chǔ)課，所涵蓋的內(nèi)容主要有：離散時間信號與系統(tǒng)的基本概念及描述方法，離散傅立葉變換及快速傅立葉變換，數(shù)字濾波器結(jié)構(gòu)及設(shè)計等。因此，加強(qiáng)該課程的建設(shè)具有重要的意義。另外我個人觀點是大學(xué)主要是培養(yǎng)自己的自學(xué)能力，老師只是個引導(dǎo)者，所以學(xué)習(xí)效果如何關(guān)鍵看自己的對學(xué)習(xí)的態(tài)度和付出程度。結(jié)合專業(yè)現(xiàn)狀，將課堂教學(xué)一部分變?yōu)槎嗝襟w教學(xué)，盡量將一些理論分析用圖形手段展示出來，以增強(qiáng)我們的感性認(rèn)識。二．設(shè)計內(nèi)容錄制一段個人自己的語音信號，并對錄制的信號進(jìn)行采樣；畫出采樣后語音信號的時域波形和頻譜圖；給定濾波器的性能指標(biāo)，采用窗函數(shù)法和雙線性變換法設(shè)計濾波器，并畫出濾波器的頻率響應(yīng)；然后用自己設(shè)計的濾波器對采集的信號進(jìn)行濾波，畫出濾波后信號的時域波形和頻譜，并對濾波前后的信號進(jìn)行對比，分析信號的變化；回放語音信號；換一個與你性別相異的人錄制同樣一段語音內(nèi)容，分析兩段內(nèi)容相同的語音信號頻譜之間有什么特點；再錄制一段同樣長時間的背景噪聲疊加到你的語音信號中，分析疊加前后信號頻譜的變化，設(shè)計一個合適的濾波器，能夠把該噪聲濾除；三．設(shè)計原理，利用函數(shù)wavrecord()，wavwrite(),wavread(),wavplay()對語音信號進(jìn)行錄制，存儲，讀取，回放。連續(xù)信號經(jīng)理想抽樣后時域、頻域發(fā)生的變化（理想抽樣信號與連續(xù)信號頻譜之間的關(guān)系）理想抽樣信號能否代表原始信號、如何不失真地還原信號即由離散信號恢復(fù)連續(xù)信號的條件（抽樣定理）理想采樣過程描述： 時域描述：?a(t)=xa(t)dT(t)=229。165。d(tnT)16