低通
fp=1000;fc=1200;As=100;Ap=1;Fs=3000;
wc=2*pi*fc/Fs; wp=2*pi*fp/Fs;
wd=wc-wp;
beta=0。
1102*(As-8。7);
N=ceil((As-7。95)/2。286/wd);
wn= kaiser(N 1,beta);
ws=(wp wc)/2/pi;
b=fir1(N,ws,wn);
figure(1);
freqz(b,1);
x=fftfilt(b,y);
X=fft(x);
figure(2);
subplot(2,2,1);plot(abs(Y1));axis([0,2000,0,100。
0]);
title('滤波前信号频谱');
subplot(2,2,2);plot(abs(X));axis([0,2000,0,100。0]);
title('滤波后信号频谱');
subplot(2,2,3);plot(y);
title('滤波前信号波形');
subplot(2,2,4);plot(x);
title('滤波后信号波形');
sound(x,fs,bits);
带通
fp1=1200 ;fp2=3000 ;fc1=1000 ;fc2=3200 ;As=100 ;Ap=1 ;Fs=35000 ;
wp1=2*pi*fp1/Fs; wc1=2*pi*fc1/Fs; wp2=2*pi*fp2/Fs; wc2=2*pi*fc2/Fs;
wdel=wp1-wc1;
beta=0。
1102*(As-8。7);
N=ceil((As-7。95)/2。286/wdel);
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);
figure(2);
subplot(2,2,1);plot(abs(Y1));axis([0,2000,0,50。
0]);
title('滤波前信号频谱');
subplot(2,2,2);plot(abs(X));axis([0,2000,-1。0,1。0]);
title('滤波后信号频谱')
subplot(2,2,3);plot(z1);
title('滤波前信号波形');
subplot(2,2,4);plot(x);
title('滤波后信号波形');
sound(x,fs,bits);
高通
fp=4800;fc=5000;As=100;Ap=1;Fs=200000;
wc=2*pi*fc/Fs; wp=2*pi*fp/Fs;
wdel=wc-wp;
beta=0。
1102*(As-8。7);
N=ceil((As-7。95)/2。286/wdel);
wn= kaiser(N,beta);
ws=(wp wc)/2/pi;
b=fir1(N-1,ws,'high',wn);
figure(1);
freqz(b,1);
x=fftfilt(b,z1);
X=fft(x);
figure(2);
subplot(2,2,1);plot(abs(Y1));axis([0,2000,0,50。
0]);
title('滤波前信号频谱');
subplot(2,2,2);plot(abs(X));axis([0,2000,-1。0,1。
0]);
title('滤波后信号频谱');
subplot(2,2,3);plot(z1);
title('滤波前信号波形');
subplot(2,2,4);plot(x);
title('滤波后信号波形');
sound(x,fs,bits);。
全部