回転逆振り子の粒子群最適化調整バックステッピング追跡制御のラピッド制御プロトタイピングシミュレーション

サマリー

このプロトコルは、固定ステップのリアルタイムシミュレーション環境内で粒子群最適化を調整したバックステップ追跡制御を評価するための標準化された迅速制御プロトタイピング手順を確立します。

要約

このプロトコルの主な目的は、非線形制御システムにおけるパーティクルスウォーム最適化(PSO)ベースのゲインチューニングを評価するための再現可能な固定ステップシミュレーションフレームワークを提供することです。実装は、古田型振り子モデルの定式化から始まり、その後2msの固定ステップ実行環境内でバックステップコントローラの統合が進みます。この手法は、非理想的な実装制約の特性評価、多目的的な粒子群最適化探索空間の定義、自動オフラインチューニングの実行、そして標準化された軌道追跡および妨害除去シナリオ群を通じて得られるパラメータの評価という体系的な4段階プロセスを含みます。この構成は高性能な産業用ワークステーションと標準化された信号インターフェースを活用し、同じ制御アーキテクチャ内で一貫した繰り返し試験比較を可能にします。この設計は、ベースラインの手動調整バックステッピングコントローラーと、まったく同じコントロール構造を共有するPSO最適化型を比較し、ゲイン選択の影響を隔離しています。コントローラーの性能は、ステップ軌道追従、混合周波数正弦波トラッキング、擾乱除去の3つの異なる運用シナリオで評価されます。10件の繰り返し試験の統計解析では、PSOに基づく最適化によりステップトラッキングRMSEを0.065から0.050radに減少し、振り子のピーク変動を33.1%減衰させました。これらの改善はRMS制御の努力を22.1%削減したとともに達成され、最適化されたパラメータがリアプノフに基づく枠組み内でより効率的なエネルギー分配を促進したことを示しています。最終的に、この手法は、後の物理ハードウェア実装前に非線形制御戦略を評価するための構造化されたシミュレーションフレームワークを提供します。

概要

古田型回転式逆振り子は、そのオープンループ不安定性と複雑なアンダーアクチュエーテッドダイナミクス^1,2,3により、非線形制御アルゴリズムの検証における基本的なベンチマークを示しています。実装指向の評価を検討する前に、これらの理論設計には厳密なシミュレーションベースの中間試験が必要です。したがって、このプロトコルは、固定ステップのシミュレーション環境におけるバックステップ追跡コントローラに対するパーティクルスウォーム最適化(PSO)による性能変化を体系的に評価するための標準化された迅速な制御プロトタイピングシミュレーションフレームワークを確立します。

より広範な文献では、回転倒振り子の制御を改善するための構造的修正やパラメータ調整を探る研究が数多く行われています。最適化ベースの設計が広く見られます。例えば、PSOはコントローラパラメータの選択に利用され、ファジィハイブリッド制御アーキテクチャ^4,5に統合されています。他のバイオインスパイアメタヒューリスティックと比較して、PSOは低次元連続探索空間での迅速な収束性と最小限のハイパーパラメータチューニング要件から、このフレームワークに特に選ばれています。近年の文献では、多様で複雑な制御シナリオにおけるインテリジェントな最適化アルゴリズムの必要性がますます強調されています。例えば、高度な最適化技術はモデル参照適応制御(MRAC)や分数順序フレームワークと効果的に組み合わされ、非線形サーボプラントの追跡精度を向上させています^6,7。さらに、最適化ベースのチューニングは、複雑な電気機械システムの結合動力学や固有の制約を管理する上で非常に有利であることが証明されています ^8,9。

最近の研究では、PSOおよびそのハイブリッド型が太陽光発電アレイにおける最大パワーポイント追跡の有効性を大幅に向上させることが検証されており、部分遮^{蔽条件下での}堅牢なパラメータ識別が示されています。ロボティクスにおいて、PSOは並列マニピュレータの拡張線形および非線形比例微分制御設計の最適化に成功裏に活用され、軌道^追跡誤差を最小化しています。さらに、PSOと適応型バックステップスライディングモード制御の統合は、空気圧式人工筋肉作動の懸垂質量における振動抑制に極めて重要であることが証明^{されています。}基本的なパラメータ選択を超えて、現代の信号処理とロバスト最適化戦略の統合は、現実的でノイズの多い物理条件下でのクローズドループ安定性を維持するために不可欠です^13,14。

最近の研究では、不確実条件下での同時ジョイント角度トラッキングと振り子安定化にも焦点が当てられています。これは、ロバスト一般化動的反転¹⁵および適応神経推定¹⁶を用いています。さらに、バックステップ制御アーキテクチャは、建物の振動抑制のためのスライドモード設計の組み込み、電子スロットルバルブに準スライド観測器を活用し、高精度直流モーター速度調整のための非線形擾乱観測器の統合など、さまざまな機械システムにわたる複雑な擾乱に対応するために広範に進化しています 17,18,19.これらの多様な応用は、バックステップ設計の多様性と堅牢な推定や最適化戦略の組み合わせを強調しています。

現代の文献で明らかな重要な弱点は、構造制御器の変更とパラメータ調整の利点を混同していることです。多くの比較研究では全く異なる制御アーキテクチャが比較されており、性能向上が基本的なアルゴリズムによるものか、単に優れた利得選択によるものかを判断することは不可能である(^20,21)。さらに、既存のシミュレーション研究は理想的な運転条件を前提とし、基本的な安定化のみに重点を置いています。現実的な実装制約によるパフォーマンス低下に対応できないことが多いです。この方法論はこれらのギャップに直接対応しています。動的軌道追跡パラダイム内でシミュレーションされたサンプリング遅延、センサーノイズ、減衰ミスマッチを導入することで、提案されたプロトコルは固定制御構造に対するPSOの最適化効果を評価します。

従来の数値積分とは異なり、このプロトコルはパラメータ調整の効率を構造制御者の変動から切り離しつつ、固定ステップのタイミング制約を強制することで差別化されます。新しい制御アーキテクチャを導入するのではなく、この手法は固定ステップのリアルタイムシミュレーション環境内で単一のバックステップ追跡コントローラを評価するものです。手動で調整したベースラインと、まったく同じコントローラのPSO最適化バリアントと比較することで、ステップトラッキング、正弦波トラッキング、妨害除去の観測された違いを利得最適化プロセスに帰属させるよう設計されています。本研究の主な貢献は、PSOベースの利得選択の影響を構造制御者の変動から分離するベンチマーク指向のリアルタイムシミュレーションプロトコルの開発です。

具体的には、この作業は以下の通りです:(1) 実装制約をエミュレートするための2msの固定ステップ実行環境を確立すること;(2) ステップ、混合周波数正弦波、トルクパルス擾乱シナリオを含む多面的な評価スイートを統合;(3) 標準化された非理想的なシミュレーション条件下でのゲインチューニング効果を比較するための定量的なベンチマークを提供します。この手法は、固定ステップタイミング、測定ノイズ、遅延、減衰ミスマッチ、飽和制約下で非線形制御則の中間シミュレーション-評価段階を必要とする制御研究者やシステムエンジニア向けに特別に設計されています。これは、追跡精度と内部状態安定性が現実的な測定ノイズや通信遅延の下でバランスを取る必要がある、過小作動電気機械システムに特に適用されます。

プロトコル

このプロトコルは人体被験者、動物実験、臨床標本を含みません。手順は、非線形電気機械制御システムを表現する固定ステップシミュレーション環境内で完全に実行されます。本研究では、物理的な回転倒振り子実験、ハードウェアインザループの検証、物理的な展開試験は行われていません。

1. プラント建設および信号コンベンションの設立

1. 古田型回転倒振り子の高精度シミュレーションモデルを構築します。系の2つの一般化座標を定義します:回転アーム角θと、垂直方向からの振り子偏差角α。
2. 系状態ベクトルを x = [θ, α, θ̇, α̇]T と定義し、α = 0 を不安定な直立平衡とします。オイラー・ラグランジュ方程式を用いて、非線形の慣性行列、コリオリ力、そして能動回転アームと受動振り子を結合させる重力ベクトルを明示的に定義する数学的プラントモデル^{を作成します。}
3. プラント、コントローラー、ポストプロセッシングスクリプト全体で、同じエンコーダ極性を一貫して使用し、正のθと正のαを定義してください。
4. 表1に詳述された名目上の機械的および電気機械的パラメータを厳密に用いて、固定ステップのリアルタイムシミュレーションソルバーを確立します。材料表で指定されたリアルタイム実行環境内でクローズドループシステムを構成します。コントローラの更新期間を2msに設定し、データログ間隔を20msに固定してください。
5. プラントモデルに3つの特定の非理想的効果を取り入れ、現実的な実装上の制限を模倣します。アームおよび振り子フィードバックチャネルに、それぞれ標準偏差0.003ラジアドと0.004ラジットのゼロ平均ホワイトノイズを注入します。
6. 測定された振り子信号に4msの転送遅延を導入します。同時に、モデリングミスマッチをシミュレートするためにコントローラ設計で用いられる公称値に対してシミュレートされたプラント減衰係数を8%増加させます。
7. モーターコマンドの出力は厳密に±10Vに制限してください。
   注意:|θ|0.70 rad を超える場合、または |α|放射能は0.35を超えています。さらに、制御コマンドが±10V飽和制限に100ms以上維持された場合にシステムリセットをトリガーするウォッチドッグタイマーを実装してください。

2. 初期状態の条件と前置詞

1. 各テストの前にシステムを下向き静止状態に初期化します。
2. 振り子を逆さまの位置に振り上げ、安定制御装置を作動させる事前配置タスクを実行します。この前置付け段階は正式なパフォーマンス評価から除外してください。
3. 振り子偏差が少なくとも1.0秒にわたり直立平衡から±0.05ラジアド以内に持続していることを確認し、すべての評価試験で動的初期条件をほぼ同等に保つことを確認してからのみ正式なデータ収集を開始します。
4. この条件を満たしたらトライアル時間を0.0秒にリセットし、指定された設定ポイントプロファイルの実行を直ちに開始します。すべての管制官が同じ事前配置手順を使用し、同じ受理基準を満たすことを確認してください。

3. ベースラインのバックステップコントローラ実装

1. 追尾誤差 eθ(t) は、命令された腕角参照 θr(t) に従って方程式を用いて定義します。
   eθ(t) = θr(t) − θ(t) (1)
2. 式2に従って複合滑り面s(t)を構築します。この定式化の中で、λを厳密に正の誤差曲面傾き係数として定義します。
   s(t) = ėθ(t) + λeθ(t) (2)
3. 基準的なモデルベースの補償項と誤差安定化比例項を積分して参照バックステップコントローラを作成します。追跡誤差変数を数学的に z1 = eθ(t) と定義し、z2 = s(t) と定義します。
4. トラッキング誤差収束を評価するためにリアプノフ安定性解析を実施します。選択された制御リャプノフ関数を時間に対して微分し、V̇ = z1 + z2ż2 を得ます。この微分に定式化された制御則を代入し、V̇ ≤ −k1z1² − k2z2² ≤ 0であることを検証します。この負半定定条件は理論的に閉ループ系の誤差が原点に収束することを保証します。
5. このリアプノフに基づく安定性条件を保持するために選択された境界内で多次元パラメータ空間を探索するよう、次のPSOアルゴリズムを制約します。
6. 標準的な境界層バックステップ^{の実践に従い}、平衡付近で高周波のチャタリングを軽減するために、不連続符号関数を連続飽和関数sat(x/φ)に置き換えます。φを式3内の制御遷移を平滑化するために用いられる境界層の幅と定義します。
   sat(x/φ) = 1 に対して x/φ > 1;SAT(x/φ) = x/φ |x/φ|≤ 1;sat(x/φ) = −1 x/φ < −1 (3)
7. 位置測定から角速度を推定するには、カットオフ周波数25Hzの一次フィルタリング後方差分アルゴリズムを適用します。
8. 表2に記載された固定ベースラインゲイン(k1 = 3.60、k2 = 1.95、λ = 2.10、φ = 0.12)を適用します。正式なテストが始まったら、これらの値を調整しないでください。

4. 粒子群最適化(PSO)によるバックステップゲイン最適化

1. リアルタイム評価の前にオフラインチューニングを実行してください。20個のパーティクルで構成されたPSOを用いて4つのコントローラパラメータ(k1、k2、λ、φ)を最適化し、実行回数を35回²⁴回に厳密に制限します。この閾値で最適化を終了し、不要な計算オーバーヘッドを防ぐべきです。経験的なオフラインチューニングでは、群れが最初の20回の反復で安定したグローバル最適適合度値に一貫して収束することが示されていました。
2. 指定された探索境界内で粒子を一様に初期化します:k1 ∈ [2.5, 5.5]、k2 ∈ [1.2, 3.0]、λ ∈ [1.2, 2.8]、φ ∈ [0.05, 1.20]。
3. 20秒のオフライントレーニングシナリオで各粒子を評価します。この評価環境が、正式なリアルタイム試験で定められた飽和限界、センサーノイズ、遅延、減衰ミスマッチを厳密に再現すること。
4. 最小化目的関数を J = 0.55RMSE + 0.25max|α| と定式化します。+ 0.20RMSU。
5. これら3つの目標成分を固定ベースライン参照値で正規化し、スウォームサーチ開始前に実行された5回の初期ベースラインランの平均として計算します。これらの正規化定数は最適化プロセス全体を通じて厳格に維持してください。
6. 最適化の過程で慣性重みを0.90から0.40に線形に減らします。両方の加速度係数を1.50に固定してください。
7. 評価段階で安全停止を引き起こす粒子構成にはあらかじめ決められたペナルティコストを割り当てます。
8. 最終イテレーション完了後、グローバルベスト解を直接抽出し、それを最適化コントローラーとして指定します。
9. 表2に詳述された最適化されたゲインを、すべての正式なPSO最適化リアルタイムトライアル(k1 = 4.19、k2 = 2.44、λ = 1.83、φ = 0.92)に実装します。

5. リアルタイム実行プラットフォームの設定

1. 制御ロジックと数学的プラントモデルを、迅速制御プロトタイピング(RCP)固定ステップリアルタイムシミュレーション環境内で同時に実行します。リアルタイムのターゲットマシンとして十分な計算能力を持つ高性能な産業用ワークステーションを活用しましょう。制御アルゴリズムを仮想アナログ出力およびエンコーダ入力チャネルを通じてシミュレートされたプラントダイナミクスと連携させ、実装レベルの信号制約をエミュレートします。
2. リアルタイムステップサイズに対して厳密な2ms間隔を強制します。ハードウェアおよび実行設定の包括的な要約は 表3 を参照してください。
3. 各トライアルを開始する前に、すべてのコントローラーの状態をリセットしてください。
4. 参照値をゼロに設定し、直立した位置で2.0秒のベースラインデータセグメントを収集します。指定された追跡プロファイルを開始する前に、このウィンドウを点検して異常な飽和イベントやセンサードリフトがないか確認してください。
5. 各新しい試行ごとに、あらかじめ定義されたシード配列から選ばれた異なるシードを用いて測定ノイズの実現を再生成します。

6. ステップトラッキングテストの実行

1. 0.0秒から2.0秒までθr = 0弧度、2.0秒から10.0秒までθr = 0.50弧度、10.0秒から20.0秒の間にθr = 0.20弧度という20秒の区分的ステップシーケンスを適用します。
2. ベースラインコントローラで10件の有効な試験を、PSO最適化コントローラで10件の有効な試験を実施します。
3. 経過時間、基準角、測定されたアーム角、トラッキング誤差、振り子偏差、制御電圧を各試行中に継続的に記録します。
4. 0.50ラジッドステップ入力後の暫過性能を特に評価し、t = 2.0秒で発生します。
5. 上昇時間(測定されたアーム角の応答が目標設定値の10%から90%に移行するのに必要な時間)を計算します。
6. オーバーシュートは、ピーク角偏差が0.50ラジアット目標を上回る割合として計算します。
7. 安定時間を計算します。これは、応答が0.50ラジヤ目標の±2%以内に最小1.0秒の範囲内に収まる最初のインスタンスとして決定されます。
8. 10.0秒前にこの和解条件を満たさなかった裁判は未解決とマークしてください。決済時間平均計算からは厳密に除外し、他のすべての指標については保持してください。
9. 定常状態誤差は、9.0秒から10.0秒までの時間的ウィンドウで追跡誤差を平均して計算します。

7. 正弦波トラッキングテストの実行

1. 式4で制御される20秒の混合周波数正弦波基準軌道を注入します。
   θr(t) = 0.26sin(2π·0.10t) + 0.12sin(2π·0.30t + 0.40) rad (4)
2. 各コントローラー構成ごとに有効なトライアルを10回クリアしてください。
3. 各試行ごとに、全試行追跡RMSE、最大絶対トラッキング誤差、制御RMS電圧、最大絶対振り子偏差、位相遅延を抽出します。
4. 測定された腕角軌道と指令参照点を±0.50秒の指定された探索ウィンドウ内で相互相関して位相遅延を推定します。測定出力が基準信号に時間的に追従する際の遅延を正の値として記録します。

8. 撹乱拒否テストの実行

1. アーム角度の基準は20秒の間ずっと0ラジヤに厳密に保ちます。
2. ちょうど6.0秒で加算モーター側トルクパルスを注入し、12.2秒で同じパルスをもう一回加えます。
3. 両方の注入イベントで擾乱振幅を0.030 N·mに設定し、パルス幅を0.12秒に設定します。各コントローラーで有効なトライアルを10回クリアしてください。
4. 各パルス開始直後の0.50秒のウィンドウ内のみのピーク絶対トラッキング誤差を計算します。その試行における代表的な撹乱後のピーク誤差として、2つの値のうち大きい方を特定し保持します。
5. 回復時間を計算します。これは、パルス開始から最初の記録されたサンプリング瞬間(|eθ|≤ 0.010 放射能。サブサンプルの回復時間は補間しないでください。データログは厳密に20ms間隔で行われます。
6. 絶対トラッキング誤差を積分し、制御RMS電圧を20秒の試験運転で累積的に計算します。

9. データエクスポートと統計的要約

1. ステップ、正弦波、撹乱の生の時間領域データセットを別々のデータファイルにエクスポートします。
2. 繰り返し試験の指標を構造化されたマトリックス形式に整理し、各試験ごとに1行、評価された指標ごとに1つの列を割り当てます。
3. 未完了または安全停止で終了した試験は厳格に除外してください。除外されたすべての実行を、同一条件下で追加の試験を行うことで、各テストシナリオごとに各制御者あたり10件の有効な試験を含む最終データセットを保証します。
4. 集計された再試験結果を平均±標準偏差として報告します。
5. 抽出した指標の分布正規性をShapiro-Wilk検定²⁵を用いて評価します。
6. 正規分布データに対して、ベースラインと最適化されたコントローラ間の統計的差異を、両側独立標本t検定を用いて解析します。非正規分布を示すパラメータに対してマン・ホイットニーU検定を適用します。
7. 統計的有意の閾値をp < 0.05に定めます。

結果

粒子群最適化(PSO)アルゴリズムは、複合目的関数の初期急速な減少と、その後徐々に収束する期間を示しました。具体的には、最適適応度値は最初の計算サイクルで1.8757に低下し、平均群適応度は2.0573から35回目の反復で1.2518に低下しました。この収束の大部分は、最初の15〜20回の反復で起こりました。この段階を超えると、グローバル最適解の軌道が収束し、群れがランダムにさまよいながら、あらかじめ定められた探索空間内の安定した領域に到達したことを示しています。最終的なPSO由来のゲインセットは、ベースラインバックステップ法則の構造的完全性を維持しつつ、安定化項間での相対的な重み付けを再分配しました(表2)。図1に示すように、最適化アルゴリズムは適応度関数の急速な減衰を促進し、初期探索段階から安定したグローバル最適解へ35回の反復で移行しました。手動で調整された基準線と比べて、最適化されたコントローラはk1およびk2の値が大きくなり、誤差面傾きλが減少し、境界層φが広がりました。

ステップトラッキング性能

比較分析により、ステップトラッキング実験中の両コントローラ間で明確な性能差が明らかになり、このプロトコルの期待される正の結果と最適でない結果が示されています。ベースラインのバックステップ制御は典型的な最適でない結果を表します。2.0秒に注入された0.50 radコマンド後も全体的な安定性を維持しましたが、初期上昇率が急激で目標軌道から早期逸脱し、持続的な残留振動が発生し、収束に至るまで長時間の時間が必要となりました。一方、PSO最適化されたコントローラは、コマンド飽和をトリガーせずに設定値に迅速かつ減衰したアプローチを取るという代表的な良好な結果を示しています。10回連続の試験で、最適化されたコントローラはステップトラッキングRMSEを0.065±0.002ラジットから0.050±0.002ラジットに低下させました(Welch t-test, p < 0.001, 表4)。過渡性指標も改善され、オーバーシュートは4.644±0.604%から3.627±0.620%に減少しました(Mann–Whitney U, p = 0.009)、上昇時間は0.788±0.025秒から0.025秒±0.548に短縮されました(Mann–Whitney U, p < 0.001)。これらの性能差は図2および図3に示されており、最適化された応答により減衰が改善され、設定値収束が速くなります。この精度向上とともに、平方根制御電圧が0.669±0.004Vから0.521±0.003Vへと22.1%削減されました。同時に、最大振り子の揺れは33.1%減衰し、5.812°から3.889°に低下しました(表4参照)。さらに、10.0秒のウィンドウ内で1.0秒の沈殿基準を満たすベースライン試験はなかったものの、最適化されたコントローラーは70%の試行で成功し、平均落ち着時間は1.057±0.348秒を記録しました。

正弦波トラッキング性能

両制御戦略とも、20秒の試験期間中、混合周波数正弦波基準の下で有界追跡を維持しました。パフォーマンスギャップはステップ応答テストほど視覚的に顕著ではありませんでしたが、繰り返しの試験を通じた統計結果は最適化されたアプローチでパフォーマンスが向上したことを示しました。ベースライン制御信号はトラッキングエラー範囲が広く、ターゲット波形に比べて明らかに遅れていました。対照的に、最適化されたコントローラは低周波および高周波の変動において設定値をより正確に追跡し、基準軌道に沿ったピーク誤差の分布をより厳密に実現しました。その結果、最適化後のフルトライアルRMSEは0.061±0.003radから0.040±0.003radに減少しました。 表5に詳述されているように、最適化されたコントローラは制御RMS電圧を0.557 ± 0.002 Vから0.402 ±0.002 Vに低下させ、最大振り子変位を2.156 ± 0.079°から1.464 ±0.061°に制限し、位相同期とエネルギー効率の向上を示しました(p < 0.001、 表5; 図4)。

撹乱の除去とクロスシナリオの堅牢性

6.0秒および12.2秒の加算トルクパルス注入により、ベースラインコントローラーの応答にアーム角度の偏差が大きくなり、振り子の動きが広がりました。最適化されたコントローラーも摂動を即座に認識しましたが、摂動後の変動は低くなっていました。ゼロ基準均衡への回復過程は質的により速い傾向を示しました(0.014秒から0.000秒へ)。しかし、この特定の指標は統計的有意性を達成せず(p = 0.078)、主に20msのデータログ解像度が測定可能なフロア効果を誘発したためです。繰り返しの試行を経て、最適化されたパラメータ化により摂動後のピーク誤差が減少しました。20秒の持続時間で計算された絶対誤差(IAE)の積分は0.377±0.012から0.248±0.015に減少し、制御RMS電圧も0.656±0.003Vから0.508±0.005Vへと低下しました。振り子振動の最大振幅も同様に11.456±0.117°から7.561±0.063°に減衰しました(表6)。コントローラーのランキングは3つのテストプロファイルすべてで一貫していました。PSO最適化されたバックステップ戦略は、過渡誤差の低減、定常状態の不正確さの低減、受動振り子の移動の最小化を実現しました(p < 0.001、 表6; 図5)全体のアクチュエーター負荷増加を要求することなく。一次評価指標の試験間分布は 図6にまとめられています。

データの利用可能性:

代表的なデータ、試験レベルの指標、あらかじめ定義されたランダムシード、PSOチューニングレコードを含む図や表の作成に使用されたサポートデータセットはZenodoに保存されており、DOI 10.5281/zenodo.20395966のもとで一般公開されています。

図1:コントローラチューニングにおける粒子群最適化の収束。(A) 35回の反復における群れの最適適応度値;(B) 各時間点の平均適応度値。複合物対物レンズは初期探査段階で急速に衰減しましたが、PSO最適化されたバックステッピング制御則の安定した領域に収束するにつれて徐々に減速しました。 この図の拡大版はこちらをクリックしてご覧ください。

図2:ベースラインおよびPSO最適化バックステーピングコントローラの代表的なステップトラッキング応答。 (A) 20秒の全ステップシーケンスにおける指令アーム角度基準およびアーム角応答の測定;(B) 0.50ラッドステップ直後に2.0秒で応答を拡大表示;(C) 2回目の0.20ラッド指令セグメント後の応答拡大表示は10.0秒から始まる。最適化後、PSO最適化コントローラはベースラインコントローラに比べて上昇時間の速さ、ピーク誤差の低さ、優れた減衰性能を示しました。 この図の拡大版はこちらをクリックしてご覧ください。

図3:2つのコントローラ下でのアーム角トラッキング誤差と制御挙動。 (A) ステップトラッキングテスト中のアーム角トラッキング誤差;(B) 同じ試行のための制御電圧信号;(C) ステップ応答時の振り子偏差角。最適化されたコントローラーは誤差の変動が小さく、振動範囲が狭く、低いRMSコマンドレベルで比較的滑らかな動きを示しました。 この図の拡大版はこちらをクリックしてご覧ください。

図4:混合周波数基準下での代表的な正弦波トラッキング性能。 (A) 20秒トラッキング実験内で正弦波参照を指令し、腕角応答を測定した;(B) 同じ時間範囲内での追跡誤差;(C) 正弦波追尾過程における振り子の偏角。PSO最適化されたコントローラは、ベースラインコントローラよりもコンポジットコマンドをより正確に追跡し、位相遅延の低減と連続サイクルでの誤差の範囲の縮小を示しました。 この図の拡大版はこちらをクリックしてご覧ください。

図5:ゼロ参照調節下での代表的な撹乱-拒絶応答。 (A) 擾乱除去試験中のアーム角追尾誤差(6.0秒および12.2秒でのトルクパルス擾乱を含む);(b) 同じ試験中の振り子偏差角;(C) 擾乱除去中の電圧信号を制御する。ベースラインコントローラと比較して、PSO最適化制御は撹乱後の最大偏差が低く、設定値への回復も速いことが示されました。 この図の拡大版はこちらをクリックしてご覧ください。

図6:3つのテスト状況すべてにおける一次評価指標のクロストライアル分布。 (A) ステップトラッキングRMSE;(B) ステップトラッキングの上昇時間;(C) 正弦波追跡RMSE;(D) 正弦波位相遅延;(E) 撹乱後のトラッキングエラーのピーク;(F) 擾乱・拒否条件下での最大振り子変位。各パネルは各制御条件下での10件の有効な試験の分布を示しています。ベースラインおよびPSO最適化されたバックステップ制御システムの中心傾向と分散が直接比較されます。 この図の拡大版はこちらをクリックしてご覧ください。

パラメータ	シンボル	価値	部隊
ロータリーアーム長	Lr	0.215	m
振り子重心長	Lp	0.168	m
ロータリーアーム質量	MR	0.254	KG
振り子質量	mp	0.097	KG
アーム粘性減衰	Br	0.0031	N·m·s/rad
振り子粘性減衰	Bp	0.0018	N·m·s/rad
モータートルク定数	Kt	0.053	N·m/A
エンコーダの解像度	Nエンc	4096	カウント/レヴ
コントローラー更新期間	Ts	0.002	s

表1:非線形回転倒振り子モデルの公称系パラメータ。 リアルタイム非線形プラントモデルおよび制御システム開発のための標準的な幾何学的、慣性的、減衰、作動、センサーのパラメータリスト。正式な試験適用時に追加される減衰ミスマッチは独立であり、これらの名目数値から除外されます。

コントローラー	ゲイン	価値	概要
ベースラインバックステップ	K1	3.6	手作業で調整された初期安定化ゲイン
ベースラインバックステップ	K2	1.95	手作業で調整された仮想制御ゲイン
ベースラインバックステップ	λ	2.1	誤差面傾き
ベースラインバックステップ	φ	0.12	境界層幅
PSO最適化されたバックステッピング	K1	4.19	35回の反復後の全域最良粒子
PSO最適化されたバックステッピング	K2	2.44	35回の反復後の全域最良粒子
PSO最適化されたバックステッピング	λ	1.83	35回の反復後の全域最良粒子
PSO最適化されたバックステッピング	φ	0.92	35回の反復後の全域最良粒子

表2:ベースラインのバックステッピングコントローラーゲインおよびPSO最適化後のバックステッピングコントローラーゲイン。 表は両コントローラの動作状態に対する利得の組み合わせを示しています。基準値は手動調整時に固定され、PSO最適化された利得とは、35回の最適化反復後に保持された最終的なグローバル最適粒子を指します。

アイテム	仕様
リアルタイムターゲットコンピュータ	産業用PC、Intel i7クラスCPU、16GB RAM
実行モード	固定ステップリアルタイムシミュレーション
コントローラー更新期間	2ミリ秒
信号輸出間隔	20 ms
信号インターフェース	エンコーダー入力およびアナログ出力、±10V
プラント構成	非線形古田型回転倒振り子モデル
注入された非理想的効果	測定ノイズ、振り子チャネル遅延、減衰ミスマッチ
測定ノイズ	アーム角に0.003の放射線;振り子角での0.004放射度
振り子チャネル遅延	4ミリ秒
プラント減衰ミスマッチ	名目上のコントローラモデルに対して+8%
コマンド飽和	±10 V
安全停止基準	|θ|> 0.70 rad, |α|> 0.35 rad、または飽和度 > 100 ms
事前配置受理	|α|≤ トライアル開始の少なくとも1.0秒前から0.05の放射線を保持してください
ステップ参照プロファイル	0 rad(0.0-2.0秒)、0.50 rad(2.0-10.0秒)、0.20 rad(10.0-20.0秒)
正弦波基準プロファイル	0.26 sin(2π·0.20t) + 0.12 sin(2π·0.30t + 0.40) rad
妨害プロトコル	6.0秒と12.2秒の2つの加算トルクパルス
擾乱振幅と幅	0.030 N·m、0.12秒
コントローラーごとの有効なトライアルとシナリオ	10

表3:リアルタイムシミュレーション環境、信号設定、試験資格ルール。 実行プラットフォームの概要、コントローラのリフレッシュレート、信号の出力間隔、非理想的な効果の注入、安全停止条件、参照軌道、妨害設定、有効な試験基準。

メートル法	ベースラインバックステップ n=10	PSO最適化バックステップ n=10	p値	統計的検定
RMSE(放射能動員)の追跡	0.065 ± 0.002	0.050 ± 0.002	<0.001	ウェルチt検定
定常状態誤差、9.0–10.0秒(rad)	-0.007 ± 0.001	-0.004 ± 0.001	<0.001	マン・ホイットニー大学
オーバーシュート(%)	4.644 ± 0.604	3.627 ± 0.620	0.009	マン・ホイットニー大学
立ち上がり時間(複数)	0.788 ± 0.025	0.548 ± 0.025	<0.001	マン・ホイットニー大学
落ち着時間(s)	有効な合着裁判なし(0/10)	1.057 ± 0.348(7/10の有効な和解裁判)	NA	比較しない
制御RMS(V)	0.669 ± 0.004	0.521 ± 0.003	<0.001	ウェルチt検定
最大振り子偏差(デグ)	5.812 ± 0.156	3.889 ± 0.109	<0.001	ウェルチt検定

表4:1コントローラーあたり10件の有効な試行におけるステップトラッキングパフォーマンス指標。 この表は、RMSE、定常状態誤差、オーバーシュート、上昇時間、設定値収束速度、制御RMS電圧、最大振り子偏差の追跡を含むステップトラッキング実験の繰り返し試験の要約統計を示しています。値は特に記載がない限り標準偏差±平均として報告されます。

PSO最適化バックステップ n=10	p値	統計的検定
0.040 ± 0.003	<0.001	ウェルチt検定
0.088 ± 0.006	<0.001	ウェルチt検定
0.232 ± 0.019	<0.001	ウェルチt検定
0.402 ± 0.002	<0.001	ウェルチt検定
1.464 ± 0.061	<0.001	ウェルチt検定

表5:10件の有効な試験における正弦波追跡性能指標(1制御者あたり)。 この表は、正弦波トラッキングテストの繰り返し試行要約を示しており、全試行平均平方根誤差(RMSE)、最大絶対追跡誤差、位相遅延、制御RMS電圧が含まれます。最大の振り子偏差も提供されています。正の位相遅延は、測定された出力が設定値に追随していることを示します。

p値	統計的検定
<0.001	ウェルチt検定
0.078	マン・ホイットニー大学
<0.001	ウェルチt検定
<0.001	ウェルチt検定
<0.001	ウェルチt検定

表6:10件の有効な試験における制御者ごとの妨害拒否性能指標。 擾乱除去実験の繰り返し試験データには、摂動後のトラッキング誤差のピーク、回復時間、絶対誤差積分(IAE)、制御RMS電圧、最大振り子偏差が含まれます。信号のエクスポート設定により、測定可能な回復時間は20msのログ期間によって制限されます。

ディスカッション

本研究の主な目的は、単に孤立した条件下であるコントローラが優れていると宣言するだけでなく、ここで評価したPSOベースの利得最適化フレームワークが、定義されたシミュレーションされた非理想的シナリオにおいてバックステップコントローラの性能を向上させることを示すことです。回転倒振り子は、その高度に非線形で非最小位相、そして作動不足^{の特性を持つ}ため、この評価の優れたベンチマークとなります。パフォーマンス指標を超えて、このプロトコルの技術的整合性はリアルタイムシミュレーション環境の意図的な特性評価にかかっています。重要な実装ステップとして、8%の減衰ミスマッチと4msのセンサー遅延を制御注入することが挙げられます。これらの特定の非理想的効果により、理想的条件だけでなく形式評価で用いられる摂動と同じ摂動下でPSOアルゴリズムがゲインを最適化できるからです。

ステップトラッキング、正弦波トラッキング、擾乱除去プロトコルを統一固定ステップ実行設定でテストした結果、シミュレーション結果は一貫したパフォーマンス傾向を示しました。特に、ステップおよび正弦波応答の両方で観察された強化、特に追尾誤差の低減、オーバーシュートの最小化、位相遅延の減衰が、RMS制御電圧の測定可能な低下と同時に達成されました。これは、最適化プロセスが単にドライブユニットを飽和させることで追跡精度を向上させたのではなく、むしろ、既存のリャプノフに基づくバックステップの枠組み内で制御力を再分配したのです。このパラメータ再分配の一貫性は、プロトコルStep 2.4で定義された前置受け入れ基準によって大きく影響されます。振り子が指定された平衡ウィンドウ内に長期間保持されることを厳密に保証することが不可欠です。そうでなければ、不規則な初期過渡現象が、その後のオフライン調整フェーズにおける早期のコスト関数停滞や非収束群行動を引き起こす可能性があります。過小作動機械システムにおいて、この再分配はアクティブ座標(回転アーム)と受動座標(振り子)の結合管理に制御努力を割り当てる方法を決定するため重要です²⁷。

さらに、アーム角度追尾誤差と内部振り子の移動が同時に減少することは、提案されたチューニングプロトコルの利点を示しています。欠作動システムの制御において、揮発性の高い内部状態を犠牲にして正確な座標追跡を達成することは、一般的に最適とは言えないか失敗した設計と見なされます。最適化された制御則により、作動していない振り子の反作用振幅を33.1%削減しつつ、追尾精度も向上しました。この複合的な改良により、PSOアルゴリズムは高利得フィードバックに依存しない寄生振動よりも、アクティブ座標と受動座標間の非線形結合をよりバランスよく調整する利得の組み合わせを特定したことが示唆されます。このバランスの取れた性能は、純粋な軌道追跡よりも制御効率の指標として有用であり、エネルギー^{変換過程中に}内部状態の動態が十分に境界が確保されていたことを示しています。同様に、撹乱除去結果は、シミュレーションされた試験条件下での最適化された利得の実用的な堅牢性を支持しています。20msのロギング間隔は正確なサブサンプル回復時間を曖昧にするフロア効果を誘発しましたが、同時にピーク後の擾乱誤差、累積絶対誤差、最大振り子変位の減少は、最適化されたパラメータが外部パルス注入のシステムダイナミクス伝播度を制限していること^{を示しています。}

このプロトコルの成功裏の実行には、一般的な実装ボトルネックへの対処が必要です。例えば、群れの軌道が35回の反復で安定しない場合、研究者はまずプロトコルステップ4.2で定められた探索境界が安定リャプノフ領域を十分に含んでいるかを確認する必要があります。さらに、リアルタイムターゲットマシンが妨害除去フェーズ中に頻繁に安全停止をトリガーする場合、プラットフォームがサポートしている場合はコントローラ更新期間を短縮したり、逆差フィルタのカットオフ周波数を改良することで、最適化ロジックを損なうことなくクローズドループの整合性を維持することができます。

これらの有望な結果にもかかわらず、この手法はシミュレーション環境に内在する抽象性に制約されています。実装されたプロトコルは、サンプリング遅延、コマンド飽和、減衰ミスマッチ、測定ノイズなど、純粋なオフライン数値積分ではしばしば無視される要素を考慮していました。しかし、本研究には物理的な回転倒振り子実験、物理的なハードウェアインザループ検証、物理的な展開試験は含まれていません。したがって、これらの結果は、指定された固定ステップおよび非理想的なシミュレーション条件下でのPSOベースのゲインチューニングのパフォーマンス効果を示すシミュレーションベースの証拠として解釈されるべきです。将来の研究では、シミュレーション内で摩擦、アクチュエータ、コンプライアンス、量子化効果を組み込むことで、モデル化された不確実性セットを拡張し、ハードウェア指向の個別研究^{が検討される}前に検討される^{可能性があります。}

本研究は、バックステッピング制御のためのPSOベースの利得最適化を評価するための迅速制御プロトタイピングシミュレーションプロトコルを確立します。10回の繰り返し試験の定量的結果により、最適化されたコントローラはステップトラッキングRMSEを0.065から0.050radに低下させ、手動調整ベースラインと比べてRMS制御負荷を22.1%削減することが示されています。振り子の移動における33.1%の減衰は、模擬された非理想的条件下での内部状態安定性の向上を示しています。このシミュレーション環境はノイズ、遅延、減衰ミスマッチ、飽和を含みますが、物理システムで見られる複雑なクーロン摩擦、モーターデッドバンド、構造コンプライアンス、センサー量子化を完全に再現することはできません。今後の研究は、シミュレーションモデルを拡張して追加の非理想効果を含めることや、同じ順位付けパターンがより広範なシミュレーション不確実性条件で安定しているかどうかの検証に焦点を当てる予定です。

材料

この記事で使用された材料の一覧

名前	会社	カタログ番号	コメント
High-performance Workstation (Windows 11 Pro)	Various / Custom Build	N/A	固定ステップのリアルタイムシミュレーションと後処理のためのホストマシン。
MATLAB (Version R2024a)	MathWorks	https://www.mathworks.com/products/matlab.html	非線形プラントモデリング、コントローラーコーディング、データエクスポート、パラメータ管理。
Simulink (Version R2024a)	MathWorks	https://www.mathworks.com/products/simulink.html	Furuta型振り子のブロック図モデル構築とコントローラー実行。
Simulink Desktop Real-Time (Version R2024a)	MathWorks	https://www.mathworks.com/products/simulink-desktop-real-time.html	コントロールモデルのデスクトップ実行のための固定ステップリアルタイムカーネル。
Python (Version 3.11)	Python Software Foundation	https://www.python.org/	二次的データ処理、統計処理、図の準備。
NumPy (Version 1.26)	NumPy Developers	https://numpy.org/	エクスポートされた試行データの数値配列操作。
pandas (Version 2.2)	pandas Developers	https://pandas.pydata.org/	繰り返し試行データの整理とサマリーテーブルの生成。
SciPy (Version 1.13)	SciPy Developers	https://scipy.org/	統計テストと信号分析ユーティリティ。
Matplotlib (Version 3.8)	Matplotlib Developers	https://matplotlib.org/	収束、トラッキング、分布図の生成。