資源簡介
ilter濾波器稱為一維數(shù)字濾波器。filtfilt濾波器稱為零相位數(shù)字濾波。其濾波算法是基于filter而來的。只是filtfilt實(shí)現(xiàn)了零相位。其基本實(shí)現(xiàn)過程為先讓信號(hào)用filter濾波,再將信號(hào)時(shí)域反轉(zhuǎn)再次通過filter濾波,這樣兩次濾波后相位實(shí)現(xiàn)了零相位。filter濾波有明顯的延遲,filtfilt濾波延時(shí)比filter小很多。
代碼片段和文件信息
%濾波器比較
%%
%產(chǎn)生原始信號(hào)
close?all
f1=10;
f2=20;
f3=30;
Fs=150;
t=0:1/Fs:1;
y1=5*sin(2*pi*f1*t)+10*sin(2*pi*f2*t)+8*sin(2*pi*f3*t);?
%%
%濾波
fp=25;%濾掉頻率超過25Hz的信號(hào)
fs=40;
wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60;???%DF指標(biāo)(低通濾波器的通、阻帶邊界頻)
[Nwp]=ellipord(wpwsrprs);?%調(diào)用ellipord計(jì)算橢圓DF階數(shù)N和通帶截止頻率wp
[BA]=ellip(Nrprswp);??????%調(diào)用ellip計(jì)算橢圓帶通DF系統(tǒng)函數(shù)系數(shù)向量B和A
y2=filter(BAy1);?%濾波
y3=filtfilt(BAy1);
%%
%FFT頻譜分析,分析濾波前后的頻率成分
N=length(y1);
nfft=2^nextpow2(N);
fft_y1=fft(y1nfft)/nfft;
fft_y2=fft(y2nfft)/nfft;
fft_y3=fft(y3nfft)/nfft;
f=Fs/2*
評(píng)論
共有 條評(píng)論