【正文】 (5)。%采樣間隔 n1=1:N/2。%播放疊加后信號 N=length(yy)。%對噪聲信號進行采樣yy=x+x1。 C: 39。)。[x,fs,bits]=wavread(39。)。title(39。幅值39。plot(f,abs(y0(n1)))。y0=fft(x1)。幅值/A39。grid。噪聲信號的時域波形39。plot(x1)。%頻帶寬度 figure(5)。%采樣間隔n1=1:N/2。%播放噪聲信號 N=length(x1)。)。%存儲語音信號 對噪聲進行頻譜的分析[x1,fs,bits]=wavread(39。C:39。%暫停命令 wavplay(x,Fs)。)。%文字提示 fprintf(39。錄音結(jié)束39。)。x=wavrecord(time*Fs,Fs,39。錄音中......39。%文字提示 pause。按Enter鍵錄音%ds39。%抽樣頻率 time=3。%顯示縱坐標表示的意義 grid。幅值 39。)。%加標題 xlabel(39。原始男生信號的頻譜圖39。%在四行兩列的第三個窗口顯示圖形 plot(f,abs(y0(n1)))。%加網(wǎng)格y0=fft(y)。)。%加標題 ylabel(39。原始男生信號的時域波形39。plot(y)。f=[(n11)*(2*pi/N)]/pi。%返回采樣點數(shù) df=fs/N。% 對語音信號進行采樣wavplay(y,fs)。 C: 39。%顯示縱坐標表示的意義 grid。幅值 39。)。%加標題 xlabel(39。原始女生信號的頻譜圖39。plot(f,abs(y0(n1)))。%加網(wǎng)格y0=fft(x)。)。%加標題 ylabel(39。原始女生信號的時域波形39。plot(x)。%頻帶寬度 figure(1)。%采樣間隔 n1=1:N/2。%播放原始信號 N=length(x)。)。%存儲語音信號 分析男女聲音的頻譜x=wavread(39。C:39。%暫停命令 wavplay(x,Fs)。)。%文字提示 fprintf(39。錄音結(jié)束39。)。x=wavrecord(time*Fs,Fs,39。錄音中......39。%文字提示 pause。按Enter鍵錄音%ds39。%抽樣頻率 time=3。sound(xdaitong, Fs,bits)。%在Matlab中，函數(shù)sound可以對聲音進行回放，其調(diào)用格式: sound(xhigh, Fs,bits)。%播放濾波后信號 圖形如下：原始信號的時域波形幅值/A01012x 10原始信號的頻譜圖34幅值/分析：三個濾波器濾波后的聲音與原來的聲音都發(fā)生了變化。%加標題 grid。帶通濾波器濾波后信號的頻譜圖39。)。ylabel(39。頻率w/pi39。%繪制模值 axis([0 1 0 50])。%對濾波后的信號進行快速傅立葉變換 subplot(4,2,8)。%顯示縱坐標表示的意義 grid。幅值/A39。)。%加標題 xlabel(39。帶通濾波后信號的時域波形39。%在四行兩列的第六個窗口顯示圖形 plot(ys1)。%將信號送入高通濾波器濾波 figure(4)。[h1,w]=freqz(h)。band39。n=ceil(*pi)。wn2=(ws2+wc2)/2。width=ws1wc1。ws1=2*pi*1200/fs。wc1=2*pi*1000/fs。%播放濾波后信號3．3帶通濾波器 fs=8000。%標題 grid。高通濾波器濾波后信號的頻譜圖39。)。ylabel(39。頻率w/pi39。%繪制模值 axis([0 1 0 50])。%將濾波后的語音信號進行快速傅立葉變換 subplot(4,2,7)。%顯示縱坐標的表示意義 grid。幅值/A39。)。%顯示標表示的意義 title(39。幅值/A39。)。%信號的時域波形 xlabel(39。%將信號送入高通濾波器濾波 subplot(4,2,5)。[h1,w]=freqz(h)。high39。n=ceil(*pi/width)。width=wswc。ws=2*5000/fs。%播放濾波后信號3．2高通濾波器濾波 fs=8000。%標題 grid。低通濾波器濾波后信號的頻譜圖39。)。ylabel(39。頻率w/pi39。%在四行兩列的第四個窗口顯示圖形 plot(f, abs(fftwave(n1)))。%顯示標表示的意義 grid。幅值/A39。)。%加標題 xlabel(39。低通濾波后信號的時域波形39。%在四行兩列的第二個窗口顯示圖形 plot(ys)。%將濾波后的語音信號進行快速傅立葉變換 figure(4)。ys=fftfilt(h,x)。,kaiser(n+1,beta))。h=fir1(n,wn/pi,39。wn=(ws+wc)/2。ws=2*pi*1200/fs。beta=。%加標題 grid。原始信號的頻譜圖39。)。%顯示橫坐標表示的意義 ylabel(39。頻率w/pi39。)。%原始信號的頻譜圖 title(39。%快速傅立葉變換 subplot(4,2,3)。%顯示縱坐標的表示意義 grid。幅值/A39。)。%錄制信號的時域波形title(39。subplot(4,2,1)。f=[(n11)*(2*pi/N)]/pi。%返回采樣點數(shù) df=fs/N。%播放原始信號wavplay(x,fs)。C:39。grid on。耗損（dB）39。)。xlabel(39。巴特沃斯帶通濾波器39。plot(w,abs(h))。figure(4)。)。[B,A]=butter(N,wc,39。s39。was=2*tan(ws/2)。Fs1=1。ws=2*pi*fs/Fs。rs=100。Fs=8000。%帶通濾波fp=[1200,3000]。grid on。耗損（dB）39。)。xlabel(39。巴特沃斯高通濾波器39。plot(w,abs(h))。figure(1)。)。,39。[B,A]=butter(N,wc,39。s39。was=2*tan(ws/2)。Fs1=1。ws=2*pi*fs/Fs。rs=100。Fs=22050。%高通濾波fp=4800。gridon。耗損（dB）39。)。xlabel(39。巴特沃斯低通濾波器39。plot(w,abs(h))。figure(1)。)。[B,A]=butter(N,wc,39。s39。was=2*tan(ws/2)。Fs1=1。wp=2*pi*fp/Fs。rp=1。fs=1200。%畫網(wǎng)格 去除無效點MATLAB程序如下。)。ylabel(39。頻率39。)。%畫出截取后的時域圖形 title(39。, 7[4000,24000])。x=wavread(39。%加網(wǎng)格 原始信號的時域與頻譜圖針對實際發(fā)出聲音落后錄制動作半拍的現(xiàn)象，如何拔除對無效點的采樣的問題: 出現(xiàn)這種現(xiàn)象的原因主要是錄音開始時，人的反應慢了半拍，導致出現(xiàn)了一些無效點，而后而出現(xiàn)的無效的點，主要是已經(jīng)沒有聲音的動作，先讀取聲音出來，將原始語音信號時域波形圖畫出來，根據(jù)己得到的信號，可以在第二次讀取聲音的后面設定采樣點，取好有效點，畫出濾除無效點后的語音信號時域波形圖，對比可以看出。)。%顯示縱坐標表示的意義 title(39。幅值 39。)。%加標題 xlabel(39。原始信號的頻譜圖39。plot(f,abs(y0(n1)))。%快速傅立葉變換 figure(2)。%顯示縱坐標的表示意義 grid。幅值/A39。)。%錄制信號的時域波形 title(39。subplot(2,1,1)。f=[(n11)*(2*pi/N)]/pi。%返回采樣點數(shù) df=fs/N。)。%按任意鍵播放語音信號wavwrite(x,Fs,39。pause。按Enter鍵回放錄音39。)。%錄制語音信號 fprintf(39。double39。)。%暫停命令 fprintf(39。,time)。%錄音時間 fprintf(39。程序如下：Fs=8000。2Wh 或 fs179。Ws/2即連續(xù)信號是帶限的，且信號最高頻率不超過抽樣頻率的二分之一，則可不失真恢復。238。Xa(jW)=237。如果：236。TTk=165。Xa(jWjk)=229。2p1165。xa(t)DT(jW)171。?a(t)Xa(jW)171。229。165。n=165。xa(t)d(tnT)=229。，F(xiàn)IR濾波器利用函數(shù)fftfilt對信號進行濾波。=fft(x)對采集的信號做快速傅立葉變換，并用[h1,w]=freqz(h)進行DTFT變換。第一篇：數(shù)字信號處理課程設計..課程設計報告課程名稱： 數(shù)字信號處理 課題名稱： 語音信號的處理與濾波姓 名： 學 號： 院 系： 專業(yè)班級： 指導教師： 完成日期： 2013年7月2日目錄第1部分 課程設計報告………………………………………3 一．設計目的……………………………………………3 二．設計內(nèi)容……………………………………………3 三．設計原理……………………………………………3 四．具體實現(xiàn)……………………………………………5 …………………………………5 ………………………8 …………………13 ………………………………19 ……………………19 ……………………………22 五． 結(jié)果分析……………………………………………27 第2部分 課程設計總結(jié)………………………………28 一． 參考文獻……………………………………………28第1部分 課程設計報告一．設計目的綜合運用本課程的理論知識進行頻譜分析以及濾波器設計，通過理論推導得出相應結(jié)論，并利用MATLAB作為工具進行實現(xiàn)，從而復習鞏固課堂所學的理論知識，提高對所學知識的綜合應用能力，并從實踐上初步實現(xiàn)對數(shù)字信號的處理。二．設計內(nèi)容錄制一段個人自己的語音信號，并對錄制的信號進行采樣；畫出采樣后語音信號的時域波形和頻譜圖；給定濾波器的性能指標，采用窗函數(shù)法和雙線性變換法設計濾波器，并畫出濾波器的頻率響應；然后用自己設計的濾波器對采集的信號進行濾波，畫出濾波后信號的時域波形和頻譜，并對濾波前后的信號進行對比，分析信號的變化；回放語音信號；換一個與你性別相異的人錄制同樣一段語音內(nèi)容，分析兩段內(nèi)容相同的語音信號頻譜之間有什么特點；再錄制一段同樣長時間的背景噪聲疊加到你的語音信號中，分析疊加前后信號頻譜的變化，設計一個合適的濾波器，能夠把該噪聲濾除；三．設計原理，利用函數(shù)wavrecord()，wavwrite(),wavread(),wavplay()對語音信號進行錄制，存儲，讀取，回放。 DF線性相位的概念，即線性相位對h(n)、H(w)及零點的約束，了解四種FIR DF的頻響特點。連續(xù)信號經(jīng)理想抽樣后時域、頻域發(fā)生的變化（理想抽樣信號與連續(xù)信號頻譜之間的關系）理想抽樣信號能否代表原始信號、如何不失真地還原信號即由離散信號恢復連續(xù)信號的條件（抽樣定理）理想采樣過程描述： 時域描述：?a(t)=xa(t)dT(t)=229。xa(nT)d(tnT)xn=165。165。dT(t)=頻域描述：利用傅氏變換的性質(zhì)，時域相乘頻域卷積，若n=165。d(tnT)165。xXa(jW)171。dT(t)則有?(jW)=1X(jW)*D(jW)XaaT2p1165。?Xa(jW)=229。Xa(jWjkWs)Tk=165。?(jW)與X(jW)的關系：理想抽樣信號的頻譜是連續(xù)信號頻譜的Xaa周期延拓,重復周期為Ws(采樣角頻率)。X(jW)239。a239。0WWs/2W179。奈奎斯特采樣定理：要使實信號采樣后能夠不失真還原，采樣頻率必須大于信號最高頻率的兩倍：Ws179。2fh四．具體實現(xiàn)在MATLAB中用wavrecord、wavread、wavplay、wavwrite對聲音進行錄制、讀取、回放、存儲。%抽樣頻率 time=3。按Enter鍵錄音%ds39。%文字提示 pause。錄音中......39。x=wavrecord(time*Fs,Fs,39。)。錄音結(jié)束39。%文字提示 fprintf(39。)。%暫停命令wavplay(x,Fs)。C:39。%存儲語音信號N=length(x)。%采樣間隔 n1=1:N/2。%頻帶寬度 figure(2)。plot(x)。原始信號的時域波形39。%加標題 ylabel(39。)。%加網(wǎng)格y0=fft(x)。subplot(2,1,2)。%原始信號的頻譜圖 title(39。)。頻率w/pi39。%顯示橫坐標表示的意義 ylabel(39。)。原始信號的頻譜圖39。%加標題grid。這樣就可以解決這個問題。C:39。%從4000點截取到24000結(jié)束 plot(x)。截取后的聲音時域圖形39。%標題 xlabel(39。)。振幅39。grid。%低通濾波fp=1000。Fs=22050。rs=100。ws=2*pi*fs/Fs。wap=2*tan(wp/2)。[N,wc]=buttord(wap,was,rp,rs,39。)。s39。[Bz,Az]=bilinear(B,A,Fs1)。[h,w]=freqz(Bz,Az,512,Fs1*22050)。title(39。)。頻率（HZ）39。ylabel(39。)。9 巴特沃思低通濾波器MATLAB程序如下。fs=5000。rp=1。wp=2*pi*fp/Fs。T=1。wap=2*tan(wp/2)。10[N,wc]=buttord(wap,was,rp,rs,39。)。high39。s39。[Bz,Az]=bilinear(B,A,Fs1)。[h,w]=freqz(Bz,Az,512,Fs1*22050)。title(39。)。頻率（HZ）39。ylabel(39。)