數(shù)字信號(hào)處理matlab濾波器課程設(shè)計(jì)(編輯修改稿)
2024-07-12 00:46 本頁面
【文章內(nèi)容簡(jiǎn)介】 91。 plot(n,Y1)。 圖像輸出如圖 2： 圖 2 語音信號(hào)頻譜分析圖 3. 設(shè)計(jì)數(shù)字濾波器和對(duì)信號(hào)濾波 （ 1）窗函數(shù)設(shè)計(jì)低通濾波器 程序設(shè)計(jì)如下： clear。close all [z1,fs,bits]=wavread(39。d:\39。) y1=z1(1:8192)。 Y1=fft(y1)。 fp=1000。fc=1200。As=100。Ap=1。Fs=8000。 wc=2*pi*fc/Fs。 wp=2*pi*fp/Fs。 wdel=wcwp。 beta=*()。 N=ceil((As8)/)。 wn= kaiser(N+1,beta)。 ws=(wp+wc)/2/pi。 b=fir1(N,ws,wn)。 figure(1)。 freqz(b,1)。 x=fftfilt(b,z1)。 X=fft(x,8192)。 figure(2)。 subplot(2,2,1)。plot(abs(Y1))。axis([0,1000,0,])。 title(39。濾波前信號(hào)頻譜 39。)。 subplot(2,2,2)。plot(abs(X))。axis([0,1000,0,])。 title(39。濾波后信號(hào)頻譜 39。)。 subplot(2,2,3)。plot(z1)。 title(39。濾波前信號(hào)波形 39。)。 subplot(2,2,4)。plot(x)。 title(39。濾波前信號(hào)波 形 39。)。 sound(x,fs,bits)。 圖形分析如下： （ 2）窗函數(shù)設(shè)計(jì)高通濾波器 程序設(shè)計(jì)如下： clear。close all [z1,fs,bits]=wavread(39。d:\39。) y1=z1(1:8192)。 Y1=fft(y1)。 fp=4800。fc=5000。As=100。Ap=1。Fs=8000。 wc=2*pi*fc/Fs。 wp=2*pi*fp/Fs。 wdel=wcwp。 beta=*()。 N=ceil((As8)/)。 wn= kaiser(N,beta)。 ws=(wp+wc)/2/pi。 b=fir1(N1,ws,39。high39。,wn)。 figure(1)。 freqz(b,1)。 x=fftfilt(b,z1)。 X=fft(x,8192)。 figure(2)。 subplot(2,2,1)。plot(abs(Y1))。axis([0,1000,0,])。 title(39。濾波前信號(hào)頻譜 39。)。 subplot(2,2,2)。plot(abs(X))。axis([0,1000,0,])。 title(39。濾波后信號(hào)頻譜 39。)。 subplot(2,2,3)。plot(z1)。 title(39。濾波前信號(hào)波形 39。)。 subplot(2,2,4)。plot(x)。 title(39。濾波前信號(hào)波形 39。)。 sound(x,fs,bits)。 圖形分析如下： （ 3）窗函數(shù)設(shè)計(jì)帶通濾波器 程序設(shè)計(jì)如下： clear。close all [z1,fs,bits]=wavread(39。d:\39。) y1=z1(1:8192)。 Y1=fft(y1)。 fp1=1200 。fp2=3000 。fc1=1000 。fc2=3200 。As=100 。Ap=1 。Fs=8000 。 wp1=2*pi*fp1/Fs。 wc1=2*pi*fc1/Fs。 wp2=2*pi*fp2/Fs。 wc2=2*pi*fc2/Fs。 wdel=wp1wc1。 beta=*()。 N=ceil((As8)/)。 ws =[(wp1+wc1)/2/pi,(wp2+wc2)/2/pi]。 wn= kaiser(N+1,beta)。 b=fir1(N,ws,wn)。 figure(1)。 freqz(b,1) x=fftfilt(b,z1)。 X=fft(x,8192)。 figure(2)。 subplot(2,2,1)。plot(abs(Y1))。axis([0,1000,0,])。 title(39。濾波前信號(hào)頻譜 39。)。 subplot(2,2,2)。plot(abs(X))。axis([0,2021,0,])。 title(39。濾波后信號(hào)頻譜 39。) subplot(2,2,3)。plot(z1)。 title(39。濾波前信號(hào)波形 39。)。 subplot(2,2,4)。plot(x)。 title(39。濾波前信號(hào)波形 39。)。 sound(x,fs,bits)。 圖形分析如下： （ 4）雙線性變換法設(shè)計(jì)低通濾波器 ①選用 butter 程序設(shè)計(jì)如下： clear。close all [z1,fs,bits]=wavread(39。d:\39。) y1=z1(1:8192)。 Y1=fft(y1)。 fp=1000。fc=1200。As=100。Ap=1。Fs=8000。 wc=2*fc/Fs。 wp=2*fp/Fs。 [N,ws]=buttord(wc,wp,Ap,As)。 [b,a]=butter(N,ws)。 figure(1)。 freqz(b,a,512,Fs)。 x=filter(b,a,z1)。 X=fft(
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1