有人能解释一下为什么这段代码在Python和Matlab中给出了不同的结果吗?
Matlab软件
z = [0.97, 0.85, 0];
p = [0.97, 0.94, 0];
k = 0.2;
disp(zp2sos(z,p,k))
输出量:
0.2000 -0.1700 0 1.0000 -0.9400 0
1.0000 -0.9700 0 1.0000 -0.9700 0
巨蟒
from scipy.signal import zpk2sos
z = [0.97, 0.85, 0]
p = [0.97, 0.94, 0]
k = 0.2
print(zpk2sos(z,p,k))
输出量:
[[ 0.2 0. 0. 1. 0. 0. ]
[ 1. -1.82 0.8245 1. -1.91 0.9118]]
1条答案
按热度按时间rryofs0p1#
我在github上找到了这段代码,这段代码定义了一个zpk2sos版本,它的输出和MATLAB的zp2sos相同。你可能需要从代码中删除一些过时的包,但除此之外,它对我来说很好用。
https://gist.github.com/endolith/4525003/8ac88d203f874fd2b97498d50da818bbf5fac0f8