過去問一覧第二種電気主任技術者2024年度 機械・制御直並列制御系の伝達関数と時間応答・周波数応答の導出
一覧に戻る

直並列制御系の伝達関数と時間応答・周波数応答の導出

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

伝達関数G1(s) G_{1} (s) ,G2(s) G_{2} (s) に対して図1のように直列結合した制御系と図2のように並列結合した制御系が与えられている。U(s) U (s) は入力量,Y1(s) Y_{1} (s) は直列結合した場合の出力量,Y2(s) Y_{2} (s) は並列結合した場合の出力量をあらわしている。U(s) U (s) ,Y1(s) Y_{1} (s) ,Y2(s) Y_{2} (s) は時間信号u(t) u (t) ,y1(t) y_{1} (t) ,y2(t) y_{2} (t) をそれぞれラプラス変換したものである。

伝達関数G1(s) G_{1} (s) ,G2(s) G_{2} (s) が,
G1(s)=3s+2G2(s)=5s2+5s+6\begin{aligned}G_{1} (s) &=\frac {3}{s+2} \\G_{2} (s) &=\frac {5}{s^{2}+5s+6} \end{aligned} のように与えられているとき,以下の問に答えよ。

問題図
(1) 図1の直列結合された制御系のU(s) U (s) からY1(s) Y_{1} (s) の伝達関数を求めよ。なお,伝達関数は,一つの有理関数で表すとし,分母及び分子はs s の多項式で示すこと。
(2) (1)で求めた直列結合された伝達関数に単位ステップ信号を加えたときの出力量Y1(s) Y_{1} (s) の定常値を求めよ。
(3) 図2の並列結合された制御系のインパルス信号を加えたときのY2(s) Y_{2} (s) の時間応答を求めよ。
(4) (3)で求めた並列結合された伝達関数にu(t)=2e−t u(t)=2\mathrm {e}^{-t} を加えたときのY2(s) Y_{2} (s) の時間応答を求めよ。
(5) (3)で求めた並列結合された伝達関数の周波数応答において,周波数を十分に大きくしたときの位相を求めよ。

解説

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

本問は、古典制御理論の根幹をなすブロック線図の結合(直列・並列)、最終値の定理、逆ラプラス変換による時間応答解析、および周波数伝達関数の位相特性を問う総合問題です。

  • 直列結合と並列結合: 図1の直列結合伝達関数は積 W1(s)=G1(s)G2(s)W_1(s) = G_1(s)G_2(s)、 図2の並列結合伝達関数は和 W2(s)=G1(s)+G2(s)W_2(s) = G_1(s) + G_2(s) となります。
  • 最終値の定理: 系が安定(sY(s)sY(s) の極の実部がすべて負)であるとき、出力の定常値は次式で与えられます。 lim⁡t→∞y(t)=lim⁡s→0sY(s)\lim_{t \to \infty} y(t) = \lim_{s \to 0} s Y(s)
  • ラプラス変換の対応関係: L[δ(t)]=1,L[u(t)]=1s,L[e−at]=1s+a\mathcal{L}[\delta(t)] = 1, \quad \mathcal{L}[u(t)] = \frac{1}{s}, \quad \mathcal{L}[\mathrm{e}^{-at}] = \frac{1}{s+a}
  • 部分分数分解(ヘビサイドの展開定理): 単極のみをもつ伝達関数 F(s)=P(s)(s−p1)(s−p2)…F(s) = \frac{P(s)}{(s-p_1)(s-p_2)\dots} の各係数は、留数法 Ck=lim⁡s→pk(s−pk)F(s)C_k = \lim_{s \to p_k} (s-p_k) F(s) により短時間で計算可能です。

小問(1)の解説:直列結合された伝達関数の導出

図1

問題図
において、G1(s)G_1(s) の出力が直接 G2(s)G_2(s) の入力となっているため、全体の伝達関数 W1(s)=Y1(s)U(s)\displaystyle W_1(s) = \frac{Y_1(s)}{U(s)} は各ブロックの積となります。

W1(s)=G1(s)G2(s)=3s+2⋅5s2+5s+6=15(s+2)(s2+5s+6)W_1(s) = G_1(s)G_2(s) = \frac{3}{s+2} \cdot \frac{5}{s^2+5s+6} = \frac{15}{(s+2)(s^2+5s+6)}

問題文の指定「分母及び分子は ss の多項式で示すこと」に従い、分母を展開します。

(s+2)(s2+5s+6)=s(s2+5s+6)+2(s2+5s+6)=s3+5s2+6s+2s2+10s+12=s3+7s2+16s+12\begin{aligned} (s+2)(s^2+5s+6) &= s(s^2+5s+6) + 2(s^2+5s+6) \\[5pt] &= s^3 + 5s^2 + 6s + 2s^2 + 10s + 12 \\[5pt] &= s^3 + 7s^2 + 16s + 12 \end{aligned}

したがって、求める伝達関数は次のようになります。

Y1(s)U(s)=15s3+7s2+16s+12\frac{Y_1(s)}{U(s)} = \frac{15}{s^3+7s^2+16s+12}

小問(2)の解説:単位ステップ応答の定常値

単位ステップ信号のラプラス変換は U(s)=1sU(s) = \displaystyle \frac{1}{s} であるため、出力量 Y1(s)Y_1(s) は次式となります。

Y1(s)=W1(s)U(s)=15s3+7s2+16s+12⋅1s=15(s+2)2(s+3)sY_1(s) = W_1(s) U(s) = \frac{15}{s^3+7s^2+16s+12} \cdot \frac{1}{s} = \frac{15}{(s+2)^2(s+3)s}

sY1(s)sY_1(s) の極は s=−2s = -2(重極)および s=−3s = -3 であり、すべて実部が負(複素平面の左半面)に存在するため、系は安定であり最終値の定理が適用できます。

lim⁡t→∞y1(t)=lim⁡s→0sY1(s)=lim⁡s→0(s⋅15s3+7s2+16s+12⋅1s)=lim⁡s→015s3+7s2+16s+12=1512=54=1.25\begin{aligned} \lim_{t \to \infty} y_1(t) &= \lim_{s \to 0} s Y_1(s) \\[5pt] &= \lim_{s \to 0} \left( s \cdot \frac{15}{s^3+7s^2+16s+12} \cdot \frac{1}{s} \right) \\[5pt] &= \lim_{s \to 0} \frac{15}{s^3+7s^2+16s+12} \\[5pt] &= \frac{15}{12} = \frac{5}{4} = 1.25 \end{aligned}

1.25(または 54)1.25 \quad \left( \text{または } \frac{5}{4} \right)

小問(3)の解説:並列結合のインパルス応答

図2

問題図
より、並列結合された制御系の伝達関数 W2(s)W_2(s) は各要素の和になります。

W2(s)=Y2(s)U(s)=G1(s)+G2(s)=3s+2+5s2+5s+6W_2(s) = \frac{Y_2(s)}{U(s)} = G_1(s) + G_2(s) = \frac{3}{s+2} + \frac{5}{s^2+5s+6}

ここで分母を因数分解すると s2+5s+6=(s+2)(s+3)s^2+5s+6 = (s+2)(s+3) です。通分して整理します。

W2(s)=3(s+3)(s+2)(s+3)+5(s+2)(s+3)=3s+14(s+2)(s+3)W_2(s) = \frac{3(s+3)}{(s+2)(s+3)} + \frac{5}{(s+2)(s+3)} = \frac{3s+14}{(s+2)(s+3)}

入力が単位インパルス信号であるとき U(s)=1U(s) = 1 であるため、出力量は Y2(s)=W2(s)Y_2(s) = W_2(s) と一致します。これを逆ラプラス変換するために部分分数分解します。

Y2(s)=3s+14(s+2)(s+3)=As+2+Bs+3Y_2(s) = \frac{3s+14}{(s+2)(s+3)} = \frac{A}{s+2} + \frac{B}{s+3}

ヘビサイドの展開定理(留数法)を用いると、係数 A,BA, B は直ちに求まります。

A=(s+2)Y2(s)∣s=−2=3s+14s+3∣s=−2=3(−2)+14−2+3=8B=(s+3)Y2(s)∣s=−3=3s+14s+2∣s=−3=3(−3)+14−3+2=−5\begin{aligned} A &= \left. (s+2) Y_2(s) \right|_{s=-2} = \left. \frac{3s+14}{s+3} \right|_{s=-2} = \frac{3(-2)+14}{-2+3} = 8 \\[5pt] B &= \left. (s+3) Y_2(s) \right|_{s=-3} = \left. \frac{3s+14}{s+2} \right|_{s=-3} = \frac{3(-3)+14}{-3+2} = -5 \end{aligned}

したがって、Y2(s)=8s+2−5s+3Y_2(s) = \displaystyle \frac{8}{s+2} - \frac{5}{s+3} となります。逆ラプラス変換を行うと、時間応答 y2(t)y_2(t) は次のようになります。

y2(t)=8e−2t−5e−3t(t≥0)y_2(t) = 8\mathrm{e}^{-2t} - 5\mathrm{e}^{-3t} \quad (t \ge 0)

【別解(スマートな部分分数分解)】
G2(s)=5(s+2)(s+3)=5s+2−5s+3G_2(s) = \frac{5}{(s+2)(s+3)} = \frac{5}{s+2} - \frac{5}{s+3} と先に分解しておけば、 W2(s)=G1(s)+G2(s)=3s+2+(5s+2−5s+3)=8s+2−5s+3W_2(s) = G_1(s) + G_2(s) = \frac{3}{s+2} + \left(\frac{5}{s+2} - \frac{5}{s+3}\right) = \frac{8}{s+2} - \frac{5}{s+3} と通分することなく瞬時に分解形を得ることができます。

小問(4)の解説:指数関数入力に対する時間応答

入力 u(t)=2e−tu(t) = 2\mathrm{e}^{-t} のラプラス変換は U(s)=2s+1U(s) = \displaystyle \frac{2}{s+1} です。小問(3)で求めた W2(s)W_2(s) に乗じることで、出力 Y2(s)Y_2(s) を導きます。

Y2(s)=W2(s)U(s)=3s+14(s+2)(s+3)⋅2s+1=6s+28(s+1)(s+2)(s+3)Y_2(s) = W_2(s) U(s) = \frac{3s+14}{(s+2)(s+3)} \cdot \frac{2}{s+1} = \frac{6s+28}{(s+1)(s+2)(s+3)}

これを部分分数に展開します。

Y2(s)=Cs+1+Ds+2+Es+3Y_2(s) = \frac{C}{s+1} + \frac{D}{s+2} + \frac{E}{s+3}

ヘビサイドの展開定理を用いて各定数を求めます。

C=6s+28(s+2)(s+3)∣s=−1=6(−1)+28(−1+2)(−1+3)=221⋅2=11D=6s+28(s+1)(s+3)∣s=−2=6(−2)+28(−2+1)(−2+3)=16(−1)⋅1=−16E=6s+28(s+1)(s+2)∣s=−3=6(−3)+28(−3+1)(−3+2)=10(−2)⋅(−1)=5\begin{aligned} C &= \left. \frac{6s+28}{(s+2)(s+3)} \right|_{s=-1} = \frac{6(-1)+28}{(-1+2)(-1+3)} = \frac{22}{1 \cdot 2} = 11 \\[5pt] D &= \left. \frac{6s+28}{(s+1)(s+3)} \right|_{s=-2} = \frac{6(-2)+28}{(-2+1)(-2+3)} = \frac{16}{(-1) \cdot 1} = -16 \\[5pt] E &= \left. \frac{6s+28}{(s+1)(s+2)} \right|_{s=-3} = \frac{6(-3)+28}{(-3+1)(-3+2)} = \frac{10}{(-2) \cdot (-1)} = 5 \end{aligned}

よって、Y2(s)Y_2(s) は以下のように分解されます。

Y2(s)=11s+1−16s+2+5s+3Y_2(s) = \frac{11}{s+1} - \frac{16}{s+2} + \frac{5}{s+3}

これを逆ラプラス変換することにより、時間応答 y2(t)y_2(t) が求まります。

y2(t)=11e−t−16e−2t+5e−3t(t≥0)y_2(t) = 11\mathrm{e}^{-t} - 16\mathrm{e}^{-2t} + 5\mathrm{e}^{-3t} \quad (t \ge 0)

※検算: t=0t=0 を代入すると y2(0)=11−16+5=0y_2(0) = 11 - 16 + 5 = 0 となり、初期値定理 lim⁡s→∞sY2(s)=0\lim_{s\to\infty} sY_2(s) = 0 と正しく一致します。

小問(5)の解説:高周波領域における周波数応答の位相

小問(3)で求めた並列結合伝達関数 W2(s)=3s+14(s+2)(s+3)W_2(s) = \displaystyle \frac{3s+14}{(s+2)(s+3)} に s=jωs = \mathrm{j}\omega を代入すると、周波数伝達関数 W2(jω)W_2(\mathrm{j}\omega) が得られます。

W2(jω)=3jω+14(2+jω)(3+jω)W_2(\mathrm{j}\omega) = \frac{3\mathrm{j}\omega + 14}{(2+\mathrm{j}\omega)(3+\mathrm{j}\omega)}

複素数の偏角(位相)の性質 arg⁡(ABC)=arg⁡(A)−arg⁡(B)−arg⁡(C)\arg\left(\frac{A}{BC}\right) = \arg(A) - \arg(B) - \arg(C) より、位相 ϕ(ω)\phi(\omega) は次のように表せます。

ϕ(ω)=tan⁡−1(3ω14)−tan⁡−1(ω2)−tan⁡−1(ω3)\phi(\omega) = \tan^{-1}\left(\frac{3\omega}{14}\right) - \tan^{-1}\left(\frac{\omega}{2}\right) - \tan^{-1}\left(\frac{\omega}{3}\right)

ここで、周波数を十分に大きくする(ω→∞\omega \to \infty)と、各逆正接項はすべて tan⁡−1(∞)=90∘\tan^{-1}(\infty) = 90^\circ(π2 rad\frac{\pi}{2}\,\mathrm{rad})に収束します。

lim⁡ω→∞ϕ(ω)=90∘−90∘−90∘=−90∘\lim_{\omega \to \infty} \phi(\omega) = 90^\circ - 90^\circ - 90^\circ = -90^\circ

(なお、高周波領域において分母・分子の最高次数の項に注目すると、 W2(jω)≈3jω(jω)2=3jω=−j3ωW_2(\mathrm{j}\omega) \approx \frac{3\mathrm{j}\omega}{(\mathrm{j}\omega)^2} = \frac{3}{\mathrm{j}\omega} = -\mathrm{j}\frac{3}{\omega} となり、純虚数かつ係数が負(−j-\mathrm{j})であることから、位相が −90∘-90^\circ となることが視覚的にも容易に確認できます。)

−90∘(または −π2 rad)-90^\circ \quad \left( \text{または } -\frac{\pi}{2} \ \mathrm{rad} \right)

解答のまとめ

  • (1) Y1(s)U(s)=15s3+7s2+16s+12\displaystyle \frac{Y_1(s)}{U(s)} = \frac{15}{s^3+7s^2+16s+12}
  • (2) 1.25(54)\displaystyle 1.25 \quad \left( \frac{5}{4} \right)
  • (3) y2(t)=8e−2t−5e−3t(t≥0)\displaystyle y_2(t) = 8\mathrm{e}^{-2t} - 5\mathrm{e}^{-3t} \quad (t \ge 0)
  • (4) y2(t)=11e−t−16e−2t+5e−3t(t≥0)\displaystyle y_2(t) = 11\mathrm{e}^{-t} - 16\mathrm{e}^{-2t} + 5\mathrm{e}^{-3t} \quad (t \ge 0)
  • (5) −90∘(−π2 rad)\displaystyle -90^\circ \quad \left( -\frac{\pi}{2} \ \mathrm{rad} \right)