PythonのScipy - 離散化の問題

Question

これはmath.stackexchangeに適していると思ったら教えてください！

オープンソースソフトウェアをより独立させようとしている人の中で、私はPythonには初めてです。私は制御システムを扱うためにScipyパッケージを使用しています。私がPythonをデバッグするのに少し便利になるまでには時間がかかるでしょう。今、コンテキストが設定されていることを、ここで私が直面しています問題は次のとおりです。

私はストレートフォワード伝達関数を離散化し、ボード線図を作成したい：

P_d = sp.signal.cont2discrete(([1],[32,80,80,40,10,1]),0.1,method='zoh'); 
w_d, mag_d, phase_d = sp.signal.bode(P_d); 

私が手にエラーが

語ります

location "location/Controller.py" 
Traceback (most recent call last): 

File "location/Controller.py", line 17, in <module> 
w_d, mag_d, phase_d     = sp.signal.bode(P_d); 

File "location\Anaconda3\lib\site-packages\scipy\signal\ltisys.py", line 2398, in bode 
w, y = freqresp(system, w=w, n=n) 

File "location\Anaconda3\lib\site-packages\scipy\signal\ltisys.py", line 2472, in freqresp 
raise ValueError("freqresp() requires a SISO (single input, single " 

ValueError: freqresp() requires a SISO (single input, single output) system. 

実際の場所を単語自体に置き換えました。私は離散化プロセスを理解して何から

はをSISOないにSISOシステムを変換しているようです。どうして？ Pythonで離散時間システムのボード線図を作るための別のコマンドがありますか？

必要に応じてお気軽にご相談ください。

Answer 1

scipy.signalの一部の機能は、連続システムおよび離散システムで動作するものがあります。の離散関数は、連続関数と同じ名前をもち、接頭辞はdです（概要はhereを参照）。 これらの関数は、システムをタプルとして表し、以下のように解釈します。

離散系：

1. （dltiのインスタンス）
2. （NUM、DEN、DT）
3. （零点、極、ゲイン、DT）
4. （A、B、 C、D、DT）

連続系：

1. （LTIのインスタンス）
2. （NUM、DEN）
3. （零点、極、ゲイン）
4. （A、B、C、D）

cont2discreteは、ケース2（num、den）として([1], [32,80,80,40,10,1])と解釈します。形式（num、den、dt）の一致する離散システムを返します。このシステムを継続的なシステムを必要とするbodeに渡すと、予測できない結果が（ゼロ、ポール、ゲイン）として解釈されます。

使用する正しい機能は、discrete（num、den、dt）システムで正しく動作するdbodeです。