Polynomial Chaos-based Input Shaper Design under Time-Varying Uncertainty

技術要約：時変不確実性下における多項式展開を用いたインプットシェイパー設計

問題提起
不確実性の定量化（UQ）は、堅牢な工学システム、特にパラメータが時間とともに動的に変化するシステムを設計する上で極めて重要である。古典的な確率論的アプローチは、固定された分布や時不変の仮定に依存することが多く、これは無人航空機（UAV）のペイロード制御や宇宙機の姿勢制御といった複雑なシナリオにおいては限界がある。多項式展開（PCE）は、不確実性伝播のための効率的なスペクトル・フレームワークを提供するが、標準的なPCE手法は、確率密度関数が進化する長期の積分において性能が低下する場合があり、固定された直交多項式基底が最適ではなくなる可能性がある。さらに、既存のロバスト・インプットシェイピング技術は、計算コストの高いモンテカルロ（MC）シミュレーションに頼ることなく、時変不確実性を効率的に処理することに苦慮することが多い。本研究は、時変不確地性にさらされる動的システムの精密な運動制御器（具体的にはインプットシェイパー）を設計するという課題に取り組むものである。

手法
本研究では、ロバストなインプットシェイパーを設計するために、時間遅れフィルタ（TDF）を統合した**侵入型多項式展開（intrusive PCE）**フレームワークを提案する。この手法は、バネ定数（およびそれに伴う固有振動数 $\omega_n$）が時間とともに変化する区間不確実性を伴う、減衰のないバネ・質量系を用いて検証される。

1. システムモデリング： システムは、$[0, t_1]$ および $[t_1, t_2]$ の2つの時間区間にモデル化される。第1区間では、固有振動数 $\omega_n$ は一様分布に従う。$t_1$ において、システムは、固有振動数 $\omega_m$ が異なる一様分布に従う第2区間へと遷移する。極めて重要な点は、$t_1$ におけるシステムの状態（位置と速度）が第2区間の初期条件として機能することであり、これは、第2区間の解が初期の不確実性（$\zeta_1$）と新しい不確実性（$\zeta_2$）の両方に依存することを意味する。
2. 侵入型PCE定式化： 著者らは、ガラーキン投影を用いた侵入型アプローチを採用している。システム方程式は、一様分布に対するルジャンドル多項式の直交級数へと展開される。これらの多項式の係数は、結果として得られる微分方程式系の求解によって決定される。
   • 時変性を扱うために、$t_1$ におけるPCE係数は、直交条件を用いて単変数展開（区間1）から多変数展開（区間2）へとマッピングされる。これにより、遷移時における状態とその導関数の連続性が確保される。
3. インプットシェイパー設計： 本フレームワークは、残留エネルギーを最小化するために時間遅れフィルタ（TDF）を設計する。以下の3つの設計を比較する：
   • 非ロバストTDF： 公称極の位置に1組の零点を配置する。
   • ロバストTDF： 公称位置に2組の零点を配置する。
   • グローバル感度解析（GSA）TDF： 特定の目的関数 $\Xi$ を最小化するように、インパルス振幅（$A_i$）と時間遅れ（$T_i$）を最適化する。目的関数は、期待残留エネルギー（$E[V_{res}]$）または残留エネルギーの分散（$Var(V_{res})$）として定義される。

主な貢献

• 時変不確実性の処理： 本論文は、不確実性パラメータが時間とともに変化するシステムに対して、侵入型PCEを適用するための具体的な手法を導入している。これは、異なる不確実性分布間でのPCE係数の転送手順、および解の連続性を保証する詳細な数学的手順を記述している。
• 効率的なロバスト制御設計： 本研究は、侵入型PCEが、モンテカルロ・シミュレーションよりも効率的に時変不確実性を考慮したインプットシェイパーを設計できることを示している。
• グローバル感度解析の統合： 本研究は、グローバル感度解析をTDF設計プロセスに統合し、システムの期待性能または残留エネルギーのロバスト性（分散）のいずれかを標的として、シェイパー・パラメータを最適化することを可能にしている。

結果
数値シミュレーションは、質量1 kgの、時変バネ定数を持つバネ・質量系に対して行われた。

• 収束性と精度： 位置と速度の期待値および分散に関するPCEの結果は、10,000サンプルのモンテカルロシミュレーションの結果と一致した。多項式次数 $P=30$ が、MCの精度に一致するために十分であることが判明した。期待値については、次数10と30の間で有意な差は見られなかった。
• 計算効率： PCEフレームワークは、10,000サンプルのモンテカルロシミュレーションよりも約12倍高速であり、計算負荷の高い不確実性伝播における利点が強調された。
• シェイパーの性能：
  • すべてのTDF設計は、非ロバスト設計と比較して、期待残留エネルギーおよびその分散を有意に減少させた。
  • GSA TDFは、標準的なロバストTDFよりも優れた性能を示した。具体的には、分散を最適化したGSA TDF（$\Xi_2$）は、標準的なロバストTDFと比較して、期待残留エネルギーを2分の1に、分散を5分の1に減少させた。
  • 分散を最適化したGSA TDFは、異なる不確実性ペアに対してもほぼ一定の残留エネルギーを示し、システム変動に対する高い無感覚性を示した。
  • GSA TDFは一般に優れた性能を提供したが、第2区間の周波数が平均に近い特定の組み合わせにおいては、標準的なロバストTDFがわずかに低い残留エネルギーを示す場合があることが著者らによって指摘された。ただし、その差は、顕著な改善が見られた境界領域のケースと比較すると、極めて軽微であった。

意義と主張
本論文は、提案された侵入型PCEフレームワークが、時変不確実性を伴うインプットシェイパーを設計するための、計算コストが低く正確なモンテカルロ法の代替手段となることを主張している。不確実性展開をシステム方程式に直接組み込むことで、システム応答および目的関数に対するより深い解析的洞察が得られる。

著者らは、このアプローチが、モンテカルロシミュレーションが計算的に困難となる非線形システムにおける振動抑制および残留エネルギー削減において特に有利であると述べている。結論として、現在の研究は線形バネ・質量系に焦点を当てているが、本手法は非線形および高忠実度モデルへの将来的な拡張のための基礎を築くものであるとしている。また、今後の課題として、より幅広い汎用性のための非侵入型PCE手法の探索、および空中ペイロード輸送アプリケーションに向けたこれらの技術の実験的検証を挙げている。