數(shù)字信號(hào)處理課程設(shè)計(jì)---數(shù)字濾波器設(shè)計(jì)及在語(yǔ)音信號(hào)(編輯修改稿)
2025-02-12 14:13 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 （n，Wn，‘ftype’）； 說明：butter函數(shù)可以設(shè)計(jì)低通，帶通，高通和帶阻數(shù)字濾波器，其特性可使通帶內(nèi)的幅度響應(yīng)最大限度地平坦，但會(huì)損失截止頻率處的下降斜度，使幅度相應(yīng)衰減較慢。 [b，a]=butter（n，Wn）可以設(shè)計(jì)截止頻率為Wn的n階低通Butterworth濾波器，其中截止頻率Wn應(yīng)滿足0= Wn =1,Wn=（采樣頻率）。使用butter函數(shù)設(shè)計(jì)濾波器，可以使通帶內(nèi)的幅度響應(yīng)最大限度地平坦，但會(huì)損失截止頻率處的下降斜度，因此，butter函數(shù)主要用于設(shè)計(jì)通帶平坦的數(shù)字濾波器3 課程設(shè)計(jì)設(shè)計(jì)步驟及結(jié)果分析[x1,fs,bits]=wavread(39。d:/39。)。sound(x1,fs,bits)。y1=fft(x1,32768)。figure(1)subplot(2,1,1)plot(x1)。title(39。原始語(yǔ)音信號(hào)時(shí)域波形39。)。axis([0 230000 1 1])subplot(2,1,2)plot(abs(y1))。title(39。原始語(yǔ)音信號(hào)頻譜39。)axis([0 32768 0 1500])[x1,fs,bits]=wavread(39。d:/39。)。y1=fft(x1,32768)。x2=awgn(x1,10)。%sound(x2,fs,bits)y2=fft(x2,32768)。f=0:fs/32768:fs/32768*32767figure(2)subplot(3,2,1)plot(x2)。title(39。加單頻正弦信號(hào)的語(yǔ)音信號(hào)時(shí)域波形39。)axis([0 230000 2 2])y3=fft(x3,32768)。subplot(3,2,4)plot(f,abs(y3))。xlabel(39。頻率（ Hz）39。)。 ylabel(39。幅度（dB）39。)。title(39。加單頻正弦信號(hào)的語(yǔ)音信號(hào)頻譜39。)axis([0 23000 0 1500])N=length(x1)1。t=0:1/fs:N/fs。 a=[*sin(2*pi*2000*t)]39。 b=[*sin(2*pi*3000*t)]39。 c=[*sin(2*pi*3500*t)]39。 x4=x1+a+b+c。%sound(x4,fs,bits)y4=fft(x4,32768)。 %對(duì)加噪音后的信號(hào)取1024個(gè)點(diǎn)的fft變換subplot(3,2,5)。plot(x4) title(39。加多個(gè)正弦干擾信號(hào)后的信號(hào)時(shí)域圖39。)。axis([0 230000 2 2])ylabel(39。幅值39。)。subplot(3,2,6)。plot(f,abs(y4))。title(39。加多個(gè)正弦干擾信號(hào)后的信號(hào)頻譜圖39。)axis([0 23000 0 1500])ylabel(39。幅值39。)。xlabel(39。頻率（Hz）39。)。[y,fs,bit]=wavread(39。d:/39。)。fr=1000。wr=2*pi*fr。Window=boxcar(8)。b=fir1(7,fr/(fs/2),Window)。[h,w]=freqz(b,1)。subplot(221)。plot(w*fs/(2*pi),abs(h))。title(39。低通濾波39。)。y1=filter(b,1,y)。y2=fftfilt(b,y)。subplot(222)。plot(y)。title(39。原始語(yǔ)音信號(hào)39。)。subplot(223)。plot(y1)。title(39。IIR濾波后語(yǔ)音信號(hào)39。)。subplot(224)。plot(y2)。title(39。FIR濾波后語(yǔ)音信號(hào)39。)。sound(y1,fs,bit)。[y,fs,bit]=wavread(39。d:/39。)。fr=1000。fp=3000wr=2*pi*fr。Window=blackman(16)。b=fir1(15,[fr/(fs/2),fp/(fs/2)],Window)。[h,w]=freqz(b,1)。subplot(221)。plot(w*fs/(2*pi),abs(h))。title(39。帶通濾波39。)。y1=filter(b,1,y)。y2=fft
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1