What indices of the original sequence are retained in the output of the decimate function?

Question

Gaëtan Garcia 2022-11-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1866583-what-indices-of-the-original-sequence-are-retained-in-the-output-of-the-decimate-function

评论： Gaëtan Garcia 2022-11-29

I am using the decimate function y = decimate(x,r) and I need to know to which index in x the value y(i) corresponds.

This is not explicitly addressed in the documentation. In the first example the number of points is k*r+1 and according to the code it seems y(1) corresponds to index 1 in x, while and y(length(y)) correponds to index length(x) in x.

In the second example of the documentation, where length(x) is not of the form k*r+1, it seems y(length(y)) corresponds to the last index of the original vector x.

Is that always the case? Would that depend on the type of filter applied (e.g. Chebyshev vs FIR) ?

Thanks for any input on the matter.

Regards.

G. Garcia

For convenience, here is the second example of the documentation:

r = 13;
n = 16:365;
lx = length(n);
x = sin(2*pi*n/153) + cos(2*pi*n/127);
plot(0:lx-1,x,'o')
hold on
y = decimate(x,r,5);
stem(lx-1:-r:0,fliplr(y),'ro','filled','markersize',4)
legend('Original','Decimated','Location','south')
xlabel('Sample number')
ylabel('Signal')

And the associated graph, showing correspondance between last element of x and last element of y:

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

请先登录，再进行评论。

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

Answer 1

Voss 2022-11-29

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1866583-what-indices-of-the-original-sequence-are-retained-in-the-output-of-the-decimate-function#answer_1115428

From item #3 in the Algorithms section of the documentation for decimate:

"In the resampled sequence (y), y(end) matches x(end) when the IIR filter is used, and y(1) matches x(1) when the FIR filter is used."

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

Voss 2022-11-29

By default, decimate uses a lowpass Chebyshev Type I infinite impulse response (IIR) filter of order 8.

Gaëtan Garcia 2022-11-29

Okay, now I get it. I did not know that the Chebyshev filter was an IIR filter, so I thought the information in the documentation was incomplete. Now to be honest, I should have realized that a filter was either FIR or IIR...

Now I really have all I need. Thank you again for being patient enough to connect the dots...

Regards.

G. Garcia

请先登录，再进行评论。