線形制御系設計

    線形システムの表現
    線形システムの特性解析
    離散時間システム
    制御器設計
    追従制御 (FF)
    外乱抑圧制御 (FB)
    観測器設計
    むだ時間を含む制御系

    モーションコントロール

    加速度制御
    剛体の姿勢制御

    モータドライブ

    駆動制御法
    磁極位置推定法

    システム同定

    ホワイトボックスモデル
    パラメトリックモデル

伝達関数の振幅・位相特性

線形システムは重ね合わせの原理より入出力伝達特性を周波数領域で議論することができる。

周波数伝達関数

重ね合わせの原理より,線形システムでは各周波数について入出力の振幅増幅率と位相差が与えられる。すなわち,入出力伝達特性は周波数を基準として議論することができる。伝達関数 G(s)G(s) に対して以下の極限が存在する場合にこれを周波数伝達関数 F(ω)F(\omega) と呼ぶ。
F(ω)=lim⁡s→jωG(s) \begin{align} F(\omega) = \lim_{s \rightarrow j\omega} G(s) \end{align}
今,G(s)G(s) が重複度 nin_{\rm i} の安定極 pip_{\rm i} (i=1...m)({\rm i}=1...{\rm m}) を持つシステムとする。また,G(s)G(s) の極およびその重複度を px, nx\bm{p}_{x},\ \bm{n}_{x} として以下のように記述する。
px∈Rmnx∈Rm \begin{align} \bm{p}_{x} &\in \mathbb{R}^{\rm m} \\ \bm{n}_{x} &\in \mathbb{R}^{\rm m} \\ \end{align}
このシステムに対して振幅 CC の周波数 ωa\omega_{\rm a} の正弦波 u(t)u(t) を入力することを考える。
u(t)=Csin⁡(ωat) \begin{align} u(t) &= C\sin(\omega_{\rm a}t) \end{align}
このとき,システムの出力は次のように記述される。
Y(s)=G(s)U(s)U(s)=Cωas2+ωa2 \begin{align} Y(s) &= G(s)U(s)\\ U(s) &= \frac{C\omega_{\rm a}}{s^2+\omega_{\rm a}^2} \end{align}
U(s)U(s) の極 ±jωa\pm j\omega_{\rm a} は虚軸上に存在し,G(s)G(s) の安定極 pip_{\rm i} と重なることはない。ここで,Y(s)Y(s) の極およびその重複度を py, ny\bm{p}_{y},\ \bm{n}_{y} とすると,次のように記述される。
py=[px−jωa+jωa]∈Rm+2ny=[nx11]∈Rm+2 \begin{align} \bm{p}_{y} &= \begin{bmatrix} \bm{p}_{x} & -j\omega_{\rm a} & +j\omega_{\rm a} \end{bmatrix} \in \mathbb{R}^{\rm m+2} \\ \bm{n}_{y} &= \begin{bmatrix} \bm{n}_{x} & 1 & 1 \end{bmatrix} \in \mathbb{R}^{\rm m+2} \\ \end{align}
このとき,Y(s)Y(s) は次のように部分分数分解表現される。
Y(s)=∑i=1m+2∑j=1ny,iAi,j(s−py,i)j=∑i=1m∑j=1nx,iAi,j(s−px,i)j+Am+1,1s+jωa−Am+2,1s−jωaAi,j≡1(ny,i−j)!lim⁡s→py,idny,i−jdsny,i−j[(s−py,i)ny,iY(s)]\begin{align} Y(s) &= \sum^{\rm m+2}_{i=1}\sum^{n_{y, i}}_{j=1} \frac{A_{i,j}}{(s-p_{y, i})^{j}} \\ &= \sum^{\rm m}_{i=1}\sum^{n_{x, i}}_{j=1} \frac{A_{i,j}}{(s-p_{x, i})^{j}} + \frac{A_{\rm m+1, 1}}{s+j\omega_{\rm a}} - \frac{A_{\rm m+2, 1}}{s-j\omega_{\rm a}} \\ A_{i,j} &\equiv \frac{1}{(n_{y, i} - j)!} \lim _{s\rightarrow p_{y,i}} \frac{d^{n_{y, i}-j}}{ds^{n_{y, i}-j}}[(s - p_{y, i})^{n_{y, i}}Y(s)] \end{align}
また,Am+1,1, Am+2,1A_{\rm m+1, 1},\ A_{\rm m+2, 1} は以下のように記述される。
Am+1,1=lim⁡s→−jωa(s+jωa)Y(s)=lim⁡s→−jωaG(s)Cωas−jωa=−C2jG(−jωa)Am+2,1=lim⁡s→+jωa(s−jωa)Y(s)=lim⁡s→+jωaG(s)Cωas+jωa=+C2jG(jωa)\begin{align} A_{\rm m+1, 1} &= \lim _{s\rightarrow -j\omega_{\rm a}} (s + j\omega_{\rm a})Y(s) \\ &= \lim _{s\rightarrow -j\omega_{\rm a}} G(s)\frac{C\omega_{\rm a}}{s-j\omega_{\rm a}} \\ &= -\frac{C}{2j}G(-j\omega_{\rm a}) \\ A_{\rm m+2, 1} &= \lim _{s\rightarrow +j\omega_{\rm a}} (s - j\omega_{\rm a})Y(s) \\ &= \lim _{s\rightarrow +j\omega_{\rm a}} G(s)\frac{C\omega_{\rm a}}{s+j\omega_{\rm a}} \\ &= +\frac{C}{2j}G(j\omega_{\rm a}) \end{align}
ここで,逆ラプラス変換から原関数 y(t)y(t) を以下のように計算することができる。
y(t)=∑i=1m(∑j=1nx,iAi,j(j−1)!tj−1)epx,it+Am+1,1e−jωat+Am+2,1ejωat=∑i=1m(∑j=1nx,iAi,j(j−1)!tj−1)epx,it+C2j(ejωatG(jωa)−e−jωatG(−jωa))\begin{align} y(t) &= \sum^{\rm m}_{i=1} \left( \sum^{n_{x, i}}_{j=1}\frac{A_{i,j}}{(j-1)!}t^{j-1} \right) e^{p_{x,i}t} + A_{\rm m+1, 1}e^{-j\omega_{\rm a}t} + A_{\rm m+2, 1}e^{j\omega_{\rm a}t} \\ &= \sum^{\rm m}_{i=1} \left( \sum^{n_{x, i}}_{j=1}\frac{A_{i,j}}{(j-1)!}t^{j-1} \right) e^{p_{x,i}t} + \frac{C}{2j} \left(e^{j\omega_{\rm a}t}G(j\omega_{\rm a}) -e^{-j\omega_{\rm a}t}G(-j\omega_{\rm a}) \right) \end{align}
システムの極が全て安定であるため,px,i<0p_{x, i} < 0 であり右辺第一項は時間経過とともに 00 に収束する。すなわち,十分に時間が経過した際には以下の式が成立する。
y(t)=C2j(ejωatG(jωa)−e−jωatG(−jωa))\begin{align} y(t) &= \frac{C}{2j} \left(e^{j\omega_{\rm a}t}G(j\omega_{\rm a}) -e^{-j\omega_{\rm a}t}G(-j\omega_{\rm a}) \right) \end{align}
ここで,G(jωa)G(j\omega_{\rm a}) と G(−jωa)G(-j\omega_{\rm a}) は共役な複素数であるため,次の性質を持つ。
∣G(jωa)∣=∣G(−jωa)∣∠G(jωa)=−∠G(−jωa)\begin{align} |G(j\omega_{\rm a})| &= |G(-j\omega_{\rm a})| \\ \angle G(j\omega_{\rm a}) &= - \angle G(-j\omega_{\rm a}) \end{align}
G(jωa)G(j\omega_{\rm a}) の振幅および位相を gG, ϕGg_{\rm G},\ \phi_{\rm G} と表現すれば,G(jωa), G(−jωa)G(j\omega_{\rm a}),\ G(-j\omega_{\rm a}) は次のように記述される。
gG≡∣G(jωa)∣ϕG≡∠G(jωa)⇔{G(jωa)=gGejϕGG(−jωa)=gGe−jϕG\begin{align} &g_{\rm G} \equiv |G(j\omega_{\rm a})|\\ &\phi_{\rm G} \equiv \angle G(j\omega_{\rm a})\\ \Leftrightarrow & \left\{ \begin{matrix} G(j\omega_{\rm a}) = g_{\rm G}e^{j\phi_{\rm G}} \\ G(-j\omega_{\rm a}) = g_{\rm G}e^{-j\phi_{\rm G}} \end{matrix} \right. \end{align}
したがって,y(t)y(t) は次のように表現される。
y(t)=gGC12j(ej(ωa+ϕG)t−e−j(ωa+ϕG)t)=gGCsin⁡(ωat+ϕG)\begin{align} y(t) &= g_{\rm G}C\frac{1}{2j}\left( e^{j(\omega_{\rm a} + \phi_{\rm G})t} -e^{-j(\omega_{\rm a}+\phi_{\rm G})t} \right) \\ &= g_{\rm G}C \sin (\omega_{\rm a}t + \phi_{\rm G}) \end{align}
すなわち,入出力信号は等しい周波数を持ち,入力は振幅増幅 gGg_{\rm G} および 位相変動 ϕG\phi_{\rm G} を伴って出力となる。この振幅・位相変化特性は周波数伝達関数を用いて記述することができる。
gG=∣F(ωa)∣ϕG=F(ωa)\begin{align} g_{\rm G} &= |F(\omega_{\rm a})|\\ \phi_{\rm G} &= F(\omega_{\rm a}) \end{align}

振幅・位相特性の計算

有理関数形式の周波数伝達関数 F(ω)F(\omega) が分子多項式および分母多項式が共通因子を持たない実係数複素多項式 N(ω),D(ω)N(\omega), D(\omega) によって表現される場合について考える。
F(ω)≡N(jω)D(jω) \begin{align} F(\omega) &\equiv \frac{N(j\omega)}{D(j\omega)} \end{align}
このとき,振幅・位相特性は以下のように記述される。
∣F(ω)∣=∣N(jω)∣∣D(jω)∣∠F(ω)=∠N(jω)−∠D(jω) \begin{align} |F(\omega)| &= \frac{|N(j\omega)|}{|D(j\omega)|} \\ \angle F(\omega) &= \angle N(j\omega) - \angle D(j\omega) \end{align}

周波数伝達関数の合成

複素数 cc が実数 a,ba, b を用いて以下のように直交座標表示される場合について考える。
c≡a+jb \begin{align} c \equiv a + jb \end{align}
このとき,複素数 cc は振幅 ∣c∣|c| と位相 ∠c\angle c を用いて以下のように極座標表示することができ,振幅成分と偏角成分の積として表現することができる。
c=∣c∣ej∠c∣c∣=a2+b2∠c=tan⁡−1(ba) \begin{align} c &= |c| e^{j\angle c} \\ |c| &= \sqrt{a^2 + b^2} \\ \angle c &= \tan^{-1}\left(\frac{b}{a}\right) \end{align}
すなわち,任意の複素数 di (i=1,...,n)d_{i}\ (i=1,...,n) の総乗 ee について,以下の式が成立する。
e≡∏i=1ndi∣e∣=∏i=1n∣di∣∠e=∑i=1n∠di \begin{align} e &\equiv \prod^{n}_{i=1} d_{i}\\ |e| &= \prod^{n}_{i=1} |d_{i}| \\ \angle e &= \sum^{n}_{i=1} \angle d_{i} \end{align}
周波数伝達関数の積についても同様の規則が成立し,合成された伝達関数の振幅は各伝達関数の振幅の総乗,位相は各伝達関数の位相の総和となる。