Time-domain frequency filter

16 vues (au cours des 30 derniers jours)
Henry
Henry le 16 Avr 2012
I have a set of time-domain data that I wish to filter. I need to use some sort of notch filter so as no delay occurs as a result of the filtering. I only want to use the data with frequencies between 0 and 3 hz to analyse afterwards. The end result must be in the time-domain - I have tried converting to the frequency domain, selecting the appropriate frequencies and then transferring back into the time-domain, which it does not let me do.
How do I perform this procedure totally in the time-domain?

Réponse acceptée

Wayne King
Wayne King le 16 Avr 2012
An example (requires Signal Processing Toolbox)
d = fdesign.lowpass('Fp,Fst,Ap,Ast',3,4,0.5,50,10);
Hd = design(d);
output = filtfilt(Hd.Numerator,1,inputdata);

Plus de réponses (3)

Wayne King
Wayne King le 16 Avr 2012
Sounds like you just need a lowpass filter since you want to only analyze data below 3 Hz. What is your sampling frequency?
You can do what you want by designing your filter using fdesign.lowpass and then using filtfilt() to filter the data.
filtfilt() implements zero-phase filtering and operates in the time domain and returns time-domain data.
  1 commentaire
Henry
Henry le 16 Avr 2012
my sampling rate is 10Hz, thanks for your help

Connectez-vous pour commenter.


Henry
Henry le 16 Avr 2012
What is the 0.5 and 50 stand for in the design?

Eshwar Sachi
Eshwar Sachi le 3 Mai 2012
As indicated in the function, 0.5 is the pass band ripple and 50 is the stop band attenuation, in units dB. You may find the following link useful for designing filters in MATLAB: http://www.mathworks.com/discovery/filter-design.html

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by