數(shù)字信號處理matlab濾波器課程設(shè)計(編輯修改稿)
2025-07-12 00:46 本頁面
【文章內(nèi)容簡介】 91。 plot(n,Y1)。 圖像輸出如圖 2： 圖 2 語音信號頻譜分析圖 3. 設(shè)計數(shù)字濾波器和對信號濾波 （ 1）窗函數(shù)設(shè)計低通濾波器 程序設(shè)計如下： 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。濾波前信號頻譜 39。)。 subplot(2,2,2)。plot(abs(X))。axis([0,1000,0,])。 title(39。濾波后信號頻譜 39。)。 subplot(2,2,3)。plot(z1)。 title(39。濾波前信號波形 39。)。 subplot(2,2,4)。plot(x)。 title(39。濾波前信號波 形 39。)。 sound(x,fs,bits)。 圖形分析如下： （ 2）窗函數(shù)設(shè)計高通濾波器 程序設(shè)計如下： 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。濾波前信號頻譜 39。)。 subplot(2,2,2)。plot(abs(X))。axis([0,1000,0,])。 title(39。濾波后信號頻譜 39。)。 subplot(2,2,3)。plot(z1)。 title(39。濾波前信號波形 39。)。 subplot(2,2,4)。plot(x)。 title(39。濾波前信號波形 39。)。 sound(x,fs,bits)。 圖形分析如下： （ 3）窗函數(shù)設(shè)計帶通濾波器 程序設(shè)計如下： 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。濾波前信號頻譜 39。)。 subplot(2,2,2)。plot(abs(X))。axis([0,2021,0,])。 title(39。濾波后信號頻譜 39。) subplot(2,2,3)。plot(z1)。 title(39。濾波前信號波形 39。)。 subplot(2,2,4)。plot(x)。 title(39。濾波前信號波形 39。)。 sound(x,fs,bits)。 圖形分析如下： （ 4）雙線性變換法設(shè)計低通濾波器 ①選用 butter 程序設(shè)計如下： 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è)計相關(guān)推薦
鄂ICP備17016276號-1