【正文】 plot(w,abs(H))。r.39。fs=4000。)。subplot(2,3,5)。)。在xc(n)和xd(n)末尾補(bǔ)零，用N=32點(diǎn)FFT分析這兩個(gè)信號(hào)的幅頻特性，觀察幅頻特性發(fā)生了什么變化？?jī)煞N情況的FFT頻譜還有相同之處嗎？這些變化說明了什么？(5)用FFT分別實(shí)現(xiàn)xa(n)（p＝8，q＝2）和xb(n)（a＝，f＝）的16點(diǎn)循環(huán)卷積和線性卷積。0238。二、實(shí)驗(yàn)內(nèi)容實(shí)驗(yàn)中用到的信號(hào)序列 a)高斯序列236。),title(39。a=(+3*i)*n。（2）學(xué)會(huì)簡(jiǎn)單的矩陣輸入和數(shù)據(jù)讀寫。 修改前程序圖 修改前程序圖 頻譜分析 左聲道及右聲道波形 實(shí)驗(yàn)分析所以，從本實(shí)驗(yàn)可知當(dāng)已知一個(gè)數(shù)字音源后，可以利用計(jì)算機(jī)的處理能力，用數(shù)字的方式通過計(jì)算模擬回聲效應(yīng)。收發(fā)時(shí)鐘信號(hào)CLKX1和CLKR1由A IC23的串行數(shù)據(jù)輸入時(shí)鐘BCLK提供, 并由A IC23的幀同步信號(hào)LRCIN、LRCOUT啟動(dòng)串口數(shù)據(jù)傳輸。當(dāng)已知一個(gè)數(shù)字音源后，可以利用計(jì)算機(jī)的處理能力，用數(shù)字的方式通過計(jì)算模擬回聲效應(yīng)。同時(shí)也是眾構(gòu)成思想交流和感情溝通的最主要的途徑。從上次的信號(hào)系統(tǒng)實(shí)驗(yàn)的初步使用到這一次的深入了解，有了更深刻的認(rèn)識(shí)。ylabel(39。i=i+(i==0)*eps。freqz(b,1,512)。Window=hamming(N+1)。頻率39。[Z,P,K]=buttap(N)。幅度39。s39。%阻帶截止頻率 rp=1。本實(shí)驗(yàn)通過使用MATLAB函數(shù)對(duì)數(shù)字濾波器進(jìn)行設(shè)計(jì)和和實(shí)現(xiàn)，要求掌握IIR數(shù)字巴特沃斯濾波器、數(shù)字切比雪夫?yàn)V波器的設(shè)計(jì)原理、設(shè)計(jì)方法和設(shè)計(jì)步驟；能根據(jù)給定的濾波器指標(biāo)進(jìn)行濾波器設(shè)計(jì)；同時(shí)也加深學(xué)生對(duì)數(shù)字濾波器的常用指標(biāo)和設(shè)計(jì)過程的理解。ylabel(39。x1=+*exp(1*j*w)+*exp(2*j*w)+*exp(3*j*w)+*exp(4*j*w)+*exp(5*j*w)。二、實(shí)驗(yàn)原理離散時(shí)間系統(tǒng)的系統(tǒng)函數(shù)定義為系統(tǒng)零狀態(tài)響應(yīng)的Z變化與激勵(lì)的Z變化之比。32點(diǎn)39。imaginary39。frequency in pi unit39。plot(w/pi,angX)。realX=real(H)。)。(1)N=10。grid on。stem(n,x)。A=1。x=[zeros(1,N1)ones(1,N)]。title(39。y=[zeros(1,NM+1)1 zeros(1,NM1)]。n=0：N1 x=exp（（a.+j*w）*n）stem(n,real(x))或stem(n,imag(x))三、實(shí)驗(yàn)內(nèi)容及分析236。二、實(shí)驗(yàn)原理離散時(shí)間信號(hào)是指在離散時(shí)刻才有定義的信號(hào)，簡(jiǎn)稱離散信號(hào)，或者序列。：時(shí)域圖的振幅大小與性別無關(guān)，只與說話人音量大小有關(guān)，音量越大，振幅越大。頻率w/pi39。%將濾波后的語音信號(hào)進(jìn)行快速傅立葉變換 figure(5)。頻率w/pi39。疊加信號(hào)的時(shí)域波形39。)。grid。C:39。%文字提示 pause。%在四行兩列的第三個(gè)窗口顯示圖形 plot(f,abs(y0(n1)))。% 對(duì)語音信號(hào)進(jìn)行采樣wavplay(y,fs)。%加網(wǎng)格y0=fft(x)。)。x=wavrecord(time*Fs,Fs,39。%加標(biāo)題 grid。幅值/A39。n=ceil(*pi)。)。%顯示標(biāo)表示的意義 title(39。width=wswc。%在四行兩列的第四個(gè)窗口顯示圖形 plot(f, abs(fftwave(n1)))。ys=fftfilt(h,x)。)。)。耗損（dB）39。s39。耗損（dB）39。[B,A]=butter(N,wc,39。gridon。[B,A]=butter(N,wc,39。)。)。%顯示縱坐標(biāo)的表示意義 grid。%按任意鍵播放語音信號(hào)wavwrite(x,Fs,39。,time)。TTk=165。xa(t)d(tnT)=229。此外，有些軟件總是會(huì)存在很多的漏洞，像是圖像處理完后下一次打開原圖像部分還會(huì)出現(xiàn)上一次處理的圖像之類的種種問題。在老師的指導(dǎo)和幫助下，我們也都將之逐一地完善了。xa(nT)d(tnT)xn=165。?(jW)與X(jW)的關(guān)系：理想抽樣信號(hào)的頻譜是連續(xù)信號(hào)頻譜的Xaa周期延拓,重復(fù)周期為Ws(采樣角頻率)。%文字提示 pause。C:39。%加網(wǎng)格y0=fft(x)。%加標(biāo)題grid。grid。s39。9 巴特沃思低通濾波器MATLAB程序如下。high39。)。)。)。%加標(biāo)題 ylabel(39。%顯示縱坐標(biāo)表示的意義 title(39。%信號(hào)送入濾波器濾波，ys為輸出 fftwave=fft(ys)。%繪制模值 xlabel(39。wn=(ws+wc)/2。高通濾波后信號(hào)的時(shí)域波形39。%顯示標(biāo)表示的意義title(39。h=fir1(n,wn/pi,39。)。%網(wǎng)格wavplay(ys1,8000)。double39。%播放原始信號(hào),解決落后半拍wavplay(x,fs)。%快速傅立葉變換 subplot(2,2,2)。%播放原始信號(hào) N=length(y)。%原始信號(hào)的頻譜圖 title(39。%暫停命令 fprintf(39。)。ylabel(39。%對(duì)錄入信號(hào)進(jìn)行采樣[x1,fs,bits]=wavread(39。)。)。subplot(3,2,5)。)。頻率圖中，女生高 27頻成分較多。離散序列通常用x（n）來表示，自變量必須是整數(shù)。1n=0編制程序產(chǎn)生單位沖激序列d(n)=237。stem(n,y)。單位階躍序列39。y=[zeros(1,NM+1)ones(1,NM)]。w=2*pi。title(39。分析:該序列具有周期性,、編制程序產(chǎn)生復(fù)正弦序列x(n)=e(2+j學(xué)號(hào)后兩位）n，并繪出其圖形。n=0:N1。grid on。imagX=imag(H)。grid。)。)。)。在MATLAB中系統(tǒng)函數(shù)的零極點(diǎn)可通過函數(shù)roots得到，也可借助函數(shù)tf2zp得到，tf2zp的語句格式為[Z,P,K]=tf2zp(B,A)其中，B與A分別表示H（z）的分子與分母多項(xiàng)式的系數(shù)向量。x2=*exp(1*j*w)+*exp(2*j*w)+*exp(3*j*w)+*exp(4*j*w)+*exp(5*j*w)。振幅39。二、實(shí)驗(yàn)原理在IIR濾波器的設(shè)計(jì)中，常用的方法是：先根據(jù)設(shè)計(jì)要求尋找一個(gè)合適的模擬原型濾波器，然后根據(jù)一定的準(zhǔn)則將此模擬原型濾波器轉(zhuǎn)換為數(shù)字濾波器。%通帶最大衰減 rs=35。)。)。[Bap,Aap]=zp2tf(Z,P,K)。)。w1=。title(39。h=(1).^n.*sin((i).*Wc)./((i).*pi)。h39。對(duì)這種語言環(huán)境也有了新的了解。 DSP系統(tǒng)特點(diǎn)和設(shè)計(jì)基本原則DSP（digital signal processor）是一種獨(dú)特的微處理器，是以數(shù)字信號(hào)來處理大量信息的器件。簡(jiǎn)單的講，可以在原聲音流中疊加延遲一段時(shí)間后的聲流，實(shí)現(xiàn)回聲效果。DX1和DR1分別與A IC23 的D IN 和DOUT 相連, 從而完成VC5509與AIC23間的數(shù)字信號(hào)通信。簡(jiǎn)單的講，可以在原聲音流中疊加延遲一段時(shí)間后的聲流，實(shí)現(xiàn)回聲效果。（3）掌握簡(jiǎn)單的繪圖命令。stem(x2)c)x(n)=3cos(+)+2sin(+)0163。cos39。(np)q239。0163。n=0:15。title(39。stem(abs(Xb),39。[num1,den1]=impinvar(B,A,fs)。OmegaC=2*pi*fc。)。(2)用矩形窗設(shè)計(jì)一個(gè)21階的線性相位低通FIR數(shù)字濾波器，截止頻率Wc＝，求出濾波器系數(shù)，并繪出濾波器的幅頻特性。stem(n,win)subplot(2,1,2)。plot(f,abs(h1),39。同一圖中畫兩條曲線： fc=1000。s39。Xb=fft(xb,16)。.39。(3)觀察三角波和反三角波序列的時(shí)域和幅頻特性，用N=8點(diǎn)FFT分析信號(hào)序列xc(n)和觀察兩者的序列形狀和頻譜曲線有什么異同？繪出兩序列及其幅頻特性xd(n)的幅頻特性，曲線。8n239。(4)應(yīng)用FFT實(shí)現(xiàn)序列的線性卷積和相關(guān)。x(t)39。x2=exp(a)。第五篇：數(shù)字信號(hào)處理實(shí)驗(yàn)報(bào)告南京郵電大學(xué)實(shí) 驗(yàn) 報(bào) 告實(shí)驗(yàn)名稱_____熟悉MATLAB環(huán)境 ___ 快速傅里葉變換及其應(yīng)用 ____IIR數(shù)字濾波器的設(shè)計(jì)_ FIR數(shù)字濾波器的設(shè)計(jì)課程名稱 數(shù)字信號(hào)處理A班級(jí)學(xué)號(hào)_______09002111___________ 姓 名 王都超開課時(shí)間 2011/2012學(xué)年，第 二 學(xué)期實(shí)驗(yàn)一熟悉MATLAB環(huán)境一、實(shí)驗(yàn)?zāi)康模?）熟悉MATLAB的主要操作命令。當(dāng)uDelay和uEffect的數(shù)值增大時(shí)，數(shù)字回聲的效果就會(huì)越加的明顯。SCLK和SDIN是AIC23控制端口的移位時(shí)鐘和數(shù)據(jù)輸入端,分別與VC5509的I2C模塊端口SCL和SDA相連。另外，生活中的回聲的成分比較復(fù)雜，有反射、漫反射、折射，還有回聲的多次反射、折射效果。語言是人類特有的功能，聲音是人類常用的工具，是相互傳遞信息的主要手段。實(shí)驗(yàn)心得體會(huì)：這次實(shí)驗(yàn)使我進(jìn)一步加深了對(duì)MATLAB軟件的使用。)。i=nalpha。b=fir1(N,ws/pi,Window)。程序1:海明窗設(shè)計(jì)N=48。xlabel(39。)。ylabel(39。1000hz[N,Wn]=cheb1ord(wp,ws,rp,rs,39。%通帶邊界頻率ws=*pi。數(shù)字濾波器你包括IIR（無限沖激響應(yīng)）和FIR（有限沖激響應(yīng)）型，在設(shè)計(jì)時(shí)通常采用不同的方法。)。w=[0:20:500]*pi/500。學(xué)會(huì)運(yùn)用MATLAB分析離散時(shí)間系統(tǒng)的系統(tǒng)函數(shù)的零極點(diǎn)；學(xué)會(huì)運(yùn)用MATLAB分析系統(tǒng)函數(shù)的零極點(diǎn)分布與其時(shí)域特性的關(guān)系；學(xué)會(huì)運(yùn)用MATLAB進(jìn)行離散時(shí)間系統(tǒng)的頻率特性分析。title(39。ylabel(39。xlabel(39。subplot(223)。angX=angle(H)。指數(shù)序列39。其中a=學(xué)號(hào)后兩位、a=1/“學(xué)號(hào)后兩位”。)。x=A*cos(w*n+phi)。程序:(1)N=5。%學(xué)號(hào)01 n=(N1):N1。stem(n,x)。%學(xué)號(hào)01 n=(N1):N1。與連續(xù)負(fù)指數(shù)信號(hào)一樣，我們將負(fù)指數(shù)序列實(shí)部和虛部的波形分開討論，得到如下結(jié)論：1）當(dāng)a0時(shí)，負(fù)指數(shù)序列x（n）的實(shí)部和虛部分別是按指數(shù)規(guī)律增長(zhǎng)的正弦振蕩序列；2)當(dāng)a3)當(dāng)a=0時(shí)，負(fù)指數(shù)序列x（n）即為虛指數(shù)序列，其實(shí)部和虛部分別是等幅的正弦振蕩序列；長(zhǎng)度為N的實(shí)指數(shù)序列x（n）可以通過下面的MATLAB命令實(shí)現(xiàn)。通過本實(shí)驗(yàn)，學(xué)生將學(xué)習(xí)如何用MATLAB產(chǎn)生一些常見的離散時(shí)間信號(hào)，實(shí)現(xiàn)信號(hào)的卷積運(yùn)算，并通過MATLAB中的繪圖工具對(duì)產(chǎn)生的信號(hào)進(jìn)行觀察，加深對(duì)常用離散信號(hào)和信號(hào)卷積和運(yùn)算的理解。通過帶通濾波器后，聲音比較適中。%標(biāo)題 xlabel(39。%信號(hào)送入濾波器濾波，ys為輸出 fftwave=fft(ys)。xlabel(39。%信號(hào)的時(shí)域波形title(39。 C: 39。)。%按任意鍵播放語音信號(hào) wavwrite(x,Fs,39。,time)。%快速傅立葉變換subplot(2,2,4)。)。%顯示縱坐標(biāo)的表示意義 grid。 C: 39。)。)。ylabel(39。wn=[wn1 wn2]。幅值/A39。)。wc=2*4800/fs。%網(wǎng)格subplot(4,2,4)。[h1,w]=freqz(h)。幅值 39。原始信號(hào)的時(shí)域波形39。ylabel(39。[N,wc]=buttord(wap,was,rp,rs,39。ylabel(39。)。)。)。振幅39。原始信號(hào)的頻譜圖39。)。%暫停命令wavplay(x,Fs)。按Enter鍵錄音%ds39。Xa(jWjkWs)Tk=165。連續(xù)信號(hào)經(jīng)理想抽樣后時(shí)域、頻域發(fā)生的變化（理想抽樣信號(hào)與連續(xù)信號(hào)頻譜之間的關(guān)系）理想抽樣信號(hào)能否代表原始信號(hào)、如何不失真地還原信號(hào)即由離散信號(hào)恢復(fù)連續(xù)信號(hào)的條件（抽樣定理）理想采樣過程描述： 時(shí)域描述：?a(t)=xa(t)dT(t)=229。剛開始的時(shí)候由于我對(duì)MATLAB語言的不熟練，所以還不能很容易地就把每個(gè)函數(shù)建立起來，不過經(jīng)過了一段實(shí)際操作后，只要明白了圖像處理的方法做起來都很順手了！我們完成了數(shù)字圖像處理軟件的功能后，并沒有滿足，還拓展了一些功能，比如說圖像的打開和保存，對(duì)于一些有影響到處理效果的參數(shù)還設(shè)置了輸入框。在實(shí)習(xí)過程中，我們并沒有很嚴(yán)格的計(jì)劃工作時(shí)間，但整體效率還是很高的，相對(duì)于別的組來說，我們完成得算比較快的。n=165。如果：236。%暫停命令 fprintf(39。)。%快速傅立葉變換 figure(2)。%加網(wǎng)格 原始信號(hào)的時(shí)域與頻譜圖針對(duì)實(shí)際發(fā)出聲音落后錄制動(dòng)作半拍的現(xiàn)象，如何拔除對(duì)無效點(diǎn)的采樣的問題: 出現(xiàn)這種現(xiàn)象的原因主要是錄音開始時(shí)，人的反應(yīng)慢了半拍，導(dǎo)致出現(xiàn)了一些無效點(diǎn)，而后而出現(xiàn)的無效的點(diǎn)，主要是已經(jīng)沒有聲音的動(dòng)作，先讀取聲音出來，將原始語音信號(hào)時(shí)域波形圖畫出來，根據(jù)己得到的信號(hào)，可以在第二次讀取聲音的后面設(shè)定采樣點(diǎn)，取好有效點(diǎn)，畫出濾除無效點(diǎn)后的語音信號(hào)時(shí)域波形圖，對(duì)比可以看出。%畫網(wǎng)格 去除無效點(diǎn)MATLAB程序如下。)。%高通濾波fp=4800。,39。grid on。[B,A]=butter(N,wc,39。grid on。幅值/A39。原始信號(hào)的頻譜圖39。%將濾波后的語音信號(hào)進(jìn)行快速傅立葉變換 figure(4)。頻率w/pi39。n=ceil(*pi/width)。)。高通濾波器濾波后信號(hào)的頻譜圖39。band39。%顯示縱坐標(biāo)表示的意義 grid。%播放濾波后信號(hào) 圖形如下：原始信號(hào)的時(shí)域波形幅值/A01012x 10原始信號(hào)的頻譜圖34幅值/分析：三個(gè)濾波器濾波后的聲音與原來的聲音都發(fā)生了變化。)。%播放原始信號(hào) N=length(x)。plot(f,abs(y0(n1)))。%返回采樣點(diǎn)數(shù) df=fs/N。原始男生信號(hào)的頻譜圖39。錄音中......39。%存儲(chǔ)語音信號(hào) 對(duì)噪聲進(jìn)行頻譜的分析[x1,fs,bits]=wavread(39。幅值/A39。 C: 39。xlabel(39。ylabel(39。plot(ys)。y