查找数组的峰值并创建位序列（不使用 findpeaks）答案

【问题标题】：Finding peaks of an array and creating a sequence of bits (without using findpeaks)查找数组的峰值并创建位序列（不使用 findpeaks）
【发布时间】：2018-03-24 14:56:56
【问题描述】：

给定一个数组 V，它是一个微分方程组的数值解，我想确定这个数组的局部正峰值。每次我有一个值大于 0 的峰值时，我都会将其编码为 1，每次我有一个负值的峰值时，我都会将其编码为 0（我不能使用 findpeaks）。

我尝试了以下方法：

function Suite_bits = message_digital(V)
j = 1;
bit(1) = 0;
for i=2:numel(V)-1
    if V(i-1)<V(i) && V(i+1)<V(i)
        if V(i)>0
            bit(j) = 1;
            j = j + 1; 
        else
            bit(j) = 0;
            j = j + 1;
        end 
    end
    Suite_bits = bit;
end

但是由于向量 V 中的一些小异常，步骤之间的值波动（记住 V 是步长 h = 0.1 的数值求解 ode 系统的解），我得到的峰值比图中实际出现的峰值要多.有什么想法可以解决这个问题吗？

【问题讨论】：

您可能会发现以下内容很有用。 stackoverflow.com/questions/22583391/…

标签： arrays algorithm matlab signal-processing

【解决方案1】：

使用 MATLAB 提供的“平滑”函数可以避免 V 向量的波动。默认情况下，此函数使用 5 点移动平均线平滑数据。您可以更改此函数的参数，使移动平均线包含多于或少于 5 个点。因此，根据波动的幅度，可以安排这个参数，以便算法平滑和忽略波动。这是一个可以工作的代码示例：

function Suite_bits = message_digital(V)
j = 1;
V = smooth(V, 15);
bit(1) = 0;
for i=2:numel(V)-1
    if V(i)>V(i-1) && V(i)>V(i+1) 
        if V(i)>0
            bit(j) = 1;
            j = j + 1; 
        else
            bit(j) = 0;
            j = j + 1;
        end 
    end
    Suite_bits = bit;
end
end

【讨论】：