?? fir_b.m
字號(hào):
x1=wavread('音樂.wav');fs=22050;wp1=2*pi*1200/fs;wp2=2*pi*3000/fs;ws1=2*pi*1000/fs;ws2=2*pi*3200/fs;wp=(wp1+ws1)/2;ws=(wp2+ws2)/2;wdelta=wp1-ws1;N=ceil(8*pi/wdelta); %取整wn=[wp ws];beta=3.39;window=kaiser(N+1,beta); [b,a]=fir1(N,wn/pi,'high',window); figure(1)freqz(b,a,512);title('FIR帶通濾波器');f2=filter(b,a,x1);figure(2)subplot(2,1,1)plot(x1)title('FIR帶通濾波器濾波前的時(shí)域波形');subplot(2,1,2)plot(f2);title('FIR帶通濾波器濾波后的時(shí)域波形');F0=fft(f2,1024);f=fs*(0:511)/1024;figure(3)y2=fft(x1,1024);subplot(2,1,1);plot(f,abs(y2(1:512)));title('FIR帶通濾波器濾波前的頻譜')xlabel('頻率/Hz');ylabel('幅值');subplot(2,1,2)F2=plot(f,abs(F0(1:512)));title('FIR帶通濾波器濾波后的頻譜')xlabel('頻率/Hz');ylabel('幅值'); sound(f2,44100); %濾波后的語音信號(hào)?? 快捷鍵說明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號(hào)
Ctrl + =
減小字號(hào)
Ctrl + -