【文章內(nèi)容簡介】 *pi*1200/fs。ws2=2*pi*3000/fs。width=ws1wc1。wn1=(ws1+wc1)/2。wn2=(ws2+wc2)/2。wn=[wn1 wn2]。n=ceil(*pi)。h=fir1(n,wn/pi,39。band39。,kaiser(n+1,beta))。[h1,w]=freqz(h)。ys1= fftfilt(h,x)。%將信號送入高通濾波器濾波 figure(4)。subplot(4,2,6)。%在四行兩列的第六個(gè)窗口顯示圖形 plot(ys1)。%繪制后信號的時(shí)域的圖形 title(39。帶通濾波后信號的時(shí)域波形39。)。%加標(biāo)題 xlabel(39。頻率w/pi39。)。ylabel(39。幅值/A39。)。%顯示縱坐標(biāo)表示的意義 grid。%網(wǎng)格fftwave=fft(ys1)。%對濾波后的信號進(jìn)行快速傅立葉變換 subplot(4,2,8)。%在四行兩列的第八個(gè)窗口顯示圖形plot(f, abs(fftwave(n1)))。%繪制模值 axis([0 1 0 50])。xlabel(39。頻率w/pi39。)。ylabel(39。幅值/A39。)。%顯示標(biāo)表示的意義 title(39。帶通濾波器濾波后信號的頻譜圖39。)。%加標(biāo)題 grid。%網(wǎng)格wavplay(ys1,8000)。%播放濾波后信號 圖形如下：原始信號的時(shí)域波形幅值/A01012x 10原始信號的頻譜圖34幅值/分析：三個(gè)濾波器濾波后的聲音與原來的聲音都發(fā)生了變化。其中低通的濾波后與原來聲音沒有很大的變化，其它兩個(gè)都又明顯的變化sound(xlow,Fs,bits)。%在Matlab中，函數(shù)sound可以對聲音進(jìn)行回放，其調(diào)用格式: sound(xhigh, Fs,bits)。%sound(x, Fs, bits)。sound(xdaitong, Fs,bits)。 錄制一段異性的聲音進(jìn)行頻譜分析Fs=8000。%抽樣頻率 time=3。%錄音時(shí)間 fprintf(39。按Enter鍵錄音%ds39。,time)。%文字提示 pause。%暫停命令 fprintf(39。錄音中......39。)。x=wavrecord(time*Fs,Fs,39。double39。)。%錄制語音信號 fprintf(39。錄音結(jié)束39。)。%文字提示 fprintf(39。按Enter鍵回放錄音39。)。pause。%暫停命令 wavplay(x,Fs)。%按任意鍵播放語音信號wavwrite(x,Fs,39。C:39。)。%存儲語音信號 分析男女聲音的頻譜x=wavread(39。 C: 39。)。%播放原始信號,解決落后半拍wavplay(x,fs)。%播放原始信號 N=length(x)。%返回采樣點(diǎn)數(shù) df=fs/N。%采樣間隔 n1=1:N/2。f=[(n11)*(2*pi/N)]/pi。%頻帶寬度 figure(1)。subplot(2,2,1)。plot(x)。%錄制信號的時(shí)域波形title(39。原始女生信號的時(shí)域波形39。)。%加標(biāo)題 ylabel(39。幅值/A39。)。%顯示縱坐標(biāo)的表示意義 grid。%加網(wǎng)格y0=fft(x)。%快速傅立葉變換 subplot(2,2,2)。plot(f,abs(y0(n1)))。%原始信號的頻譜圖 title(39。原始女生信號的頻譜圖39。)。%加標(biāo)題 xlabel(39。頻率w/pi39。)。%顯示橫坐標(biāo)表示的意義 ylabel(39。幅值 39。)。%顯示縱坐標(biāo)表示的意義 grid。%加網(wǎng)格[y,fs,bits]=wavread(39。 C: 39。)。% 對語音信號進(jìn)行采樣wavplay(y,fs)。%播放原始信號 N=length(y)。%返回采樣點(diǎn)數(shù) df=fs/N。%采樣間隔 n1=1:N/2。f=[(n11)*(2*pi/N)]/pi。%頻帶寬度 subplot(2,2,3)。plot(y)。%錄制信號的時(shí)域波形title(39。原始男生信號的時(shí)域波形39。)。%加標(biāo)題 ylabel(39。幅值/A39。)。%顯示縱坐標(biāo)的表示意義 grid。%加網(wǎng)格y0=fft(y)。%快速傅立葉變換subplot(2,2,4)。%在四行兩列的第三個(gè)窗口顯示圖形 plot(f,abs(y0(n1)))。%原始信號的頻譜圖 title(39。原始男生信號的頻譜圖39。)。%加標(biāo)題 xlabel(39。頻率w/pi39。)。%顯示橫坐標(biāo)表示的意義 ylabel(39。幅值 39。)。%顯示縱坐標(biāo)表示的意義 grid。%加網(wǎng)格分析：就時(shí)域圖看，男生的時(shí)域圖中振幅比女生的高，對于頻譜圖女生的高頻成分比較多Fs=8000。%抽樣頻率 time=3。%錄音時(shí)間 fprintf(39。按Enter鍵錄音%ds39。,time)。%文字提示 pause。%暫停命令 fprintf(39。錄音中......39。)。x=wavrecord(time*Fs,Fs,39。double39。)。%錄制語音信號fprintf(39。錄音結(jié)束39。)。%文字提示 fprintf(39。按Enter鍵回放錄音39。)。pause。%暫停命令 wavplay(x,Fs)。%按任意鍵播放語音信號 wavwrite(x,Fs,39。C:39。)。%存儲語音信號 對噪聲進(jìn)行頻譜的分析[x1,fs,bits]=wavread(39。 C: 39。)。%對語音信號進(jìn)行采樣wavplay(x1,fs)。%播放噪聲信號 N=length(x1)。%返回采樣點(diǎn)數(shù) df=fs/N。%采樣間隔n1=1:N/2。f=[(n11)*(2*pi/N)]/pi。%頻帶寬度 figure(5)。subplot(3,2,1)。plot(x1)。%信號的時(shí)域波形 title(39。噪聲信號的時(shí)域波形39。)。grid。ylabel(39。幅值/A39。)。y0=fft(x1)。%快速傅立葉變換subplot(3,2,2)。plot(f,abs(y0(n1)))。%噪聲信號的頻譜圖 ylabel(39。幅值39。)。title(39。噪聲信號的頻譜圖39。)。fs=8000。[x,fs,bits]=wavread(39。 C: 39。)。%對錄入信號進(jìn)行采樣[x1,fs,bits]=wavread(39。 C: 39。)。%對噪聲信號進(jìn)行采樣yy=x+x1。%將兩個(gè)聲音疊加：wavplay(yy,fs)。%播放疊加后信號 N=length(yy)。%返回采樣點(diǎn)數(shù) df=fs/N。%采樣間隔 n1=1:N/2。f=[(n11)*(2*pi/N)]/pi。%頻帶寬度 figure(5)。subplot(3,2,3)。plot(yy,39。LineWidth39。,2)。%信號的時(shí)域波形title(39。疊加信號的時(shí)域波形39。)。xlabel(39。時(shí)間/t39。)。ylabel(39。幅值/A39。)。grid。y0=fft(yy)。%快速傅立葉變換 subplot(3,2,4)。plot(f,abs(y0(n1)))。%疊加信號的頻譜圖 title(39。疊加信號的頻譜圖39。)。xlabel(39。頻率w/pi39。)。ylabel(39。幅值/db39。)。grid。 設(shè)計(jì)一個(gè)合適的濾波器將噪聲濾除 fs=18000。%采樣頻率 Wp=2*1000/fs。%通帶截至頻率 Ws=2*2000/fs。%阻帶截至頻率 Rp=1。%最大衰減 Rs=100。%最小衰減[N,Wn]=buttord(Wp,Ws,Rp,Rs)。%buttord函數(shù)（n為階數(shù)，Wn為截至頻率）[num,den]=butter(N,Wn)。%butter函數(shù)（num為分子系數(shù)den為分母系數(shù)）[h,w]=freqz(num,den)。%DTFT變換ys=filter(num,den,yy)。%信號送入濾波器濾波，ys為輸出 fftwave=fft(ys)。%將濾波后的語音信號進(jìn)行快速傅立葉變換 figure(5)。subplot(3,2,5)。plot(ys)。%信號的時(shí)域波形title(39。低通濾波后信號的時(shí)域波形39。)。%加標(biāo)題 ylabel(39。幅值/A39。)。%顯示標(biāo)表示的意義 grid。%網(wǎng)格 subplot(3,2,6)。plot(f, abs(fftwave(n1)))。%繪制模值 title(39。低通濾波器濾波后信號的頻譜圖39。)。%標(biāo)題 xlabel(39。頻率w/pi39。)。ylabel(39。幅值/A39。)。%顯示標(biāo)表示的意義 grid。%加網(wǎng)格wavplay(ys,8000)。%播放濾波后信號 grid。圖形如下：噪聲信號的時(shí)域波形1100噪聲信號的頻譜圖幅值/ 10疊加信號的時(shí)域波形10101時(shí)間/t2200幅值/db34幅值/，常會出現(xiàn)實(shí)際發(fā)出聲音落后錄制動作半拍，可在[x,fs,bits]=wavread(39。d:39。)加 窗[x,fs,bits]=wavread(39。d:39。，[100 10000])，窗的長度可根據(jù)需要定義。，把高頻濾除，聲音變得比較低沉。當(dāng)通過高通濾波器后，把低頻濾除，聲音變得比較就尖銳。通過帶通濾波器后，聲音比較適中。：時(shí)域圖的振幅大小與性別無關(guān)，只與說話人音量大小有關(guān)，音量越大，振幅越大。頻率圖中，女生高 27頻成分較多。，噪聲與原信號明顯區(qū)分,但通過低通濾波器后，噪聲沒有濾除，信號產(chǎn)生失真。原因可能為噪聲與信號頻率相近無法濾除。第2部分 課程設(shè)計(jì)總結(jié)通過本次課程設(shè)計(jì)，使我們對數(shù)字信號處理相關(guān)知識有了更深刻的理解，尤其是對各種濾波器的設(shè)計(jì)。在設(shè)計(jì)的過程中遇到了很多問題，剛剛開始時(shí)曾天真的認(rèn)為只要把以前的程序改了參數(shù)就可以用了，可是問題沒有我想象中的那么簡單，單純的搬程序是不能解決問題的。通過查閱資料和請教同學(xué)收獲了很多以前不懂的理論知識。再利用所學(xué)的操作，發(fā)現(xiàn)所寫的程序還是沒有能夠運(yùn)行，通過不斷地調(diào)試，運(yùn)行，最終得出了需要的結(jié)果。整個(gè)過程中學(xué)到了很多新的知識，特別是對Matlab的使用終于有些了解。在以后的學(xué)習(xí)中還需要深入了解這方面的內(nèi)容。在這次的課程設(shè)計(jì)中讓我體會最深的是：知識來不得半點(diǎn)的馬虎。也認(rèn)識到自己的不足，以后要進(jìn)一步學(xué)習(xí)。八．參考文獻(xiàn)[1]數(shù)字信號處理教程(第三版)程佩青 清華大學(xué)出版社 [2]MATLAB信號處理 劉波 文忠 電子工業(yè)出版社 [3] 王宏 清華大學(xué)出版社[4]MATLAB基礎(chǔ)與編程入門 張威 西安電子科技大學(xué)出版社[5] 數(shù)字信號處理及其MATLAB實(shí)驗(yàn) 趙紅怡 張常 化學(xué)工業(yè)出版社[6]MATLAB信號處理詳解 陳亞勇等 人民郵電出版社 [7] 數(shù)字信號處理錢同惠 機(jī)械工業(yè)出版社 29第四篇：數(shù)字信號處理課程設(shè)計(jì)目 錄摘要...........................................................................................................................................1 1 緒論..............................................................................................................................................2 DSP系統(tǒng)特點(diǎn)和設(shè)計(jì)基本原則......................................................................................2 國內(nèi)外研究動態(tài).............................................................................................................2 2系統(tǒng)設(shè)計(jì)........................................................................................................................................3 3硬件設(shè)計(jì)........................................................................................................................................5 硬件結(jié)構(gòu)...........................................................................................................................5 硬件電路設(shè)計(jì)...................................................................................................................7 總輸入電路.....................................................