【正文】 脈沖響應不變法： fc=1000。)。)。n163。n163。h=filter(a,b,x1)。)。a=(+3*i)*n。另外，導師嚴謹求實的治學態(tài)度、一絲不茍的工作作風和高尚的人格魅力，都給了學生很大感觸，使學生終生受益。這時從耳機中能聽到麥克風中的輸入語音放送。接口的原理框圖, 如下圖所示。因此，語音編碼一直是通信和信號處理的研究熱點，并其取得了驚人的進展，目前在PC機上的語音編碼已經趨于成熟，而如何在嵌入式系統(tǒng)中實時實現(xiàn)語音壓縮編碼則是近些年來語音信號處理領域的研究熱點之一。關鍵詞：語音處理；DSP；TMS320VC5509。axis([0,2,0,])。)。: N=9。Window=blackman(N+1)。常用的窗函數(shù)有矩形窗、海明窗和布萊克曼窗等。幅度39。wp2=2*fs*tan(wp1/2)。grid on。194。（2）chebycheby2。ylabel(39。subplot(2,2,2)。另一種形式為[H，w]= freqz（B，A，N，‘whole’）與第一種方式不同之處在于角頻率的范圍由[0，π]擴展到[0，2π]。title(39。subplot(2,1,1)。grid。axis([0 ])。title(39。x=[,]。a=1。虛部39。x=A*sin(w*n+phi)。A=1。pn247。)。u(n)u(n“學號后兩位”)，并繪出其圖形。)。長度為N的正弦序列x（n）可以通過下面的MATLAB命令獲得。也認識到自己的不足，以后要進一步學習。d:39。plot(f, abs(fftwave(n1)))。%最小衰減[N,Wn]=buttord(Wp,Ws,Rp,Rs)。plot(f,abs(y0(n1)))。subplot(3,2,3)。噪聲信號的頻譜圖39。subplot(3,2,1)。按Enter鍵回放錄音39。%加網(wǎng)格分析：就時域圖看，男生的時域圖中振幅比女生的高，對于頻譜圖女生的高頻成分比較多Fs=8000。幅值/A39。)。)。%按任意鍵播放語音信號wavwrite(x,Fs,39。,time)。幅值/A39。)。ws2=2*pi*3000/fs。xlabel(39。頻率w/pi39。%加網(wǎng)格wavplay(ys,8000)。ylabel(39。n=ceil(*pi /width)。%加標題 xlabel(39。%頻帶寬度 figure(4)。)。T=1。)。wap=2*tan(wp/2)。頻率（HZ）39。wap=2*tan(wp/2)。%標題 xlabel(39。%顯示橫坐標表示的意義 ylabel(39。原始信號的時域波形39。%文字提示 fprintf(39。2fh四．具體實現(xiàn)在MATLAB中用wavrecord、wavread、wavplay、wavwrite對聲音進行錄制、讀取、回放、存儲。dT(t)則有?(jW)=1X(jW)*D(jW)XaaT2p1165。二．設計內容錄制一段個人自己的語音信號，并對錄制的信號進行采樣；畫出采樣后語音信號的時域波形和頻譜圖；給定濾波器的性能指標，采用窗函數(shù)法和雙線性變換法設計濾波器，并畫出濾波器的頻率響應；然后用自己設計的濾波器對采集的信號進行濾波，畫出濾波后信號的時域波形和頻譜，并對濾波前后的信號進行對比，分析信號的變化；回放語音信號；換一個與你性別相異的人錄制同樣一段語音內容，分析兩段內容相同的語音信號頻譜之間有什么特點；再錄制一段同樣長時間的背景噪聲疊加到你的語音信號中，分析疊加前后信號頻譜的變化，設計一個合適的濾波器，能夠把該噪聲濾除；三．設計原理，利用函數(shù)wavrecord()，wavwrite(),wavread(),wavplay()對語音信號進行錄制，存儲，讀取，回放。數(shù)字信號處理課程屬于專業(yè)基礎課，所涵蓋的內容主要有：離散時間信號與系統(tǒng)的基本概念及描述方法，離散傅立葉變換及快速傅立葉變換，數(shù)字濾波器結構及設計等。數(shù)字信號處理技術正飛速發(fā)展，它不但自成一門學科，更是以不同形式影響和滲透到其他學科：它與國民經濟息息相關，與國防建設緊密相連；它影響或改變著我們的生產、生活方式，因此受到人們普遍的關注。可以從事硬件方面，比如說PCB，別小看這門技術，平時我們在試驗時制作的簡單，這一技術難點就在于板的層數(shù)越多，要做的越穩(wěn)定就越難，這可是非常有難度的，如果學好了學精了，也是非常好找工作的。實驗課主要是以MATLAB為平臺，充分利用MATLAB的數(shù)字信號處理各種功能讓學生親自動手將課堂所學進行仿真實現(xiàn)。?a(t)Xa(jW)171。Ws/2即連續(xù)信號是帶限的，且信號最高頻率不超過抽樣頻率的二分之一，則可不失真恢復。%錄制語音信號 fprintf(39。subplot(2,1,1)。%加標題 xlabel(39。%畫出截取后的時域圖形 title(39。wp=2*pi*fp/Fs。巴特沃斯低通濾波器39。ws=2*pi*fs/Fs。plot(w,abs(h))。rs=100。plot(w,abs(h))。%返回采樣點數(shù) df=fs/N。%原始信號的頻譜圖 title(39。ws=2*pi*1200/fs。%加標題 xlabel(39。低通濾波器濾波后信號的頻譜圖39。%將信號送入高通濾波器濾波 subplot(4,2,5)。%將濾波后的語音信號進行快速傅立葉變換 subplot(4,2,7)。wc1=2*pi*1000/fs。%在四行兩列的第六個窗口顯示圖形 plot(ys1)。頻率w/pi39。%抽樣頻率 time=3。)。plot(x)。)。原始男生信號的時域波形39。幅值 39。錄音結束39。%采樣間隔n1=1:N/2。幅值39。%采樣間隔 n1=1:N/2。grid。%通帶截至頻率 Ws=2*2000/fs。)。圖形如下：噪聲信號的時域波形1100噪聲信號的頻譜圖幅值/ 10疊加信號的時域波形10101時間/t2200幅值/db34幅值/，常會出現(xiàn)實際發(fā)出聲音落后錄制動作半拍，可在[x,fs,bits]=wavread(39。整個過程中學到了很多新的知識，特別是對Matlab的使用終于有些了解。（4）正弦序列x（n）這里，正弦序列的參數(shù)都是實數(shù)。stem(n,x)。1n179。stem(n,x)。2230。)。w=2*pi。)。grid on。三、實驗內容及分析（１）編程計算并畫出下面DTFT的實部、虛部、幅度和相位譜。ylabel(39。title(39。axis([0 11 1])。xn=cos(5*pi*n/16)。stem(n,xn)。利用系統(tǒng)的頻率響應，可以分析系統(tǒng)對各種頻率成分的響應特性，并推出系統(tǒng)的特性（高通、低通、帶通、帶阻等）。title(39。title(39。所以從模擬傳遞函數(shù)可直接通過代數(shù)置換得到數(shù)字濾波器的傳遞函數(shù)。233。[P,Q]=freqs(H,W)。ts=1/fs。xlabel(39。但是將h（n）作為濾波器的脈沖響應有兩個問題：一是它是無限長的；二是它是非因果的。)。w163。axis([0,8,])。ylabel(39。在高度信息化的今天，語音信號處理是信息高速公路、多媒體技術、辦公自動化、現(xiàn)代通信及智能系統(tǒng)等新興領域應用的核心技術之一。一方面，對聲學語音學統(tǒng)計模型的研究逐漸深入，魯棒的語音識別、基于語音段的建模方法及隱馬爾可夫模型與人工神經網(wǎng)絡的結合成為研究的熱點。AIC23通過麥克風輸入或者立體聲音頻輸入采集模擬信號, 并把模擬信號轉化為數(shù)字信號, 存儲到DSP的內部RAM中,以便DSP處理。用戶可以根據(jù)電阻阻值調節(jié)增益的大小，使語音輸入輸出達到最佳效果。參考文獻[1][M].北京:中國水利水電出版社,2004.[2]王安民,陳明欣,[M].北京:清華大學出版社,2007 [3]彭啟琮,[M].成都:電子科技大學出版社,1997 [4]李宏偉,[J].解放軍理工大學學報,2001(1):41~44 [5][M].梁曉雯,裴小平,李玉虎,:清華大學出版社,2006 [6][M].北京:機械工業(yè)出版社,2003比較圖4和圖5,可以看到1200Hz以上的頻譜明顯得到了抑制。n163。x=cos(100*pi*t).*sin(pi*t)。b=[1 ]。eansin(2pfn)xb(n)=237。n4239。)。subplot(2,3,2)。(2)觀察雙線性變換及脈沖響應不變法設計的濾波器的頻域特性，了解雙線性變換法及脈沖響應不變法的特點。[B,A]=butter(3, OmegaC,39。)。stem(abs(Xa),39。7 其他(1)觀察高斯序列的時域和幅頻特性，固定信號xa(n)中參數(shù)p=8，改變q的值，使q分別等于2，4，8，觀察它們的時域和幅頻特性，了解當q取不同值時，對信號序列的時域幅頻特性的影響；固定q=8，改變p，使p分別等于8，13，14，觀察參數(shù)p變化對信號序列的時域及幅頻特性的影響，觀察p等于多少時，會發(fā)生明顯的泄漏現(xiàn)象，混疊是否也隨之出現(xiàn)？記錄實驗中觀察到的現(xiàn)象，繪出相應的時域序列和幅頻特性曲線。xc(n)=237。(3)了解應用FFT進行信號頻譜分析過程中可能出現(xiàn)的問題，以便在實際中正確應用FFT。),ylabel(39。15n=0:1:15。我們始終是一個團結、友愛、積極向上的集體。分別調整uDelay和uEffect的取值，使他們保持在01023范圍內，同時聽聽耳機中的輸出有何變化。CS接數(shù)字地, 定義了I2 C總線上AIC23的外設地址, 通過將CS接到高電平或低電平, 可以選擇A IC23作為從設備在I2 C總線上的地址。同時，反射聲音的物體對聲波的反射能力，決定了聽到的回聲的強弱和質量。另外通過語音來傳遞信息是人類最重要的、最有效、最常用的交換信息的形式。)。n39。n=(0:8)。ws=[w1,w2]。三、實驗內容及分析分別用海明窗和布萊克曼窗設計一個48階的FIR帶通濾波器，通帶為Wn=[]。grid on。s39。)。229。程序:wp=*pi。（二）波形圖圖41 幅頻、相頻響應曲線、零極點分布圖實驗六 IIR數(shù)字濾波器的設計一、實驗目的從理論上講，任何的線性是不變（LTI）離散時間系統(tǒng)都可以看做一個數(shù)字濾波器，因此設計數(shù)字濾波器實際就是設計離散時間系統(tǒng)。幅度部分39。a=[1,]。（１）（２）實驗四 離散時間LTI系統(tǒng)的Z域分析一、實驗目的本實驗通過使用MATLAB函數(shù)對離散時間系統(tǒng)的一些特性進行仿真分析，以加深對離散時間系統(tǒng)的零極點、穩(wěn)定性，頻率響應等概念的理解。stem(k,abs(Xk))。)。grid。axis([0 0 ])。magX=abs(H)。title(39。編制程序產生指數(shù)序列x(n)=an，并繪出其圖形。正弦信號39。n=0::N1。x(n)=sin(2n)并繪出其圖形。M=1。x=[zeros(1,N1)ones(1,N)]。M=1。當|a|1時，單邊指數(shù)序列發(fā)散；當|a|0時，該序列均取正值；當a（6）負指數(shù)序列x（n）當a=0時，得到虛指數(shù)序列x（n）。使用MATLAB軟件可以很方便地產生各種常見的離散時間信號，而且它還具有強大繪圖功能，便于用戶直觀地處理輸出結果。當通過高通濾波器后，把低頻濾除，聲音變得比較就尖銳。)。%DTFT變換ys=filter(num,den,yy)。)。,2)。[x,fs,bits]=wavread(39。噪聲信號的時域波形39。%暫停命令 wavplay(x,Fs)。按Enter鍵錄音%ds39。%加網(wǎng)格y0=fft(y)。 C: 39。)。%存儲語音信號 分析男女聲音的頻譜x=wavread(39。錄音中......39。帶通濾波器濾波后信號的頻譜圖39。)。wn2=(ws2+wc2)/2。ylabel(39。幅值/A39。ws=2*5000/fs。%顯示標表示的意義 grid。,kaiser(n+1,beta))。%顯示橫坐標表示的意義 ylabel(39。%錄制信號的時域波形title(39。)。was=2*tan(ws/2)。)。s39。耗損（dB）39。s39。ylabel(39。%顯示縱坐標表示的意義 title(39。幅值/A39。pause。%錄音時間 fprintf(39。Xa(jWjk)=229。，F(xiàn)IR濾波器利用函數(shù)fftfilt對信號進行濾波。隨著計算機技術、微電子技術、數(shù)字信號處理理論和方法的發(fā)展，半個世紀以來，尤其是最近的三十來年里，數(shù)字信號處理的方法和應用得到了飛躍式的發(fā)展，數(shù)字信號處理的地位和作用變得越來越重要。第二單元的課程我們理解了時域離散信號（序列）的傅立葉變換，時域離散信號Z變換，時域離散系統(tǒng)的頻域分析。對于我們通信專業(yè)，我覺得是個很好的專業(yè)，現(xiàn)在這個專業(yè)很熱門，這個專業(yè)以后就業(yè)的方向也很多，就業(yè)面很廣。近年來，國外及國內有些學校對一般電類專業(yè)該課程的教學主要強調應用性學習，主要介紹數(shù)字信號處理的用途和用法，而對其深奧的理論推導僅做一般介紹，并給學生提供進行實驗的機會，以激發(fā)學生對該課程的興趣和學習主動性。dT(t)=頻域描述：利用傅氏變換的性質，時域相乘頻域卷積，若n=165。a239。x=wavrecord(time*Fs,Fs,39。%采樣間隔 n1=1:N/2。%原始信號的頻譜圖 title(39。C:39。Fs=22050。[h,w]=freqz(Bz,Az,512,Fs1*22050)。rp=1。[Bz,Az]=bilinear(B,A,Fs1)。fs=[1000,3200]。[Bz,Az]=bilinear(B,A,Fs1)。)。%加網(wǎng)格y0=fft(x)。%加網(wǎng)格3．1低通濾波器濾波 fs=8000。%信號的時域波形title(39。幅值/A39。,kaiser(n+2,beta))。)。%加網(wǎng)格wavplay(ys,8000)。ys1= fftfilt(h,x)。%在四行兩列的第八個窗口顯示圖形plot(f, abs(fftwave(n1)))。%sound(x, Fs, bits)。)。f=[(n11)*(2*pi/N)]/pi。)。%頻帶寬度 subplot(2,2,3)。頻率w/pi39。double39。%對語音信號進行采樣wavplay(x1,fs)。%快速傅立葉變換subplot(3,2,2)。%將兩個聲音疊加：wavplay(yy,fs)。ylabel(39。grid。)。%顯示標表示的意義 grid。在設計的過程中遇到了很