Downsampled data exceeds the input data

2022 11 18
2 个回答
回答已采纳
更新时间：2022 11 18
2 次查看（30 天）

0 个投票

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
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

1 个投票

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

0 个投票

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更早的评论 隐藏 1更早的评论

Thank you.
Both 'interp1' and 'griddedInterpolant' function expects specific data points 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.
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.
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.
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.

请先登录，再进行评论。

类别

在 帮助中心 和 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