[論文レビュー] Many physical laws are ridge functions

Paul G. Constantine, Zachary del Rosario+1 more|arXiv (Cornell University)|May 25, 2016

Sparse and Compressive Sensing Techniques参考文献 11被引用数 10

ひとこと要約

この論文は、古典的次元解析（Buckingham Pi定理を介して）を低次元構造と結びつけることで、多くの物理法則が特定の方向に沿って定数である「リッジ関数」であることを示している。次元解析で特定された部分空間が、データ駆動型次元削減から得られるアクティブ部分空間と正確に一致することを示しており、物理法則が本質的にリッジ構造を示していることを検証している。

ABSTRACT

A ridge function is a function of several variables that is constant along certain directions in its domain. Using classical dimensional analysis, we show that many physical laws are ridge functions; this fact yields insight into the structure of physical laws and motivates further study into ridge functions and their properties. We also connect dimensional analysis to modern subspace-based techniques for dimension reduction, including active subspaces in deterministic approximation and sufficient dimension reduction in statistical regression.

研究の動機と目的

• 古典的次元解析と現代の次元削減技術との間の理論的関係を確立すること。
• 物理法則が本質的に低次元構造、特にリッジ関数形式を示しているかどうかを調査すること。
• 次元解析から得られる部分空間が、数値積分および勾配ベースの手法から導かれるアクティブ部分空間を含んでいることを検証すること。
• 次元不変性を通じて、リッジ関数が物理系において自然に出現することを示すこと。
• 古典的物理学の原則と、次元削減における現代の計算手法を統合すること。

提案手法

• m個の物理変数に対して次元ベクトルを用いて線形代数的に次元量を表現し、次元ベクトル行列Dを構築する。
• Buckingham Pi定理を適用して無次元群を導出し、独立な次元のk次元部分空間（k ≤ m）を形成する。
• f(x) = g(Aᵀx) としてリッジ関数を定義する。ここでAは次元行列Dᵀの核空間を張る。この定義により、fはAに直交する方向に沿って定数となる。
• 数値積分と有限差分を用いて、C = ∫∇f(x)∇f(x)ᵀp(x)dx の共分散行列Cを推定し、アクティブ部分空間を計算する。
• 各アクティブ部分空間の基底ベクトルが、次元解析の部分空間の基底によって線形的に表現可能かどうかをチェックすることで、部分空間の包含関係をテストする。
• ガウス求積法と有限差分スキームを用い、残差ノルムの収束を誤差解析とともに検証する。

実験結果

リサーチクエスチョン

• RQ1 次元不変性によって、物理法則が本質的にリッジ関数構造を示すのか？
• RQ2 次元解析で特定された部分空間が、勾配ベースの次元削減から得られるアクティブ部分空間を含むのか？
• RQ3 物理モデルの数値積分によるアクティブ部分空間が、次元解析の部分空間とどの程度一致するのか？
• RQ4 微分および積分における数値誤差が、部分空間の包含関係の検証にどのように影響するのか？
• RQ5 古典的次元解析と、アクティブ部分空間や十分次元削減などの現代的技法との関係は何か？

主な発見

• 次元解析から導かれる物理法則は、数学的にリッジ関数と同等であり、関数は入力変数の低次元射影にのみ依存する。
• 次元行列Dᵀの核空間から得られる次元解析の部分空間は、物理モデルの勾配の数値積分により特定されたアクティブ部分空間を含んでいる。
• 数値実験により、部分空間包含誤差を測る残差ノルムr²が、有限差分ステップサイズhを小さくするにつれてO(h)でゼロに収束することが示され、理論的部分空間への収束が確認された。
• 層流および乱流の両ケースにおいて、アクティブ部分空間は次元解析の部分空間の真の部分空間であり、広範な入力パrameter値の範囲で検証された。
• 結果から、リッジ構造がモデリングの副産物ではなく、次元不変性に起因する物理法則の根本的性質であることが確認された。
• 本研究は、古典的物理学（次元解析）と現代の計算手法（アクティブ部分空間）の間の厳密な接続を確立し、物理法則が本質的に低次元構造を持つことを示した。