问题描述

如何在逐日流量数据中提取洪峰流量数据，记录出现洪峰流量的日期？
洪峰事件：两侧连续两日流量低于该日流量即为洪峰事件。

函数用法

我的软件版本：matlab R2016b
findpeaks
函数要求的最低版本：R2007b
findpeaks官方帮助文档

用法
pks = findpeaks(data)
[pks,locs] = findpeaks(data)
[pks,locs,w,p] = findpeaks(data)
[] = findpeaks(data,x)
[] = findpeaks(data,Fs)
[] = findpeaks(___,Name,Value)
findpeaks(___)

（1）pks = findpeaks(data)
查找data的局部最大值，data一般一维数组。pks为局部最大值组成的数组。
局部最大值的定义为，相邻的两侧数值比该数值更小。

（2）[pks,locs] = findpeaks(data)
用法同（1），locs为pks在data中的索引。

（3）[pks,locs,w,p] = findpeaks(data)
用法同（2），w为每个峰的宽度，p为每个峰的高度。

（4）其他用法待后续详细挖掘

计算案例

计算思路：利用findpeaks查找局部最大值序列（极值），然后对极值序列进行筛选，去除不满足的极值。

流量数据（30个数据）
data = [0.7 1.7 1.7 2.0 1.3 1.5 2.0 19.0 30.0 42.0 29.0 35.0 32.0 34.0 31.0 21.0 35.0 24.0 21.0 23.0 18.0 15.0 12.0 10.4 6.9 6.7 5.2 4.6 4.0 4.6];

% pks为洪峰流量，p_ind为洪峰流量出现的索引
[pks, p_ind] = findpeaks(data);
% 洪峰事件的次数
m = length(p_ind);
% 预定义新的洪峰流量序列、洪峰事件的索引
pks_2 = [];
p_ind_2 = [];
% 按照本文的定义，洪峰事件：两侧连续两日流量低于该日流量即为洪峰事件。
% 对原洪峰流量序列进行筛选，满足前后连续两日流量低于该日流量的数据为新的洪峰流量
% 从第三个洪峰流量数据开始检查，直到倒数第三个洪峰流量数据
for i = 3:m-2
if pks(i) > data(p_ind(i)-1) && pks(i) > data(p_ind(i)-2) &&…
pks(i) > data(p_ind(i)+1) && pks(i) > data(p_ind(i)+2)
pks_2 = [pks_2, pks(i)];
p_ind_2 = [p_ind_2, p_ind(i)];
end
end
% 组成新的洪峰流量序列
pks_2 = [pks(1:2), pks_2, pks(end-1:end)];
p_ind_2 = [p_ind(1:2), p_ind_2, p_ind(end-1:end)];

% tips 有点小问题是两侧边界的四个洪峰流量数据未经检测，需要后续人工检验。