【正文】 a=(N1)/2。*w)。line([0,10],[0,0])。plot(w/pi,magH)。title(39。做這兩個(gè)實(shí)驗(yàn)的時(shí)候程序有點(diǎn)困難，但經(jīng)過(guò)細(xì)心的改寫(xiě)圖形最終出來(lái)了。TLV320AIC23B1 緒論語(yǔ)音是人類(lèi)相互間所進(jìn)行的通信的最自然和最簡(jiǎn)潔方便的形式，語(yǔ)音通信是一種理想的人機(jī)通信方式。它不僅具有可編程性，而且其實(shí)時(shí)運(yùn)行速度可達(dá)每秒數(shù)以千萬(wàn)條復(fù)雜指令程序，遠(yuǎn)遠(yuǎn)超過(guò)通用微處理器，是數(shù)字化電子世界中日益重要的電腦芯片。2系統(tǒng)設(shè)計(jì)在實(shí)際生活中，當(dāng)聲源遇到物體時(shí)會(huì)發(fā)生反射，反射的聲波和聲源聲波一起傳輸，聽(tīng)者會(huì)發(fā)現(xiàn)反射聲波部分比聲源聲波慢一些，類(lèi)似人們面對(duì)山體高聲呼喊后可以在過(guò)一會(huì)兒聽(tīng)到回聲的現(xiàn)象。本次實(shí)驗(yàn)的程序流程圖如下： 程序流程圖本次實(shí)驗(yàn)的系統(tǒng)框圖如下： 系統(tǒng)框圖3硬件設(shè)計(jì) 硬件結(jié)構(gòu), 系統(tǒng)主要包括VC5509和A IC23 兩個(gè)模塊。 AIC23與VC5509接口原理圖系統(tǒng)中A IC23的主時(shí)鐘12 MHz直接由外部的晶振提供。 語(yǔ)音輸入電路 語(yǔ)音輸出電路 語(yǔ)音輸出電路語(yǔ)音信號(hào)通道包括模擬輸入和模擬輸出兩個(gè)部分。將觀察窗口中bEcho的取值改成非0值。5 總結(jié)與心得體會(huì)通過(guò)本次課程設(shè)計(jì)，我明白了細(xì)節(jié)決定成敗這句話(huà)的道理，在實(shí)驗(yàn)中，有很多注意的地方，都被忽視了，導(dǎo)致再花費(fèi)更多的時(shí)間去修改，這嚴(yán)重影響了試驗(yàn)的進(jìn)度。在此，學(xué)生謹(jǐn)向老師致以最真摯的感激和最崇高的敬佩之情。二、實(shí)驗(yàn)內(nèi)容(1)數(shù)組的加、減、乘、除和乘方運(yùn)算。stem(x1)b)x(n)=e(+3j)n 0163。t163。xlabel(39。(8)求以下差分方程所描述系統(tǒng)的單位脈沖響應(yīng)h(n), 0163。stem(n,h)axis([1 ])實(shí)驗(yàn)二快速傅里葉變換及其應(yīng)用一、實(shí)驗(yàn)?zāi)康?1)在理論學(xué)習(xí)的基礎(chǔ)上，通過(guò)本實(shí)驗(yàn)，加深對(duì)FFT的理解，熟悉MATLAB中的有關(guān)函數(shù)。0238。15其他c)三角波序列 236。n163。34163。xa=exp((np).^2/q)。Xa=fft(xa,16)。A=1。title(39。title(39。fs=4000。[h1,w]=freqz(num1,den1)。二、實(shí)驗(yàn)內(nèi)容（1）P162 設(shè)采樣周期T=250ms(采樣頻率fs=4kHz)，分別用脈沖響應(yīng)不變法和雙線(xiàn)性變換法設(shè)計(jì)一個(gè)三階巴特沃思低通濾波器，其3dB邊界頻率為fc=1kHz。.39。.39。Xa(k)=FFT[xa(n)]的波形 39。xa波形39。p=8。0163。n163。0163。xa(n)=237。n=0:1:N1。)(6)給定一因果系統(tǒng)H(z)=(1+頻響應(yīng)和相頻響應(yīng)。r39。n163。x1=.^n。（4）用MATLAB編程并學(xué)會(huì)創(chuàng)建函數(shù)。在整個(gè)課程設(shè)計(jì)的選題、研究和撰寫(xiě)過(guò)程中，老師都給了我精心的指導(dǎo)、熱忱的鼓勵(lì)和支持，他的精心點(diǎn)撥為我開(kāi)拓了研究視野，修正了寫(xiě)作思路，對(duì)課程設(shè)計(jì)的完善和質(zhì)量的提高起到了關(guān)鍵性的作用。當(dāng)然通過(guò)復(fù)雜運(yùn)算，可以計(jì)算各種效應(yīng)的混響效果。4 實(shí)驗(yàn)結(jié)果及分析 實(shí)驗(yàn)結(jié)果按“F5”鍵運(yùn)行，注意觀察窗口中的bEcho=0，表示數(shù)字回聲功能沒(méi)有激活。 硬件電路設(shè)計(jì) 總輸入電路 總輸入電路從左到右各部分電路為：話(huà)筒，開(kāi)關(guān)，語(yǔ)音輸入電路，UA741高增益放大電路，有源二階帶 通濾波器。AIC23能夠?qū)崿F(xiàn)與VC5509 DSP的McBSP端口的無(wú)縫連接, 使系統(tǒng)設(shè)計(jì)更加簡(jiǎn)單。當(dāng)然通過(guò)復(fù)雜運(yùn)算，可以計(jì)算各種效應(yīng)的混響效果。在通信越來(lái)越發(fā)達(dá)的當(dāng)今世界，尤其最近幾十年，語(yǔ)音壓縮編碼技術(shù)在移動(dòng) 通信、IP電話(huà)通信、保密通信、衛(wèi)星通信以及語(yǔ)音存儲(chǔ)等很多方面得到了廣泛的應(yīng)用。其工作原理是接收模擬信號(hào)，轉(zhuǎn)換為0或1的數(shù)字信號(hào)。本文設(shè)計(jì)了一個(gè)基于TMS320VC5509定點(diǎn)的語(yǔ)音信號(hào)處理系統(tǒng)，實(shí)現(xiàn)對(duì)語(yǔ)音信號(hào)的采集、處理與回放等功能，為今后復(fù)雜的語(yǔ)音信號(hào)處理算法的研究和實(shí)時(shí)實(shí)現(xiàn)提供一個(gè)通用平臺(tái)。在實(shí)驗(yàn)的過(guò)程中，我對(duì)數(shù)字濾波器的整個(gè)過(guò)程有了很好的理解和掌握。)。)。矩形窗設(shè)計(jì)h(n)39。%矩形窗函數(shù)設(shè)計(jì)的系統(tǒng)脈沖響應(yīng) w=(0:1:500)*2*pi/500。w163。布萊克曼窗39。程序2:萊克曼窗設(shè)計(jì)N=48。w2=。設(shè)計(jì)時(shí)，要根據(jù)阻帶的最小衰減和過(guò)渡帶寬度來(lái)選擇恰當(dāng)?shù)拇昂瘮?shù)類(lèi)型和窗口長(zhǎng)度N。ylabel(39。)ylabel(39。[b,a]=lp2lp(Bap,Aap,Wn)。ws1=ws*ts?；贐utterworth型模擬濾波器原型使用雙線(xiàn)性變換方法設(shè)計(jì)數(shù)字濾波器的，要求參數(shù)為截止頻率Wp=；通帶最大衰減Ap=1dB；阻帶截止頻率Ws=；阻帶最小衰減AS=40dB。plot(Q*Fs/(2*pi),abs(P))。[Z,P,K]=cheby1(N,rp,Wn,39。249。%阻帶最小衰減Fs=1000。（1）butter。IIR濾波器的階數(shù)就等于所選的模擬原型濾波器的階數(shù)，所以其階數(shù)確定主要是在模擬原型濾波器中進(jìn)行的。)。)。)。x22=x2+(x2==0)*eps。一種形式為[H,w]= reqz(B,A,N)其中，B與A分別表示H（z）分子和分母多項(xiàng)式的系數(shù)向量；N為正整數(shù)，默認(rèn)值為512；返回值w包含[0，π]范圍內(nèi)的N個(gè)頻率等分點(diǎn)；返回值H則是離散時(shí)間系統(tǒng)頻率響應(yīng)在0~π范圍內(nèi)N個(gè)頻率處的值。它的作用是將H（z）的有理分式表示式轉(zhuǎn)換為零極點(diǎn)增益形式。stem(k,abs(Xk))。figure N=64。Xk=fft(xn,N)。title(39。plot(w/pi,imagX)。ylabel(39。)。xlabel(39。)。subplot(221)。３－１clear。實(shí)驗(yàn)三 離散時(shí)間信號(hào)的頻域分析一、實(shí)驗(yàn)?zāi)康男盘?hào)的頻域分析是信號(hào)處理中一種有效的工具。n=0:N1。a=1。subplot(2,1,2)stem(n,imag(x)),title(39。N=3。n=0::N1。余弦信號(hào)39。分析:該序列具有周期性,且輸出為余弦信號(hào).(2)N=5。phi=0。編制程序產(chǎn)生正弦序列x(n)=cos(2pn)、x(n)=cos231。z=xy。單位階躍序列39。)。、u(n“學(xué)號(hào)后兩位”)及0np0238。title(39。單位沖激序列39?！?并繪出其圖及d(n”學(xué)號(hào)后兩位0n185。可以證明，只有當(dāng)2π/w為有理數(shù)時(shí)，正弦序列具有周期性。常見(jiàn)的離散信號(hào)如下：（1）單位沖激序列δ（n）如果δ（n）在時(shí)間軸上延遲了k個(gè)單位，得到δ（nk），即長(zhǎng)度為N的單位沖激序列δ（n）可以通過(guò)下面的MATLAB命令獲得。在這次的課程設(shè)計(jì)中讓我體會(huì)最深的是：知識(shí)來(lái)不得半點(diǎn)的馬虎。，噪聲與原信號(hào)明顯區(qū)分,但通過(guò)低通濾波器后，噪聲沒(méi)有濾除，信號(hào)產(chǎn)生失真。)加 窗[x,fs,bits]=wavread(39。ylabel(39。%網(wǎng)格 subplot(3,2,6)。plot(ys)。%最大衰減 Rs=100。ylabel(39。%快速傅立葉變換 subplot(3,2,4)。xlabel(39。%頻帶寬度 figure(5)。 C: 39。title(39。幅值/A39。%頻帶寬度 figure(5)。%存儲(chǔ)語(yǔ)音信號(hào) 對(duì)噪聲進(jìn)行頻譜的分析[x1,fs,bits]=wavread(39。%文字提示 fprintf(39。錄音中......39。%顯示縱坐標(biāo)表示的意義 grid。原始男生信號(hào)的頻譜圖39。%加標(biāo)題 ylabel(39。%返回采樣點(diǎn)數(shù) df=fs/N。幅值 39。plot(f,abs(y0(n1)))。原始女生信號(hào)的時(shí)域波形39。%播放原始信號(hào) N=length(x)。%暫停命令 wavplay(x,Fs)。)。按Enter鍵錄音%ds39。%播放濾波后信號(hào) 圖形如下：原始信號(hào)的時(shí)域波形幅值/A01012x 10原始信號(hào)的頻譜圖34幅值/分析：三個(gè)濾波器濾波后的聲音與原來(lái)的聲音都發(fā)生了變化。ylabel(39。%顯示縱坐標(biāo)表示的意義 grid。帶通濾波后信號(hào)的時(shí)域波形39。band39。ws1=2*pi*1200/fs。高通濾波器濾波后信號(hào)的頻譜圖39。%繪制模值 axis([0 1 0 50])。)。%信號(hào)的時(shí)域波形 xlabel(39。n=ceil(*pi/width)。%標(biāo)題 grid。頻率w/pi39。)。%將濾波后的語(yǔ)音信號(hào)進(jìn)行快速傅立葉變換 figure(4)。wn=(ws+wc)/2。原始信號(hào)的頻譜圖39。)。幅值/A39。f=[(n11)*(2*pi/N)]/pi。grid on。巴特沃斯帶通濾波器39。[B,A]=butter(N,wc,39。ws=2*pi*fs/Fs。grid on。巴特沃斯高通濾波器39。,39。Fs1=1。%高通濾波fp=4800。xlabel(39。)。Fs1=1。%畫(huà)網(wǎng)格 去除無(wú)效點(diǎn)MATLAB程序如下。)。%加網(wǎng)格 原始信號(hào)的時(shí)域與頻譜圖針對(duì)實(shí)際發(fā)出聲音落后錄制動(dòng)作半拍的現(xiàn)象，如何拔除對(duì)無(wú)效點(diǎn)的采樣的問(wèn)題: 出現(xiàn)這種現(xiàn)象的原因主要是錄音開(kāi)始時(shí)，人的反應(yīng)慢了半拍，導(dǎo)致出現(xiàn)了一些無(wú)效點(diǎn)，而后而出現(xiàn)的無(wú)效的點(diǎn)，主要是已經(jīng)沒(méi)有聲音的動(dòng)作，先讀取聲音出來(lái)，將原始語(yǔ)音信號(hào)時(shí)域波形圖畫(huà)出來(lái)，根據(jù)己得到的信號(hào)，可以在第二次讀取聲音的后面設(shè)定采樣點(diǎn)，取好有效點(diǎn)，畫(huà)出濾除無(wú)效點(diǎn)后的語(yǔ)音信號(hào)時(shí)域波形圖，對(duì)比可以看出。)。%快速傅立葉變換 figure(2)。%錄制信號(hào)的時(shí)域波形 title(39。)。)。%暫停命令 fprintf(39。2Wh 或 fs179。如果：236。xa(t)DT(jW)171。n=165。第二篇：數(shù)字信號(hào)處理課程設(shè)計(jì)..課程設(shè)計(jì)報(bào)告課程名稱(chēng)： 數(shù)字信號(hào)處理 課題名稱(chēng)： 語(yǔ)音信號(hào)的處理與濾波姓 名： 學(xué) 號(hào)： 院 系： 專(zhuān)業(yè)班級(jí)： 指導(dǎo)教師： 完成日期： 2013年7月2日目錄第1部分 課程設(shè)計(jì)報(bào)告………………………………………3 一．設(shè)計(jì)目的……………………………………………3 二．設(shè)計(jì)內(nèi)容……………………………………………3 三．設(shè)計(jì)原理……………………………………………3 四．具體實(shí)現(xiàn)……………………………………………5 …………………………………5 ………………………8 …………………13 ………………………………19 ……………………19 ……………………………22 五． 結(jié)果分析……………………………………………27 第2部分 課程設(shè)計(jì)總結(jié)………………………………28 一． 參考文獻(xiàn)……………………………………………28第1部分 課程設(shè)計(jì)報(bào)告一．設(shè)計(jì)目的綜合運(yùn)用本課程的理論知識(shí)進(jìn)行頻譜分析以及濾波器設(shè)計(jì)，通過(guò)理論推導(dǎo)得出相應(yīng)結(jié)論，并利用MATLAB作為工具進(jìn)行實(shí)現(xiàn)，從而復(fù)習(xí)鞏固課堂所學(xué)的理論知識(shí)，提高對(duì)所學(xué)知識(shí)的綜合應(yīng)用能力，并從實(shí)踐上初步實(shí)現(xiàn)對(duì)數(shù)字信號(hào)的處理。因此，我們班的同學(xué)在這一個(gè)學(xué)期的學(xué)習(xí)中，我們都感覺(jué)比較輕松。我選擇了這個(gè)專(zhuān)業(yè)，在這里讀了三年關(guān)于通信知識(shí)的書(shū)，我還是想以后畢業(yè)能夠從事這個(gè)方面的工作，現(xiàn)在學(xué)了通信原理、數(shù)字信號(hào)處理這些很有用的專(zhuān)業(yè)課，所以，我在以后的學(xué)習(xí)中，我會(huì)把這些方面的知識(shí)學(xué)扎實(shí)，從事技術(shù)這一塊要能吃苦，我也做好了準(zhǔn)備，現(xiàn)在還很年輕，年輕的時(shí)候多吃點(diǎn)苦沒(méi)什么，為了我自己美好的將來(lái)，我會(huì)努力學(xué)好這個(gè)專(zhuān)業(yè)的。第六單元的課程我們理解數(shù)字濾波器的基本概念，模擬濾波器的設(shè)計(jì),巴特沃斯濾波器的設(shè)計(jì),切比雪夫?yàn)V波器的設(shè)計(jì)，脈沖響應(yīng)不變法設(shè)計(jì)無(wú)限脈沖響應(yīng)字?jǐn)?shù)字濾波器，雙線(xiàn)性變換法設(shè)計(jì)無(wú)限脈沖響應(yīng)字?jǐn)?shù)字濾波器，數(shù)字高通﹑帶通﹑帶阻濾波器的設(shè)計(jì)。第一篇：數(shù)字信號(hào)處理學(xué)習(xí)心得體會(huì)數(shù)字信號(hào)處理學(xué)習(xí)心得一、課程認(rèn)識(shí)和內(nèi)容理解《數(shù)字信號(hào)處理》是我們通信工程和電子類(lèi)專(zhuān)業(yè)的一門(mén)重要的專(zhuān)業(yè)基礎(chǔ)課程，主要任務(wù)是研究數(shù)字信號(hào)處理理論的基本概念和基本分析方法，通過(guò)建立數(shù)學(xué)模型和適當(dāng)?shù)臄?shù)學(xué)分析處理，來(lái)展示這些理論和方法的實(shí)際應(yīng)用。第五單元的課程我們學(xué)了網(wǎng)絡(luò)結(jié)構(gòu)的表示方法——信號(hào)流圖，無(wú)限脈沖響應(yīng)基本網(wǎng)絡(luò)結(jié)構(gòu)，有限脈沖響應(yīng)基本網(wǎng)絡(luò)結(jié)構(gòu)，時(shí)域離散系統(tǒng)狀態(tài)變量分析法。也可以從事軟件方面，這實(shí)際上要我們具備比較好的模電和數(shù)電的基礎(chǔ)知識(shí)。利用實(shí)驗(yàn)課讓我們來(lái)編程做仿真，體會(huì)信號(hào)處理課程的樂(lè)趣，這樣子激發(fā)了學(xué)生的興趣、提高了教學(xué)的效果。實(shí)驗(yàn)課還可以通過(guò)用DSP試驗(yàn)箱實(shí)現(xiàn)數(shù)字信號(hào)處理的功能向?qū)W生進(jìn)行演示。xa(nT)d(tnT)xn=165。xXa(jW)171。?(jW)與X(jW)的關(guān)系：理想抽樣信號(hào)的頻譜是連續(xù)信號(hào)頻譜的Xaa周期延拓,重復(fù)周期為Ws(采樣角頻率)。奈奎斯特采樣定理：要使實(shí)信號(hào)采樣后能夠不失真還原，采樣頻率必須大于信號(hào)最高頻率的兩倍：Ws179。%文字提示 pause。錄音結(jié)束39。C:39。plot(x)。%加網(wǎng)格y0=fft(x)。頻率w/pi39。%加標(biāo)題grid。截取后的聲音時(shí)域圖形39。grid。ws=2*pi*fs/Fs。s39。)。9 巴特沃思低通濾波器MATLAB程序如下。T=1。high39。title(39。)。wp=2*pi*fp/Fs。)。title(39。)。%采樣間隔 n1=1:N/2。%加標(biāo)題 ylabel(39。原始信號(hào)的頻譜圖39。%顯示縱坐標(biāo)表示的意義 title(39。width=wswc。%信號(hào)送入濾波器濾波，ys為輸出 fftwave=fft(ys)。頻率w/pi39。%繪制模值 xlabel(39。)。wn=(ws+wc)/2。%在四行兩列的第五個(gè)窗口顯示圖形 plot(ys)。高通濾波后信號(hào)的時(shí)域波形39。%在四行兩列的第七個(gè)窗口顯示圖形 plot(f,abs(fftwave(n1)))。%顯示標(biāo)表示的意義title(39。wc2=2*pi*3200/fs。h=fir1(n,wn/pi,39。%繪制后信號(hào)的時(shí)域的圖形 title(39。)。)。%網(wǎng)格wavplay(ys1,8000)。%錄音時(shí)間 fprintf(39。double39。pause。%播放原始信號(hào),解決落后半拍wavplay(x,fs)。%錄制信號(hào)的時(shí)域波形title(39。%快速傅立葉變換 subplot(2,2,2)。%顯示橫坐標(biāo)表示的意義 ylabel(39。%播放原始信號(hào) N=length(y)。)。%原始信號(hào)的頻譜圖 title(39。)。%暫停命令 fprintf(39。)。)。f=[(n11)*(2*pi/N)]/pi。ylabel(39。)。%對(duì)錄入信號(hào)進(jìn)行采樣[x1,fs,bits]=wavread(39。f=[(n