むだ時間のパデ近似シミュレーター 戻る
制御工学

むだ時間のパデ近似シミュレーター

純粋なむだ時間 e^(−sT) を有理伝達関数で近似する「パデ近似」を体験するツールです。近似次数やむだ時間を変えると、位相特性・有効帯域、そして右半面零点が生み出すステップ応答の初期アンダーシュートがリアルタイムで分かります。

パラメータ設定
むだ時間 T
s
入力が出力に現れるまでの純粋な時間遅れ
パデ近似の次数 n
[n/n] 近似。大きいほど正確だが右半面零点が増える
プラント時定数 τ
s
むだ時間とカスケード接続する1次プラントの時定数
プラントゲイン K
1次プラントの定常ゲイン(ステップ応答の最終値)
計算結果
近似次数 n
右半面零点の数
位相誤差 @ω=1/T (deg)
有効帯域 誤差<5° (rad/s)
ステップ応答の二乗平均誤差
初期アンダーシュート
ステップ応答アニメーション — 真のむだ時間 vs パデ近似

青の破線が真の遅れ応答(t=T まで平坦→指数的に立ち上がり)、オレンジがパデ近似の応答。立ち上がり前の「へこみ」が右半面零点による初期アンダーシュートです。

位相特性 — 真のむだ時間とパデ近似の位相
ステップ応答 — 真の遅れ応答とパデ近似応答
理論・主要公式

$$e^{-sT}\approx\frac{N_n(s)}{D_n(s)},\qquad e^{-sT}\approx\frac{1-sT/2}{1+sT/2}\ \ (n=1)$$

むだ時間 e^(−sT) の [n/n] パデ近似。分子 Nₙ(s) と分母 Dₙ(s) はともに n 次多項式。n=1 では 1次/1次の有理関数になる。

$$e^{-sT}\approx\frac{1-sT/2+(sT)^2/12}{1+sT/2+(sT)^2/12}\ \ (n=2)$$

2次パデ近似。係数は c_k = (2n−k)!·n! / ((2n)!·k!·(n−k)!) で与えられ、Dₙ(s)=Σ c_k(sT)^k、Nₙ(s)=Σ c_k(−sT)^k。

$$\left|\frac{N_n(j\omega)}{D_n(j\omega)}\right|=1,\qquad \angle e^{-j\omega T}=-\omega T$$

パデ近似はオールパス(全域でゲイン1)であり、n 個の右半面零点を寄与する。誤差はすべて位相に現れ、周波数が高いほど真値 −ωT から離れる。

むだ時間のパデ近似とは

🙋
「むだ時間」って、入力してから出力が出るまでの遅れのことですよね。これって伝達関数で書けないんですか?
🎓
書けるよ。むだ時間 T の伝達関数はちゃんと e^(−sT) という形になる。問題は、これが「指数関数」だってことなんだ。例えば配管を流れる水の温度を測るとき、ヒーターから温度計まで水が届くのに数秒かかる。その遅れがまさに e^(−sT)。でもこれは分子も分母も多項式じゃない「無理関数」だから、極や零点が無限個あって、根軌跡もボード線図設計も状態空間モデルもそのままじゃ使えないんだよ。
🙋
なるほど…じゃあパデ近似は、その e^(−sT) を無理やり多項式の比に直すってことですか?
🎓
そのとおり。パデ近似は e^(−sT) を「分子 Nₙ(s) ÷ 分母 Dₙ(s)」という有理関数で近似する。両方とも n 次多項式で、n=1 なら一番有名な (1−sT/2)/(1+sT/2)。左のスライダーで「次数 n」を1にすると、まさにこの式が分母に現れる。これなら極も零点も有限個になるから、普通のコントローラ設計のなかにむだ時間を組み込めるんだ。
🙋
上のステップ応答グラフを見ると、パデ近似のほうが立ち上がる前に一瞬下に「へこむ」んですけど、これバグじゃないんですか?真のむだ時間はただ平らなのに。
🎓
いいところに気づいたね。それはバグじゃなくて、パデ近似の本質的な性質なんだ。[n/n] パデ近似は分子に n 個の「右半面零点」を持つ。右半面零点を持つ系は「非最小位相系」と呼ばれて、ステップを加えた瞬間に出力が目標と逆へ動く「逆応答」を必ず示す。本物のむだ時間は t
🙋
じゃあ次数 n を上げれば、その「へこみ」は消えて本物に近づくんですか?
🎓
そう、n を上げるとアンダーシュートは小さく・短くなって、ステップ応答も位相特性も真のむだ時間にぐっと近づく。左で n を1から5へ動かすと、上の「有効帯域」の数字(位相誤差が5°を超える周波数)が単調に広がるのが見えるはずだ。ただし n を上げた分だけ右半面零点も増えて系が高次になる。だから実務では、必要な制御帯域でむだ時間の位相が十分正確になる「最小の次数」を選ぶ。PID設計なら1次か2次で足りることが多いよ。
🙋
位相のグラフだけ大きくずれて、振幅は出てこないのはなぜですか?
🎓
むだ時間は「振幅は変えず、位相だけ遅らせる」要素だからだよ。e^(−jωT) は大きさが常に1で、角度だけ −ωT 回る。パデ近似もこの性質をきちんと受け継いでいて、|Nₙ/Dₙ|=1 が厳密に成り立つ「オールパス(全域通過)」になっている。だから近似の誤差は振幅には一切出ず、すべて位相の側に現れる。位相のグラフだけ見ればいい、というのはそういう理由なんだ。

よくある質問

純粋なむだ時間 e^(−sT) は s の指数関数で、分母・分子が多項式である有理伝達関数ではありません。根や極を持たない無理関数なので、根軌跡法・ボード線図設計・状態空間モデルといった「有理伝達関数を前提とした」制御設計手法をそのまま適用できません。パデ近似は e^(−sT) を分子 Nₙ(s) と分母 Dₙ(s) の比という有理関数で近似します。これにより極・零点が定義でき、通常のコントローラ設計や数値シミュレーションにむだ時間を組み込めるようになります。
次数 n を上げると、真のむだ時間の位相 −ωT を正確に再現できる周波数帯域が広がります。本ツールの「有効帯域」は位相誤差が初めて5°を超える周波数で、n を1から5へ上げると単調に広がります。一方で [n/n] パデ近似は n 個の右半面(不安定)零点を導入するため、次数を上げるほどステップ応答の初期アンダーシュート(逆応答)の形が複雑になります。実務では、対象とする制御帯域でむだ時間の位相が十分正確になる最小の次数を選びます。多くのPID設計では1次か2次で足ります。
[n/n] パデ近似は分子に n 個の右半面零点(s 平面の右側にある零点)を持ちます。右半面零点を持つ系は「非最小位相系」と呼ばれ、ステップ入力を加えた直後に出力が目標と逆方向へ動く「逆応答(アンダーシュート)」を示します。本来のむだ時間は t
オールパスとは、すべての周波数 ω でゲイン(振幅比)が1である性質です。[n/n] パデ近似 Nₙ(jω)/Dₙ(jω) は、分子と分母が共役の関係にあるため |Nₙ/Dₙ|=1 が厳密に成り立ちます。これは真のむだ時間 e^(−jωT) が |e^(−jωT)|=1 で振幅を変えず純粋な位相遅れだけを与える性質と一致します。つまりパデ近似はむだ時間の「振幅は変えず位相だけ遅らせる」という本質を正しく捉えており、近似による誤差はすべて位相の側に現れます。

実世界での応用

プロセス制御(化学プラント・温度制御):配管の輸送遅れ、熱交換器の応答遅れ、分析計のサンプリング遅れなど、プロセス系にはむだ時間が必ず付きまといます。PIDチューニングやモデル予測制御の設計では、このむだ時間をパデ近似で有理化してから安定余裕・ゲイン交差周波数を評価します。むだ時間が大きいプラントほど位相余裕を食いつぶすため、近似精度がそのまま制御性能の見積もり精度に直結します。

サーボ系・通信遅延を含む制御:ネットワーク越しのリモート制御や、センサ・アクチュエータの計算遅れを含むディジタル制御系では、ループ内のむだ時間が発振の原因になります。設計段階でむだ時間をパデ近似に置き換え、ナイキスト線図やボード線図で「むだ時間がどれだけ位相を食うか」を定量化し、制御帯域を安全側に決めます。

スミス予測器など、むだ時間補償の理解:むだ時間補償器(スミス予測器)を学ぶには、まず「むだ時間が有理伝達関数として扱えるとどうなるか」を体感するのが近道です。パデ近似は補償器の内部モデルとしても使われ、近似次数を変えたときに閉ループ挙動がどう変わるかを把握しておくことは、補償器の設計・調整に直結します。

制御工学の教育・CAEの前段検討:むだ時間を含む系の根軌跡・周波数応答は、解析ソフトでも「まずパデ近似に置換してから」計算されることが多くあります。本ツールのように次数と誤差の関係を可視化しておくと、解析ツールが内部で何をしているか、なぜ高次にすると逆応答が複雑になるかを理解した上で結果を読めます。

よくある誤解と注意点

まず多い誤解が、「次数 n を上げれば上げるほど良い」というものです。確かに位相の正確な帯域は広がりますが、[n/n] パデ近似は n 個の右半面零点を必ず導入します。これらを含んだモデルでコントローラを設計すると、近似が生む非最小位相のクセまで補償しようとして、かえって設計が難しくなることがあります。さらにモデル次数が上がると計算負荷も増えます。「対象の制御帯域でむだ時間の位相が十分正確になる最小次数」を選ぶのが定石で、やみくもな高次化は禁物です。

次に、「パデ近似はむだ時間と完全に同じ」と思い込むこと。パデ近似はあくまで有理関数による近似であり、真のむだ時間とは特に高周波で位相が大きくずれます。本ツールの位相特性グラフでも、ある周波数を超えるとパデ近似の位相が真値 −ωT から離れていくのが見えます。近似が信頼できるのは「有効帯域」の内側だけです。制御帯域がこの有効帯域を超えていると、安定余裕の見積もりが楽観的になり、実機で予想外の発振を招きます。

最後に、初期アンダーシュート(逆応答)を「数値誤差」や「シミュレーションのバグ」と誤認することです。ステップ応答が立ち上がる前に逆へ動くのは、右半面零点を持つ非最小位相系の正しい挙動です。これはパデ近似が「むだ時間の待ち時間」を有理系で表現するために避けられない代償であり、次数を上げると小さくはなりますが完全には消えません。逆応答を見て近似を疑うのではなく、それが非最小位相の現れだと理解しておくことが、むだ時間系の解析では重要です。

使い方ガイド

  1. むだ時間T(秒)と近似次数n(1~6)をスライダで設定します。例:T=0.5秒、n=3次
  2. ゲインkと時定数τを調整し、閉ループシステムの安定性を確認します。k=2、τ=0.1秒の場合、右半面零点の生成を観察できます
  3. 「シミュレート」ボタンでパデ近似の周波数応答とステップ応答を描画。位相誤差が5°以下の有効帯域幅がrad/s単位で表示されます

具体的な計算例

むだ時間T=1.0秒をn=4次パデ近似する場合、伝達関数e^(-sT)≈(s⁴-20s³+180s²-840s+1680)/(s⁴+20s³+180s²+840s+1680)となります。ω=1rad/sで位相誤差は約2.8°、有効帯域は0~0.4rad/sです。同じT=1.0秒をn=2次で近似すると位相誤差は12.5°に拡大し、制御系の応答性が大幅に低下します。プロセス制御でT=2秒のバルブむだ時間にn=5次パデ近似を適用すればステップ応答の二乗平均誤差は0.015以下に抑制可能です

実務での注意点

  1. n≥3で近似次数を増やすと右半面零点が出現し、制御ゲイン設定時のアンダーシュート(-15~-25%)を招きやすくなります。PID制御では微分項の追加が必須です
  2. 有効帯域が狭い場合(T=2秒、n=2次で0.2rad/s以下)は、スミス予測器やむだ時間補償器の導入を検討してください
  3. ステップ応答の二乗平均誤差が0.05を超える場合、実装用には高次パデ近似(n=5~6)またはモデル予測制御への切り替えが必要です