【ステップバイステップ解法】
(1) 目標値 R(s) から制御量 Y(s) までの閉ループ伝達関数 W(s)
ブロック線図
を見ると、内側の加算点に速度帰還
K2sY(s)、外側の加算点に主フィードバック
Y(s) がそれぞれ負帰還されています。各加算点の出力信号を定義して立式します。
① 外側の減算点の出力(制御偏差)を E(s) とすると、
E(s)=R(s)−Y(s)
② これに比例ゲイン K1 が掛かり、さらに内側の減算点で帰還信号 K2sY(s) が減算されます。制御対象 s(s+2)(s+4)1 への入力信号を U(s) とおくと、
U(s)=K1E(s)−K2sY(s)=K1{R(s)−Y(s)}−K2sY(s)
③ 出力 Y(s) は、入力 U(s) に制御対象の伝達関数を乗じたものであるため、
Y(s)=s(s+2)(s+4)1U(s)=s(s+2)(s+4)K1R(s)−(K1+K2s)Y(s)
両辺に分母 s(s+2)(s+4) を掛けます。
s(s+2)(s+4)Y(s)=K1R(s)−(K1+K2s)Y(s)
Y(s) を含む項をすべて左辺に移項して整理します。
[s(s+2)(s+4)+K2s+K1]Y(s)=K1R(s)
ここで、多項式を展開します。
s(s+2)(s+4)=s(s2+6s+8)=s3+6s2+8s
したがって、括弧内は
s3+6s2+(8+K2)s+K1
となり、閉ループ伝達関数 W(s)=R(s)Y(s) は次のようになります。
W(s)=s3+6s2+(8+K2)s+K1K1
(2) 特性根を −1,−2 とするゲイン K1,K2 およびその他の特性根
閉ループ系の特性方程式は、伝達関数 W(s) の分母多項式を 0 と置いたものです。
s3+6s2+(8+K2)s+K1=0
この3次方程式の3つの根を −1、−2、そして残りのもう一つの根を s3 とおきます。
【解法:根と係数の関係の利用】
3次の代数方程式 s3+a2s2+a1s+a0=0 において、根の和は −a2 に等しくなります。
本問の特性方程式における s2 の係数は 6 ですから、
(−1)+(−2)+s3=−6
−3+s3=−6⟹s3=−3
これにより、直ちにその他の特性根が −3 であることが求まります。
したがって、特性多項式は次のように因数分解されます。
(s+1)(s+2)(s+3)=0
この左辺を展開します。
(s2+3s+2)(s+3)=s3+6s2+11s+6=0
特性方程式 s3+6s2+(8+K2)s+K1=0 と各次数の係数を比較すると、
- s の係数:8+K2=11⟹K2=3
- 定数項:K1=6
K1=6,K2=3,その他の特性根: −3
(3) 単位インパルス応答 y(t)
小問(2)で得られた K1=6、K2=3 を代入すると、閉ループ伝達関数は次のようになります。
W(s)=(s+1)(s+2)(s+3)6
入力が単位インパルス信号 r(t)=δ(t) のとき、そのラプラス変換は R(s)=1 です。したがって、出力 Y(s) は伝達関数そのものに一致します。
Y(s)=W(s)R(s)=(s+1)(s+2)(s+3)6
この有理関数を次のように部分分数分解します。
Y(s)=s+1A+s+2B+s+3C
各係数 A,B,C はヘビサイドの展開定理(留数計算)を用いて効率的に求められます。
- A=s→−1lim(s+1)Y(s)=(−1+2)(−1+3)6=(1)(2)6=3
- B=s→−2lim(s+2)Y(s)=(−2+1)(−2+3)6=(−1)(1)6=−6
- C=s→−3lim(s+3)Y(s)=(−3+1)(−3+2)6=(−2)(−1)6=3
よって、Y(s) は次のように展開されます。
Y(s)=s+13−s+26+s+33
公式 L−1[s+a1]=e−at を用いて各項をラプラス逆変換すると、単位インパルス応答 y(t) が求まります。
y(t)=3e−t−6e−2t+3e−3t(t≥0)
【検算・物理的妥当性の確認】
t=0 を代入すると
y(0)=3(1)−6(1)+3(1)=0 となり、制御対象(3次遅れ系)の初期応答として整合します。
また、
t→∞ で
y(t)→0 となり、インパルス入力に対する安定な系の応答として完全に成立していることが確認できます。
【最終結論】
(1) 閉ループ伝達関数 W(s)
W(s)=s3+6s2+(8+K2)s+K1K1
(2) ゲイン K1,K2 およびその他の特性根
K1=6,K2=3,その他の特性根: −3
(3) 単位インパルス応答 y(t)
y(t)=3e−t−6e−2t+3e−3t(t≥0)