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