python波形分割_用Python裁剪波形顶部

最简单的方法是创建一个引入“软剪裁”的映射，值是1:1线性的，直到接近阈值，然后压缩，直到达到某个更大的阈值。在

编辑：这个想法是合理的，但是实现还有一点不尽如人意。我又想了想，换了个更好的。在

目标是使从压缩开始到硬剪裁限制的平滑过渡。线性变换的斜率为1，截断信号的斜率为0。从斜率1转换为斜率0的最干净的方法是使用四分之一的正弦波，因此更新后的函数基于sin。在from math import *

hard_limit = 32767

linear_limit = 23197 # -3 dB

clip_limit = linear_limit + int(pi/2 * (hard_limit - linear_limit))

def soft_clip(n):

amplitude, sign = abs(n), 1 if n >= 0 else -1

if amplitude <= linear_limit:

return n

if amplitude >= clip_limit:

return hard_limit * sign

scale = hard_limit - linear_limit

compression = scale * sin(float(amplitude - linear_limit) / scale)

return (linear_limit + int(compression)) * sign

下面是一个在振幅为40000的正弦波上使用上述方法的示例，背景是未夹持的版本。在