線形制御系設計

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

    モーションコントロール

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

    モータドライブ

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

    システム同定

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

ブロック線図を用いたシステムの表現

ブロック線図は信号の流れを視覚的に確認するために使用される。

ブロック線図

入力 rr から出力 yy までの伝達関数が GG と表現される場合、次のように表現する。
image is here
システム全体を表現する場合には複数のブロックを接続することで表現する。 そのような場合にはシステムの線形性等を利用してブロック線図を等価変換し、視認性の高い表現に変形することができる。

ブロック線図の合成

システムが線形であるとき、線形写像の和および合成も線形写像となることから、並列接続されたブロックの伝達関数は次のようになる。
image is here
Gry(s)=i=1nGi(s)\begin{align} G_{\rm ry}(s) = \sum^n_{i=1} G_{i}(s) \end{align}
また、直列接続されたブロックの伝達関数は次のようになる。
image is here
Gry(s)=i=1nGi(s)\begin{align} G_{\rm ry}(s) = \prod^n_{i=1} G_{i}(s) \end{align}
線形写像の合成は可換であるため、直列接続されたブロックの順序を入れ替えた場合にも入出力伝達関数は変化しない。

フィードバック系の表現

入出力の開ループ伝達特性が Gop(s)G_{\rm op}(s), 帰還伝達特性が H(s)H(s)の伝達経路から構成されるフィードバック系について考える。
image is here
この入出力特性は、次の代数方程式を解くことで得られる。
y=Gop(rHy)y=Gop1+GHr\begin{align} y &= G_{\rm op}(r - Hy) \\ \therefore y &= \frac{G_{\rm op}}{1 + GH}r \end{align}
伝達関数の積 GHGH は経路を一巡し帰還するまでの伝達特性を表し、一巡伝達関数と呼ばれる。フィードバック制御系の閉ループ伝達関数 GclG_{cl} は、開ループ伝達特性伝達特性 GopG_{\rm op}、一巡伝達特性 LL を用いて次のように記述される。
Gcl=Gop1+L\begin{align} G_{cl} = \frac{G_{\rm op}}{1 + L} \end{align}

波動現象

フィードバックシステムにおける開ループ伝達と一巡伝達の作用を理解するために、一次元波動方程式で表現される波動現象について確認する。波動現象は波の伝達によって記述することが可能であり、波の進行と反射からなるフィードバック系として説明することができる。ここでは有界領域 x[0,l]x \in [0, l] において定義される以下の一次元波動方程式について考える。ただし、t,x,ct, x, c はそれぞれ時刻、座標、波速を表し、y(x,t)y(x,t)は時刻 tt、座標 xx における波の変位を表す。
22ty(x,t)=c222xy(x,t)\begin{align} \frac{\partial^2}{\partial^2 t}y(x, t) = c^2 \frac{\partial^2}{\partial^2 x}y(x, t) \end{align}
このシステムの期状態は以下のように記述されるものとする。
y(x,0)=0\begin{align} y(x, 0) &= 0 \end{align}
このシステムの始点は強制変位 u(t)u(t) で与えられ、また終点は自由端であるとして境界条件を以下のように記述する。
y(0,t)=u(t)xy(l,t)=0\begin{align} y(0, t) &= u(t) \\ \frac{\partial}{\partial x} y(l, t) &= 0 \end{align}
この条件下で偏微分方程式を解くことにより、次の解を得る。
y(x,s)=excs+e2lxcs1+e2lcsu(s)\begin{align} y(x, s) = \frac{e^{-\frac{x}{c}s} + e^{-\frac{2l - x}{c}s}}{1 + e^{-\frac{2l}{c}s}} u(s) \end{align}
ここで、任意の座標y(x,s)y(x, s)y(x,t)y(x, t) のラプラス変換を表す。始点の変位 u(t)u(t) から座標 xx の変位 y(x,t)y(x, t) までの伝達関数は時間遅れ要素を用いて表現される。この遅れは波の伝播の遅れと見ることができる。伝達関数の分子に存在する二つの時間遅れ要素は、始点から発生した波が座標 xx に到達するまでの時間 x/cx/c と、終点で反射した後に座標 xx に到達するまでの時間 (2lx)/c(2l - x)/c を表す。すなわち、分子は始点から発生した波の座標 xx における進行波と反射波の和となり、開ループ伝達特性を表すことがわかる。また、伝達関数の分母に存在する時間遅れ要素は始点から発生した波の往復時間 2l/c2l/c と一致し、一巡伝達関数であることが確認できる。
image is here

ブロック線図の分解

実数 ai (i=1,...,n)a_{i}\ (i=1,...,n) および bj (j=1,...,m)b_{j}\ (j=1,...,m) により伝達関数 GG が次のように表現される場合について考える。ただし、nmn \le m とする。
G=k=0nakskk=0mbksk\begin{align} G &= \frac{\sum^{n}_{k=0} a_{k}s^k}{\sum^{m}_{k=0} b_{k}s^k} \end{align}
この伝達関数に関して次のように変形する。
G=1bmsm1+1bmsmk=0m1bkskk=0naksk\begin{align} G &= \frac{\frac{1}{b_{m}s^m}}{1 + \frac{1}{b_{m}s^m}\sum^{m-1}_{k=0} b_{k}s^k} \sum^{n}_{k=0} a_{k}s^k \end{align}
ここで、次のブロック線図は上記の伝達特性を実現するものとなる。
image is here

ダイナミクスの図示

ブロック線図は微分方程式で記述されるシステムに対する視認性を向上し、システムの特性に対する理解を補助することができる。 高次のダイナミクスを持つシステムは複雑な入出力伝達関数を有するが、微分方程式からブロック線図を構築することで物理的に明快な表現を与える。 また、この表現ではシステム内部の信号の流れが確認できることから、物理的知見からの設計を促進する。

二慣性共振系

振動現象を記述するモデルとして、多慣性共振系が挙げられる。ここでは簡単のため二慣性共振系の伝達関数を計算することを目的とし、以下の連立微分方程式で記述されるシステムを考える。ただし、m1,m2m_{1}, m_{2} および x1,x2x_{1}, x_{2} は慣性体 1,21, 2 の質量および座標とし、k,fk, f は復元力定数および慣性体 11 に作用する外力とする。
m1x¨1=k(x1x2)+fm2x¨2=k(x1x2)\begin{align} m_{1}\ddot{x}_{1} &= -k (x_{1} - x_{2}) + f \\ m_{2}\ddot{x}_{2} &= k (x_{1} - x_{2}) \end{align}
この連立微分方程式からブロック線図を作成すると、次の表現を得る。
image is here
ここで、ブロック線図より次の関係式を見つけることができる。
x2=km2s2+kx1x1=1s2m2s2+km1m2s2+(m1+m2)kfx2=1s2km1m2s2+(m1+m2)kf\begin{align} x_{2} &= \frac{k}{m_{2}s^2 + k} x_{1} \\ x_{1} &= \frac{1}{s^2}\frac{m_{2}s^2 + k}{m_{1}m_{2}s^2 + (m_{1} + m_{2})k}f \\ \therefore x_{2} &= \frac{1}{s^2}\frac{k}{m_{1}m_{2}s^2 + (m_{1} + m_{2})k}f \end{align}