【正文】 5。dT(t)則有?(jW)=1X(jW)*D(jW)XaaT2p1165。Xa(jWjkWs)Tk=165。X(jW)239。0WWs/2W179。2fh四．具體實(shí)現(xiàn)在MATLAB中用wavrecord、wavread、wavplay、wavwrite對聲音進(jìn)行錄制、讀取、回放、存儲。按Enter鍵錄音%ds39。錄音中......39。)。%文字提示 fprintf(39。%暫停命令wavplay(x,Fs)。%存儲語音信號N=length(x)。%頻帶寬度 figure(2)。原始信號的時(shí)域波形39。)。subplot(2,1,2)。)。%顯示橫坐標(biāo)表示的意義 ylabel(39。原始信號的頻譜圖39。這樣就可以解決這個(gè)問題。%從4000點(diǎn)截取到24000結(jié)束 plot(x)。%標(biāo)題 xlabel(39。振幅39。%低通濾波fp=1000。rs=100。wap=2*tan(wp/2)。)。[Bz,Az]=bilinear(B,A,Fs1)。title(39。頻率（HZ）39。)。fs=5000。wp=2*pi*fp/Fs。wap=2*tan(wp/2)。)。s39。[h,w]=freqz(Bz,Az,512,Fs1*22050)。)。ylabel(39。MATLAB程序如下。rp=1。T=1。[N,wc]=buttord(wap,was,rp,rs,39。s39。[h,w]=freqz(Bz,Az,512,Fs1*1000)。)。ylabel(39。123．將聲音信號送入濾波器濾波x=wavread(39。%播放原始信號 N=length(x)。%頻帶寬度 figure(4)。原始信號的時(shí)域波形39。)。plot(f,abs(y0(n1)))。%加標(biāo)題 xlabel(39。幅值 39。)。wc=2*pi*1000/fs。n=ceil(*pi /width)。[h1,w]=freqz(h)。subplot(4,2,2)。)。ylabel(39。%網(wǎng)格subplot(4,2,4)。)。%顯示標(biāo)表示的意義title(39。%加網(wǎng)格wavplay(ys,8000)。wc=2*4800/fs。h=fir1(n,wn/pi, 39。ys=fftfilt(h,x)。頻率w/pi39。)。%標(biāo)題 ylabel(39。%網(wǎng)格fftwave=fft(ys)。xlabel(39。幅值/A39。)。beta=。ws2=2*pi*3000/fs。wn=[wn1 wn2]。,kaiser(n+1,beta))。subplot(4,2,6)。)。ylabel(39。%網(wǎng)格fftwave=fft(ys1)。xlabel(39。幅值/A39。)。其中低通的濾波后與原來聲音沒有很大的變化，其它兩個(gè)都又明顯的變化sound(xlow,Fs,bits)。 錄制一段異性的聲音進(jìn)行頻譜分析Fs=8000。,time)。)。%錄制語音信號 fprintf(39。按Enter鍵回放錄音39。%按任意鍵播放語音信號wavwrite(x,Fs,39。 C: 39。%返回采樣點(diǎn)數(shù) df=fs/N。subplot(2,2,1)。)。%顯示縱坐標(biāo)的表示意義 grid。%原始信號的頻譜圖 title(39。頻率w/pi39。)。)。%采樣間隔 n1=1:N/2。%錄制信號的時(shí)域波形title(39。幅值/A39。%快速傅立葉變換subplot(2,2,4)。)。%顯示橫坐標(biāo)表示的意義 ylabel(39。%加網(wǎng)格分析：就時(shí)域圖看，男生的時(shí)域圖中振幅比女生的高，對于頻譜圖女生的高頻成分比較多Fs=8000。,time)。)。%錄制語音信號fprintf(39。按Enter鍵回放錄音39。%按任意鍵播放語音信號 wavwrite(x,Fs,39。 C: 39。%返回采樣點(diǎn)數(shù) df=fs/N。subplot(3,2,1)。)。)。%噪聲信號的頻譜圖 ylabel(39。噪聲信號的頻譜圖39。 C: 39。)。%返回采樣點(diǎn)數(shù) df=fs/N。subplot(3,2,3)。%信號的時(shí)域波形title(39。時(shí)間/t39。)。plot(f,abs(y0(n1)))。xlabel(39。幅值/db39。%采樣頻率 Wp=2*1000/fs。%最小衰減[N,Wn]=buttord(Wp,Ws,Rp,Rs)。%信號送入濾波器濾波，ys為輸出 fftwave=fft(ys)。%信號的時(shí)域波形title(39。幅值/A39。plot(f, abs(fftwave(n1)))。%標(biāo)題 xlabel(39。幅值/A39。%播放濾波后信號 grid。d:39。通過帶通濾波器后，聲音比較適中。原因可能為噪聲與信號頻率相近無法濾除。再利用所學(xué)的操作，發(fā)現(xiàn)所寫的程序還是沒有能夠運(yùn)行，通過不斷地調(diào)試，運(yùn)行，最終得出了需要的結(jié)果。也認(rèn)識到自己的不足，以后要進(jìn)一步學(xué)習(xí)。通過本實(shí)驗(yàn)，學(xué)生將學(xué)習(xí)如何用MATLAB產(chǎn)生一些常見的離散時(shí)間信號，實(shí)現(xiàn)信號的卷積運(yùn)算，并通過MATLAB中的繪圖工具對產(chǎn)生的信號進(jìn)行觀察，加深對常用離散信號和信號卷積和運(yùn)算的理解。n=（N1）：N1 x=[zeros(1,N1)1 zeros(1,N1)]； stem(n,x)延遲K個(gè)采樣點(diǎn)的長度為N的單位沖激序列δ（nk）（kn=0：N1 y=[zeros(1,M)1 zeros(1,NM1)]； stem（n，y）（2）單位階躍序列u（n）如果u（n）在時(shí)間軸上延遲了k個(gè)單位，得到u（nk），即長度為N的單位階躍序列u（n）可以通過下面的MATLAB命令獲得。因此，用MATLAB表示矩形序列可利用上面的單位階躍序列組合而成。長度為N的正弦序列x（n）可以通過下面的MATLAB命令獲得。與連續(xù)負(fù)指數(shù)信號一樣，我們將負(fù)指數(shù)序列實(shí)部和虛部的波形分開討論，得到如下結(jié)論：1）當(dāng)a0時(shí)，負(fù)指數(shù)序列x（n）的實(shí)部和虛部分別是按指數(shù)規(guī)律增長的正弦振蕩序列；2)當(dāng)a3)當(dāng)a=0時(shí)，負(fù)指數(shù)序列x（n）即為虛指數(shù)序列，其實(shí)部和虛部分別是等幅的正弦振蕩序列；長度為N的實(shí)指數(shù)序列x（n）可以通過下面的MATLAB命令實(shí)現(xiàn)。0238。x=[zeros(1,N1)1 zeros(1,N1)]。)。%學(xué)號01 n=(N1):N1。單位沖激序列39。236。u(n)u(n“學(xué)號后兩位”)，并繪出其圖形。stem(n,x)。grid on。x=[zeros(1,NM+1)ones(1,NM)]。)。%學(xué)號01 n=(N1):N1。stem(n,z)。grid on。pn247。程序:(1)N=5。n=0::N1。余弦信號39。A=1。x=A*cos(w*n+phi)。)。A=1。x=A*sin(w*n+phi)。)。n=0::N1。實(shí)部39。虛部39。其中a=學(xué)號后兩位、a=1/“學(xué)號后兩位”。%學(xué)號01 x=a.^n。)。a=1。指數(shù)序列39。在離散信號的頻域分析中，通常將信號表示成單位采樣序列的線性組合，而在頻域中，將信號表示成復(fù)變量或的線性組合。（２）使用到的MATLAB命令有基于DTFT離散時(shí)間信號分析函數(shù)以及求解序列的DFT函數(shù)。x=[,]。angX=angle(H)。plot(w/pi,magX)。)。title(39。subplot(223)。frequency in pi unit39。)。axis([0 ])。xlabel(39。real part39。)。grid。ylabel(39。虛部 part39。n=0:N1。subplot(2,1,1)。title(39。n=0:N1。subplot(2,1,1)。title(39。學(xué)會運(yùn)用MATLAB分析離散時(shí)間系統(tǒng)的系統(tǒng)函數(shù)的零極點(diǎn)；學(xué)會運(yùn)用MATLAB分析系統(tǒng)函數(shù)的零極點(diǎn)分布與其時(shí)域特性的關(guān)系；學(xué)會運(yùn)用MATLAB進(jìn)行離散時(shí)間系統(tǒng)的頻率特性分析。若要獲得系統(tǒng)函數(shù)H（z）的零極點(diǎn)分布圖，可直接應(yīng)用zplane函數(shù)，其語句格式為Zplane（B，A）其中，B與A分別表示H（z）的分子和分母多項(xiàng)式的系數(shù)向量。可根據(jù)系統(tǒng)的傳遞函數(shù)H（z）求單位沖激響應(yīng)h（n）的函數(shù)impz、filter等。另一種形式為[H，w]= freqz（B，A，N，‘whole’）與第一種方式不同之處在于角頻率的范圍由[0，π]擴(kuò)展到[0，2π]。w=[0:20:500]*pi/500。x=x1./x22。zplane(b,a)。subplot(2,2,2)。)。subplot(2,2,4)。)。ylabel(39。數(shù)字濾波器你包括IIR（無限沖激響應(yīng)）和FIR（有限沖激響應(yīng)）型，在設(shè)計(jì)時(shí)通常采用不同的方法。IIR數(shù)字濾波器的設(shè)計(jì)方法如下：（1）沖激響應(yīng)不變法。與沖激響應(yīng)不變法比較，雙線性變化的主要優(yōu)點(diǎn)是靠頻率的非線性關(guān)系得到s平面與z平面的單值一一對應(yīng)關(guān)系，整個(gè)值對應(yīng)于單位圓一周。（2）chebycheby2。%通帶邊界頻率ws=*pi。%188。179。194。1000hz[N,Wn]=cheb1ord(wp,ws,rp,rs,39。s39。freqs(H,W)。grid on。ylabel(39。程序: wp=*pi。fs=500。wp2=2*fs*tan(wp1/2)。)。[bz,az]=bilinear(b,a,fs)。grid on。幅度39。xlabel(39。幅度(dB)39。二、實(shí)驗(yàn)原理如果系統(tǒng)的沖激響應(yīng)h（n）為已知，則系統(tǒng)的輸入輸出關(guān)系為y（n）=x（n）*h（n）對于低通濾波器，只要設(shè)計(jì)出低通濾波器的沖激響應(yīng)函數(shù)，就可以由式得到系統(tǒng)的輸出了。常用的窗函數(shù)有矩形窗、海明窗和布萊克曼窗等。程序1:海明窗設(shè)計(jì)N=48。ws=[w1,w2]。海明窗39。Window=blackman(N+1)。b=fir1(N,ws/pi,Window)。)。ej(wp)=163。: N=9。i=nalpha。H=h*exp(j*n39。stem(n,h)。)。)。subplot(212)。)。axis([0,2,0,])。實(shí)驗(yàn)心得體會：這次實(shí)驗(yàn)使我進(jìn)一步加深了對MATLAB軟件的使用。IIR數(shù)字濾波器的設(shè)計(jì)讓我知道了巴特沃思濾波器和切比雪夫?yàn)V波器的頻率特性，還有雙線性變換及脈沖響應(yīng)不變法設(shè)計(jì)的濾波器的頻率特性。語音是人類相互之間進(jìn)行交流時(shí)使用最多、最自然、最基本也是最重要的信息載體。關(guān)鍵詞：語音處理；DSP；TMS320VC5509。語言是人類特有的功能，聲音是人類常用的工具，是相互傳遞信息的主要手段。再對數(shù)字信號進(jìn)行修改、刪除、強(qiáng)化，并在其他系統(tǒng)芯片中把數(shù)字?jǐn)?shù)據(jù)解譯回模擬數(shù)據(jù)或?qū)嶋H環(huán)境格式。但是它的快速發(fā)展可以說是從1940年前后Dudley的聲碼器和Potter等人的可見語音開始的；20世紀(jì)60年代中期形成的一系列數(shù)字信號處理的理念和技術(shù)基礎(chǔ)；到了80年代，由于矢量量化、隱馬爾可夫模型和人工神經(jīng)網(wǎng)絡(luò)等相繼被應(yīng)用于語音信號處理，并經(jīng)過不斷改進(jìn)與完善，使得語音信號處理技術(shù)產(chǎn)生了突破性的進(jìn)展。因此，語音編碼一直是通信和信號處理的研究熱點(diǎn)，并其取得了驚人的進(jìn)展，目前在PC機(jī)上的語音編碼已經(jīng)趨于成熟，而如何在嵌入式系統(tǒng)中實(shí)時(shí)實(shí)現(xiàn)語音壓縮編碼則是近些年來語音信號處理領(lǐng)域的研究熱點(diǎn)之一。另外，生活中的回聲的成分比較復(fù)雜，有反射、漫反射、折射，還有回聲的多次反射、折射效果。如此產(chǎn)生的回聲，我們稱之為數(shù)字回聲?？刂仆ǖ乐辉谂渲肁IC23 的內(nèi)部寄存器時(shí)工作, 而當(dāng)傳輸音頻數(shù)據(jù)時(shí)則處于閑置狀態(tài)。接口的原理框圖, 如下圖所示。SCLK和SDIN是AIC23控制端口的移位時(shí)鐘和數(shù)據(jù)輸入端,分別與VC5509的I2C模塊端口SCL和SDA相連。 總輸出電路 總輸出電路從左到右各部分電路為：LM386高頻功率放大器及其外圍器件連接電路，語音輸出電路，開關(guān)，揚(yáng)聲器。下圖中LLINEOUT 為左聲道輸出，RLINEOUT為右聲道輸出。這時(shí)從耳機(jī)中能聽到麥克風(fēng)中的輸入語音放送。當(dāng)uDelay和uEffect的數(shù)值增大時(shí)，數(shù)字回聲的效果就會越加的明顯。聲音放送可以加入數(shù)字回聲，數(shù)字回聲的強(qiáng)弱和與原聲的延遲均可在程序中設(shè)定和調(diào)整。因此，本次課程設(shè)計(jì)讓我受益匪淺。另外，導(dǎo)師嚴(yán)謹(jǐn)求實(shí)的治學(xué)態(tài)度、一絲不茍的工作作風(fēng)和高尚的人格魅力，都給了學(xué)生很大感觸，使學(xué)生終生受益。第五篇：數(shù)字信號處理實(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ù)字信號處理A班級學(xué)號_______09002111___________ 姓 名 王都超開課時(shí)間 2011/2012學(xué)年，第 二 學(xué)期實(shí)驗(yàn)一熟悉MATLAB環(huán)境一、實(shí)驗(yàn)?zāi)康模?）熟悉MATLAB的主要操作命令。（5）觀察離散系統(tǒng)的頻率響應(yīng)。D =2 E =F = G =2434096(2)用MATLAB實(shí)現(xiàn)下列序列： a)x(n)= 0163。a=(+3*i)*n。x2=exp(a)。15(4)繪出下列時(shí)間函數(shù)的圖形，對x軸、y軸以及圖形上方均須加上適當(dāng)?shù)臉?biāo)注： a)x(t)=sin(2pt)0163。4s t=0::4。)。x(t)39。2z1+z2)/(+)，求出并繪制H(z)的幅(7)計(jì)算序列{821 2 3}和序列{2 313}的離散卷積，并作圖表示卷積結(jié)果。a=[12]。h=filter(a,b,x1)。(4)應(yīng)用FFT實(shí)現(xiàn)序列的線性卷積和相關(guān)。e239。15 其他b)衰減正弦序列236。n163。8n239。34163。xd(n)=237。n163。(3)觀察三角波和反三角波序列的時(shí)域和幅頻特性，用N=8點(diǎn)FFT分析信號序列xc(n)和觀察兩者的序列形狀和頻譜曲線有什么異同？繪出兩序列及其幅頻特性xd(n)的幅頻特性，曲線。q=2。.39。)。.39。)。xb=exp(a*n).*sin(2*pi*f*n)。)。Xb=fft(xb,16)。)。實(shí)驗(yàn)過程與結(jié)果（含實(shí)驗(yàn)程序、運(yùn)行的數(shù)據(jù)結(jié)果和圖形）；實(shí)驗(yàn)三 IIR數(shù)字濾波器的設(shè)計(jì)一、實(shí)驗(yàn)?zāi)康?1)掌握雙線性變換法及脈沖響應(yīng)不變法設(shè)計(jì)IIR數(shù)字濾波器的具體設(shè)計(jì)方法及其原理，熟悉用雙線性變換法及脈沖響應(yīng)不變法設(shè)計(jì)低通、高通和帶通IIR數(shù)字濾波器的計(jì)算機(jī)編程。脈沖響應(yīng)不變法： fc=1000。s39。f =