清空筛选条件
清空筛选条件

how to find time interval of peak value in a timeseries ?

3 次查看（过去 30 天）
Harish Dhanasekaran Velayutha Rajan
评论： Star Strider 2017-8-4
i have a load profile data set.. Time vs Load Demand as a timeseries object. If I have a threshold for peak as 800KW How do I find at what time my value goes above 800 kW ? The time series object is for every 30 mins interval over 24 hrs.
I need to find the highlighted time interval in the picture. Any idea?

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

采纳的回答

Star Strider
Star Strider 2017-8-4
If you want to interpolate to find the ‘exact’ time, this works:
t = 0 : 30 : 1440;
LD = 500 + 400*sin(2*pi*t/2880);
zci = @(v) find(v(:).*circshift(v(:), [-1 0]) <= 0); % Returns Approximate Zero-Crossing Indices Of Argument Vector
X800_idx = zci(LD - 800); % Subtract 800, Find Approximate Zero-Crossing
LD800_Rng = LD(X800_idx-1 : X800_idx+1)'; % Define Range For Interpolation
t800_Rng = t(X800_idx-1 : X800_idx+1)'; % Define Range For Interpolation
t800 = interp1(LD800_Rng, t800_Rng, 800, 'linear', 'extrap'); % Find ‘Exact’ Time (Minutes)
figure(1)
plot(t, LD)
hold on
plot(t800, 800, '*r')
hold off
grid
  1 个评论
Star Strider
Star Strider 2017-8-4
A slightly more robust version that finds the points where ‘LD’ crosses 800:
t = 0 : 30 : 1440; % Create Data
LD = 500 + 400*sin(2*pi*t/2880); % Create Data
zci = @(v) find(v(:).*circshift(v(:), [-1 0]) <= 0); % Returns Approximate Zero-Crossing Indices Of Argument Vector
X800_idx = zci(LD - 800); % Subtract 800, Find Approximate Zero-Crossing
for k1 = 1:numel(X800_idx)
LD800_Rng = LD(X800_idx(k1)-1 : X800_idx(k1)+1)'; % Define Range For Interpolation
t800_Rng = t(X800_idx(k1)-1 : X800_idx(k1)+1)'; % Define Range For Interpolation
t800(k1) = interp1(LD800_Rng, t800_Rng, 800, 'linear', 'extrap'); % Find ‘Exact’ Time (Minutes)
end
figure(1)
plot(t, LD) % Plot Load Curve
hold on
plot(t800, [1 1]*800, '*r') % Plot ‘Load = 800’ Points
hold off
grid

请先登录，再进行评论。

更多回答（1 个）

fbaillon
fbaillon 2017-8-4
You can use the find function:
T= ... % Time
LD= .... % Load Demand
timeYouWant=T(find(LD>800e3,1,'first))
Something like that...

类别

在 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