scipy.signal.besselap

scipy.signal.besselap(N, norm='phase')[源代码]

返回N阶贝塞尔过滤模拟原型的(z,p,k)。

参数
N集成

过滤勋章。

norm{‘阶段’,‘延迟’,‘料盒’},可选

频率归一化:

phase

过滤被归一化,使得相位响应在角度(例如,rad/s)截止频率为1时到达其中点。低通和高通滤波器都会发生这种情况,因此这是“相位匹配”的情况。 [6]

震级响应渐近线与同阶的巴特沃斯过滤相同,截止值为 Wn

这是默认设置,与MATLAB的实现相匹配。

delay

过滤被规格化,使得通带中的群延迟为1(例如,1秒)。这是通过求解贝塞尔多项式得到的“自然”型。

mag

过滤经过归一化处理,在角频率1处的增益幅度为-3dB。邦德称之为“频率归一化”。 [1]

0.18.0 新版功能.

退货
zndarray

传递函数的零点。始终为空数组。

pndarray

传递函数的极点。

k标量

传递函数的增益。对于相位规格化,该值始终为1。

参见

bessel

使用此原型设计过滤功能

注意事项

为了找到极点位置,需要生成近似的起始点 [2] 关于普通贝塞尔多项式的零点 [3], 然后是Aberth-Ehrlich方法 [4] [5] 在Kv(X)Bessel函数上使用,以计算更精确的零,然后将这些位置围绕单位圆进行反转。

参考文献

1

C.R.邦德,“贝塞尔·过滤常数”,http://www.crbond.com/papers/bsf.pdf

2

Campos和Calderon,“贝塞尔多项式零点的近似闭式公式”, arXiv:1105.0957

3

Thomson,W.E.,“具有最大平坦频率特性的延迟网络”,电气工程师学会论文集,第三部分,1949年11月,第96卷,第44期,第487-490页。

4

ABERTH,“同时寻找多项式的所有零点的迭代方法”,“计算数学”,第27卷,第122期,1973年4月。

5

Ehrlich,“一种修正的牛顿多项式方法”,“ACM通讯”,第10卷,第2期,第107-108页,1967年2月, DOI:10.1145/363067.363115

6

米勒和博恩,“贝塞尔过滤交叉,及其与他人的关系”,Rane Note 147,1998年,https://www.ranecommercial.com/legacy/note147.html