問題のポイント・使用する基本公式
本問は、フィードフォワード制御および局所フィードバックを含む2自由度制御系の構造解析、安定判別、定常特性を総合的に問う問題です。
- インパルス応答:伝達関数 G(s) の逆ラプラス変換により得られます。
公式:L[sinωt]=s2+ω2ω⟺L−1[s2+ω21]=ω1sinωt
- ブロック線図の基本関係式(図 参照):
E(s)U(s)Y(s)=R(s)−Y(s)=F(s)R(s)+C1(s)E(s)−C2(s)Y(s)=G(s)U(s)
- ラウス・フルビッツの安定判別法:3次の特性方程式 a3s3+a2s2+a1s+a0=0(全係数正)において、系が安定である条件は a2a1−a3a0>0、安定限界となる条件は a2a1−a3a0=0 です。
- 最終値の定理(定常偏差):閉ループ系が安定であるとき、定常偏差 ess は次式で求まります。
ess=t→∞lime(t)=s→0limsE(s)
(1) 制御対象の単位インパルス応答
単位インパルス入力 δ(t) のラプラス変換は 1 です。したがって、制御対象単体の出力 Y(s) は伝達関数そのものとなります。
Y(s)=G(s)⋅1=s2+41=s2+221
正弦波の逆ラプラス変換の公式を適用できるよう、係数を調整します。
Y(s)=21⋅s2+222
両辺を逆ラプラス変換すると、求めるインパルス応答 g(t) は次のようになります。
g(t)=21sin2t(t≥0)
(2) F(s)=0 のときの伝達関数 R(s)E(s)
図
より、ブロック線図の基本関係式を立てます。
E(s)U(s)Y(s)=R(s)−Y(s)=C1(s)E(s)−C2(s)Y(s)=G(s)U(s)⋯①⋯②⋯③
③式を②式に代入して U(s) について整理します。
U(s){1+C2(s)G(s)}U(s)U(s)=C1(s)E(s)−C2(s)G(s)U(s)=C1(s)E(s)=1+C2(s)G(s)C1(s)E(s)
これを Y(s)=G(s)U(s) に代入すると、Y(s)=1+C2(s)G(s)C1(s)G(s)E(s) となり、①式に代入します。
E(s)R(s)R(s)E(s)=R(s)−1+C2(s)G(s)C1(s)G(s)E(s)=[1+1+C2(s)G(s)C1(s)G(s)]E(s)=1+C2(s)G(s)1+C1(s)G(s)+C2(s)G(s)E(s)=1+C1(s)G(s)+C2(s)G(s)1+C2(s)G(s)
与えられた各関数 G(s)=s2+41、C1(s)=K1(1+s1)=sK1(s+1)、C2(s)=K2s を代入します。
R(s)E(s)=1+sK1(s+1)⋅s2+41+K2s⋅s2+411+K2s⋅s2+41
分子・分母に s(s2+4) を乗じて整理します。
R(s)E(s)=s(s2+4)+K1(s+1)+K2s2s(s2+4)+K2s2=s3+4s+K1s+K1+K2s2s3+4s+K2s2
降べきの順に整理すると、求める伝達関数が得られます。
R(s)E(s)=s3+K2s2+(K1+4)s+K1s3+K2s2+4s(=s3+K2s2+(K1+4)s+K1s(s2+K2s+4))
(3) 制御系が安定限界となるときの条件と応答
閉ループ系の特性方程式は、(2)で求めた伝達関数の分母多項式を 0 と置いたものです。
D(s)=s3+K2s2+(K1+4)s+K1=0
ラウスの安定判別法を適用するためにラウス表を作成します。
s31K1+4 s2K2K1 s1K2K2(K1+4)−K10 s0K1
K1>0、K2>0 であるため、系が安定限界となるのは s1 行の第1列要素が 0 になるときです。
K2K2(K1+4)−K1=0⟹K2(K1+4)=K1(または K2=K1+4K1)
【安定限界における応答】
この条件 K1=K2(K1+4) を特性方程式に代入して因数分解すると、
D(s)=s3+K2s2+(K1+4)s+K2(K1+4)=(s+K2)(s2+K1+4)=0
特性根は s=−K2 および純虚数根 s=±jK1+4 となります。実部が負の根 s=−K2 による成分は減衰しますが、虚軸上の根 s=±jK1+4 による成分は減衰も発散もしません。したがって、制御系は振幅が減衰しない持続振動(等振幅の周期振動)を示します。
関係:K2(K1+4)=K1(または K1K2+4K2−K1=0)
応答:減衰しない一定振幅の持続振動(角周波数 ω=K1+4 [rad/s] の振動)を示す。
(4) 空欄 (a)~(f) の導出
図
のブロック線図から関係式を立てます。
操作量 U(s) の合流点に注目すると、
U(s)=F(s)R(s)+C1(s)E(s)−C2(s)Y(s)
ここで Y(s)=G(s)U(s) を代入すると、
U(s)=F(s)R(s)+C1(s)E(s)−C2(s)G(s)U(s)
また、偏差 E(s) の合流点に注目すると、
E(s)=R(s)−Y(s)=R(s)−G(s)U(s)
問題文の式形式と係数を比較することで、空欄 (a)~(d) が求まります。
- (a):F(s)
- (b):C1(s)
- (c):C2(s)G(s)
- (d):G(s)
続いて、U(s) を消去して伝達関数 R(s)E(s) を求めます。
U(s) の式をまとめると、
{1+C2(s)G(s)}U(s)=F(s)R(s)+C1(s)E(s)⟹U(s)=1+C2(s)G(s)F(s)R(s)+C1(s)E(s)
これを E(s)=R(s)−G(s)U(s) に代入します。
E(s){1+C2(s)G(s)}E(s){1+C1(s)G(s)+C2(s)G(s)}E(s)=R(s)−G(s)1+C2(s)G(s)F(s)R(s)+C1(s)E(s)={1+C2(s)G(s)}R(s)−F(s)G(s)R(s)−C1(s)G(s)E(s)={1+C2(s)G(s)−F(s)G(s)}R(s)
両辺を整理して比を取ると、
R(s)E(s)=1+{C1(s)G(s)+C2(s)G(s)}1+{C2(s)G(s)−F(s)G(s)}
問題文の形式 1+ (f) 1+ (e) と比較して、空欄 (e) および (f) が定まります。
(a) F(s)
(b) C1(s)
(c) C2(s)G(s)
(d) G(s)
(e) C2(s)G(s)−F(s)G(s) (または {C2(s)−F(s)}G(s))
(f) C1(s)G(s)+C2(s)G(s) (または {C1(s)+C2(s)}G(s))
(5) ランプ状目標値 r(t)=t に対する定常偏差
(4)で得た伝達関数に、F(s)=4、C1(s)=K1(1+s1)=sK1(s+1)、C2(s)=K2s、G(s)=s2+41 を代入します。
R(s)E(s)=1+s(s2+4)K1(s+1)+s2+4K2s1+K2s⋅s2+41−4⋅s2+41=s(s2+4)+K1(s+1)+K2s2s(s2+4)+K2s2−4s
分子の 4s と −4s が相殺されるため、
R(s)E(s)=s3+K2s2+(K1+4)s+K1s3+K2s2=s3+K2s2+(K1+4)s+K1s2(s+K2)
ランプ状の目標値 r(t)=t のラプラス変換は R(s)=s21 です。したがって、偏差 E(s) は次のようになります。
E(s)=R(s)E(s)R(s)=s3+K2s2+(K1+4)s+K1s2(s+K2)⋅s21=s3+K2s2+(K1+4)s+K1s+K2
系が安定である(K2(K1+4)>K1)とき、最終値の定理を適用して定常偏差 ess を求めます。
ess=t→∞lime(t)=s→0limsE(s)=s→0lim[s⋅s3+K2s2+(K1+4)s+K1s+K2]=0⋅K1K2=0
0
解答のまとめ
- (1) 制御対象のインパルス応答:g(t)=21sin2t(t≥0)
- (2) 伝達関数:R(s)E(s)=s3+K2s2+(K1+4)s+K1s3+K2s2+4s
- (3) 安定限界条件と応答:K2(K1+4)=K1、応答は減衰しない持続振動となる。
- (4) 空欄:(a) F(s)、(b) C1(s)、(c) C2(s)G(s)、(d) G(s)、(e) C2(s)G(s)−F(s)G(s)、(f) C1(s)G(s)+C2(s)G(s)
- (5) ランプ入力に対する定常偏差:0