PI制御チューニング : Ziegler Nicholsの限界感度法の続きです.

0. むだ時間 - 導入

システムの入力に対する出力の時間遅れの1種で、時間領域で以下のように表現でき、
$$y(t)=x(t-T_d)$$ 伝達関数では、
$$G(s)=e^{-s{T}_d}$$ と表される.

このような性質が付与されたシステムを制御する時、
制御コントローラにしてみれば、制御してから結果に反映されるまで時間がかかる状態になるので、
ゲインを高く設定しすぎるとハンチングを起こすなどそわそわしている状態になります。

こういった性質は制御分野以外でも最近とりあげられていて、
「世界はシステムで動く」というドネラさんの本でも、
例えば在庫管理的な場面でむだ時間が生じると現象が振動的になるという話が出ていました。
(制御工学で社会分析している印象を受ける本です. 応用としては面白いかも.)

• 世界はシステムで動く ―― いま起きていることの本質をつかむ考え方
• ドネラ・H・メドウズ

⇒ むだ時間は実世界では必ず存在し、その影響が制御設計やチューニングを難しくします

1. 1次遅れ系とむだ時間

むだ時間が付与された1次遅れ系のベクトル軌跡を見てみる.
1次遅れ系の時定数を時定数1(sec), むだ時間を1(sec)として伝達関数は、
$$G_2(s)=\frac{1}{s+1}{e}^{-s}$$

これを複素平面で表現するため極形式にすると、
$$G_2(\omega i)=\frac{1}{\omega i+1}{e}^{-\omega i}$$
から
$$z=|G_2(\omega i)|\cdot exp(i\cdot arg(G_2(\omega i)))=\frac{1}{\sqrt{{\omega }^2+1}}exp(i(-arctan(\omega )+(-\omega )))$$
これを図で表すと分かる通り、時間遅れのない状態(水色)に時間遅れをつけるとオレンジの線のようにぐるぐるができ、不安定点(Re=-1)に近づき、より不安定になっていることがわかります.

水色（$\frac{1}{s+1}$)、オレンジ($\frac{1}{s+1}{e}^{-s}$)
Note: お絵かきにはエクセルを使い、ω = [0.01, 100] (rad/s)の範囲で計算しました.

$\omega$	r	$\theta$	Real	Im
0.01 ～ 100	= 1/sqrt(A2^2+1)	=- atan2(1, A2) - A2	=B2 * cos(C2)	=B2 * sin(C2)

$\omega$

r

$\theta$

Real

Im

0.01 ～ 100

= 1/sqrt(A2^2+1)

=- atan2(1, A2) - A2

=B2 * cos(C2)

=B2 * sin(C2)

2. むだ時間1次遅れ系のPI制御へジーグラ・ニコルス限界感度法を適用

先の$P(s)=\frac{1}{\tau s+1}$ として、プラント $P_2(s)=P(s)e^{-Ls}$ にPI制御 $C(s)=K_p+\frac{K_i}{s}$ を施すと、
$$G(s)=\frac{P\cdot C\cdot e^{-Ls}}{1+P\cdot C\cdot e^{-Ls}}=\frac{K_p{e}^{-Ls}s+K_i{e}^{-Ls}}{\tau s^2+(1+K_p{e}^{-Ls})s+K_i{e}^{-Ls}}$$
$K_u$を求めるために、$K_i=0$とすると、
$$G(s)=\frac{K_p{e}^{-Ls}}{\tau s+1+K_p{e}^{-Ls}}$$
これに対し、当初、$K_u$を求めるために$e^{-Ls}$のパデ近似を用いて代数的に範囲を求めようと考えたのですが、文献[2]で検証している通り、特に1次,2次などの低次のパデ近似では良い解にたどりつけないようです.

$K_u$を求めるにあたってナイキスト線図を使って幾何的に考えることにします.

まず、$K_p=1$として
$$G(s)=\frac{e^{-Ls}}{\tau s+1+K_p{e}^{-Ls}}$$

のナイキスト線図を描きます．

この時

$$z=|G(\omega i)|\cdot exp(i\cdot arg(G(\omega i)))=\frac{1}{\sqrt{(\tau \omega {)}^2+1}}exp(i(-arctan(\tau \omega )+(-L\omega )))$$

を描くことになるが、以下図中の緑の点の座標を求めます．

$\omega$が大きくなるほどゼロ点に近づくので、
zが複素平面左半面の実軸上と交わる正で最小の${\omega }_0$を求めればよい.

従って 、$-arctan(\tau {\omega }_0)+(-L{\omega }_0)=-\pi$で、最小のものが求める周波数です．
※ この時、$z({\omega }_0)=-\frac{1}{\sqrt{(\tau {\omega }_0{)}^2+1}}$

ここで $K_p=\sqrt{(\tau {\omega }_0{)}^2+1}$とするとベクトル線図は以下の図の水色の線のようになり、これが $K_u$となります．

拡大した図

まとめ:
以上の検討から$K_p$の限界である$K_u$は、
$$K_u=\sqrt{(\tau {\omega }_0{)}^2+1}$$ ただし、
$$arctan(\tau {\omega }_0)=\pi -(L{\omega }_0)$$

($\tau \cdot {\omega }_0$は、(0,$\pi /2$)の範囲 )

また、この時、振動周期(s)は$T_u=\frac{2\pi }{{\omega }_0}$である.

${\omega }_0$は数値解析で求めるほかありませんが、二分法など使うと良いでしょう.

二分法 - 高精度計算サイト http://keisan.casio.jp/exec/system/1180917668

2.1. お試し計算

例えば、$$G(s)=\frac{1}{s+1}{e}^{-s}$$として$K_u$を求めると、$K_u\approx 2.26({\omega }_0\approx 2.029;T_u\approx 3.097sec)$と求まります．

シミュレーションで、この値を使ってステップ応答を求めると以下の結果になります．

ジーグラ・ニコルスの限界感度法に従って、

制御の種類	$K_p$	$K_d$	$K_i$
P制御	$0.5K_u$	0	0
PI制御	$0.45K_u$	0	$\frac{0.45K_u}{0.83T_u}$
PID制御	$0.6K_u$	$0.075T_u{K}_u$	$\frac{0.6K_u}{0.5T_u}$

から$K_p=1.018$, $K_i=0.396$としてステップ応答を求めたものがこちら.

やっとZNの限界感度法がうまく行きました.

参考

[1] The Time Delay
http://lpsa.swarthmore.edu/BackGround/TimeDelay/TimeDelay.html

[2] PID Controllers for Systems with Time-Delay
http://msc.berkeley.edu/PID/modernPID3-delay.pdf