過去問一覧第二種電気主任技術者2013年度 機械・制御フィードフォワード補償を含む制御系の伝達関数と安定判別
一覧に戻る

フィードフォワード補償を含む制御系の伝達関数と安定判別

PID
第二種電気主任技術者 2013年度 機械・制御 問4

図の制御系において,R(s) R(s) ,E(s) E(s) ,U(s) U(s) 及びY(s) Y(s) は,目標値r(t) r(t) ,偏差e(t) e(t) ,操作量u(t) u(t) 及び出力y(t) y(t) をそれぞれラプラス変換したものである。以下では,K1>0 K_{1}>0 ,K2>0 K_{2}>0 とする。次の問に答えよ。

問題図
(1) 図の制御対象の伝達関数がG(s)=1s2+4 \displaystyle G(s)=\frac {1}{s^{2}+4} で与えられるとき,制御対象の単位インパルス応答を求めよ。
(2) 上記(1)において,さらにF(s)=0 F(s)=0 ,C1(s)=K1(1+1s) \displaystyle C_{1}(s)=K_{1}\left( 1+\frac {1}{s}\right) ,C2(s)=K2s C_{2}(s)=K_{2}s のとき,R(s) R(s) からE(s) E(s) までの伝達関数をK1 K_{1} とK2 K_{2} を用いて表せ。
(3) 上記(2)において,図の制御系が安定限界となるとき,K1 K_{1} とK2 K_{2} が満たすべき関係を求めよ。またこのとき制御系が示す応答について簡単に述べよ。
(4) 次に,F(s)≠0 F(s)≠0 の場合を考える。次の空欄(a)~(f)に入る数式をF(s) F(s) ,C1(s) C_{1}(s) ,C2(s) C_{2}(s) ,G(s) G(s) を用いてできるだけ簡略化して表せ。

図から
U(s)= (a) ⋅R(s)+ (b) ⋅E(s)− (c) ⋅U(s)E(s)=R(s)− (d) ⋅U(s)\begin{aligned}U(s)&=\fbox { (a) } \cdot R(s) +\fbox { (b) } \cdot E(s) -\fbox { (c) } \cdot U(s) \\[ 5pt ] E(s)&=R(s) -\fbox { (d) } \cdot U(s) \\[ 5pt ] \end{aligned} が成り立つ。これらから,R(s) R(s) からE(s) E(s) までの伝達関数を求めると
E(s)R(s)=1+ (e) 1+ (f) \begin{aligned}\frac {E(s)}{R(s)}&=\frac {1+\fbox { (e) }}{1+\fbox { (f) }} \\[ 5pt ]\end{aligned} となる。
(5) 上記(4)において,F(s)=4 F(s)=4 ,C1(s)=K1(1+1s) \displaystyle C_{1}(s)=K_{1}\left( 1+\frac {1}{s}\right) ,C2(s)=K2s C_{2}(s)=K_{2}s ,G(s)=1s2+4 \displaystyle G(s)=\frac {1}{s^{2}+4} を代入したとき,ランプ状の目標値r(t)=t r(t)=t に対する定常偏差を求めよ。

解説

問題のポイント・使用する基本公式

本問は、フィードフォワード制御および局所フィードバックを含む2自由度制御系の構造解析、安定判別、定常特性を総合的に問う問題です。

  • インパルス応答:伝達関数 G(s)G(s) の逆ラプラス変換により得られます。
    公式:L[sin⁡ωt]=ωs2+ω2  ⟺  L−1[1s2+ω2]=1ωsin⁡ωt\mathcal{L}[\sin \omega t] = \displaystyle\frac{\omega}{s^2+\omega^2} \iff \mathcal{L}^{-1}\left[\frac{1}{s^2+\omega^2}\right] = \frac{1}{\omega}\sin \omega t
  • ブロック線図の基本関係式(図
    問題図
    参照): E(s)=R(s)−Y(s)U(s)=F(s)R(s)+C1(s)E(s)−C2(s)Y(s)Y(s)=G(s)U(s)\begin{aligned} E(s) &= R(s) - Y(s) \\[4pt] U(s) &= F(s)R(s) + C_1(s)E(s) - C_2(s)Y(s) \\[4pt] Y(s) &= G(s)U(s) \end{aligned}
  • ラウス・フルビッツの安定判別法:3次の特性方程式 a3s3+a2s2+a1s+a0=0a_3 s^3 + a_2 s^2 + a_1 s + a_0 = 0(全係数正)において、系が安定である条件は a2a1−a3a0>0a_2 a_1 - a_3 a_0 > 0、安定限界となる条件は a2a1−a3a0=0a_2 a_1 - a_3 a_0 = 0 です。
  • 最終値の定理(定常偏差):閉ループ系が安定であるとき、定常偏差 esse_{\mathrm{ss}} は次式で求まります。 ess=lim⁡t→∞e(t)=lim⁡s→0sE(s)e_{\mathrm{ss}} = \lim_{t \to \infty} e(t) = \lim_{s \to 0} s E(s)

(1) 制御対象の単位インパルス応答

単位インパルス入力 δ(t)\delta(t) のラプラス変換は 11 です。したがって、制御対象単体の出力 Y(s)Y(s) は伝達関数そのものとなります。

Y(s)=G(s)⋅1=1s2+4=1s2+22Y(s) = G(s) \cdot 1 = \frac{1}{s^2 + 4} = \frac{1}{s^2 + 2^2}

正弦波の逆ラプラス変換の公式を適用できるよう、係数を調整します。

Y(s)=12⋅2s2+22Y(s) = \frac{1}{2} \cdot \frac{2}{s^2 + 2^2}

両辺を逆ラプラス変換すると、求めるインパルス応答 g(t)g(t) は次のようになります。

g(t)=12sin⁡2t(t≥0)g(t) = \displaystyle\frac{1}{2}\sin 2t \quad (t \ge 0)

(2) F(s)=0F(s)=0 のときの伝達関数 E(s)R(s)\displaystyle\frac{E(s)}{R(s)}

図

問題図
より、ブロック線図の基本関係式を立てます。

E(s)=R(s)−Y(s)⋯①U(s)=C1(s)E(s)−C2(s)Y(s)⋯②Y(s)=G(s)U(s)⋯③\begin{aligned} E(s) &= R(s) - Y(s) & \cdots ① \\[5pt] U(s) &= C_1(s)E(s) - C_2(s)Y(s) & \cdots ② \\[5pt] Y(s) &= G(s)U(s) & \cdots ③ \end{aligned}

③式を②式に代入して U(s)U(s) について整理します。

U(s)=C1(s)E(s)−C2(s)G(s)U(s){1+C2(s)G(s)}U(s)=C1(s)E(s)U(s)=C1(s)1+C2(s)G(s)E(s)\begin{aligned} U(s) &= C_1(s)E(s) - C_2(s)G(s)U(s) \\[5pt] \{1 + C_2(s)G(s)\}U(s) &= C_1(s)E(s) \\[5pt] U(s) &= \frac{C_1(s)}{1 + C_2(s)G(s)}E(s) \end{aligned}

これを Y(s)=G(s)U(s)Y(s) = G(s)U(s) に代入すると、Y(s)=C1(s)G(s)1+C2(s)G(s)E(s)Y(s) = \displaystyle\frac{C_1(s)G(s)}{1 + C_2(s)G(s)}E(s) となり、①式に代入します。

E(s)=R(s)−C1(s)G(s)1+C2(s)G(s)E(s)R(s)=[1+C1(s)G(s)1+C2(s)G(s)]E(s)=1+C1(s)G(s)+C2(s)G(s)1+C2(s)G(s)E(s)E(s)R(s)=1+C2(s)G(s)1+C1(s)G(s)+C2(s)G(s)\begin{aligned} E(s) &= R(s) - \frac{C_1(s)G(s)}{1 + C_2(s)G(s)}E(s) \\[5pt] R(s) &= \left[ 1 + \frac{C_1(s)G(s)}{1 + C_2(s)G(s)} \right] E(s) = \frac{1 + C_1(s)G(s) + C_2(s)G(s)}{1 + C_2(s)G(s)} E(s) \\[5pt] \frac{E(s)}{R(s)} &= \frac{1 + C_2(s)G(s)}{1 + C_1(s)G(s) + C_2(s)G(s)} \end{aligned}

与えられた各関数 G(s)=1s2+4G(s) = \displaystyle\frac{1}{s^2+4}、C1(s)=K1(1+1s)=K1(s+1)sC_1(s) = K_1\left(1 + \frac{1}{s}\right) = \frac{K_1(s+1)}{s}、C2(s)=K2sC_2(s) = K_2 s を代入します。

E(s)R(s)=1+K2s⋅1s2+41+K1(s+1)s⋅1s2+4+K2s⋅1s2+4\frac{E(s)}{R(s)} = \frac{1 + K_2 s \cdot \displaystyle\frac{1}{s^2+4}}{1 + \frac{K_1(s+1)}{s} \cdot \displaystyle\frac{1}{s^2+4} + K_2 s \cdot \displaystyle\frac{1}{s^2+4}}

分子・分母に s(s2+4)s(s^2+4) を乗じて整理します。

E(s)R(s)=s(s2+4)+K2s2s(s2+4)+K1(s+1)+K2s2=s3+4s+K2s2s3+4s+K1s+K1+K2s2\begin{aligned} \frac{E(s)}{R(s)} &= \frac{s(s^2+4) + K_2 s^2}{s(s^2+4) + K_1(s+1) + K_2 s^2} \\[5pt] &= \frac{s^3 + 4s + K_2 s^2}{s^3 + 4s + K_1 s + K_1 + K_2 s^2} \end{aligned}

降べきの順に整理すると、求める伝達関数が得られます。

E(s)R(s)=s3+K2s2+4ss3+K2s2+(K1+4)s+K1(=s(s2+K2s+4)s3+K2s2+(K1+4)s+K1)\displaystyle\frac{E(s)}{R(s)} = \frac{s^3 + K_2 s^2 + 4s}{s^3 + K_2 s^2 + (K_1 + 4)s + K_1} \quad \left( = \frac{s(s^2 + K_2 s + 4)}{s^3 + K_2 s^2 + (K_1 + 4)s + K_1} \right)

(3) 制御系が安定限界となるときの条件と応答

閉ループ系の特性方程式は、(2)で求めた伝達関数の分母多項式を 00 と置いたものです。

D(s)=s3+K2s2+(K1+4)s+K1=0D(s) = s^3 + K_2 s^2 + (K_1 + 4)s + K_1 = 0

ラウスの安定判別法を適用するためにラウス表を作成します。

s31K1+4 s2K2K1 s1K2(K1+4)−K1K20 s0K1\begin{array}{c|cc} s^3 & 1 & K_1 + 4 \ s^2 & K_2 & K_1 \ s^1 & \displaystyle\frac{K_2(K_1+4) - K_1}{K_2} & 0 \ s^0 & K_1 & \end{array}

K1>0K_1 > 0、K2>0K_2 > 0 であるため、系が安定限界となるのは s1s^1 行の第1列要素が 00 になるときです。

K2(K1+4)−K1K2=0  ⟹  K2(K1+4)=K1(または K2=K1K1+4)\frac{K_2(K_1 + 4) - K_1}{K_2} = 0 \implies K_2(K_1 + 4) = K_1 \quad \left( \text{または } K_2 = \frac{K_1}{K_1+4} \right)

【安定限界における応答】
この条件 K1=K2(K1+4)K_1 = K_2(K_1+4) を特性方程式に代入して因数分解すると、

D(s)=s3+K2s2+(K1+4)s+K2(K1+4)=(s+K2)(s2+K1+4)=0D(s) = s^3 + K_2 s^2 + (K_1 + 4)s + K_2(K_1 + 4) = (s + K_2)(s^2 + K_1 + 4) = 0

特性根は s=−K2s = -K_2 および純虚数根 s=±jK1+4s = \pm j\sqrt{K_1+4} となります。実部が負の根 s=−K2s = -K_2 による成分は減衰しますが、虚軸上の根 s=±jK1+4s = \pm j\sqrt{K_1+4} による成分は減衰も発散もしません。したがって、制御系は振幅が減衰しない持続振動(等振幅の周期振動)を示します。

関係:K2(K1+4)=K1K_2(K_1 + 4) = K_1(または K1K2+4K2−K1=0K_1 K_2 + 4K_2 - K_1 = 0)
応答:減衰しない一定振幅の持続振動(角周波数 ω=K1+4 [rad/s]\omega = \sqrt{K_1 + 4}\ [\mathrm{rad/s}] の振動)を示す。

(4) 空欄 (a)~(f) の導出

図

問題図
のブロック線図から関係式を立てます。

操作量 U(s)U(s) の合流点に注目すると、

U(s)=F(s)R(s)+C1(s)E(s)−C2(s)Y(s)U(s) = F(s)R(s) + C_1(s)E(s) - C_2(s)Y(s)

ここで Y(s)=G(s)U(s)Y(s) = G(s)U(s) を代入すると、

U(s)=F(s)R(s)+C1(s)E(s)−C2(s)G(s)U(s)U(s) = F(s)R(s) + C_1(s)E(s) - C_2(s)G(s)U(s)

また、偏差 E(s)E(s) の合流点に注目すると、

E(s)=R(s)−Y(s)=R(s)−G(s)U(s)E(s) = R(s) - Y(s) = R(s) - G(s)U(s)

問題文の式形式と係数を比較することで、空欄 (a)~(d) が求まります。

  • (a):F(s)F(s)
  • (b):C1(s)C_1(s)
  • (c):C2(s)G(s)C_2(s)G(s)
  • (d):G(s)G(s)

続いて、U(s)U(s) を消去して伝達関数 E(s)R(s)\displaystyle\frac{E(s)}{R(s)} を求めます。
U(s)U(s) の式をまとめると、

{1+C2(s)G(s)}U(s)=F(s)R(s)+C1(s)E(s)  ⟹  U(s)=F(s)R(s)+C1(s)E(s)1+C2(s)G(s)\{1 + C_2(s)G(s)\}U(s) = F(s)R(s) + C_1(s)E(s) \implies U(s) = \frac{F(s)R(s) + C_1(s)E(s)}{1 + C_2(s)G(s)}

これを E(s)=R(s)−G(s)U(s)E(s) = R(s) - G(s)U(s) に代入します。

E(s)=R(s)−G(s)F(s)R(s)+C1(s)E(s)1+C2(s)G(s){1+C2(s)G(s)}E(s)={1+C2(s)G(s)}R(s)−F(s)G(s)R(s)−C1(s)G(s)E(s){1+C1(s)G(s)+C2(s)G(s)}E(s)={1+C2(s)G(s)−F(s)G(s)}R(s)\begin{aligned} E(s) &= R(s) - G(s)\frac{F(s)R(s) + C_1(s)E(s)}{1 + C_2(s)G(s)} \\[5pt] \{1 + C_2(s)G(s)\}E(s) &= \{1 + C_2(s)G(s)\}R(s) - F(s)G(s)R(s) - C_1(s)G(s)E(s) \\[5pt] \{1 + C_1(s)G(s) + C_2(s)G(s)\}E(s) &= \{1 + C_2(s)G(s) - F(s)G(s)\}R(s) \end{aligned}

両辺を整理して比を取ると、

E(s)R(s)=1+{C2(s)G(s)−F(s)G(s)}1+{C1(s)G(s)+C2(s)G(s)}\frac{E(s)}{R(s)} = \frac{1 + \{C_2(s)G(s) - F(s)G(s)\}}{1 + \{C_1(s)G(s) + C_2(s)G(s)\}}

問題文の形式 1+ (e) 1+ (f) \displaystyle\frac{1 + \fbox{ (e) }}{1 + \fbox{ (f) }} と比較して、空欄 (e) および (f) が定まります。

(a) F(s)F(s)
(b) C1(s)C_1(s)
(c) C2(s)G(s)C_2(s)G(s)
(d) G(s)G(s)
(e) C2(s)G(s)−F(s)G(s)C_2(s)G(s) - F(s)G(s) (または {C2(s)−F(s)}G(s)\{C_2(s)-F(s)\}G(s))
(f) C1(s)G(s)+C2(s)G(s)C_1(s)G(s) + C_2(s)G(s) (または {C1(s)+C2(s)}G(s)\{C_1(s)+C_2(s)\}G(s))

(5) ランプ状目標値 r(t)=tr(t)=t に対する定常偏差

(4)で得た伝達関数に、F(s)=4F(s)=4、C1(s)=K1(1+1s)=K1(s+1)sC_1(s)=K_1\left(1+\displaystyle\frac{1}{s}\right) = \frac{K_1(s+1)}{s}、C2(s)=K2sC_2(s)=K_2 s、G(s)=1s2+4G(s)=\displaystyle\frac{1}{s^2+4} を代入します。

E(s)R(s)=1+K2s⋅1s2+4−4⋅1s2+41+K1(s+1)s(s2+4)+K2ss2+4=s(s2+4)+K2s2−4ss(s2+4)+K1(s+1)+K2s2\begin{aligned} \frac{E(s)}{R(s)} &= \frac{1 + K_2 s \cdot \frac{1}{s^2+4} - 4 \cdot \frac{1}{s^2+4}}{1 + \frac{K_1(s+1)}{s(s^2+4)} + \frac{K_2 s}{s^2+4}} \\[5pt] &= \frac{s(s^2+4) + K_2 s^2 - 4s}{s(s^2+4) + K_1(s+1) + K_2 s^2} \end{aligned}

分子の 4s4s と −4s-4s が相殺されるため、

E(s)R(s)=s3+K2s2s3+K2s2+(K1+4)s+K1=s2(s+K2)s3+K2s2+(K1+4)s+K1\frac{E(s)}{R(s)} = \frac{s^3 + K_2 s^2}{s^3 + K_2 s^2 + (K_1+4)s + K_1} = \frac{s^2(s + K_2)}{s^3 + K_2 s^2 + (K_1+4)s + K_1}

ランプ状の目標値 r(t)=tr(t)=t のラプラス変換は R(s)=1s2R(s) = \displaystyle\frac{1}{s^2} です。したがって、偏差 E(s)E(s) は次のようになります。

E(s)=E(s)R(s)R(s)=s2(s+K2)s3+K2s2+(K1+4)s+K1⋅1s2=s+K2s3+K2s2+(K1+4)s+K1E(s) = \frac{E(s)}{R(s)} R(s) = \frac{s^2(s + K_2)}{s^3 + K_2 s^2 + (K_1+4)s + K_1} \cdot \frac{1}{s^2} = \frac{s + K_2}{s^3 + K_2 s^2 + (K_1+4)s + K_1}

系が安定である(K2(K1+4)>K1K_2(K_1+4) > K_1)とき、最終値の定理を適用して定常偏差 esse_{\mathrm{ss}} を求めます。

ess=lim⁡t→∞e(t)=lim⁡s→0sE(s)=lim⁡s→0[s⋅s+K2s3+K2s2+(K1+4)s+K1]=0⋅K2K1=0e_{\mathrm{ss}} = \lim_{t \to \infty} e(t) = \lim_{s \to 0} s E(s) = \lim_{s \to 0} \left[ s \cdot \frac{s + K_2}{s^3 + K_2 s^2 + (K_1+4)s + K_1} \right] = 0 \cdot \frac{K_2}{K_1} = 0

00

解答のまとめ

  • (1) 制御対象のインパルス応答:g(t)=12sin⁡2t(t≥0)g(t) = \displaystyle\frac{1}{2}\sin 2t \quad (t \ge 0)
  • (2) 伝達関数:E(s)R(s)=s3+K2s2+4ss3+K2s2+(K1+4)s+K1\displaystyle\frac{E(s)}{R(s)} = \frac{s^3 + K_2 s^2 + 4s}{s^3 + K_2 s^2 + (K_1 + 4)s + K_1}
  • (3) 安定限界条件と応答:K2(K1+4)=K1K_2(K_1 + 4) = K_1、応答は減衰しない持続振動となる。
  • (4) 空欄:(a) F(s)F(s)、(b) C1(s)C_1(s)、(c) C2(s)G(s)C_2(s)G(s)、(d) G(s)G(s)、(e) C2(s)G(s)−F(s)G(s)C_2(s)G(s) - F(s)G(s)、(f) C1(s)G(s)+C2(s)G(s)C_1(s)G(s) + C_2(s)G(s)
  • (5) ランプ入力に対する定常偏差:00