Downsampled data exceeds the input data

2 次查看（过去 30 天）
Chandramouli Jambunathan
评论： Chandramouli Jambunathan 2022-11-18
I have an input data sampled at 64 Hz and I want it converted to 60Hz. I tried using 'resample' function and 'downsample' function within matlab,but in both the cases the output data exceeds the input data at some sample points. For example, the input data contains angles in degrees which is under -180 deg to +180 deg, but the downsampled data contains 220 deg,which is wrong. how to avoid these?
In the attached picture, you can find all the input and output are overlayed. Can someone please help me?
  2 个评论
Mathieu NOE
Mathieu NOE 2022-11-18
hello
for the resampling task , use interp1 with a new time vector defined at rate = 60 Hz
to avoid data above / below +/- 180 ° use wrap
Chandramouli Jambunathan
Chandramouli Jambunathan 2022-11-18
Thank you.. interp1 expects the coordinates of the query points (i.e. xq) and I dont have specific co-ordinates or sample points for my 60Hz data output. All I wanted is to convert an array of data which is at 64Hz to 60Hz.
Now I am using 'interp' and 'decimate' functions to upscale and downsample respectively. And the Wrap function seems to help me to limit the angles.
Thank you.
But is it a bug in the 'interp' function as it gives output outside the input data range(as you can see in the graph).Ideally it should interpolate within the input data range.

请先登录，再进行评论。

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

采纳的回答

Mathieu NOE
Mathieu NOE 2022-11-18
a demo is maybe more efficient here
so the original data gets resampled at 60 Hz by using only the start / end points of the time vector of the original data (sampled at 64 Hz)
the two data sets do match visually on the graph even though the points are not on the same grid (of course)
samples=1000;
% original data at 64 Hz rate
dt1 = 1/64;
t1= 1 + (0:samples-1)*dt1;
y1 = 180*(square(t1));
plot(t1,y1,'b-*')
% resample the data at 60Hz
dt2 = 1/60;
t2= t1(1):dt2:t1(end); % NB wee use only the first and last value of t1 to create t2
y2 = interp1(t1,y1,t2);
plot(t1,y1,'b',t2,y2,'r')

更多回答（1 个）

Jan
Jan 2022-11-18
A linear interpolation avoids output points outside the range of input point. Use interp1 or the modern and faster griddedInterpolant .
A fast C-mex function for this job: https://www.mathworks.com/matlabcentral/fileexchange/25463-scaletime
  3 个评论
显示 1更早的评论
Jan
Jan 2022-11-18
I do not understand the meaning of "Both 'interp1' and 'griddedInterpolant' function expects specific data points". But Mathieu has show already, how to use these interpolation methods.
Chandramouli Jambunathan
Chandramouli Jambunathan 2022-11-18
I misunderstood that I need to supply the sample points for 60Hz.. that's why I mentioned earlier to your post that those function expects specific data points.. but I was wrong.. Mathieu has already shown with an example and that worked. Thanks for your comments too.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Multirate Signal Processing 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by