How can I add missing hours (rows) in a time series matrix?

Question

Katerina F 2014-8-15

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/151004-how-can-i-add-missing-hours-rows-in-a-time-series-matrix

评论： Star Strider 2014-8-15

example time series.txt

An example from an excel file with month-day-hour-solar irradiance is attached as a text file with data for two days. I have to work with many files like this with data for each year and I do not know where the missing values are. What I would like to do is to identify the missing hours and insert the missing rows with the correct month-day-hour, and the indication NaN or something else in the missing data.

Could you please let me know which code I could use to solve this?

Thanks, K.

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

请先登录，再进行评论。

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

Answer 1

Star Strider 2014-8-15

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/151004-how-can-i-add-missing-hours-rows-in-a-time-series-matrix#answer_148756

No file attached.

For each file, I would use the beginning and end year-month-day-hour values for each day, converted to date numbers with the datenum function, and create a continuous vector of hours between them.

Then use those vectors with the interp1 function to interpolate the missing data.

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

Evan 2014-8-15

These might help:

Both are available toolbox free, I believe. There also is the interp function, but that does require a toolbox.

Star Strider 2014-8-15

在 MATLAB Online 中打开

Actually, you don’t need datenum with your data format. (It is a built-in core MATLAB function, as is interp1.)

Does this do what you want?

D = dlmread('example time series.txt');             % Read Data
DIdx = find(D(:,2) == 1, 1, 'last');                % Find Where Day#1 Ends
Hq = [0:23]';                                       % Create 24-Hr Interpolation Vector
% Interpolate:
Irad(1:24,:)  = interp1(D(1:DIdx,3), D(1:DIdx,4), Hq, 'linear','extrap');
Irad(25:48,:) = interp1(D(DIdx+1:end,3), D(DIdx+1:end,4), Hq, 'linear','extrap');
% Create Output Matrix (Match Input Format):
TimeSeries = [ones(48,1)  [ones(24,1)*1; ones(24,1)*2]  [Hq; Hq]  Irad];

The TimeSeries matrix matches the input data in case you need all that. If the first column is not always a series of 1, and you need it, we will have to work on a way to reproduce it correctly. If you don’t need it, discard it. The other columns create the day number 1 or 2, the hour vector for each day [0:23], and the interpolated or extrapolated irradiance values.

请先登录，再进行评论。