I have a datafile contains data sets (size 200* 15) .how can I apply on all the data a lowpass butterworth filter of 2-order with cutoff fre 10 HZ

조회 수: 1(최근 30일)
Rahul Pandey
Rahul Pandey 2016년 12월 29일
댓글: David J. Mack 2017년 1월 11일
How can I call filter function like [b,a] = butter(n,Wn) or [b,a] = butter(n,Wn,ftype) for whole data file and collect them in another for further use.

채택된 답변

David J. Mack
David J. Mack 2016년 12월 29일
편집: David J. Mack 2016년 12월 29일
Hey Rahul,
use FILTER with the coefficients created with BUTTER, e.g.:
%Assuming X is your 200x15 data matrix and fSInHz is the sampling rate of the data.
filtFCInHz = 10; %Cut-off frequency.
filtOrder = 2; %Order.
[b,a] = butter(filtOrder,filtFCInHz/(fSInHz/2)); %Normalize against Nyquist frequency.
Y = filter(b,a,X); %Y is the filtered 200x15 output.
If you want to apply zero-phase filtering (which is recommended for time series), use FILTFILT instead of FILTER.
Greetings, David
  댓글 수: 4
David J. Mack
David J. Mack 2017년 1월 11일
Do you have a time vector for your data? If so, use:
fSInHz = mode(diff(t,[],1),1);
Assuming t is a 200xk column-oriented array with the time for each sample in seconds. If k==1 (uniform sampling rate) fSInHz is a scalar, if k==15 (different sampling rates for each column in X) it is a 1x15 row vector.

댓글을 달려면 로그인하십시오.

추가 답변(0개)

Community Treasure Hunt

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

Start Hunting!

Translated by