數(shù)字信號(hào)課程設(shè)計(jì)-語音信號(hào)的處理與濾波-文庫吧
2025-01-01 16:58 本頁面
【正文】 n21,Kaiser(n21+1,beta))。[h,w]=freqz(b21,1)。plot(w/pi,abs(h))。title(39。圖——6 FIR低通濾波器39。)。頻率響應(yīng)圖為：圖——6FIR高通濾波器設(shè)計(jì)程序：ft=25000。fpts=[4800 5000]。dev=[ ]。mag=[1 0]。[n22,wn22,beta,ftype]=kaiserord(fpts,mag,dev,ft)。b22=fir1(n22,wn22,39。high39。,Kaiser(n22+1,beta))。[h,w]=freqz(b22,1)。plot(w/pi,abs(h))。title(39。圖——7 FI高通濾波器39。)。頻率響應(yīng)圖為：圖——7FIR帶通濾波器的設(shè)計(jì)程序：fp1=1200。fp2=3000。fs1=1000。fs2=3200。ft=10000。ws1=fs1/(ft/2)*pi。ws2=fs2/(ft/2)*pi。wp1=fp1/(ft/2)*pi。wp2=fp2/(ft/2)*pi。wp=[wp1,wp2]。ws=[ws1,ws2]。deltaw=wp1ws1。N0=ceil(11*pi/deltaw)。N=N0+mod(N0+1,2)。windows=(blackman(N))39。wc1=(ws1+wp1)/2。wc2=(ws2+wp2)/2。hd=ideallp(wc2,N)ideallp(wc1,N)。b=hd.*windows。[h,w]=freqz(b,1)。plot(w/pi,abs(h))。 title(39。圖——8 FIR帶通濾波器39。)。 頻率響應(yīng)圖為：圖——84. 對(duì)“”語音信號(hào)進(jìn)行濾波 分別采用設(shè)計(jì)好的低通，高通，帶通濾波器對(duì)“”語音信號(hào)進(jìn)行濾波，用MATLAB畫出濾波后的時(shí)頻域圖，并回放錄播后的語音信號(hào)。（1） FIR低通濾波程序：ft=5500。fpts=[1000 1200]。dev=[ ]。mag=[1 0]。[n21,wn21,beta,ftype]=kaiserord(fpts,mag,dev,ft)。b21=fir1(n21,wn21,Kaiser(n21+1,beta))。[h,w]=freqz(b21,1)。[y,fs,nbits]=wavread(‘39。)。%加載語音信號(hào)f=filter(b21,1,y)。wavplay(f)。f1=fft(f)。subplot(2,1,1)plot(f)。title(39。低通濾波后的波形39。)。subplot(2,1,2)plot(abs(f1))。title(39。低通濾波后的頻譜圖39。)運(yùn)行結(jié)果圖為附件：圖——9（2） FIR高通濾波程序：ft=25000。fpts=[4800 5000]。dev=[ ]。mag=[1 0]。[n22,wn22,beta,ftype]=kaiserord(fpts,mag,dev,ft)。b22=fir1(n22,wn22,39。high39。,Kaiser(n22+1,beta))。[h,w]=freqz(b22,1)。[y,fs,nbits]=wavread(39。39。)。%加載語音信號(hào)f=filter(b22,1,y)。wavplay(f)。f1=fft(f)。subplot(2,1,1)plot(f)。title(39。高通濾波后的波形39。)。subplot(2,1,2)plot(abs(f1))。title(39。高通濾波后的頻譜圖39。)運(yùn)行結(jié)果圖為附件：圖——10（3） FIR帶通濾波程序：fp1=1200。fp2=3000。fs1=1000。fs2=3200。ft=10000。ws1=fs1/(ft/2)*pi。ws2=fs2/(ft/2)*pi。wp1=fp1/(ft/2)*pi。wp2=fp2/(ft/2)*pi。wp=[wp1,wp2]。ws=[ws1,ws2]。deltaw=wp1ws1。N0=ceil(11*pi/deltaw)。N=N0+mod(N0+1,2)。windows=(blackman(N))39。wc1=(ws1+wp1)/2。wc2=(ws2+wp2)/2。hd=ideallp(wc2,N)ideallp(wc1,N)。b=hd.*windows。[h,w]=freqz(b,1)。 [y,fs,nbits]=wavread(39。39。)。%加載語音信號(hào) f=filter(b,1,y)。 wavplay(f)。 f1=fft(f)。 subplot(2,1,1) plot(f)。 title(39。帶通濾波后的波形39。)。 subplot(2,1,2) plot(abs(f1))。 title(39。帶通濾波后的頻譜圖39。) 運(yùn)行結(jié)果圖為附件：圖——11（4） IIR低通濾波程序：fp=1000。 ft=5500。%取抽樣頻率 fs=1200。 wp=2*pi*fp/ft。%通帶頻率 ws=2*pi*fs/ft。%截止頻率 Omgp=2*ft*tan(wp/2)。 Omgs=2*ft*tan(ws/2)。 [n,Omgc]=buttord(Omgp,Omgs,1,100,39。s39。)。 [z,p,k]=buttap(n)。 B=k*real(poly(z))。 A=real(poly(p))。 [b,a]=lp2lp(B,A,Omgc)。 [ba,aa]=bilinear(b,a,ft)。[Ha,w]=freqz(ba,aa)。[y,fs,nbits]=wavread(39。39。)。%加載語音信號(hào)f=filter(ba,aa,y)。wavplay(f)。f1=fft(f)。subplot(2,1,1)plot(f)。title(39。低通濾波后的波形39。)。subplot(2,1,2)plot(abs(f1))。title(39。低通濾波后的頻譜圖39。)運(yùn)行結(jié)果圖為附件：圖——9——1（5） IIR高通濾波程序： fp=5000。 ft=25000。%取抽樣頻率 fs=4800。 Rp=1。As=100。 wp=2*pi*fp/ft。%通帶頻率 ws=2*pi*fs/ft。%截止頻率 Omgp=2*ft*tan(wp/2)。 Omgs=2*ft*tan(ws/2)。 [n,Omgc]=ellipord(Omgp,Omgs,Rp,As,39。s39。)。 %計(jì)算階數(shù)n和截止頻率 [z,p,k]=ellipap(n,Rp,As)。 B=k*real(poly(z))。 A=real(poly(p))。 [b,a]=lp2hp(B,A,Omgc)。 [ba,aa]=bilinear(b,a,ft)。[Ha,w]=freqz(ba,aa)。[y,fs,nbits]=wavread(39。39。)。%加載語音信號(hào)f=filter(ba,aa,y)。wavplay(f)。f1=fft(f)。subplot(2,1,1)plot(f)。title(39。IIR高通濾波后的波形39。)。subplot(2,1,2)plot(abs(f1))。title(39。IIR高通濾波
規(guī)章制度相關(guān)推薦
鄂ICP備17016276號(hào)-1