ValueError ("분모 다항식은 랭크 -1 배열이어야합니다.")

Question

lti transient response analysis using Python(numpy, scipy, matplotlib)에 다음 코드가 있습니다. 나는 파이썬으로 새롭다. 나는 플롯해야하는 전달 행렬을 가지고있다.

나는 mathwork: tf을 발견했습니다.

from numpy import min, max 
from scipy import linspace 
from scipy.signal import lti, step, impulse 

num00 = [0.0] 
den00 = [0.0] 

num01 = [-2383.3] 
den01 = [1.0,160.3460,-1962.0,-314598.852] 

num10 = [1.0] 
den10 = [1.0] 

num11 = [31.9361,0,111320.0] 
den11 = [1.0,160.3460,-1962.0,-314598.852] 

num = [[num00,num01],[num10,num11]] 
den = [[den00,den01],[den10,den11]] 

tf = lti(num,den) 

t = 0  
s = 0 

# get t = time, s = unit-step response 
t , s = step(tf) 

t , s = step(tf, T = linspace(min(t), t[-1], 1000)) 

t , i = impulse(tf, T = linspace(min(t), t[-1], 1000)) 

from matplotlib import pyplot as plt 

plt.plot(t, s, t, i) 

plt.title('Transient-Response Analysis') 
plt.xlabel('Time(sec)') 
plt.ylabel('Amplitude') 
plt.hlines(1, min(t), max(t), colors='r') 
plt.hlines(0, min(t), max(t)) 
plt.xlim(xmax=max(t)) 
plt.legend(('Unit-Step Response', 'Unit-Impulse Response'), loc=0) 
plt.grid() 
plt.show() 

을 나는 오류가 다음 무엇입니까 : 문제의

>>> tf = lti(num,den) 
Traceback (most recent call last): 
    File "<stdin>", line 1, in <module> 
    File "C:\Python26\lib\site-packages\scipy\signal\ltisys.py", line 236, in __init__self.__dict__['num'], self.__dict__['den'] = normalize(*args) 
    File "C:\Python26\lib\site-packages\scipy\signal\filter_design.py", line 276, in normalize raise ValueError("Denominator polynomial must be rank-1 array.") ValueError: Denominator polynomial must be rank-1 array. 

Answer 1

일부는 당신이 전달하는 NUM/서재가 잘 형성 매트릭스되지 않는 것입니다 다음과 같이 내가 노력하고 있어요. 당신의 코드에서, 당신은 : 지금까지 NumPy와 당신이 매트릭스를 만들려고에 관한 한, 나는 단지 1, 그것은 전달 함수 행렬의 하나의 구성 요소의 실현 있기 때문에 아주 잘 작동하지 않습니다

num01 = [-2383.3] 
den01 = [1.0,160.3460,-1962.0,-314598.852] 

요소는 분자에, 요소는 분모에 4 개 있습니다. 따라서 다음과 같은 것이 필요할 것입니다 :

num01 = [ 0,  0,  0,-2383.3] 

매우 높은 수준의 분자가 있어야합니다. 내가 걸음을 밟으 려 할 때 : 아마도 당신이 기대하는 것은 아니 겠지요. python-control 패키지를 살펴 보는 것도 좋습니다. 물론 SLICOT python package과 같은 패키지의 모든 기본 설정을 가져와야합니다. 나는 그것이 궁극적으로 당신을 잘 섬길 것이라고 믿습니다.