最简单的方法是创建一个引入“软剪裁”的映射,值是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的正弦波上使用上述方法的示例,背景是未夹持的版本。在