Resampling Data using Interp1
41 views (last 30 days)
Show older comments
Lucianne Morin
on 29 Jun 2020
Commented: Star Strider
on 1 Jul 2020
Hello!
This may be a fairly obvious answer, but I'm trying to resample some data down to 100 Hz using interp1. I load my data, and then extract the first column (which is the time series) to use in my for loop. However, I get the error "Interpolation requires at least two sample points in each dimension." I tried adjusting the code to call for a(i+1) to see if maybe I needed to start later, but I got the same error. The data itself is part of my thesis so I'm not sure if I can attach it, but I've copy and pasted my code below.
Thank you in advance for any help!
load('BNI_processed.mat')
a = Force_arr;
time_matrix =(a(1:end,1)).';
for i = 1:length(time_matrix)-2
aa = length(time_matrix(i):time_matrix(i+1));
p = 61;
Force(:,i)=interp1(a(i,2),1:aa/p:aa, 'spline');
end
1 Comment
Accepted Answer
Star Strider
on 30 Jun 2020
You are doing signal processing. For that, use the resample function. It will do the interpolation, and will also use an anti-aliasing filter to remove unwanted artifacts. I am not certain what your signal is or what you want to do with it, so I cannot provide a link to the most appropriate section of the resample documentation for that.
2 Comments
Star Strider
on 1 Jul 2020
I was hoping for ‘BNI_processed.mat’. Lacking it, the resample function can do exactly what you wantpreviously, since I have done this previously.
From the documentation (there is no specific link to this line):
- y = resample(x,tx,fs) uses a polyphase antialiasing filter to resample the signal at the uniform sample rate specified in fs.
Referring to your code, this would probably be:
[force_array, p_new] = resample(a(:,2:end), p, 100);
[shoe_array, q_new] = resample(b(:,2:end), q, 100);
I cannot determine for certain that will work because I do not have ‘BNI_processed.mat’ to work with and check. I am going only by your interp1 calls.
Without ‘BNI_processed.mat’, my only option is to post this as UNTESTED CODE. It should do what you want.
.
More Answers (1)
Walter Roberson
on 1 Jul 2020
Edited: Walter Roberson
on 1 Jul 2020
load('BNI_processed.mat')
input_forces = Force_arr(:,2);
time_matrix = Force_arr(:,1);
Fs = 100;
last_time = ceil(time_matrix(end) / Fs) * Fs; %round UP
time_to_interp = 0 : 1/Fs : last_time;
Force = [time_to_interp; interp1(time_matrix, input_forces, time_to_interp)] .;
0 Comments
See Also
Categories
Find more on Multirate Signal Processing in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!