blackman 布莱克曼 滤波器.docVIP

clear;close all [z1,fs,bits]=wavread(E:\shiyan_7_updata\录音2.wav); n=0:8191; y1=z1(1:8192); Y1=fft(y1); %滤波前 %%%%%%Blackman低通滤波%%%%%%%% fp=1000;fc=1200;As=100;Ap=1;Fs=8000; wc=2*pi*fc/Fs; wp=2*pi*fp/Fs; B=wc-wp; %过渡带宽度 M=ceil(12*pi/B)-1; figure; wp=(wp+wc)/2/pi; b=fir1(M,wp,blackman(M+1)); stem(b,.);xlabel(n);ylabel(h(n)); figure; [h1,w1]=freqz(b,1); plot(w1/(2*pi)*Fs,20*log10(abs(h1))); xlabel(频率/Hz);ylabel(幅值/db); x=fftfilt(b,z1); X=fft(x,8192); %滤波后 figure; subplot(2,2,1);plot(n,Y1); title(滤波前信号频谱); subplot(2,2,2);plot(n,X); title(滤波后信号频谱); subplot(2,2,3);plot(z1); title(滤波前信号波形); subplot(2,2,4);plot(x); title(滤波后信号波形); sound(x,fs,bits); clear;close all ; [z1,fs,bits]=wavread(E:\shiyan_7_updata\录音2.wav); n=0:16383; y1=z1(1:16384); Y1=fft(y1); %%%%%Blackman高通滤波%%%%%% fp=4800;fc=5000;As=100;Ap=1;Fs=16000; wc=2*pi*fc/Fs; wp=2*pi*fp/Fs; B=wc-wp; M=ceil(12*pi/B); M0=M+mod(M+1,2);%确保长度为奇数 ws=(wp+wc)/2/pi; b=fir1(M0-1,ws,high,blackman(M0)); figure; stem(b,.);xlabel(n);ylabel(h(n)); [h,w]=freqz(b,1); figure; plot(w/(2*pi)*Fs,20*log10(abs(h))); xlabel(频率/Hz);ylabel(幅值/db); x=fftfilt(b,z1); X=fft(x,16384); figure; subplot(2,2,1);plot(n,Y1); axis([0,16000,-100,100]); title(滤波前信号频谱); subplot(2,2,2);plot(n,X); axis([0,16000,-100,100]); title(滤波后信号频谱); subplot(2,2,3);plot(z1); title(滤波前信号波形); subplot(2,2,4);plot(x); title(滤波后信号波形); sound(x,fs,bits); clear;close all; [z1,fs,bits]=wavread(E:\shiyan_7_updata\录音2.wav); %sound(z1,fs,bits); n=(0:8191); y1=z1(1:8192); Y1=fft(y1); %%%%%%Blackman带通滤波%%%%%%% 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; B=wp1-wc1; N=ceil(12*pi/B)-1; ws =[(wp1+wc1)/2/pi,(wp2+wc2)/2/pi]; figure(1); b=fir1(N,ws,blackman(N+1)); stem(b,.);xlabel(n);ylabel(h(n)); figure(2); [h1,w1]=freqz(b,1); plot(w1/(2*pi)*Fs,20*log10(abs(h1))); xlabel(频率/Hz);ylabel(幅值/db); title(160409108 李楷); x=fftfilt(b,z1); X=fft(x,8192); figure(3); subplot(2,2,1);plot(n,Y1); axis([0,8000,-100,1