基于matlab的語音信號fir數(shù)字濾波處理畢業(yè)論文(編輯修改稿)
2025-07-15 16:00 本頁面
【文章內(nèi)容簡介】 時域波形用plot(y)就可以畫出。MATLAB提供了快速傅里葉變換算法FFT計算DFT的函數(shù)fft,其調(diào)用格式如下：Xk=fft(xn,N)參數(shù)xn為被變換的時域序列向量，N是DFT變換區(qū)間長度，當(dāng)N大于xn的長度時，fft函數(shù)自動在xn后面補(bǔ)零。，當(dāng)N小于xn的長度時，fft函數(shù)計算xn的前N個元素，忽略其后面的元素。在本次設(shè)計中，我們利用fft對語音信號進(jìn)行快速傅里葉變換，就可以得到信號的頻譜特性。其程序如下：clc。clear all。[x1,Fs,bits]=wavread(39。F:\學(xué)習(xí)\新建文件夾\39。)。sound(x1,Fs)。figure(1)。subplot(2,1,1)。plot(x1)。 %做原始語音信號的時域圖形title(39。原始語音信號時域波形39。)。xlabel(39。time n39。)。ylabel(39。yinliang n39。)。axis([0,200000,])。y1=fft(x1)。 %做length(x1)點的FFTy1=fftshift(y1)。%平移，是頻率中心為0derta_Fs = Fs/length(x1)。%設(shè)置頻譜的間隔，分辨率subplot(2,1,2)。plot([Fs/2:derta_Fs: Fs/2derta_Fs],abs(y1))。%畫出原始語音信號的頻譜圖title(39。原始語音信號的頻譜39。)。axis([25000,25000,0,400])。grid on。波形為：根據(jù)語音信號的頻譜特點，大致在300~3400之間，因此考慮用帶通濾波器對噪聲部分進(jìn)行過濾。MATLAB中設(shè)計FIR濾波器的方法有窗函數(shù)法、多帶方法、最小二乘法、任意響應(yīng)法、余弦法等。本設(shè)計從比較熟知的窗函數(shù)法入手，通過不同窗函數(shù)的分析和比較進(jìn)行編程。h(n)=hd(n)w(n)W(n)所加的窗函數(shù)，通過上式將理想濾波器變?yōu)槲锢砜蓪崿F(xiàn)的濾波器。MATLAB中提供的窗函數(shù)有矩形窗、漢明窗、漢寧窗、Blackman窗、三角窗、Kaiser窗等。各種窗函數(shù)都有明顯的主瓣和旁瓣，主瓣和旁瓣的幅值衰減決定了窗函數(shù)的應(yīng)用場合。矩形窗具有最窄的主瓣，但也具有最大的旁瓣峰值。Blackman窗具有最大的旁瓣衰減，但也具有最寬的主瓣寬度。不同的窗函數(shù)在這兩方面的特點不同，因此應(yīng)根據(jù)具體的問題進(jìn)行選擇。通常來講，哈明窗和漢寧窗的主瓣具有較
范文總結(jié)相關(guān)推薦
鄂ICP備17016276號-1