ジュリーの安定判別シミュレーター — 離散時間系

ディジタル制御器の設計：マイコンやFPGAで動くPIDコントローラ、モータドライバ、電源の電流制御ループは、すべてサンプリング周期ごとに値を更新する離散時間系です。これらの閉ループ特性方程式は z の多項式になり、安定性はジュリーの判別で確認します。とくにサンプリング周期を長くすると極が単位円の外へ移動しやすくなるため、周期と制御ゲインを変数に残したまま、ジュリー条件を不等式の形で解いて安定なパラメータ範囲を求めます。

連続系のディジタル化（離散化）の検証：アナログで設計した制御器を、Tustin変換やゼロ次ホールドでディジタル実装に置き換える場面は非常に多いです。離散化すると安定領域が左半面から単位円内へ変わるため、連続系で安定だった設計が離散化後に不安定になることがあります。ジュリー判別は、離散化後の特性多項式が本当に単位円内に収まっているかを手早く検算するのに使えます。

ディジタルフィルタの安定性確認：IIRフィルタ（無限インパルス応答フィルタ）の伝達関数の分母多項式も、根が単位円内になければ出力が発散します。オーディオ処理や通信のディジタル信号処理では、係数を量子化した後でも分母の根が単位円内に留まっているかをジュリー条件で確認します。固定小数点実装で係数を丸めた結果、極が円の外へ出てしまう事故を防ぐ検算になります。

制御工学教育とCAE前の検算：離散時間系の極を数値的に求める前の「答えの当たりづけ」としても有用です。MATLABやPythonで根を計算する前にジュリー条件を見れば、安定か不安定かが即座に分かります。逆に数値計算の根がジュリー判別と矛盾していれば、係数の入力ミスや離散化のモデル化誤りを疑うサニティチェックになります。

まず最も多いのが、「連続系の感覚をそのまま離散系に持ち込む」誤解です。連続時間系では「極が原点に近いほど応答が速く安定」というイメージがありますが、離散時間系の z 平面では事情が違います。安定領域は単位円の内側であり、原点 z=0 が最も速い（実質1サンプルで収束する）応答に対応し、円周 |z|=1 に近いほど減衰が遅くなります。s 平面の「左に行くほど良い」という感覚で z 平面を見ると判断を誤ります。z 平面では「原点に近いほど良い」と読み替える必要があります。

次に、「全ての係数が正なら安定」という思い込みです。これはラウス・フルビッツでも誤解されますが、ジュリー判別では係数の符号は安定の十分条件にも必要条件にもなりません。離散時間系の安定な多項式の係数は負になることもあり（左の例 a₃=−0.8、a₁=−0.128 がそうです）、逆に係数が全て正でも根が単位円の外に出ることがあります。必ず必要条件 P(1)>0・P(−1)>0・|a₀|<1 とジュリー内部条件の両方まで確認してください。係数の符号だけを見て安心するのは危険です。

最後に、サンプリング周期と係数量子化の影響です。ジュリー判別は与えられた多項式に対しては厳密ですが、その多項式自体がサンプリング周期に強く依存します。周期を長くとると同じ連続系でも離散極が単位円の外へ移り、不安定化します。さらにディジタル実装では係数が有限ビットに丸められるため、設計上は安定でも実装後の量子化済み係数で根が円の外に出ることがあります。ジュリー条件は「ぎりぎり満たす」のではなく、根が単位円の内側に十分な余裕をもって入るように、係数とサンプリング周期に余裕を取った設計にしてください。