趣味算法-求波峰波谷最大值

给一个数组,求它的最大的波峰波谷的落差。

举例:数组  A={2, 3, 6, 5, 7, 9}, 其中 6 和 9 被看做是波峰,2和5则是波谷。D[2, 6]=4, D[6,5]=1, D=[5,9]=4. 则 Thus, MaxD(A)=4.


想法:波峰不一定是数组中的最大值,他是趋势向下的转折点,波谷也不一定是数组中的最小值,他是趋势向上的转折点。数组中最边上的两个元素比较特殊,他们必然属于转折的部分。如果A[0] < A[1]  则A[0]是波谷,反之则是波峰。

先扫描一次数组,找到其中所有的转折点,然后配对计算出他们的差,取差最大的。


测试数据:

2 3 6 5 7 9
2 3 6 5 7
2 3 6 5 7 9 10


示例程序:

#include <iostream>

using namespace std;



void scanPeak(int arr[], int length, int* pvlst)
{
    int i = 0;
    int curr_peak = 0;
    int curr_val = 0;
    bool is_continue = true;

    for (i = 1; i < length-1; i++)
    {
        if ((arr[i] < arr[i-1]) && (arr[i] < arr[i+1]))
        {
            pvlst[i] = -1;
  • 1
    点赞
  • 17
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
这段代码可能需要优化的地方包括: 1. 加载数据时可能需要检查数据是否存在异常值,以避免影响后续处理的精度。 2. 对于寻找波峰波谷的方法,可能需要使用更加精确的算法来寻找局部最大值和最小值。 下面是一种可能的优化方式: ``` % 加载数据并且检查数据 data1 = load('out3.txt'); if any(isnan(data1(:))) || any(isinf(data1(:))) error('数据存在异常值'); end data = data1(:, 1); % 对数据进行平滑处理 window_size = 5; data_smooth = smoothdata(data, 'movmean', window_size); % 对平滑后的数据进行插值 x = 1:length(data_smooth); xi = linspace(1, length(data_smooth), 10*length(data_smooth)); % 使用样条插值插值数据 data_interp = interp1(x, data_smooth, xi, 'spline'); % 寻找波峰波谷 [peaks, IndMax] = findpeaks(data_interp); [troughs, IndMin] = findpeaks(-data_interp); troughs = -troughs; % 绘制图形 figure; hold on; box on; plot(xi, data_interp); plot(xi(IndMin), data_interp(IndMin), 'r^'); plot(xi(IndMax), data_interp(IndMax), 'k*'); legend('曲线', '波谷点', '波峰点'); title('计算离散节点的波峰波谷信息', 'FontWeight', 'Bold'); ``` 其中,我们添加了对数据的检查代码,以避免加载到异常值。同时,我们使用 `findpeaks` 函数寻找波峰波谷,这个函数可以返回波峰波谷的位置和值。最后,我们还添加了对波峰波谷的可视化代码,以便观察波峰波谷的位置。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值