導入 図 1:スライディングモード制御器によって安定化されるシステムの位相平面 軌跡。初期到達段階の後、システム状態は線に沿って「スライド」する。s = 0 {\displaystyle s=0} 特定のs = 0 {\displaystyle s=0} この曲面は、この曲面に制約されたときに望ましい次数削減ダイナミクスを持つため選択されます。この場合、s = x 1 + x ˙ 1 = 0 \displaystyle s=x_{1}+{\dot {x}}_{1}=0} 表面は一次LTIシステムに対応する x ˙ 1 = − x 1 \displaystyle {\dot {x}}_{1}=-x_{1}} これは指数関数的に安定した 起源を持つ。 図 1は、スライディングモード制御下のシステムの軌跡の一例を示しています。スライディング面は次のように記述されます。s = 0 {\displaystyle s=0} 、そして、システム軌道が表面に到達してから有限時間が経過すると、表面に沿ったスライディングモードが開始されます。スライディングモードの理論的記述では、システムはスライディング面内に留まり、表面に沿って滑っているとだけ考えればよいとされています。しかし、実際のスライディングモード制御の実装では、この理論的挙動を、高周波で一般的に非決定論的なスイッチング制御信号で近似し、スライディング面の狭い範囲でシステムが「チャタリング」[ nb 1 ] を起こすようにします。チャタリングは、スライディング面の周囲にデッドバンド や境界層を設ける、あるいは他の補償方法を用いることで低減できます。システムは一般的に非線形ですが、図 1のシステムがスライディング面内に留まっている場合の理想化された (つまり、チャタリングのない) 挙動は、s = 0 {\displaystyle s=0} サーフェスは指数的に安定な 原点を持つLTI システム です。補償方法の 1 つは、 [ 2 ] [ 3 ] で提案された適応スライディングモード制御法で、推定された不確実性を使用して連続制御則を構築します。この方法では、精度を維持しながらチャタリングが除去されます (詳細は参考文献 [2] および [3] を参照)。提案された適応スライディングモードコントローラの 3 つの特徴は次のとおりです。(i) 構造化された (またはパラメトリックな) 不確実性と非構造化された不確実性 (モデル化されていないダイナミクス、未知の外部擾乱) は、集中不確実性と呼ばれる単一のタイプの不確実性項に統合されます。したがって、システムの線形にパラメータ化された動的モデルは必要なく、このアプローチの単純な構造と計算効率の良い特性により、リアルタイム制御アプリケーションに適しています。(ii) 適応スライディングモード制御スキームの設計は、最悪のシナリオ (つまり、不確実性の境界) に依存するのではなく、オンラインで推定された不確実性ベクトルに依存します。したがって、不確実性の範囲に関する事前知識は必要なく、各時点において、制御入力は存在する不確実性を補償します。(iii) スライディングモード制御理論の基本原理を用いて開発された連続制御則は、境界層アプローチでよく見られる性能と堅牢性のトレードオフなしにチャタリング現象を排除します。
直感的に言えば、スライディングモード制御は、実質的に無限のゲインを用いて、 動的システム の軌道を制限されたスライディングモード部分空間に沿ってスライドさせます。この次数削減されたスライディングモードからの軌道は、望ましい特性(例えば、システムが望ましい平衡点 に達するまで自然にスライドする)を持ちます。スライディングモード制御の主な強みは、その堅牢性 です。制御は2つの状態(例えば、「オン」/「オフ」または「前進」/「後退」)の切り替えのように単純なので、高精度である必要はなく、制御チャネルに入力されるパラメータの変動にも敏感ではありません。さらに、制御法則は連続関数ではないため、スライディングモードには 有限 時間で到達できます(つまり、漸近挙動よりも優れています)。特定の一般的な条件下では、最適性を得るためには バンバン制御 の使用が必要となります。したがって、スライディングモード制御は、幅広い動的システムに対する最適なコントローラ を記述します。
スライディングモードコントローラの応用例の 1 つは、スイッチング電力変換器によって動作する電気駆動装置の制御です。[ 4 ] : 「はじめに」これらの変換器の動作モードが不連続であるため、パルス幅変調 [ 5 ] または同様の技術[ nb 2 ] によって連続信号を離散状態しか取れない出力に適用する必要がある連続コントローラよりも、不連続スライディングモードコントローラの方が自然な実装の選択肢となります。スライディングモード制御はロボット工学で多くの応用例があります。特に、この制御アルゴリズムは、模擬荒波における無人水上船舶の追跡制御に非常に高い成功度で使用されています。[ 6 ] [ 7 ]
スライディングモード制御は、より穏やかな制御動作を持つ他の非線形制御 よりも慎重に適用する必要があります。特に、アクチュエータには遅延やその他の不完全性があるため、ハードなスライディングモード制御動作は、チャタリング、エネルギー損失、プラントの損傷、およびモデル化されていないダイナミクスの励起につながる可能性があります。[ 8 ] : 554~556 連続制御設計手法は、これらの問題の影響を受けにくく、スライディングモードコントローラを模倣することができます。[ 8 ] : 556~563
理論的基礎 以下の定理は、可変構造制御の基礎を形成する。
定理3:滑り運動 させて
∂ σ ∂ x B ( x 、 t ) {\displaystyle {\frac {\partial \sigma }{\partial {\mathbf {x} }}}B(\mathbf {x} ,t)} 非特異 である。つまり、システムには、軌道をスライディングモードに近づけることができる制御が常に存在するような制御可能性がある。そして、スライディングモードが σ ( x ) = 0 {\displaystyle \sigma (\mathbf {x} )=\mathbf {0} } が達成されると、システムはそのスライディングモードにとどまります。スライディングモードの軌道に沿って、σ ( x ) {\displaystyle \sigma (\mathbf {x} )} は一定であるため、スライディングモード軌道は微分方程式で記述される。
σ ˙ = 0 {\displaystyle {\dot {\sigma }}=\mathbf {0} } 。もしx {\displaystyle \mathbf {x} } -この微分方程式 に関して平衡が 安定して いる場合、システムは平衡に向かって滑りモード面に沿って滑ります。
スライディングモードに対する等価制御法則は、以下の式を解くことで求められる。
σ ˙ ( x ) = 0 {\displaystyle {\dot {\sigma }}(\mathbf {x} )=0} 等価制御法則の場合u ( x ) {\displaystyle \mathbf {u} (\mathbf {x} )} つまり、
∂ σ ∂ x ( f ( x 、 t ) + B ( x 、 t ) u ) ⏞ x ˙ = 0 {\displaystyle {\frac {\partial \sigma }{\partial \mathbf {x} }}\overbrace {\left(f(\mathbf {x} ,t)+B(\mathbf {x} ,t)\mathbf {u} \right)} ^{\dot {\mathbf {x} }}=\mathbf {0} } そのため同等の制御
u = − ( ∂ σ ∂ x B ( x 、 t ) ) − 1 ∂ σ ∂ x f ( x 、 t ) {\displaystyle \mathbf {u} =-\left({\frac {\partial \sigma }{\partial \mathbf {x} }}B(\mathbf {x} ,t)\right)^{-1}{\frac {\partial \sigma }{\partial \mathbf {x} }}f(\mathbf {x} ,t)} つまり、実際の制御はu {\displaystyle \mathbf {u} } 連続的 ではなく、スライディングモードを横切る急速な切り替えではσ ( x ) = 0 {\displaystyle \sigma (\mathbf {x} )=\mathbf {0} } システムは、あたかもこの継続的な制御によって駆動されているかのように動作する ように強制される。
同様に、スライディングモードにおけるシステム軌道は、
x ˙ = f ( x 、 t ) − B ( x 、 t ) ( ∂ σ ∂ x B ( x 、 t ) ) − 1 ∂ σ ∂ x f ( x 、 t ) ⏞ f ( x 、 t ) + B ( x 、 t ) u = f ( x 、 t ) ( 私 − B ( x 、 t ) ( ∂ σ ∂ x B ( x 、 t ) ) − 1 ∂ σ ∂ x ) {\displaystyle {\dot {\mathbf {x} }}=\overbrace {f(\mathbf {x} ,t)-B(\mathbf {x} ,t)\left({\frac {\partial \sigma }{\partial \mathbf {x} }}B(\mathbf {x} ,t)\right)^{-1}{\frac {\partial \sigma }{\partial \mathbf {x} }}f(\mathbf {x} ,t)} ^{f(\mathbf {x} ,t)+B(\mathbf {x} ,t)u}=f(\mathbf {x} ,t)\left(\mathbf {I} -B(\mathbf {x} ,t)\left({\frac {\partial \sigma }{\partial \mathbf {x} }}B(\mathbf {x} ,t)\right)^{-1}{\frac {\partial \sigma }{\partial \mathbf {x} }}\right)} 結果として得られるシステムは、スライディングモード微分方程式に一致する。
σ ˙ ( x ) = 0 {\displaystyle {\dot {\sigma }}(\mathbf {x} )=\mathbf {0} } スライディングモード面σ ( x ) = 0 {\displaystyle \sigma (\mathbf {x} )=\mathbf {0} } 到達段階からの軌道条件は、上記で導出されたより単純な条件に簡略化される。したがって、システムはより単純な条件に従うと想定できる。σ ˙ = 0 {\displaystyle {\dot {\sigma }}=0} システムがスライディングモードを見つける期間中の初期過渡状態後の状態。等式が成り立つと、ほぼ同じ動きが維持されます。σ ( x ) = 0 {\displaystyle \sigma (\mathbf {x} )=\mathbf {0} } あくまで近似的に成り立つ。
これらの定理から、スライディング運動は制御チャネルを通してシステムに入る十分に小さな外乱に対して不変(つまり、鈍感)であることがわかる。つまり、制御が十分に大きい限り、σ ⊺ σ ˙ < 0 {\displaystyle \sigma ^{\intercal }{\dot {\sigma }}<0} そしてσ ˙ {\displaystyle {\dot {\sigma }}} がゼロから一様に離れている場合、スライディングモードは外乱がないかのように維持されます。スライディングモード制御の特定の外乱やモデルの不確実性に対する不変性は、その最も魅力的な特徴であり、非常にロバスト です。
以下の例で説明するように、スライディングモード制御則は制約を維持できます
x ˙ + x = 0 {\displaystyle {\dot {x}}+x=0} 漸近的に安定化させるには、
x ¨ = 1 ( t 、 x 、 x ˙ ) + u {\displaystyle {\ddot {x}}=a(t,x,{\dot {x}})+u} いつ1 ( ⋅ ) {\displaystyle a(\cdot )} には有限の上限があります。この場合、スライディングモードは次のようになります。
x ˙ = − x {\displaystyle {\dot {x}}=-x} (つまり、x ˙ + x = 0 {\displaystyle {\dot {x}}+x=0} つまり、このようにシステムが制約されると、単純な安定 線形システム のように振る舞い、したがって、グローバルに指数的に安定な平衡点を持つ。( x 、 x ˙ ) = ( 0 、 0 ) {\displaystyle (x,{\dot {x}})=(0,0)} 起源。
制御設計の例 式(1 )で記述される単一入力u ( すなわち、 m = 1 {\displaystyle m=1} ) スイッチング関数は線形結合 として選択される。 重量s 私 > 0 {\displaystyle s_{i}>0} すべての人々のために1 ≤ 私 ≤ n {\displaystyle 1\leq i\leq n} 滑り面は単体で あり、σ ( x ) = 0 {\displaystyle \sigma (\mathbf {x} )=0} 軌道がこの表面に沿って滑るように強制されると、 σ ˙ ( x ) = 0 {\displaystyle {\dot {\sigma }}(\mathbf {x} )=0} など s 1 x ˙ 1 + s 2 x ˙ 2 + ⋯ + s n − 1 x ˙ n − 1 + s n x ˙ n = 0 {\displaystyle s_{1}{\dot {x}}_{1}+s_{2}{\dot {x}}_{2}+\cdots +s_{n-1}{\dot {x}}_{n-1}+s_{n}{\dot {x}}_{n}=0} これは次数削減システムです(つまり、新しいシステムは次数です)n − 1 {\displaystyle n-1} システムはこれに制約されているため( n − 1 ) {\displaystyle (n-1)} 次元スライディングモードシンプレックス)。この曲面は好ましい特性を持つ可能性があります(例えば、植物のダイナミクスがこの曲面に沿って滑るように強制されると、原点に向かって移動します)。x = 0 {\displaystyle \mathbf {x} =\mathbf {0} } 式(3 ) のリアプノフ関数 の導関数をとると、次のようになる。 V ˙ ( σ ( x ) ) = σ ( x ) T ⏞ ∂ V ∂ x σ ˙ ( x ) ⏞ d σ d t {\displaystyle {\dot {V}}(\sigma (\mathbf {x} ))=\overbrace {\sigma (\mathbf {x} )^{\text{T}}} ^{\tfrac {\partial V}{\partial \mathbf {x} }}\overbrace {{\dot {\sigma }}(\mathbf {x} )} ^{\tfrac {\operatorname {d} \sigma }{\operatorname {d} t}}} 確実にするためにV ˙ < 0 {\displaystyle {\dot {V}}<0} フィードバック制御法則u ( x ) {\displaystyle u(\mathbf {x} )} は、 { σ ˙ < 0 もし σ > 0 σ ˙ > 0 もし σ < 0 {\displaystyle {\begin{cases}{\dot {\sigma }}<0&{\text{if }}\sigma >0\\{\dot {\sigma }}>0&{\text{if }}\sigma <0\end{cases}}} したがって、製品はσ σ ˙ < 0 {\displaystyle \sigma {\dot {\sigma }}<0} なぜなら、それは負の数と正の数の積だからです。 制御法則u ( x ) {\displaystyle u(\mathbf {x} )} 選ばれるのは u ( x ) = { u + ( x ) もし σ ( x ) > 0 u − ( x ) もし σ ( x ) < 0 {\displaystyle u(\mathbf {x} )={\begin{cases}u^{+}(\mathbf {x} )&{\text{if }}\sigma (\mathbf {x} )>0\\u^{-}(\mathbf {x} )&{\text{if }}\sigma (\mathbf {x} )<0\end{cases}}} どこ u + ( x ) {\displaystyle u^{+}(\mathbf {x} )} は、式 ( 5 ) (つまり、σ ˙ {\displaystyle {\dot {\sigma }}} )は負の値 ですx {\displaystyle \mathbf {x} } u − ( x ) {\displaystyle u^{-}(\mathbf {x} )} は、式 ( 5 ) (つまり、σ ˙ {\displaystyle {\dot {\sigma }}} )は正 ですx {\displaystyle \mathbf {x} } 結果として生じる軌道は、滑り面に向かって移動するはずである。σ ( x ) = 0 {\displaystyle \sigma (\mathbf {x} )=0} 実際のシステムには遅延があるため、スライディングモード軌道はしばしばこのスライディング面に沿って前後に揺れ動きます(つまり、真の軌道は滑らかに追従しない可能性があります)。 σ ( x ) = 0 {\displaystyle \sigma (\mathbf {x} )=0} (ただし、スライドモードを解除すると必ずスライドモードに戻ります)。 x ¨ = 1 ( t 、 x 、 x ˙ ) + u {\displaystyle {\ddot {x}}=a(t,x,{\dot {x}})+u} これは2次元状態空間 で表現できる(x 1 = x {\displaystyle x_{1}=x} そしてx 2 = x ˙ {\displaystyle x_{2}={\dot {x}}} ) として { x ˙ 1 = x 2 x ˙ 2 = 1 ( t 、 x 1 、 x 2 ) + u {\displaystyle {\begin{cases}{\dot {x}}_{1}=x_{2}\\{\dot {x}}_{2}=a(t,x_{1},x_{2})+u\end{cases}}} また、以下を仮定します。すする { | 1 ( ⋅ ) | } ≤ k {\displaystyle \sup\{|a(\cdot )|\}\leq k} (つまり、| 1 | {\displaystyle |a|} 有限の上限k が既知である。このシステムでは、スイッチング関数を選択する。 σ ( x 1 、 x 2 ) = x 1 + x 2 = x + x ˙ {\displaystyle \sigma (x_{1},x_{2})=x_{1}+x_{2}=x+{\dot {x}}} 前述の例から、フィードバック制御法則を選択する必要があります。u ( x 、 x ˙ ) {\displaystyle u(x,{\dot {x}})} となることによってσ σ ˙ < 0 {\displaystyle \sigma {\dot {\sigma }}<0} 。 ここ、 σ ˙ = x ˙ 1 + x ˙ 2 = x ˙ + x ¨ = x ˙ + 1 ( t 、 x 、 x ˙ ) + u ⏞ x ¨ {\displaystyle {\dot {\sigma }}={\dot {x}}_{1}+{\dot {x}}_{2}={\dot {x}}+{\ddot {x}}={\dot {x}}\,+\,\overbrace {a(t,x,{\dot {x}})+u} ^{\ddot {x}}} いつx + x ˙ < 0 {\displaystyle x+{\dot {x}}<0} (つまり、σ < 0 {\displaystyle \sigma <0} )σ ˙ > 0 {\displaystyle {\dot {\sigma }}>0} 制御法則は次のように選択する必要があります。u > | x ˙ + 1 ( t 、 x 、 x ˙ ) | {\displaystyle u>|{\dot {x}}+a(t,x,{\dot {x}})|} いつx + x ˙ > 0 {\displaystyle x+{\dot {x}}>0} (つまり、σ > 0 {\displaystyle \sigma >0} )σ ˙ < 0 {\displaystyle {\dot {\sigma }}<0} 制御法則は次のように選択する必要があります。u < − | x ˙ + 1 ( t 、 x 、 x ˙ ) | {\displaystyle u<-|{\dot {x}}+a(t,x,{\dot {x}})|} しかし、三角不等式 により、 | x ˙ | + | 1 ( t 、 x 、 x ˙ ) | ≥ | x ˙ + 1 ( t 、 x 、 x ˙ ) | {\displaystyle |{\dot {x}}|+|a(t,x,{\dot {x}})|\geq |{\dot {x}}+a(t,x,{\dot {x}})|} そして、| 1 | {\displaystyle |a|} 、 | x ˙ | + k + 1 > | x ˙ | + | 1 ( t 、 x 、 x ˙ ) | {\displaystyle |{\dot {x}}|+k+1>|{\dot {x}}|+|a(t,x,{\dot {x}})|} したがって、制御法則によってシステムをフィードバック安定化(スライディングモードに戻る)することができる。 u ( x 、 x ˙ ) = { | x ˙ | + k + 1 もし x + x ˙ ⏟ < 0 、 − ( | x ˙ | + k + 1 ) もし x + x ˙ ⏞ σ > 0 {\displaystyle u(x,{\dot {x}})={\begin{cases}|{\dot {x}}|+k+1&{\text{if }}\underbrace {x+{\dot {x}}} <0,\\-\left(|{\dot {x}}|+k+1\right)&{\text{if }}\overbrace {x+{\dot {x}}} ^{\sigma }>0\end{cases}}} これは閉じた形 で次のように 表現できます。u ( x 、 x ˙ ) = − ( | x ˙ | + k + 1 ) サイン ( x ˙ + x ⏞ σ ) ⏟ (つまり、テスト) σ > 0 ) {\displaystyle u(x,{\dot {x}})=-(|{\dot {x}}|+k+1)\underbrace {\operatorname {sgn} (\overbrace {{\dot {x}}+x} ^{\sigma })} _{{\text{(i.e., tests }}\sigma >0{\text{)}}}} システムの軌道が強制的に移動させられると仮定すると、σ ( x ) = 0 {\displaystyle \sigma (\mathbf {x} )=0} 、 それから x ˙ = − x (つまり、 σ ( x 、 x ˙ ) = x + x ˙ = 0 ) {\displaystyle {\dot {x}}=-x\qquad {\text{(i.e., }}\sigma (x,{\dot {x}})=x+{\dot {x}}=0{\text{)}}} システムがスライディングモードに達すると、システムの2次元ダイナミクスは、グローバルに指数的に安定な平衡点を持つこの1次元システムのように振る舞います。 ( x 、 x ˙ ) = ( 0 、 0 ) {\displaystyle (x,{\dot {x}})=(0,0)} 。
自動化された設計ソリューション スライディングモード制御システムの設計にはさまざまな理論が存在するが、解析的および数値的手法で遭遇する実際的な困難のために、非常に効果的な設計手法が不足している。しかし、遺伝的アルゴリズム などの再利用可能な計算パラダイムを利用することで、最適設計の「解決不可能な問題」を、実際に解決可能な「非決定論的多項式問題」に変換することができる。これにより、スライディングモード制御のコンピュータによる自動設計が可能になる。[ 11 ]
スライディングモードオブザーバー スライディングモード制御は、状態オブザーバ の設計に使用できます。これらの非線形高ゲインオブザーバは、推定誤差ダイナミクスの座標を有限時間内にゼロにする能力があります。さらに、スイッチドモードオブザーバは、カルマンフィルタ と同様に魅力的な測定ノイズ耐性を備えています。[ 12 ] [ 13 ] 簡単にするために、この例では、 LTI システム 用のLuenberger オブザーバ の従来のスライディングモード修正を使用します。これらのスライディングモードオブザーバでは、システムがスライディングモードに入ると、オブザーバダイナミクスの次数が 1 つ減ります。この特定の例では、単一の推定状態の推定誤差が有限時間内にゼロになり、その後、他の推定誤差は指数関数的にゼロに減衰します。しかし、Drakunov によって最初に説明されたように、[ 14 ] 非線形システムのスライディングモードオブザーバは、 すべての推定状態の推定誤差を有限時間 (かつ任意に小さい) 内にゼロにするように設計されています。
ここでは、LTIシステムについて考えてみましょう。
{ x ˙ = A x + B u y = [ 1 0 0 ⋯ ] x = x 1 {\displaystyle {\begin{cases}{\dot {\mathbf {x} }}=A\mathbf {x} +B\mathbf {u} \\y={\begin{bmatrix}1&0&0&\cdots &\end{bmatrix}}\mathbf {x} =x_{1}\end{cases}}} 状態ベクトルx ≜ ( x 1 、 x 2 、 … 、 x n ) ∈ R n {\displaystyle \mathbf {x} \triangleq (x_{1},x_{2},\dots ,x_{n})\in \mathbb {R} ^{n}} 、u ≜ ( u 1 、 u 2 、 … 、 u r ) ∈ R r {\displaystyle \mathbf {u} \triangleq (u_{1},u_{2},\dots ,u_{r})\in \mathbb {R} ^{r}} は入力のベクトルであり、出力y は の最初の状態に等しいスカラーです。x {\displaystyle \mathbf {x} } 状態ベクトル。
A ≜ [ 1 11 A 12 A 21 A 22 ] {\displaystyle A\triangleq {\begin{bmatrix}a_{11}&A_{12}\\A_{21}&A_{22}\end{bmatrix}}} どこ
1 11 {\displaystyle a_{11}} は、最初の状態の影響を表すスカラーです。x 1 {\displaystyle x_{1}} それ自体で、A 21 ∈ R ( n − 1 ) {\displaystyle A_{21}\in \mathbb {R} ^{(n-1)}} は、最初の状態が他の状態に与える影響に対応する行ベクトルです。A 22 ∈ R ( n − 1 ) × ( n − 1 ) {\displaystyle A_{22}\in \mathbb {R} ^{(n-1)\times (n-1)}} は、他の状態が自身に与える影響を表す行列であり、A 12 ∈ R 1 × ( n − 1 ) {\displaystyle A_{12}\in \mathbb {R} ^{1\times (n-1)}} これは、他の状態が最初の状態に与える影響を表す列ベクトルです。目標は、状態ベクトルを推定する高ゲイン状態オブザーバを設計することである。x {\displaystyle \mathbf {x} } 測定からの情報のみを使用するy = x 1 {\displaystyle y=x_{1}} したがって、ベクトルをx ^ = ( x ^ 1 、 x ^ 2 、 … 、 x ^ n ) ∈ R n {\displaystyle {\hat {\mathbf {x} }}=({\hat {x}}_{1},{\hat {x}}_{2},\dots ,{\hat {x}}_{n})\in \mathbb {R} ^{n}} n 個の状態の推定値とする。観測者は次の形式をとる。
x ^ ˙ = A x ^ + B u + L v ( x ^ 1 − x 1 ) {\displaystyle {\dot {\hat {\mathbf {x} }}}=A{\hat {\mathbf {x} }}+B\mathbf {u} +Lv({\hat {x}}_{1}-x_{1})} どこv : R → R {\displaystyle v:\mathbb {R} \to \mathbb {R} } 推定状態間の誤差の非線形関数であるx ^ 1 {\displaystyle {\hat {x}}_{1}} そして出力y = x 1 {\displaystyle y=x_{1}} 、 そしてL ∈ R n {\displaystyle L\in \mathbb {R} ^{n}} は、典型的な線形ルーエンバーガーオブザーバ と同様の目的を果たすオブザーバゲインベクトルである。同様に、
L = [ − 1 L 2 ] {\displaystyle L={\begin{bmatrix}-1\\L_{2}\end{bmatrix}}} どこL 2 ∈ R ( n − 1 ) {\displaystyle L_{2}\in \mathbb {R} ^{(n-1)}} は列ベクトルです。さらに、e = ( e 1 、 e 2 、 … 、 e n ) ∈ R n {\displaystyle \mathbf {e} =(e_{1},e_{2},\dots ,e_{n})\in \mathbb {R} ^{n}} 状態推定誤差とする。つまり、e = x ^ − x {\displaystyle \mathbf {e} ={\hat {\mathbf {x} }}-\mathbf {x} } 誤差の動態は次のようになる。
e ˙ = x ^ ˙ − x ˙ = A x ^ + B u + L v ( x ^ 1 − x 1 ) − A x − B u = A ( x ^ − x ) + L v ( x ^ 1 − x 1 ) = A e + L v ( e 1 ) {\displaystyle {\begin{aligned}{\dot {\mathbf {e} }}&={\dot {\hat {\mathbf {x} }}}-{\dot {\mathbf {x} }}\\&=A{\hat {\mathbf {x} }}+B\mathbf {u} +Lv({\hat {x}}_{1}-x_{1})-A\mathbf {x} -B\mathbf {u} \\&=A({\hat {\mathbf {x} }}-\mathbf {x} )+Lv({\hat {x}}_{1}-x_{1})\\&=A\mathbf {e} +Lv(e_{1})\end{aligned}}} どこe 1 = x ^ 1 − x 1 {\displaystyle e_{1}={\hat {x}}_{1}-x_{1}} は、最初の状態推定値の推定誤差です。非線形制御則vは、 スライディングマニホールドを強制するように設計できます。
0 = x ^ 1 − x 1 {\displaystyle 0={\hat {x}}_{1}-x_{1}} そのため、その見積もりはx ^ 1 {\displaystyle {\hat {x}}_{1}} 実際の状況を追跡するx 1 {\displaystyle x_{1}} ある有限時間後(つまり、x ^ 1 = x 1 {\displaystyle {\hat {x}}_{1}=x_{1}} したがって、スライディングモード制御の切り替え機能は
σ ( x ^ 1 、 x ^ ) ≜ e 1 = x ^ 1 − x 1 。 {\displaystyle \sigma ({\hat {x}}_{1},{\hat {x}})\triangleq e_{1}={\hat {x}}_{1}-x_{1}.} スライディングマニホールドを実現するには、σ ˙ {\displaystyle {\dot {\sigma }}} そしてσ {\displaystyle \sigma } 常に反対の符号を持つ必要があります (つまり、σ σ ˙ < 0 {\displaystyle \sigma {\dot {\sigma }}<0} 基本的 にすべてのx {\displaystyle \mathbf {x} } )。 しかし、
σ ˙ = e ˙ 1 = 1 11 e 1 + A 12 e 2 − v ( e 1 ) = 1 11 e 1 + A 12 e 2 − v ( σ ) {\displaystyle {\dot {\sigma }}={\dot {e}}_{1}=a_{11}e_{1}+A_{12}\mathbf {e} _{2}-v(e_{1})=a_{11}e_{1}+A_{12}\mathbf {e} _{2}-v(\sigma )} どこe 2 ≜ ( e 2 、 e 3 、 … 、 e n ) ∈ R ( n − 1 ) {\displaystyle \mathbf {e} _{2}\triangleq (e_{2},e_{3},\ldots ,e_{n})\in \mathbb {R} ^{(n-1)}} これは、測定されていないすべての状態に対する推定誤差の集合です。σ σ ˙ < 0 {\displaystyle \sigma {\dot {\sigma }}<0} 、 させて
v ( σ ) = M サイン ( σ ) {\displaystyle v(\sigma )=M\operatorname {sgn} (\sigma )} どこ
M > 最大 { | 1 11 e 1 + A 12 e 2 | } 。 {\displaystyle M>\max\{|a_{11}e_{1}+A_{12}\mathbf {e} _{2}|\}.} つまり、正の定数M は、 システムの推定誤差の最大値をスケーリングした値(すなわち、初期誤差。初期誤差は有界であると仮定されるため、M を 十分に大きく選択できる。al)よりも大きくなければならない。M が十分に大きい場合、 システムは次のことを達成すると考えられる。e 1 = 0 {\displaystyle e_{1}=0} (つまり、x ^ 1 = x 1 {\displaystyle {\hat {x}}_{1}=x_{1}} )。 なぜならe 1 {\displaystyle e_{1}} この多様体に沿って一定(つまり0)であり、e ˙ 1 = 0 {\displaystyle {\dot {e}}_{1}=0} 同様である。したがって、不連続制御v ( σ ) {\displaystyle v(\sigma )} 同等の連続制御に置き換えることができるv 式 {\displaystyle v_{\text{eq}}} どこ
0 = σ ˙ = 1 11 e 1 ⏞ = 0 + A 12 e 2 − v 式 ⏞ v ( σ ) = A 12 e 2 − v 式 。 {\displaystyle 0={\dot {\sigma }}=a_{11}{\mathord {\overbrace {e_{1}} ^{{}=0}}}+A_{12}\mathbf {e} _{2}-{\mathord {\overbrace {v_{\text{eq}}} ^{v(\sigma )}}}=A_{12}\mathbf {e} _{2}-v_{\text{eq}}.} それで
v 式 ⏟ スカラー = A 12 ⏟ 1 × ( n − 1 ) ベクター e 2 ⏟ ( n − 1 ) × 1 ベクター 。 {\displaystyle {\mathord {\underbrace {v_{\text{eq}}} _{\text{scalar}}}}={\mathord {\underbrace {A_{12}} _{1\times (n-1) \atop {\text{ vector}}}}}{\mathord {\underbrace {\mathbf {e} _{2}} _{(n-1)\times 1 \atop {\text{ vector}}}}}.} この同等の制御v 式 {\displaystyle v_{\text{eq}}} 他者からの貢献を表す( n − 1 ) {\displaystyle (n-1)} 出力状態の軌跡への状態x 1 {\displaystyle x_{1}} 特に、行A 12 {\displaystyle A_{12}} エラーサブシステムの出力ベクトルのように機能します
[ e ˙ 2 e ˙ 3 ⋮ e ˙ n ] ⏞ e ˙ 2 = A 2 [ e 2 e 3 ⋮ e n ] ⏞ e 2 + L 2 v ( e 1 ) = A 2 e 2 + L 2 v 式 = A 2 e 2 + L 2 A 12 e 2 = ( A 2 + L 2 A 12 ) e 2 。 {\displaystyle {\mathord {\overbrace {\begin{bmatrix}{\dot {e}}_{2}\\{\dot {e}}_{3}\\\vdots \\{\dot {e}}_{n}\end{bmatrix}} ^{{\dot {\mathbf {e} }}_{2}}}}=A_{2}{\mathord {\overbrace {\begin{bmatrix}e_{2}\\e_{3}\\\vdots \\e_{n}\end{bmatrix}} ^{\mathbf {e} _{2}}}}+L_{2}v(e_{1})=A_{2}\mathbf {e} _{2}+L_{2}v_{\text{eq}}=A_{2}\mathbf {e} _{2}+L_{2}A_{12}\mathbf {e} _{2}=(A_{2}+L_{2}A_{12})\mathbf {e} _{2}.} したがって、推定誤差を確実にするためにe 2 {\displaystyle \mathbf {e} _{2}} 測定されていない状態については、ゼロに収束する。( n − 1 ) × 1 {\displaystyle (n-1)\times 1} ベクターL 2 {\displaystyle L_{2}} は、( n − 1 ) × ( n − 1 ) {\displaystyle (n-1)\times (n-1)} マトリックス( A 2 + L 2 A 12 ) {\displaystyle (A_{2}+L_{2}A_{12})} これはフルヴィッツ行列で ある(つまり、各固有値 の実部は負でなければならない)。したがって、これが観測可能で ある限り、e 2 {\displaystyle \mathbf {e} _{2}} システムは、典型的な線形状態オブザーバ とまったく同じ方法で安定化できます。A 12 {\displaystyle A_{12}} は出力行列(つまり「C 」)と見なされます。つまり、v 式 {\displaystyle v_{\text{eq}}} 等価制御は、測定されていない状態に関する測定情報を提供し、その推定値を漸近的にそれらに継続的に近づけることができる。一方、不連続制御はv = M サイン ( x ^ 1 − x ) {\displaystyle v=M\operatorname {sgn} ({\hat {x}}_{1}-x)} 測定状態の推定値は有限時間内に誤差がゼロになるように強制されます。さらに、白色で平均値がゼロの対称的な測定ノイズ(例えばガウスノイズ )は制御v のスイッチング周波数にのみ影響するため、ノイズは等価スライディングモード制御にはほとんど影響を与えません。v 式 {\displaystyle v_{\text{eq}}} したがって、スライディングモードオブザーバはカルマンフィルタ の ような特徴を持つ。[ 13 ]
オブザーバーの最終バージョンは次のようになる。
x ^ ˙ = A x ^ + B u + L M サイン ( x ^ 1 − x 1 ) = A x ^ + B u + [ − 1 L 2 ] M サイン ( x ^ 1 − x 1 ) = A x ^ + B u + [ − M L 2 M ] サイン ( x ^ 1 − x 1 ) = A x ^ + [ B [ − M L 2 M ] ] [ u サイン ( x ^ 1 − x 1 ) ] = A 観測 x ^ + B 観測 u 観測 {\displaystyle {\begin{aligned}{\dot {\hat {\mathbf {x} }}}&=A{\hat {\mathbf {x} }}+B\mathbf {u} +LM\operatorname {sgn} ({\hat {x}}_{1}-x_{1})\\&=A{\hat {\mathbf {x} }}+B\mathbf {u} +{\begin{bmatrix}-1\\L_{2}\end{bmatrix}}M\operatorname {sgn} ({\hat {x}}_{1}-x_{1})\\&=A{\hat {\mathbf {x} }}+B\mathbf {u} +{\begin{bmatrix}-M\\L_{2}M\end{bmatrix}}\operatorname {sgn} ({\hat {x}}_{1}-x_{1})\\&=A{\hat {\mathbf {x} }}+{\begin{bmatrix}B&{\begin{bmatrix}-M\\L_{2}M\end{bmatrix}}\end{bmatrix}}{\begin{bmatrix}\mathbf {u} \\\operatorname {sgn} ({\hat {x}}_{1}-x_{1})\end{bmatrix}}\\&=A_{\text{obs}}{\hat {\mathbf {x} }}+B_{\text{obs}}\mathbf {u} _{\text{obs}}\end{aligned}}} どこ
A 観測 ≜ A 、 {\displaystyle A_{\text{obs}}\triangleq A,} B 観測 ≜ [ B [ − M L 2 M ] ] 、 {\displaystyle B_{\text{obs}}\triangleq {\begin{bmatrix}B&{\begin{bmatrix}-M\\L_{2}M\end{bmatrix}}\end{bmatrix}},} そしてu 観測 ≜ [ u サイン ( x ^ 1 − x 1 ) ] 。 {\displaystyle u_{\text{obs}}\triangleq {\begin{bmatrix}\mathbf {u} \\\operatorname {sgn} ({\hat {x}}_{1}-x_{1})\end{bmatrix}}.} つまり、制御ベクトルを拡張することによってu {\displaystyle \mathbf {u} } スイッチング機能付きサイン ( x ^ 1 − x 1 ) {\displaystyle \operatorname {sgn} ({\hat {x}}_{1}-x_{1})} スライディングモードオブザーバはLTIシステムとして実装できます。つまり、不連続信号サイン ( x ^ 1 − x 1 ) {\displaystyle \operatorname {sgn} ({\hat {x}}_{1}-x_{1})} は、2入力LTIシステムの制御入力として見なされる。
簡略化のため、この例ではスライディングモードオブザーバが単一の状態(つまり出力)の測定値にアクセスできると仮定します。y = x 1 {\displaystyle y=x_{1}} ただし、同様の手順を使用して、状態の重み付き組み合わせのベクトル(つまり、出力が)に対するスライディングモードオブザーバを設計することもできます。y = C x {\displaystyle \mathbf {y} =C\mathbf {x} } (汎用行列C を使用)。いずれの場合も、スライディングモードは推定出力がy ^ {\displaystyle {\hat {\mathbf {y} }}} 測定された出力に従うy {\displaystyle \mathbf {y} } 誤差ゼロ(つまり、多様体)σ ( x ) ≜ y ^ − y = 0 {\displaystyle \sigma (\mathbf {x} )\triangleq {\hat {\mathbf {y} }}-\mathbf {y} =\mathbf {0} } )
注記 ↑ 「チャタリング」とは、有限の周波数と振幅を持つ振動 という望ましくない現象です。チャタリングは、制御精度の低下、可動機械部品の摩耗の増加、電力回路における熱損失の増加につながるため、有害な現象です。詳細については、Utkin, Vadim; Lee, Jason Hoon (2006年7月)、「スライディングモード制御システムにおけるチャタリング問題」 、vol. 10.1109/VSS.2006.1644542、pp. 346–350 ↑ その他のパルス型変調技術にはデルタシグマ変調 があります。
参考文献 1 2 Zinober, ASI 編 (1990).不確実システムの決定論的制御 . ロンドン: Peter Peregrinus Press. ISBN 978-0-86341-170-0 。↑ Zeinali M.; Notash L. (2010). "ロボットマニピュレータのための不確実性推定器を用いた適応スライディングモード制御". International Journal of Mechanism and Machine Theory . 45 (1): 80–90 . doi : 10.1016/j.mechmachtheory.2009.08.003 . ↑ Zeinali, Meysar (2018). 有限時間で軌道が実スライディングモードに収束するロボットマニピュレータのための一次連続適応スライディングモード制御 。第15回可変構造システムとスライディングモード制御に関する国際ワークショップ。オーストリア、グラーツ工科大学。 ↑ Utkin, Vadim I. (1993). "スライディングモード制御設計原理と電気駆動への応用". IEEE Transactions on Industrial Electronics . 40 (1): 23–36 . CiteSeerX 10.1.1.477.77 . doi : 10.1109/41.184818 . ↑ Tan, Siew-Chong (2012). スイッチング電源コンバータのスライディングモード制御:技術と実装 . Yuk-Ming Lai、Chi Kong Tse. フロリダ州ボカラトン:CRC Press. ISBN 978-1-4398-3025-3 。↑ 「模擬荒波環境下における無人船舶の自律航行と障害物回避 - ビラノバ大学」 ↑ Mahini; et al. (2013). "荒れた海における水上船舶の自律運用に関する実験装置". Robotica . 31 (5): 703– 715. doi : 10.1017/s0263574712000720 . S2CID 31903795 . 1 2 3 4 Khalil, HK (2002). 非線形システム (第3 版). アッパーサドルリバー、ニュージャージー州: プレンティスホール . ISBN 978-0-13-067389-3 。↑ フィリポフ、AF(1988)。 不連続な右辺を持つ微分方程式 。クルーワー 。ISBN 978-90-277-2699-5 。↑ ペルケッティ、W.バーボット、日本 (2002)。 工学におけるスライディング モード制御 。マルセル・デッカーのハードカバー。 ISBN 978-0-8247-0671-5 。↑ Li, Yun; et al. (1996). "遺伝的アルゴリズムによるスライディングモード制御システムの設計への自動化アプローチ" . International Journal of Control . 64 (3): 721–739 . CiteSeerX 10.1.1.43.1654 . doi : 10.1080/00207179608921865 . ↑ Utkin, Vadim; Guldner, Jürgen; Shi, Jingxin (1999). 電気機械システムにおけるスライディングモード制御 . フィラデルフィア、ペンシルベニア州: Taylor & Francis, Inc. ISBN 978-0-7484-0116-1 。1 2 Drakunov, SV (1983). 「不連続パラメータ を 持つ適応型準最適フィルタ」。 オートメーションとリモートコントロール 。44 (9): 1167–1175 。 ↑ Drakunov, SV (1992). "等価制御法に基づくスライディングモードオブザーバ" . [ 1992 ] Proceedings of the 31st IEEE Conference on Decision and Control . pp. 2368–2370 . doi : 10.1109/CDC.1992.371368 . ISBN 978-0-7803-0872-5 . S2CID 120072463 . Zinober, Alan SI 編 (1994).可変構造とリアプノフ制御 . 制御と情報科学に関する講義ノート. 第 193 巻. ロンドン: Springer-Verlag. doi : 10.1007/BFb0033675 . ISBN 978-3-540-19869-7 。 エドワーズ、C.、スパーゲン、S. (1998).スライディングモード制御。理論と応用 。ロンドン:テイラー&フランシス。ISBN 978-0-7484-0601-2 。 サバノビッチ、アシフ;フリッドマン、レオニード ;スパーゲン、サラ(2004)。可変構造システム:原理から実装まで 。ロンドン:電気技術者協会。ISBN 0-86341-350-1 。 Edwards, C.; Fossas Colet, E.; Fridman, L. (2006). Advances in Variable Structure and Sliding Mode Control . Berlin: Springer Verlag. ISBN 978-3-540-32800-1 。 Bartolini, G.; Fridman, L.; Pisano, A.; Usai, E. (2008).現代スライディングモード制御理論 ― 新しい展望と応用 ― ベルリン:Springer Verlag. ISBN 978-3-540-79016-7 。 Fridman, L.; Barbot, J.-P.; Plestan, F. (2016). Recent Trends in Sliding Mode Control . London: The Institution of Engineering and Technology. ISBN 978-1-78561-076-9 。 Shtessel, Y.; Edwards, C.; Fridman, L.; Levant, A. (2014).スライディングモード制御と観測 . バーゼル:Birkhauser. ISBN 978-0-81764-8923 。 シラ= ラミレス、ヘバート(2015)。スライディングモード制御:デルタシグマ変調アプローチ 。バーゼル:ビルクハウザー。ISBN 978-3-319-17256-9 。Fridman, L.; Moreno, J.; Bandyopadhyay, B.; Kamal Asif Chalanga, S.; Chalanga, S. (2015). Continuous Nested Algorithms: The Fifth Generation of Sliding Mode Controllers. In: Recent Advances in Sliding Modes: From Control to Intelligent Mechatronics, X. Yu, O. Efe (Eds.), pp. 5–35 . Studies in Systems, Decision and Control. Vol. 24. London: Springer-Verlag. doi : 10.1007/978-3-319-18290-2_2 . ISBN 978-3-319-18289-6 。 Li, S.; Yu, X.; Fridman, L.; Man, Z.; Wang, X. (2017).可変構造システムとスライディングモード制御の進歩 ― 理論と応用 ― ― システム、意思決定、制御に関する研究 ― 第 24巻 ― ロンドン:Springer-Verlag. ISBN 978-3-319-62895-0 。 Ferrara, A. ; Incremona, GP; Cucuzzella, M. (2019).高度最適化に基づくスライディングモード制御 . SIAM. doi : 10.1137/1.9781611975840 . ISBN 978-1611975833 . S2CID 198313071 . Steinberger, M.、Horn, M.、Fridman, L. 編 (2020).可変構造システムとスライディングモード制御 . システム、意思決定、制御に関する研究. 第 271 巻. ロンドン: Springer-Verlag. doi : 10.1007/BFb0033675 . ISBN 978-3-030-36620-9 。
さらに読む Drakunov SV、Utkin VI. (1992). 「動的システムにおけるスライディングモード制御」。International Journal of Control . 55 (4): 1029–1037 . doi : 10.1080/00207179208934270 . hdl : 10338.dmlcz/135339 .