基于matlab的語音信號采集和窗函數(shù)法濾波42784(已改無錯字)
2023-07-19 14:45:49 本頁面
【正文】 er）窗等。對于窗函數(shù)的選擇，應(yīng)考慮被分析信號的性質(zhì)與處理要求。如果僅要求精確讀出主瓣頻率，而不考慮幅值精度，則可選用主瓣寬度比較窄而便于分辨的矩形窗，例如測量物體的自振頻率等；如果分析窄帶信號，且有較強(qiáng)的干擾噪聲，則應(yīng)選用旁瓣幅度小的窗函數(shù)，如漢寧窗、三角窗等；對于隨時間按指數(shù)衰減的函數(shù)，可采用指數(shù)窗來提高信噪比。 MATLAB窗函數(shù)的實現(xiàn)利用MATLAB實現(xiàn)窗函數(shù)法設(shè)計FIR濾波器，主要是選擇合適的窗函數(shù)進(jìn)行截斷運算。先從理論上得到待逼近理想濾波器的單位脈沖響應(yīng)，再由通帶、阻帶衰減指標(biāo)確定窗函數(shù)類型，由過渡帶確定FIR濾波器階數(shù)M，最后利用MATLAB計算出窗函數(shù)的值，以及hd[k]Wn[k]的值，由此即得所設(shè)計的FIR濾波器的h[k]。MATLAB提供了許多常用的窗函數(shù)，其中部分窗函數(shù)的調(diào)用形式為 W=hanning(N) W=hamming(N) W=Blackman(N) W=Kaiser(N,beta)其中N是窗函數(shù)的長度，beta是控制kaiser窗形狀的參數(shù)。返回的變量w是一個長度為N的列向量，給出窗函數(shù)N點的取值。對于kaiser窗，MATLAB還提供了一個根據(jù)帶設(shè)計濾波器的指標(biāo)直接計算窗函數(shù)的函數(shù)，其基本調(diào)用形式為 [M,Wc,beta,ftype]=kaiserord(f,a,dev)其中：參數(shù)f表示需要設(shè)計的FIR濾波器的B個頻帶；參數(shù)a為B個頻帶的幅度值，一個通帶取1，阻帶取0；參數(shù)dev為B個頻帶中的波動值。 利用hanning窗設(shè)計低通濾波器 由題意可知：要設(shè)計濾波器濾除信號的高頻成分，即設(shè)計低通濾波器。故確定設(shè)計的低通濾波器的設(shè)計指標(biāo)為：Wp=*pi,Ws=*pi 程序如下：fs=8000。x1=wavread(39。C:\WINDOWS\Media\Windows XP 39。)。 t=(0:length(x1)1)/8000。f=fs*(0:2047)/4096。A1=。A2=。d=[A1*cos(2*pi*3600*t)+A2*sin(2*pi*3800*t)]39。x2=x1+[d d]。wp=*pi。ws=*pi。wdelta=wswp。N=ceil(*pi/wdelta)。 %取整wn=(+)*pi/2。[bz,az]=fir1(N,wn/pi,hanning(N+1))。 %選擇窗函數(shù)并歸一化截止頻率figure(1)freqz(bz,az)。grid on。axis tight。f2=filter(bz,az,x2)。figure(2)subplot(2,1,1)plot(t,x2)。grid on。axis tight。title(39。濾波前的時域波形39。)。subplot(2,1,2)plot(t,f2)。grid on。axis tight。title(39。濾波后的時域波形39。)。y3=fft(f2,4096)。f=fs*(0:2047)/4096。figure(3)y2=fft(x2,4096)。subp
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1