過去問一覧第二種電気主任技術者2023年度 機械・制御ボード線図のゲイン特性曲線からの伝達関数の同定
一覧に戻る

ボード線図のゲイン特性曲線からの伝達関数の同定

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

本問題で扱う伝達関数の全ての極と零点は複素平面上の右半平面には存在しないものと仮定する。以下の問に答えよ。ここで必要に応じて,14=20log⁡10100.7 14=20 \log _{10} 10^{0.7} ,100.7=5.0119 10^{0.7}=5.0119 を用いよ。ただし,全ての図は折れ線近似で表している。

(1) 図1に示すゲイン特性曲線から積分要素の伝達関数G1(s)=1TIs \displaystyle G_{1}\left( s\right) =\frac {1}{T_{\mathrm {I}}s} を求めよ。
問題図
(2) 図2に示すゲイン特性曲線が表す伝達関数G(s) G \left( s\right) をG(s)=G1(s)G2(s) G \left( s\right) =G_{1} \left( s\right) G_{2} \left( s\right) のように分解して考える。G1(s) G_{1} \left( s\right) を小問(1)で求めた伝達関数とするときのG2(s) G_{2} \left( s\right) を求めよ。
問題図
(3) 図2に示すゲイン特性曲線から伝達関数G(s) G \left( s\right) を求めよ。ただし,G(s)=1TI1s \displaystyle G\left( s\right) =\frac {1}{T_{\mathrm {I1}}s} とおいて,その周波数伝達関数G(jω)=1jTI1ω \displaystyle G\left( \mathrm {j}\omega \right) =\frac {1}{\mathrm {j}T_{\mathrm {I1}}\omega } のゲインがω=0.1 rad/s \omega = 0.1 \ \mathrm {rad / s} のとき40 dB 40 \ \mathrm {dB} であることを用いてTI1 T_{\mathrm {I1}} を決定せよ。
(4) 折れ線近似で表した図3に示すゲイン特性曲線から伝達関数G(s) G \left( s\right) を求めよ。ただし,G(s) G \left( s\right) を積分要素G1(s)=1TI2s \displaystyle G_{1}\left( s\right) =\frac {1}{T_{\mathrm {I2}}s} と一次遅れ要素G2(s)=11+Ts \displaystyle G_{2}\left( s\right) =\frac {1}{1+Ts} に分解して考えよ。
問題図
(5)

解説

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

本問は、折れ線近似されたボード線図(ゲイン特性曲線)から、対応する各要素(積分要素・比例要素・一次遅れ要素)のパラメータを同定し、伝達関数を導出する問題です。

  • 周波数伝達関数とゲインの定義:
    伝達関数 G(s)G(s) に対し s=jωs = \mathrm{j}\omega とおいた周波数伝達関数 G(jω)G(\mathrm{j}\omega) のゲイン g(ω) [dB]g(\omega)\ [\mathrm{dB}] は、 g(ω)=20log⁡10∣G(jω)∣g(\omega) = 20 \log_{10} |G(\mathrm{j}\omega)|
  • 積分要素 G(s)=1TIsG(s) = \dfrac{1}{T_{\mathrm{I}}s} のゲイン特性:
    g(ω)=20log⁡10∣1jωTI∣=−20log⁡10(ωTI)g(\omega) = 20 \log_{10} \left|\frac{1}{\mathrm{j}\omega T_{\mathrm{I}}}\right| = -20 \log_{10} (\omega T_{\mathrm{I}}) 角周波数 ω\omega が10倍(1 dec1\ \mathrm{dec})になるとゲインは −20 dB-20\ \mathrm{dB} 変化するため、全周波数域で傾き −20 dB/dec-20\ \mathrm{dB/dec} の直線となります。また、g(ω)=0 dBg(\omega) = 0\ \mathrm{dB} となる交点角周波数は ω=1TI\omega = \dfrac{1}{T_{\mathrm{I}}} です。
  • 一次遅れ要素 G(s)=11+TsG(s) = \dfrac{1}{1 + Ts} の折れ線近似:
    折れ点角周波数を ωc=1T\omega_{\mathrm{c}} = \dfrac{1}{T} とすると、
    • ω≤ωc\omega \le \omega_{\mathrm{c}} の低周波域:ゲイン 0 dB0\ \mathrm{dB}(傾き 0 dB/dec0\ \mathrm{dB/dec})
    • ω≥ωc\omega \ge \omega_{\mathrm{c}} の高周波域:傾き −20 dB/dec-20\ \mathrm{dB/dec} で直線的に減衰
  • 伝達関数の積とゲインの加算性:
    伝達関数が G(s)=G1(s)G2(s)G(s) = G_1(s)G_2(s) と表されるとき、その合成ゲインは各要素のゲインの和になります: 20log⁡10∣G(jω)∣=20log⁡10∣G1(jω)∣+20log⁡10∣G2(jω)∣20\log_{10}|G(\mathrm{j}\omega)| = 20\log_{10}|G_1(\mathrm{j}\omega)| + 20\log_{10}|G_2(\mathrm{j}\omega)|

【小問(1)の解説】

図1から積分要素 G1(s)=1TIsG_1(s) = \dfrac{1}{T_{\mathrm{I}}s} を決定する

周波数伝達関数は G1(jω)=1jωTIG_1(\mathrm{j}\omega) = \dfrac{1}{\mathrm{j}\omega T_{\mathrm{I}}} であり、そのゲイン g1(ω)g_1(\omega) は次式となります。

g1(ω)=20log⁡10∣1jωTI∣=−20log⁡10(ωTI) [dB]g_1(\omega) = 20 \log_{10} \left|\frac{1}{\mathrm{j}\omega T_{\mathrm{I}}}\right| = -20 \log_{10}(\omega T_{\mathrm{I}})\ [\mathrm{dB}]

問題図
より、特性曲線は ω=0.1 rad/s\omega = 0.1\ \mathrm{rad/s} のとき g1=0 dBg_1 = 0\ \mathrm{dB}(また ω=0.01\omega = 0.01 で 20 dB20\ \mathrm{dB}、ω=1\omega = 1 で −20 dB-20\ \mathrm{dB})を通っています。ω=0.1 rad/s\omega = 0.1\ \mathrm{rad/s} で g1=0 dBg_1 = 0\ \mathrm{dB} となる条件を代入すると、

−20log⁡10(0.1TI)=0-20 \log_{10}(0.1 T_{\mathrm{I}}) = 0 log⁡10(0.1TI)=0\log_{10}(0.1 T_{\mathrm{I}}) = 0 0.1TI=100=1  ⟹  TI=10 s0.1 T_{\mathrm{I}} = 10^0 = 1 \implies T_{\mathrm{I}} = 10\ \mathrm{s}

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

G1(s)=110sG_1(s) = \frac{1}{10s}

【小問(2)の解説】

G(s)=G1(s)G2(s)G(s) = G_1(s)G_2(s) と分解したときの G2(s)G_2(s) を求める

問題図
(図2)のゲイン特性は、図1と同じ傾き −20 dB/dec-20\ \mathrm{dB/dec} の直線です。同じ角周波数 ω=0.1 rad/s\omega = 0.1\ \mathrm{rad/s} におけるゲインを比較すると、

  • 図1のゲイン:0 dB0\ \mathrm{dB}
  • 図2のゲイン:40 dB40\ \mathrm{dB}

傾きが等しく一定の差を保っていることから、図2の特性は図1の特性を周波数全域にわたり +40 dB+40\ \mathrm{dB} だけ平行移動(上方シフト)させたものであることが分かります。

ゲインの加算性より、周波数伝達関数 G2(jω)G_2(\mathrm{j}\omega) のゲイン g2(ω)g_2(\omega) は、

g2(ω)=20log⁡10∣G(jω)∣−20log⁡10∣G1(jω)∣=40−0=40 [dB]g_2(\omega) = 20 \log_{10}|G(\mathrm{j}\omega)| - 20 \log_{10}|G_1(\mathrm{j}\omega)| = 40 - 0 = 40\ [\mathrm{dB}]

これは角周波数 ω\omega に依存しない一定値であるため、G2(s)G_2(s) は周波数によらない純粋な比例要素(定数ゲイン KK)です。

20log⁡10∣G2(jω)∣=40  ⟹  log⁡10∣G2(jω)∣=2  ⟹  ∣G2(jω)∣=102=10020 \log_{10}|G_2(\mathrm{j}\omega)| = 40 \implies \log_{10}|G_2(\mathrm{j}\omega)| = 2 \implies |G_2(\mathrm{j}\omega)| = 10^2 = 100

問題の仮定より極や零点は右半平面に存在しないため、G2(s)G_2(s) は正の実数定数となります。

G2(s)=100G_2(s) = 100

【小問(3)の解説】

図2のゲイン特性から直接 G(s)=1TI1sG(s) = \dfrac{1}{T_{\mathrm{I1}}s} を求める

周波数伝達関数 G(jω)=1jωTI1G(\mathrm{j}\omega) = \dfrac{1}{\mathrm{j}\omega T_{\mathrm{I1}}} のゲインは、

g(ω)=20log⁡10∣1jωTI1∣=−20log⁡10(ωTI1) [dB]g(\omega) = 20 \log_{10} \left|\frac{1}{\mathrm{j}\omega T_{\mathrm{I1}}}\right| = -20 \log_{10}(\omega T_{\mathrm{I1}})\ [\mathrm{dB}]

題意より、ω=0.1 rad/s\omega = 0.1\ \mathrm{rad/s} のとき g=40 dBg = 40\ \mathrm{dB} であるため、

−20log⁡10(0.1TI1)=40-20 \log_{10}(0.1 T_{\mathrm{I1}}) = 40 log⁡10(0.1TI1)=−2\log_{10}(0.1 T_{\mathrm{I1}}) = -2 0.1TI1=10−2=0.01  ⟹  TI1=0.1 s0.1 T_{\mathrm{I1}} = 10^{-2} = 0.01 \implies T_{\mathrm{I1}} = 0.1\ \mathrm{s}

したがって、G(s)G(s) は次のように求められます。

TI1=0.1 s,G(s)=10.1s=10sT_{\mathrm{I1}} = 0.1\ \mathrm{s},\quad G(s) = \frac{1}{0.1s} = \frac{10}{s}

※ 小問(1)と小問(2)の結果を掛け合わせると G(s)=G1(s)G2(s)=110s×100=10sG(s) = G_1(s)G_2(s) = \dfrac{1}{10s} \times 100 = \dfrac{10}{s} となり、完全に一致することが確認できます。

【小問(4)の解説】

図3の折れ線近似特性から G(s)G(s) を求める

伝達関数を積分要素 G1(s)=1TI2sG_1(s) = \dfrac{1}{T_{\mathrm{I2}}s} と一次遅れ要素 G2(s)=11+TsG_2(s) = \dfrac{1}{1+Ts} の直列結合:

G(s)=G1(s)G2(s)=1TI2s(1+Ts)G(s) = G_1(s)G_2(s) = \frac{1}{T_{\mathrm{I2}}s(1+Ts)}

として分解します。

問題図
(図3)の特性曲線を分析して各パラメータを決定します。

1. 一次遅れ要素の時定数 TT の決定

図3において、ゲインの傾きは ω=0.5 rad/s\omega = 0.5\ \mathrm{rad/s} の地点で −20 dB/dec-20\ \mathrm{dB/dec} から −40 dB/dec-40\ \mathrm{dB/dec} へと急変しています。傾きが −20 dB/dec-20\ \mathrm{dB/dec} だけ減少していることから、この折れ点角周波数が一次遅れ要素の折れ点 ωc=0.5 rad/s\omega_{\mathrm{c}} = 0.5\ \mathrm{rad/s} に対応します。

ωc=1T=0.5 rad/s  ⟹  T=10.5=2 s\omega_{\mathrm{c}} = \frac{1}{T} = 0.5\ \mathrm{rad/s} \implies T = \frac{1}{0.5} = 2\ \mathrm{s}

したがって、一次遅れ要素は G2(s)=11+2sG_2(s) = \dfrac{1}{1+2s} と決定されます。

2. 積分要素の時定数 TI2T_{\mathrm{I2}} の決定

折れ線近似では、一次遅れ要素 G2(s)G_2(s) のゲインは折れ点角周波数以下(ω≤0.5 rad/s\omega \le 0.5\ \mathrm{rad/s})において 0 dB0\ \mathrm{dB} とみなされます。

したがって、折れ点 ω=0.5 rad/s\omega = 0.5\ \mathrm{rad/s} における全体のゲイン 14 dB14\ \mathrm{dB} は、すべて積分要素 G1(s)G_1(s) が担っていることになります。

g1(0.5)=−20log⁡10(0.5TI2)=14 [dB]g_1(0.5) = -20\log_{10}(0.5 T_{\mathrm{I2}}) = 14\ [\mathrm{dB}]

問題文に与えられた条件式 14=20log⁡10100.714 = 20 \log_{10} 10^{0.7} および 100.7=5.011910^{0.7} = 5.0119 を用いて変形します。

20log⁡10(10.5TI2)=20log⁡10100.720 \log_{10}\left(\frac{1}{0.5 T_{\mathrm{I2}}}\right) = 20 \log_{10} 10^{0.7} 10.5TI2=100.7=5.0119\frac{1}{0.5 T_{\mathrm{I2}}} = 10^{0.7} = 5.0119 TI2=10.5×5.0119=25.0119 [s]T_{\mathrm{I2}} = \frac{1}{0.5 \times 5.0119} = \frac{2}{5.0119}\ [\mathrm{s}]

これより、積分要素は G1(s)=1TI2s=5.01192sG_1(s) = \dfrac{1}{T_{\mathrm{I2}}s} = \dfrac{5.0119}{2s} と求まります。

3. 伝達関数 G(s)G(s) の合成

以上を掛け合わせると、伝達関数 G(s)G(s) は次のようになります。

G(s)=G1(s)G2(s)=5.01192s⋅11+2s=5.01192s(1+2s)≒5.012s(1+2s)G(s) = G_1(s)G_2(s) = \frac{5.0119}{2s} \cdot \frac{1}{1+2s} = \frac{5.0119}{2s(1+2s)} \fallingdotseq \frac{5.01}{2s(1+2s)}

なお、分子・分母を整理して分子に比例ゲインをまとめた表記(5.01192≈2.506\dfrac{5.0119}{2} \approx 2.506)も同等に正解です。

G(s)=5.012s(1+2s)(≈2.51s(1+2s))G(s) = \frac{5.01}{2s(1+2s)}\quad \left(\approx \frac{2.51}{s(1+2s)}\right)

【正解の総括】

  • (1) G1(s)=110sG_1(s) = \dfrac{1}{10s}
  • (2) G2(s)=100G_2(s) = 100
  • (3) TI1=0.1 s,G(s)=10sT_{\mathrm{I1}} = 0.1\ \mathrm{s},\quad G(s) = \dfrac{10}{s}
  • (4) G(s)=5.012s(1+2s)(または 2.51s(1+2s))G(s) = \dfrac{5.01}{2s(1+2s)}\quad \left(\text{または } \dfrac{2.51}{s(1+2s)}\right)

本問のように「ボード線図から伝達関数を逆算する」問題では、①直線の傾きから要素の次数(積分・一次遅れなど)を見極めること、②折れ点角周波数から各要素の時定数 TT を特定すること、③折れ線近似における基準点(ωc\omega_c において一次遅れのゲインは 0 dB0\ \mathrm{dB})を利用して積分時定数やゲイン定数を決定すること、の3ステップを着実に実行することが重要です。