unexpected shift after downsampling using decimate

Question

Ray Lee 2014-11-21

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/163762-unexpected-shift-after-downsampling-using-decimate

编辑： Vivian 2024-2-20

n = 1000;
t = 1:n;
x = rand(1,n) - 0.5; 
x = sin(2*pi/100*t);
y = decimate(x,4);
figure('color','w');
ha = axes('nextplot','add','box','on');
plot(ha,t,x,'k','marker','.');
plot(ha,t(1:4:end),y,'r','marker','o'); addkeycb;

There is a shift of 3 points.

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

请先登录，再进行评论。

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

Answer 1

Star Strider 2014-11-21

3
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/163762-unexpected-shift-after-downsampling-using-decimate#answer_159802

在 MATLAB Online 中打开

Add the default filter to produce the correct decimation:

y = decimate(x,4,'fir');

3 个评论
显示 1更早的评论隐藏 1更早的评论

Star Strider 2014-11-23

My pleasure!

It may have to do with the difference between the default IIR filter (using filtfilt) and the specified FIR filter (using filter). They have different characteristics, and one may work better in some situations than the other. (I didn’t try resample to see what results it would produce. That might be worth exploring if you’re interested.)

F S 2019-8-6

This answer solves the problem but is technically wrong. The answer from JK below gives you the real explanation and solution, in case you'd rather use the default filter.

请先登录，再进行评论。

Answer 2

Jonathan Kohler 2017-6-27

3
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/163762-unexpected-shift-after-downsampling-using-decimate#answer_272164

编辑：Jonathan Kohler 2017-6-28

在 MATLAB Online 中打开

This apparent time shift is due to MATLAB's choice of initial index for the down-sampled data, and only indirectly related to the choice of filter.

As per the documentation ( https://www.mathworks.com/help/signal/ref/decimate.html , under 'Algorithms'), the first point of the original data and downsampled data are chosen to match for FIR filters, and the last point are chosen to match for IIR filters. The reason for this choice eludes me, but maybe the provided reference explains it.

Because of the difference in handling FIR filters, there is in no shift, as pointed out by Star Strider. However, to fix this for IIR filters, you need only change your choice of time values corresponding to the downsampled data. Instead of choosing

t(1:r:end)

You should choose