[論文レビュー] Unit Consistency, Generalized Inverses, and Effective System Design Methods

Jeffrey Uhlmann|arXiv (Cornell University)|Apr 24, 2016

Product Development and Customization参考文献 15被引用数 6

ひとこと要約

本稿では、対角スケーリング変換の下で次元単位を保存するユニット整合性（UC）を持つ一般化された行列逆行列を導入し、工学的システム設計における一貫性を保証する。この手法により、ムーア・ペンローズ逆行列、特異値分解（SVD）、行列分解のUC版が導出され、線形システムおよび逆問題のユニット不変性分析を通じて、ロバストなシステム整合性の検証が可能になる。

ABSTRACT

This paper examines the potential role of unit consistency as a system design principle. Unit-consistent generalized matrix inverses and unit-invariant matrix decompositions are derived in support of this principle. Applications of the methods described are illustrated with examples relating to nonlinear system identification and robustness to multiplicative noise for image database retrieval.

研究の動機と目的

• ムーア・ペンローズ逆行列のような広く使われる行列ツールにおいて、単位の整合性が欠如していることによる、システム設計における単位依存バイアスの問題を解決すること。
• 任意の非特異な対角（ユニット）変換に対して不変である一般化された行列逆行列を構築し、状態空間表現における物理的単位を保存すること。
• このユニット整合性の原則を、SVD や固有値分解などの行列分解へと拡張し、モジュラー構造や複雑なシステムにおける一貫性のある解析を可能にすること。
• データ伝送におけるチェックサムに類似した、ユニット整合性テストを通じて、基準値が不要な状態で関数的誤りを自動検出する体系的な手法を提供すること。
• 混合単位や近似解を伴うシステムにおいて、ユニット整合性をシステム工学の基本的設計原則として推進すること。

提案手法

• 行列 $\mathbf{A}$ のユニットスケーリング版 $\mathbf{S}$ と、対数的行平均および列平均から導かれる対角スケーリング行列 $\mathbf{d}_L, \mathbf{d}_R$ を用いて、$\mathbf{A}^{\text{-U}} = \text{pinv}(\mathbf{S}) \odot (\mathbf{d}_L \mathbf{d}_R^T)$ で定義される新しい一般化逆行列 $\mathbf{A}^{\text{-U}}$ を提案する。
• 各行および各列の非ゼロ要素の幾何平均を1に正規化するための、Sinkhorn型反復アルゴリズムを導入し、対角スケーリング行列 $\mathbf{d}_L$ および $\mathbf{d}_R$ を計算する。
• スケーリングされた行列 $\mathbf{S} = \mathbf{D}_L \mathbf{A} \mathbf{D}_R$ を用いて、$\mathbf{U} \mathbf{S} \mathbf{V}^T$ の形でユニット不変（UI）な特異値分解（SVD）を定義する。ここで $\mathbf{D}_L, \mathbf{D}_R$ は反復スケーリング処理から得られる。
• 幾何平均正規化のためのサイズ関数 $s_{\times}(\mathbf{u}) = \left(\prod_{i} |u_i|\right)^{1/b}$ を適用し、スケーリングに対して連続的かつ同次的であることを保証する。
• スケーリングベクトル $\mathbf{u}$ および $\mathbf{v}$ を求めるために、対数変換と反復的行・列平均の差分を用い、その後指数関数を適用して $\mathbf{d}_L = \exp(\mathbf{u})$, $\mathbf{d}_R = \exp(\mathbf{v})$ を回復する。
• 実数または複素数行列に対して、$O(mn)$ の計算量で動作する $\texttt{dscale}$, $\texttt{uinv}$, $\texttt{usvd}$, $\texttt{usv\_decomp}$ のOctave/Matlab実装を提供し、数値的安定性とユニット整合性を確保する。

実験結果

リサーチクエスチョン

• RQ1 任意の対角（ユニット）スケーリング変換に対して不変であり、システムモデルにおける物理的単位を保存する一般化された行列逆行列を構築できるか？
• RQ2 ムーア・ペンローズ逆行列およびSVDを、混合単位や任意の単位を持つ工学的システムにおいてもユニット整合性を維持する形で再定義できるか？
• RQ3 ユニット整合性テストは、基準値が不要な状態で、複雑なシステムにおける関数的誤りを効果的かつ自動的に検出するための堅牢な手法として、どの程度の有効性を示せるか？
• RQ4 特にゼロ要素やスパースな要素を含む任意の行列に対して、ユニット整合性を持つ対角スケーリング行列を効率的に計算するための反復的アルゴリズムは何か？
• RQ5 幾何平均、$L_p$-ノルム、連続的関数 $s_{a,b}$ といった異なるサイズ関数は、スケーリングされた行列の収束性および不変性特性にどのように影響を与えるか？

主な発見

• 提案されたユニット整合性逆行列 $\mathbf{A}^{\text{-U}}$ は、入力に対角ユニット変換 $\mathbf{D}$ を適用した場合、出力に対しても対応するユニット変換が適用され、次元的整合性が保たれることを保証する。
• 反復的スケーリングアルゴリズムは高速に収束し、各行および各列の非ゼロ要素の積の大きさが1になる一意なスケーリング行列 $\mathbf{S} = \mathbf{D}_L \mathbf{A} \mathbf{D}_R$ を生成する。
• UC-SVD および UC-分解手法は対角スケーリングに対して不変であり、入力変数の単位選択にかかわらず、線形システムの解析を一貫して行える。
• 幾何平均サイズ関数 $s_{\times}(\mathbf{u})$ は、連続的かつ同次的であり、ユニット不変な正規化を提供する。これは、要素行列における補題 IV.4 および V.3 の結果と一致する。
• $\texttt{dscale}$ の実装は $O(mn)$ の計算量を達成し、実数または複素数行列に対して対角スケーリング行列を正しく計算する。対数変換と反復的精錬により、数値的安定性が保証される。
• 本手法により、ユニット整合性テストが関数的整合性のチェックとして利用可能になる。もしシステムが対角スケーリングの下でもユニット整合性を保てなければ、高確率で欠陥が検出され、チェックサム検証に類似した動作が実現される。