关于提取峰谷值的程序

[xiaokang5280]

void Rainflow::find_peaks(vector<float> * peaks) {
peaks->push_back(_points[0]);


unsigned i;
for (i = 1; i < (_points.size()-1); i++) {
double slope1 = _points[i] - _points[i - 1];
double slope2 = _points[i + 1] - _points[i];

if ((slope1 * slope2) <= 0. && fabs(slope1) > 0.) {
peaks->push_back(_points[i]);
}
}
peaks->push_back(_points.back());
}
以上是一段提取峰谷值的函数，峰谷值就是满足相邻的点不相等，每三个点中的第二个点小于旁边两个点或者大于旁边两个点。
我这段程序提取出来的点多了，有大神帮我看看吗

[Isnis-fallen]

if ((slope1 * slope2) <= 0. && fabs(slope1) > 0.) 改为 if ((slope1 * slope2) < 0）

[昆仑道长]

将你提取出来的波峰数据点从大到小排序，取排在前面的几个数据就是你想要的.

[昆仑道长]

我还想到一个绝妙的算法，准备鼓捣一下，发表篇论文.

[昆仑道长]

刚百度了下，前人已经做的很好了，不需要我重复造轮子了.参见 https://blog.csdn.net/boris_wang2008/article/details/81209849 https://blog.csdn.net/zone53/article/details/77942341

[xiaokang5280]

引用 1 楼 focuslight 的回复：

if ((slope1 * slope2) <= 0. && fabs(slope1) > 0.) 改为 if ((slope1 * slope2) < 0）

这样提取出来的数又少了