Digital filter with lowpass, then filtfilt: different output

Question

Riccardo Sorvillo 2020-5-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/525236-digital-filter-with-lowpass-then-filtfilt-different-output

回答： Paul 2024-6-8

采纳的回答： Star Strider

在 MATLAB Online 中打开

Hello world!

Here I am again dealing with digital filters; I am using the latest MATLAB version (R2020a) with DSP system toolbox.

I am manipulating accelerometers data that need to be filtered by a lowpass filter: satisfactory results have been achieved with function lowpass. Since next acquistions have to pass through the very same filter, I extracted the filter built in lowpass writing

[dataflt, Hd] = lowpass(data, cut_off, fs, 'ImpulseResponse', 'fir');

Subsequent steps would be to apply Hd via filtfilt function to all the other acquisitions: as stated in past MathWorks questions, lowpass inherently uses filtfilt function. In order to double-check this statement, I filtered the data used in lowpass with filter Hd writing

dataflt = filtfilt(Hd,data);

I plotted the two outputs to compare them (in the time domain, figure below): it can be concluded that discrepancies arise at the beginning and at the end of the time history.

Why does such a difference exist? Am I neglecting anything fundamental? Is there any way to solve this problem?

I would prefer to use filtfilt (or similar) instead of lowpass because it seems to me the faster way to filter data, given that the filter is precalculated: is this statement true?

Thank you in advance for your wothwhile help.

Riccardo

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Star Strider 2020-5-13

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/525236-digital-filter-with-lowpass-then-filtfilt-different-output#answer_432284

Occasionally, filter transients can appear, such as those you are seeing. The way I deal with those (to my satisfaction at least) is to add a vector of ones multiplied by the initial value of the vector to the beginning (and similarly at the end if necessary) to the vector. I then remove these elements from the filtered signal. The length of the added elements can vary, I usually use 10 samples. Experiment to get the result you want.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Star Strider 2020-5-13

As always, my pleasure!

To use the procedure I discussed, it would not be necessary to remove any elements of the data. Add the ones vectors multiplied by the beginning (and end if necessary) values of your data, then remove those elements from the filtered data. The data itself remains intact. That is simply a way of dealing with the transients.

I believe lowpass (and related functions) use filtfilt, although they may also incorporate a way of setting the initial conditions of filtfilt to deal with the transients. (The leading and trailing vectors approach is how I deal with them.) The problem with running lowpass (or the others) in the loop as posted is that lowpass designs and implements the filter again each time it is called, and that takes time. I would design the filter once, and then use it subsequently for each signal.

The filtfilt function can also filter arrays, the constraint being that the signal lengths all have to be equal in order to create the arrays. That eliminates the need for the loop. (See the documentationn section on x for details.)

There are also ways to design filters using command-line functions. I like the filters that lowpass and the others produce, and I have learned from those designs. I do not always use them, preferring to design my own filters in some situations.

Riccardo Sorvillo 2020-5-15

Thank you very much for your valuable help: your explaination was clear and solved my problem.

Star Strider 2020-5-15

As always, my pleasure!

Thank you! I do my best!

请先登录，再进行评论。

Answer 2

Paul 2024-6-8

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/525236-digital-filter-with-lowpass-then-filtfilt-different-output#answer_1469126

在 MATLAB Online 中打开

lowpass calls filtfilt (or executes code equivalent to filtfilt) only if the underlying filter is IIR. If it's FIR, as in the question by @Riccardo Sorvillo, then lowpass "compensates for the delay introduced by the filter" by using the known delay of a FIR filter that depends on its order.

Example:

rng(100);

fs = 1e3;

t = 0:1/fs:1;

x = [1 2]*sin(2*pi*[50 250]'.*t) + randn(size(t))/10;

% lowpass() with IIR filter, same result as filtfilt()

[ylp,d] = lowpass(x,200,fs,'ImpulseResponse','iir','Steepness',0.5);

yff = filtfilt(d,x);

figure

plot(t,ylp-yff)

% lowpass() with FIR filter, different result than filtfilt()

[ylp,d] = lowpass(x,150,fs);

yff = filtfilt(d,x);

figure

plot(t,ylp-yff),grid

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Digital filter with lowpass, then filtfilt: different output

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

4 个评论
显示 2更早的评论隐藏 2更早的评论

更多回答（1 个）

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

另请参阅

类别

标签

Community Treasure Hunt

Digital filter with lowpass, then filtfilt: different output

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

4 个评论 显示 2更早的评论隐藏 2更早的评论

更多回答（1 个）

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

另请参阅

类别

标签

Community Treasure Hunt

WeChat

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

4 个评论
显示 2更早的评论隐藏 2更早的评论

0 个评论
显示 -2更早的评论隐藏 -2更早的评论