增强计算斜率符号变化 (SSC) 的代码

Question

我正在尝试编写这个坡度符号变化 (SSC) 公式：

我在下面的尝试似乎遗漏了一些东西，因为它生成-2 而不是1。我真的陷入了困境。

signal = np.sin(np.linspace(1, 10)) + np.random.normal(scale=0.1, size=len(np.linspace(1, 10)))
ssc = np.diff(np.sign(np.diff(signal)))

非常感谢您的帮助。

Answer 1

• 您可以像以前一样使用np.diff 计算时间差异：

  • x_{n} - x_{n-1} 和 np.diff(x, prepend=1)[1:-1]

    或切片：x[1:-1] - x[:-2].

  • x_{n} - x_{n+1} 与 -np.diff(x)[1:]

    或切片：x[1:-1] - x[2:]

• 您可以使用给定阈值delta 的掩码实现f：

  f = lambda x: (x >= delta).astype(float)
  

SSC 的结果表达式为：

SSC = sum(f(-np.diff(x, prepend=1)[1:-1]*np.diff(x)[1:]))