【正文】 信。當DSP完成對音頻數(shù)據(jù)的處理以后, AIC23再把數(shù)字信號轉(zhuǎn)化為模擬信號, 這樣就能夠在立體聲輸出端或者耳機輸出端聽到聲音。簡單的講，可以在原聲音流中疊加延遲一段時間后的聲流，實現(xiàn)回聲效果。另一方面，為了語音識別實用化的需要，講者自適應(yīng)、聽覺模型、快速搜索識別算法以及進一步的語言模型的研究等課題倍受關(guān)注。 DSP系統(tǒng)特點和設(shè)計基本原則DSP（digital signal processor）是一種獨特的微處理器，是以數(shù)字信號來處理大量信息的器件。通常這些信號處理的過程要滿足實時且快速高效的要求，隨著DSP技術(shù)的發(fā)展，以DSP為內(nèi)核的設(shè)備越來越多，為語音信號的處理提供了良好的平臺。對這種語言環(huán)境也有了新的了解。H振幅39。h39。title(39。h=(1).^n.*sin((i).*Wc)./((i).*pi)。p0163。title(39。grid on。w1=。對此，采取兩項措施：一是將h（n）截短；二是將其右移。)。頻率39。[Bap,Aap]=zp2tf(Z,P,K)。wp1=wp*ts。)。figure(2)。)。209。%通帶最大衰減 rs=35。MATLAB提供了一組標準的數(shù)字濾波器設(shè)計函數(shù)，大大簡化了濾波器的設(shè)計工程。二、實驗原理在IIR濾波器的設(shè)計中，常用的方法是：先根據(jù)設(shè)計要求尋找一個合適的模擬原型濾波器，然后根據(jù)一定的準則將此模擬原型濾波器轉(zhuǎn)換為數(shù)字濾波器。相位部分39。振幅39。零極點圖39。x2=*exp(1*j*w)+*exp(2*j*w)+*exp(3*j*w)+*exp(4*j*w)+*exp(5*j*w)。MATLAB提供了求離散時間系統(tǒng)頻響特性的函數(shù)freqz，調(diào)用freqz的格式主要有兩種。在MATLAB中系統(tǒng)函數(shù)的零極點可通過函數(shù)roots得到，也可借助函數(shù)tf2zp得到，tf2zp的語句格式為[Z,P,K]=tf2zp(B,A)其中，B與A分別表示H（z）的分子與分母多項式的系數(shù)向量。subplot(2,1,2)。)。k=0:1:N1。)。subplot(224)。)。相位 part39。grid。magnitude39。imagX=imag(H)。X(e)=+(n)=cosjwjw++++（２）計算32點序列5pn16，0≦n≦31的32點和64點DFT，分別繪出幅度譜圖形，并繪出該序列的DTFT圖形。grid on。(2)N=10。n=0:N1。grid on。分析:該序列具有周期性,、編制程序產(chǎn)生復(fù)正弦序列x(n)=e(2+j學(xué)號后兩位）n，并繪出其圖形。phi=0。title(39。grid on。w=2*pi。246。y=[zeros(1,NM+1)ones(1,NM)]。title(39。單位階躍序列39。0編制程序產(chǎn)生單位階躍序列u(n)=237。stem(n,y)。title(39。1n=0編制程序產(chǎn)生單位沖激序列d(n)=237。與連續(xù)的正弦信號不同，正弦序列的自變量n必須為整數(shù)。離散序列通常用x（n）來表示，自變量必須是整數(shù)。在以后的學(xué)習(xí)中還需要深入了解這方面的內(nèi)容。頻率圖中，女生高 27頻成分較多。d:39。)。%顯示標表示的意義 grid。subplot(3,2,5)。%阻帶截至頻率 Rp=1。)。y0=fft(yy)。)。f=[(n11)*(2*pi/N)]/pi。%對錄入信號進行采樣[x1,fs,bits]=wavread(39。)。ylabel(39。f=[(n11)*(2*pi/N)]/pi。)。)。%暫停命令 fprintf(39。)。%原始信號的頻譜圖 title(39。)。%播放原始信號 N=length(y)。%顯示橫坐標表示的意義 ylabel(39。%快速傅立葉變換 subplot(2,2,2)。%錄制信號的時域波形title(39。%播放原始信號,解決落后半拍wavplay(x,fs)。pause。double39。%錄音時間 fprintf(39。%網(wǎng)格wavplay(ys1,8000)。)。)。%繪制后信號的時域的圖形 title(39。h=fir1(n,wn/pi,39。wc2=2*pi*3200/fs。%顯示標表示的意義title(39。%在四行兩列的第七個窗口顯示圖形 plot(f,abs(fftwave(n1)))。高通濾波后信號的時域波形39。%在四行兩列的第五個窗口顯示圖形 plot(ys)。wn=(ws+wc)/2。)。%繪制模值 xlabel(39。頻率w/pi39。%信號送入濾波器濾波，ys為輸出 fftwave=fft(ys)。width=wswc。%顯示縱坐標表示的意義 title(39。原始信號的頻譜圖39。%加標題 ylabel(39。%采樣間隔 n1=1:N/2。)。title(39。)。wp=2*pi*fp/Fs。)。title(39。high39。T=1。9 巴特沃思低通濾波器MATLAB程序如下。)。s39。ws=2*pi*fs/Fs。grid。截取后的聲音時域圖形39。%加標題grid。頻率w/pi39。%加網(wǎng)格y0=fft(x)。plot(x)。C:39。錄音結(jié)束39。%文字提示 pause。奈奎斯特采樣定理：要使實信號采樣后能夠不失真還原，采樣頻率必須大于信號最高頻率的兩倍：Ws179。?(jW)與X(jW)的關(guān)系：理想抽樣信號的頻譜是連續(xù)信號頻譜的Xaa周期延拓,重復(fù)周期為Ws(采樣角頻率)。xXa(jW)171。xa(nT)d(tnT)xn=165。希望以后每一次實習(xí)都可以學(xué)到更多更有用的知識。在老師的指導(dǎo)和幫助下，我們也都將之逐一地完善了。第一篇：數(shù)字信號處理實習(xí)感想數(shù)字圖像處理個人體會這次實習(xí)讓我學(xué)會了如何用MATLAB語句來實現(xiàn)一些功能，雖然用起來還不是很熟練，但收獲不少。此外，有些軟件總是會存在很多的漏洞，像是圖像處理完后下一次打開原圖像部分還會出現(xiàn)上一次處理的圖像之類的種種問題。或許這就是學(xué)習(xí)的精神吧！最后，我們要感謝院里給我們提供了這次實習(xí)的機會，感謝黃忠朝老師和玉罡老師指導(dǎo)我們完成了這次實習(xí)。xa(t)d(tnT)=229。?a(t)Xa(jW)171。TTk=165。Ws/2即連續(xù)信號是帶限的，且信號最高頻率不超過抽樣頻率的二分之一，則可不失真恢復(fù)。,time)。%錄制語音信號 fprintf(39。%按任意鍵播放語音信號wavwrite(x,Fs,39。subplot(2,1,1)。%顯示縱坐標的表示意義 grid。%加標題 xlabel(39。)。%畫出截取后的時域圖形 title(39。)。wp=2*pi*fp/Fs。[B,A]=butter(N,wc,39。巴特沃斯低通濾波器39。gridon。ws=2*pi*fs/Fs。[B,A]=butter(N,wc,39。plot(w,abs(h))。耗損（dB）39。rs=100。s39。plot(w,abs(h))。耗損（dB）39。%返回采樣點數(shù) df=fs/N。)。%原始信號的頻譜圖 title(39。)。ws=2*pi*1200/fs。ys=fftfilt(h,x)。%加標題 xlabel(39。%在四行兩列的第四個窗口顯示圖形 plot(f, abs(fftwave(n1)))。低通濾波器濾波后信號的頻譜圖39。width=wswc。%將信號送入高通濾波器濾波 subplot(4,2,5)。%顯示標表示的意義 title(39。%將濾波后的語音信號進行快速傅立葉變換 subplot(4,2,7)。)。wc1=2*pi*1000/fs。n=ceil(*pi)。%在四行兩列的第六個窗口顯示圖形 plot(ys1)。幅值/A39。頻率w/pi39。%加標題 grid。%抽樣頻率 time=3。x=wavrecord(time*Fs,Fs,39。)。)。plot(x)。%加網(wǎng)格y0=fft(x)。)。% 對語音信號進行采樣wavplay(y,fs)。原始男生信號的時域波形39。%在四行兩列的第三個窗口顯示圖形 plot(f,abs(y0(n1)))。幅值 39。%文字提示 pause。錄音結(jié)束39。C:39。%采樣間隔n1=1:N/2。grid。幅值39。)。%采樣間隔 n1=1:N/2。疊加信號的時域波形39。grid。頻率w/pi39。%通帶截至頻率 Ws=2*2000/fs。%將濾波后的語音信號進行快速傅立葉變換 figure(5)。)。頻率w/pi39。圖形如下：噪聲信號的時域波形1100噪聲信號的頻譜圖幅值/ 10疊加信號的時域波形10101時間/t2200幅值/db34幅值/，常會出現(xiàn)實際發(fā)出聲音落后錄制動作半拍，可在[x,fs,bits]=wavread(39。：時域圖的振幅大小與性別無關(guān)，只與說話人音量大小有關(guān)，音量越大，振幅越大。整個過程中學(xué)到了很多新的知識，特別是對Matlab的使用終于有些了解。二、實驗原理離散時間信號是指在離散時刻才有定義的信號，簡稱離散信號，或者序列。（4）正弦序列x（n）這里，正弦序列的參數(shù)都是實數(shù)。n=0：N1 x=exp（（a.+j*w）*n）stem(n,real(x))或stem(n,imag(x))三、實驗內(nèi)容及分析236。stem(n,x)。y=[zeros(1,NM+1)1 zeros(1,NM1)]。1n179。title(39。stem(n,x)。x=[zeros(1,N1)ones(1,N)]。2230。A=1。)。stem(n,x)。w=2*pi。grid on。)。(1)N=10。grid on。)。三、實驗內(nèi)容及分析（１）編程計算并畫出下面DTFT的實部、虛部、幅度和相位譜。realX=real(H)。ylabel(39。plot(w/pi,angX)。title(39。frequency in pi unit39。axis([0 11 1])。imaginary39。xn=cos(5*pi*n/16)。32點39。stem(n,xn)。二、實驗原理離散時間系統(tǒng)的系統(tǒng)函數(shù)定義為系統(tǒng)零狀態(tài)響應(yīng)的Z變化與激勵的Z變化之比。利用系統(tǒng)的頻率響應(yīng)，可以分析系統(tǒng)對各種頻率成分的響應(yīng)特性，并推出系統(tǒng)的特性（高通、低通、帶通、帶阻等）。x1=+*exp(1*j*w)+*exp(2*j*w)+*exp(3*j*w)+*exp(4*j*w)+*exp(5*j*w)。title(39。ylabel(39。title(39。本實驗通過使用MATLAB函數(shù)對數(shù)字濾波器進行設(shè)計和和實現(xiàn)，要求掌握IIR數(shù)字巴特沃斯濾波器、數(shù)字切比雪夫濾波器的設(shè)計原理、設(shè)計方法和設(shè)計步驟；能根據(jù)給定的濾波器指標進行濾波器設(shè)計；同時也加深學(xué)生對數(shù)字濾波器的常用指標和設(shè)計過程的理解。所以從模擬傳遞函數(shù)可直接通過代數(shù)置換得到數(shù)字濾波器的傳遞函數(shù)。%阻帶截止頻率 rp=1。233。s39。[P,Q]=freqs(H,W)。幅度39。ts=1/fs。[Z,P,K]=buttap(N)。xlabel(39。頻率39。但是將h（n）作為濾波器的脈沖響應(yīng)有兩個問題：一是它是無限長的；二是它是非因果的。Window=hamming(N+1)。)。freqz(b,1,512)。w163。i=i+(i==0)*eps。axis([0,8,])。ylabel(39。ylabel(39。從上次的信號系統(tǒng)實驗的初步使用到這一次的深入了解，有了更深刻的認識。在高度信息化的今天，語音信號處理是信息高速公路、多媒體技術(shù)、辦公自動化、現(xiàn)代通信及智能系統(tǒng)等新興領(lǐng)域應(yīng)用的核心技術(shù)之一。同時也是眾構(gòu)成思想交流和感情溝通的最主要的途徑。一方面，對聲學(xué)語音學(xué)統(tǒng)計模型的研究逐漸深入，魯棒的語音識別、基于語音段的建模方法及隱馬爾可夫模型與人工神經(jīng)網(wǎng)絡(luò)的結(jié)合成為研究的熱點。當已知一個數(shù)字音源后，可以利用計算機的處理能力，用數(shù)字的方式通過計算模擬回聲效應(yīng)。AIC23通過麥克風(fēng)輸入或者立體聲音頻輸入采集模擬信號, 并把模擬信號轉(zhuǎn)化為數(shù)字信號, 存儲到DSP的內(nèi)部RAM中,以便DSP處理。收發(fā)時鐘信號CLKX1和CLKR1由A IC23的串行數(shù)據(jù)輸入時鐘BCLK提供, 并由A IC23的幀同步信號LRCIN、LRCOUT啟動串口數(shù)據(jù)傳輸。用戶可以根據(jù)電阻阻值調(diào)節(jié)增益的大小，使語音輸入輸出達到最佳效果。 修改前程序圖 修改前程序圖 頻譜分析 左聲道及右聲道波形 實驗分析所以，從本實驗可知當已知一個數(shù)字音源后，可以利用計算機的處理能力，用數(shù)字的方式通過計算模擬回聲效應(yīng)。參考文獻[1][M].北京:中國水利水電出版社,2004.[2]王安民,陳明欣,[M].北京:清華大學(xué)出版社,2007 [3]彭啟琮,[M].成都:電子科技大學(xué)出版社,1997 [4]李宏偉,[J].解放軍理工大學(xué)學(xué)報,2001(1):41~44 [5][M].梁曉雯,裴小平,李玉虎,:清華大學(xué)出版社,2006 [6][M].北京:機械工業(yè)出版社,2003比較圖4和圖5,可以看到1200Hz以上的頻譜明顯得到了抑制。（2）學(xué)會簡單的矩陣輸入和數(shù)據(jù)讀寫。n163。a=(+3*i)*n。x=cos(100*pi*t).*sin(pi*t)。),title(39。b=[1 ]。二、實驗內(nèi)容實驗中用到的信號序列 a)高斯序列236。eansin(2pfn)xb(n)=237。0238。n4239。在xc(n)和xd(n)末尾補零，用N=32點FFT分析這兩個信號的幅頻特性，觀察幅頻特性發(fā)生了什么變化？兩種情況的FFT頻譜還有相同之處嗎？這些變化說明了什么？(5)用FFT分別實現(xiàn)xa(n)（p＝8，q＝2）和xb(n)（a＝，f＝）的16點循環(huán)卷積和線性卷積。)。)。subplot(2,3,2)。subplot(2,3,5)。(2)觀察雙線性變換及脈沖響應(yīng)不變法設(shè)計的濾波器的頻域特性，了解雙線性變換法及脈沖響應(yīng)不變法的特點。)。[B,A]=butter(3, OmegaC,39。fs=4000。OmegaC=2*fs*tan(pi*fc/fs)。r.39。實驗過程與結(jié)果（含實驗程序、運行的數(shù)據(jù)結(jié)果和圖形）；實驗四FIR數(shù)字濾波器的設(shè)計一、實驗?zāi)康?1)掌握用窗函數(shù)法，頻率采樣法及優(yōu)化設(shè)計法設(shè)計FIR濾波器的原理及方法，熟悉相應(yīng)的計算機編程；(2)熟悉線性相位FIR濾波器的幅頻特性和相頻特性；(3)了解各種不同窗函數(shù)對濾波器性能的影響。plot(w,abs(H))