interpolate text file data onto different time axis

3 次查看(过去 30 天)
Hello,
I have a text file which when imported creates a 226 x 2 matrix, with time being the (:,1) and SSH being (:,2). I want to compare and plot this data on some data I have created which is on a 723 x 1 axis. I have tried interp1 and I keep getting all NaN's apart from a few numbers here:
interp1(data(:,1), data(:,2), createddata) so like interp1(226x1, 226x1, 723x1)
I also tried making my own xaxis (1:723) and trying that but it didn't work...
Can someone help me?
thanks in advance! Michael
  7 个评论
Michael Haderlein
Michael Haderlein 2014-7-28
Your createddata vector is 1:723? As all numbers in the first column (=x) are outside the interval 1:723, it's clear that you only get NaNs. You need to provide times at your createddata vector.
Michael
Michael 2014-7-28
I have just tried using the time values (given in a separate file) from 1993-2012 in a 723x1 matrix but I still get NaN's. Is this what you meant I should do? maybe I have misunderstood you.
t file = 723 x 1 (1.9930:1.206)

请先登录,再进行评论。

采纳的回答

dpb
dpb 2014-7-28
As the other Michael says, "Houston, we have a problem!"
How do you want these data to line up with each other? How does the time relationship of the first point in the above data set relate to the first point in your series of length 723? Likewise for the last. If you're trying to stretch the 226 over 723 points just on a 1:1 basis between first and last, you don't need to use interp1, just set the time for the 226 points via
t1=linspace(1,723,226);
or, generically,
t1=linspace(1,length(series2),length(series1));
If there's some other way they're to correlate, need to define what that relationship is and if the time in these data is significant create a time vector between those two points similarly except use the time values instead of 1:N

更多回答(2 个)

Star Strider
Star Strider 2014-7-28
If you’re getting NaN values from interp1, you are probably extrapolating.
Change your call to interp1 to:
ynew = interp1(data(:,1), data(:,2), createddata, 'linear', 'extrap');
and see if that does what you want. (Change ‘ynew’ to your own variable. I used it here simply to complete the line correctly.)

Michael
Michael 2014-7-28
Thanks, linspace worked a treat!

类别

Help CenterFile Exchange 中查找有关 Annotations 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by