scipy.signal.besselap¶
- scipy.signal.besselap(N, norm='phase')[源代码]¶
返回N阶贝塞尔过滤模拟原型的(z,p,k)。
- 参数
- 退货
- 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