跳到主要内容

京都大学 情報学研究科 システム科学専攻 2024年8月実施 専門科目 制御工学

Author

AKIRA (小红书:94184092292), 祭音Myyura

Description

大学公表の原題

問題1

図1のフィードバック制御系に関する以下の設問に答えよ。ただし,

P(s)=1s2+2s8,F(s)=as+1P(s) = \frac{1}{s^2 + 2s - 8}, \quad F(s) = \frac{a}{s + 1}

とし,a,Ka, K は定数パラメータとする。

(1) P(s)P(s) のステップ応答を求めよ。

(2) rr から yy への伝達関数 Gyr(s)G_{yr}(s)dd から yy への伝達関数 Gyd(s)G_{yd}(s) をそれぞれ求めよ。

(3) Gyr(s),Gyd(s)G_{yr}(s), G_{yd}(s) がともに安定なとき,cc を定数パラメータとして,ステップ入力 d(t)=cd(t) = c と単位ステップ入力 r(t)=1r(t) = 1 を同時に加えると,出力 y(t)y(t) の定常値 limty(t)\lim_{t \to \infty} y(t)22 であった。このとき,cc の値を a,Ka, K を用いて表せ。

(4) a,Ka, K を実数の範囲で自由に選べるとして,Gyr(s)G_{yr}(s) のすべての極の実部が正になる a,Ka, K が存在するか否か,理由とともに述べよ。

問題2

図2に示す内部安定なフィードバック制御系について、以下の設問に答えよ。 なお、制御対象 P(s)P(s) は最小位相で、その定常ゲイン P(0)P(0) は正の定数とする。 また、C(s)C(s) は次式で与えられる PID 制御器であり、比例ゲイン KPK_P、積分時間 TIT_I、微分時間 TDT_D はすべて正の定数とする。

C(s)=KP(1+1TIs+TDs)C(s) = K_P \left( 1 + \frac{1}{T_I s} + T_D s \right)

(1) PID 制御器 C(s)C(s) を設計したところ、一巡伝達関数のボード線図が図3のようになった。ゲイン余裕および位相余裕を図から読みとって答えよ。

(2) 設問 (1) の PID 制御器の比例ゲイン KPK_P のみを調整して、ゲイン余裕を 10 dB だけ大きくしたい。KPK_P を何倍にすればよいか答えよ。

(3) 設問 (1) の PID 制御器を用いたところ、設定値 rr をステップ状に変更した直後に操作量 uu が急激に変動することが判明した。 この変動を抑制するために、図4に示す I-PD 制御系を採用することにした。 ただし、KPK_P, TIT_I, TDT_D の値は設問 (1) の PID 制御器と同じ値とする。 PID 制御と比較して、操作量の急激な変動を抑制するために I-PD 制御が有効であると期待できる理由を簡潔に説明せよ。

(4) 図4の I-PD 制御系について、外乱 dd から制御量 yy への閉ループ伝達関数を求めよ。さらに、PID 制御と I-PD 制御で同じ制御パラメータ KPK_P, TIT_I, TDT_D を用いる場合について、求めた閉ループ伝達関数が、図2の PID 制御系の外乱 dd から制御量 yy への閉ループ伝達関数と同じになるかどうかを理由とともに答えよ。

题目描述

  1. 对图 1 所示反馈控制系统,设

    P(s)=1s2+2s8,F(s)=as+1,P(s)=\frac1{s^2+2s-8},\qquad F(s)=\frac a{s+1},

    其中 a,Ka,K 为常数参数。

    (1)求 P(s)P(s) 的阶跃响应。

    (2)分别求从 rryy 的传递函数 Gyr(s)G_{yr}(s),以及从扰动 ddyy 的传递函数 Gyd(s)G_{yd}(s)

    (3)假设 Gyr(s),Gyd(s)G_{yr}(s),G_{yd}(s) 都稳定。同时施加阶跃扰动 d(t)=cd(t)=c 和单位阶跃参考输入 r(t)=1r(t)=1,其中 cc 为常数参数;若输出稳态值满足

    limty(t)=2,\lim_{t\to\infty}y(t)=2,

    a,Ka,K 表示 cc

    (4)允许在实数范围内自由选择 a,Ka,K。判断是否存在某组 a,Ka,K,使 Gyr(s)G_{yr}(s) 的所有极点实部都为正,并说明理由。

  2. 考虑图 2 所示内部稳定的反馈控制系统。对象 P(s)P(s) 为最小相位系统,且其稳态增益 P(0)P(0) 为正常数。控制器为 PID 控制器

    C(s)=KP(1+1TIs+TDs),C(s)=K_P\left(1+\frac1{T_Is}+T_Ds\right),

    其中比例增益 KPK_P、积分时间 TIT_I、微分时间 TDT_D 均为正常数。

    (1)设计 PID 控制器后,开环传递函数的 Bode 图如图 3。由图读出增益裕度和相位裕度。

    (2)仅调整第(1)问 PID 控制器的比例增益 KPK_P,使增益裕度增加 10dB10\,\mathrm{dB}。求 KPK_P 应变为原来的多少倍。

    (3)使用第(1)问的 PID 控制器时,发现参考值 rr 发生阶跃变化后,操作量 uu 会立即剧烈变化。为抑制该变化,改用图 4 的 I-PD 控制结构,且 KP,TI,TDK_P,T_I,T_D 与原 PID 控制器相同。简要说明相较 PID 控制,为什么 I-PD 有望抑制操作量的突然变化。

    (4)对图 4 的 I-PD 系统,求从扰动 dd 到被控量 yy 的闭环传递函数。若 PID 与 I-PD 使用相同的 KP,TI,TDK_P,T_I,T_D,判断该传递函数是否与图 2 PID 系统从 ddyy 的闭环传递函数相同,并说明理由。

Kai

問題1

(1)

P(s)s=1s(s+4)(s2)=18s+124(s+4)+112(s2).\frac{P(s)}s=\frac1{s(s+4)(s-2)} =-\frac1{8s}+\frac1{24(s+4)}+\frac1{12(s-2)}.

従ってステップ応答は

y(t)=18+124e4t+112e2t,t0.\boxed{y(t)=-\frac18+\frac1{24}e^{-4t}+\frac1{12}e^{2t}},\qquad t\ge0.

(2)

ブロック線図より y=(F+KP)(ry)+dy=(F+KP)(r-y)+d となる。従って

Gyr(s)=F(s)+KP(s)1+F(s)+KP(s),Gyd(s)=11+F(s)+KP(s).\boxed{G_{yr}(s)=\frac{F(s)+KP(s)}{1+F(s)+KP(s)},\qquad G_{yd}(s)=\frac1{1+F(s)+KP(s)}}.

Δ(s)=(s+1)(s2+2s8)+a(s2+2s8)+K(s+1)\Delta(s)=(s+1)(s^2+2s-8)+a(s^2+2s-8)+K(s+1) とおけば、

Gyr(s)=a(s2+2s8)+K(s+1)Δ(s),Gyd(s)=(s+1)(s2+2s8)Δ(s).G_{yr}(s)=\frac{a(s^2+2s-8)+K(s+1)}{\Delta(s)},\qquad G_{yd}(s)=\frac{(s+1)(s^2+2s-8)}{\Delta(s)}.

(3)

L0=F(0)+KP(0)=aK/8L_0=F(0)+KP(0)=a-K/8 とおく。安定性から最終値定理を適用すると

2=Gyr(0)+cGyd(0)=L0+c1+L0.2=G_{yr}(0)+cG_{yd}(0)=\frac{L_0+c}{1+L_0}.

よって

c=2+aK8.\boxed{c=2+a-\frac K8}.

(4)

例えば a=2,K=0a=-2,K=0 とすれば、

Gyr(s)=F(s)1+F(s)=2s1.G_{yr}(s)=\frac{F(s)}{1+F(s)}=-\frac2{s-1}.

唯一の極は s=1s=1 であり、その実部は正である。従って、そのような a,Ka,K は存在する。

問題2

(1)

位相交差周波数でのゲインは約 7dB-7\,\mathrm{dB}、ゲイン交差周波数での位相は約 135-135^\circ と読み取れる。従って

GM7dB,PM45.\boxed{GM\simeq7\,\mathrm{dB},\qquad PM\simeq45^\circ}.

(2)

KPK_Pk>0k>0 倍すると、位相線図は変わらず、ゲイン線図が 20log10k20\log_{10}k だけ上下する。ゲイン余裕を 10dB10\,\mathrm{dB} 増やす条件は

20log10k=10,k=110.20\log_{10}k=-10, \qquad\boxed{k=\frac1{\sqrt{10}}}.

(3)

PID 制御では u=KP(ry)+(KP/TI)(ry)dt+KPTD(r˙y˙)u=K_P(r-y)+(K_P/T_I)\int(r-y)\,dt+K_PT_D(\dot r-\dot y) となり、設定値のステップ変化が比例項の跳躍と微分項のインパルスを直接生じさせる。I-PD 制御では

u=KPTI(ry)dtKPyKPTDy˙u=\frac{K_P}{T_I}\int(r-y)\,dt-K_Py-K_PT_D\dot y

となり、設定値 rr は積分項だけに入る。従って設定値の変化による比例キックと微分キックを抑制できる。

(4)

r=0r=0 とすると、I-PD 制御でも u=C(s)yu=-C(s)y となる。y=P(s)u+dy=P(s)u+d より

Gyd(s)=11+C(s)P(s).\boxed{G_{yd}(s)=\frac1{1+C(s)P(s)}}.

PID 制御も r=0r=0 で同じ関係を満たすため、同一パラメータのとき外乱から制御量への閉ループ伝達関数は等しい。