Barrier Certificate Enforcement

制御動作を変更し、バリア証明書の制約およびアクション範囲を満たす

R2022a 以降

• 

ライブラリ:
Simulink Control Design / Constraint Control

説明

Barrier Certificate Enforcement ブロックは、バリア証明書の制約とアクション範囲に従って、指定された制御動作の最も近くにある変更された制御動作を計算します。

このブロックは、二次計画法 (QP) ソルバーを使用して、関数 $${\left|u-u_0\right|}^2$$ を最小化する制御動作 u を検出します。ここで、u₀ は変更されていない制御動作です。

ソルバーは、以下の制約を最適化問題に適用します。

$$\begin{array}{l}{q}_x{f}_x+q_x{g}_{x}u+\gamma h_x^{\beta }\ge 0\\ u_{\min }\le u\le u_{\max }\end{array}$$

ここで、以下となります。

• f_x と g_x は、プラント ダイナミクス $$\dot{x}=f(x)+g(x)u$$ によって定義される関数。

• h_x は制御バリア関数。

• q_x は状態 x に対する制御バリア関数の偏導関数。

• γ は制約係数。

• β は制約指数。

• u_min は制御動作の下限。

• u_max は制御動作の上限。

Barrier Certificate Enforcement ブロックには、Optimization Toolbox™ ソフトウェアが必要です。

バリア証明書の適用の詳細については、制御設計へのバリア証明書の適用を参照してください。

例

PID コントローラーに対するバリア証明書の制約の適用

Barrier Certificate Enforcement ブロックを使用して、バリア証明書の制約を PID 制御アプリケーションに適用する。

• ライブ スクリプトを開く

アダプティブ クルーズ コントロールに対するバリア証明書の制約の適用

アダプティブ クルーズ コントロールに対するバリア証明書の制約を適用して速度と車間距離を維持する。

• ライブ スクリプトを開く

衝突のないロボットに対するバリア証明書の制約の適用

2 台のロボットがそれぞれのターゲット位置に衝突なしで到達するためのバリア証明書の制約を適用する。

• ライブ スクリプトを開く

衝突のないマルチロボット システムに対するバリア証明書の制約の適用

3 台のロボットがそれぞれのターゲット位置に衝突なしで到達するためのバリア証明書の制約を適用する。

• ライブ スクリプトを開く

端子

入力

u0 — 制御動作

スカラー | ベクトル

変更されていない制御動作。スカラーまたはベクトルとして指定します。

[アクション数] パラメーターが 1 の場合、u0 をスカラー信号に接続します。そうでない場合、u0 を長さが [アクション数] に等しいベクトル信号に接続します。

fx — ステート関数

スカラー | ベクトル

次のプラント ダイナミクスの方程式のステート関数 f(x)。

$$\dot{x}=f(x)+g(x)u$$

[fx] を N_x 行 1 列の信号に接続します。ここで、N_x は [状態数] パラメーターに等しくなります。

gx — 入力関数

スカラー | ベクトル | 行列

次のプラント ダイナミクスの方程式の入力関数 g(x)。

$$\dot{x}=f(x)+g(x)u$$

[gx] を N_x 行 N_u 列の信号に接続します。ここで、N_x は [状態数] パラメーターに等しく、N_u は [アクション数] パラメーターに等しくなります。

hx — 制御バリア関数

スカラー | ベクトル

制御バリア関数。以下のプラント状態の安全セットとして定義されます。

$$\{x:h(x)\ge 0\}$$

[hx] を N_c 行 1 列の信号に接続します。ここで、N_c は [バリア証明書の数] パラメーターに等しくなります。

qx — 制御バリア関数の偏導関数

スカラー | ベクトル | 行列

プラント状態に対する制御バリア関数の偏導関数。

$$q(x)=\frac{\partial h}{\partial x}$$

[qx] を N_c 行 N_x 列の信号に接続します。ここで、N_c は [バリア証明書の数] パラメーターに等しく、N_x は [状態数] パラメーターに等しくなります。

umax — アクション信号の上限

スカラー | ベクトル

アクション信号の実行時の上限を指定するには、この入力端子を有効にします。この端子が無効化されている場合、ブロックは制御動作に上限を適用しません。

[アクション数] パラメーターが 1 の場合、umax をスカラー信号に接続します。そうでない場合、umax を長さが [アクション数] に等しいベクトル信号に接続します。

依存関係

この入力端子を有効にするには、[上限 (umax) に外部ソースを使用] パラメーターを選択します。

umin — アクション信号の下限

スカラー | ベクトル

アクション信号の実行時の下限を指定するには、この入力端子を有効にします。この端子が無効化されている場合、ブロックは制御動作に下限を適用しません。

[アクション数] パラメーターが 1 の場合、umin をスカラー信号に接続します。そうでない場合、umin を長さが [アクション数] に等しいベクトル信号に接続します。

依存関係

この入力端子を有効にするには、[下限 (umin) に外部ソースを使用] パラメーターを選択します。

出力

u* — 変更済み制御動作

スカラー | ベクトル

QP ソルバーによって返される変更済み制御動作。

ソルバーが最大反復回数に達する前に解を見つけた場合、u* はこの最適解を出力します。

ソルバーが最大反復回数に達した場合、最適化は停止し、u* は準最適解を出力します。

最初の最適化問題が実行不可能な場合、返される制御動作は、ブロックが制約またはアクション範囲を無視するように設定されているかどうかにより異なります。詳細については、[exitflag] パラメーターを参照してください。

[アクション数] パラメーターが 1 の場合、u* はスカラー信号を出力します。そうでない場合、u* は、長さが [アクション数] に等しいベクトル信号を出力します。

exitflag — 最適化ステータス

1 | 0 | 負の整数

QP ソルバーの最適化ステータス。次の表に取り得るステータス値を示します。

終了フラグ	説明
1	ソルバーは、すべての制約と範囲が有効な状態で最適解に収束しました。この場合、u* は最適な制御動作を出力します。
0	ソルバーが最大反復回数に達しました。u* の制御動作の出力は、準最適の場合があります。
負の整数	最初の最適化問題が実行不可能であり、次のいずれかのシナリオが適用されます。 アクション範囲なしで最適化を再実行しましたが、実行可能解を生成しませんでした。 制約範囲なしで最適化を再実行しましたが、実行可能解を生成しませんでした。 この場合、u* の制御動作の出力はゼロです。

依存関係

この出力端子を有効にするには、[最適化ステータス] パラメーターを選択します。

パラメーター

[パラメーター] タブ

状態数 — プラント状態の数

1 (既定値) | 正の整数

プラントの状態の数を指定します。

プログラムでの使用

ブロック パラメーター: nx

タイプ: 文字ベクトル

既定: '1'

アクション数 — 制御動作の数

1 (既定値) | 正の整数

範囲を適用し、最適化するアクションの数を指定します。

プログラムでの使用

ブロック パラメーター: nu

タイプ: 文字ベクトル

既定: '1'

バリア証明書の数 — バリア証明書の制約の数

1 (既定値) | 正の整数

適用するバリア証明書の制約の数を指定します。

プログラムでの使用

ブロック パラメーター: nc

タイプ: 文字ベクトル

既定: '1'

制約係数 — 制約係数

10 (既定値) | 正のスカラー | ベクトル

バリア証明書の制約に制約係数 γ を指定します。

[バリア証明書の数] パラメーターが 1 の場合、[制約係数] を有限の正のスカラーとして指定します。そうでない場合、[制約係数] を有限の正のスカラー値または [バリア証明書の数] と等しい長さの正のスカラーの列ベクトルとして指定できます。

プログラムでの使用

ブロック パラメーター: gamma

タイプ: 文字ベクトル

既定: '10'

制約指数 — 制約指数

1 (既定値) | 正の奇数の整数 | ベクトル

バリア証明書の制約に制約指数 β を指定します。

[バリア証明書の数] パラメーターが 1 の場合、[制約指数] を正の奇数の整数として指定します。そうでない場合、[制約指数] を正の奇数の整数または [バリア証明書の数] と等しい長さをもつ正の奇数の整数の列ベクトルとして指定できます。

プログラムでの使用

ブロック パラメーター: beta

タイプ: 文字ベクトル

既定: '1'

上限 (umax) に外部ソースを使用 — アクションの上限の入力端子を追加

off (既定値) | on

このパラメーターを選択して、アクションの外部上限の umax 入力端子を追加します。

プログラムでの使用

ブロック パラメーター: external_umax

タイプ: 文字ベクトル

値: 'off' | 'on'

既定: 'off'

下限 (umin) に外部ソースを使用 — アクションの下限の入力端子を追加

off (既定値) | on

このパラメーターを選択して、アクションの外部下限の umin 入力端子を追加します。

プログラムでの使用

ブロック パラメーター: external_umin

タイプ: 文字ベクトル

値: 'off' | 'on'

既定: 'off'

[ブロック] タブ

サンプル時間 — 最適化サンプル時間

0.1 (既定値) | 正のスカラー

最適化を実行するサンプル時間を指定します。

プログラムでの使用

ブロック パラメーター: Ts

タイプ: 文字ベクトル

既定: '0.1'

最大反復回数 — 最適化の最大反復回数

200 (既定値) | 正の整数

最適化を反復する最大回数を指定します。

プログラムでの使用

ブロック パラメーター: maxiter

タイプ: 文字ベクトル

既定: '200'

制約の許容誤差 — 制約違反の許容誤差

1e-6 (既定値) | 非負のスカラー

制約違反の許容誤差の値を指定します。

プログラムでの使用

ブロック パラメーター: tol

タイプ: 文字ベクトル

既定: '1e-6'

最適化ステータス — 終了フラグ出力端子を追加

off (既定値) | on

このパラメーターを選択して、QP ソルバーの最適化ステータスに対する [exitflag] 出力端子を追加します。

プログラムでの使用

ブロック パラメーター: exitflag

タイプ: 文字ベクトル

値: 'off' | 'on'

既定: 'off'

拡張機能

C/C++ コード生成

Simulink® Coder™ を使用して C および C++ コードを生成します。

Barrier Certificate Enforcement ブロックは、倍精度の信号のコード生成のみをサポートします。

バージョン履歴

R2022a で導入

R2023b: ライブラリの場所の変更

Barrier Certificate Enforcement ブロックは Simulink Control Design/Constraint Control サブライブラリに含まれるようになりました。

R2023b では、以前のリリースで保存されたモデルを開くと、Barrier Certificate Enforcement ブロックのライブラリ リンクが自動的に更新されます。

参考

Constraint Enforcement

トピック

• 制御設計へのバリア証明書の適用