數(shù)字信號(hào)處理matlab濾波器課程設(shè)計(jì)(編輯修改稿)
2025-02-12 16:58 本頁面
【文章內(nèi)容簡(jiǎn)介】 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,2000,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(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,10
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1