[論文レビュー] Fast acoustic source imaging using multi-frequency sparse data

Ala Alzaalig, Guanghui Hu+2 more|arXiv (Cornell University)|Dec 7, 2017

Ultrasonics and Acoustic Wave Propagation参考文献 1被引用数 10

ひとこと要約

本稿では、多周波数のスパースな遠方場または近方場データを用いた音源像の高速で反復処理を要しないサンプリング手法を提案する。行列-ベクトル乗算を用いて指標関数を計算することで、限られたデータからでも高い精度で音源のサポートを効率的に再構成でき、周波数を複数用いることで1つの観測方向からの観測でも音源の位置とサイズを回復可能であることを証明している。

ABSTRACT

We consider the acoustic source imaging problems using multiple frequency data. Using the data of one observation direction/point, we prove that some information (size and location) of the source support can be recovered. A non-iterative method is then proposed to image the source for the Helmholtz equation using multiple frequency far field data of one or several observation directions. The method is simple to implement and extremely fast since it only computes an indicator function on the interested domain using only matrix vector multiplications. Numerical examples are presented to validate the effectiveness of the method.

研究の動機と目的

• スパースで広帯域の多周波数測定値から音源サポートを再構成する課題に対処すること。
• 一様周波数データでは拡張された源を一意に特定できないという制限を克服すること。
• 直接問題の反復的解法を避ける計算効率の高い手法を開発すること。
• 多周波数データを用いて1つまたは数個の観測方向からのみでも、源サポート回復の理論的・数値的妥当性を確立すること。
• 空間的分布や周波数依存分布を示す源に対して、本手法のロバスト性を示すこと。

提案手法

• 多周波数の遠方場または近方場データを用いた、修正された因子分解法に基づく反復処理を要しないサンプリング手法を提案する。
• 遠方場パターンと観測方向を含む行列-ベクトル乗算により、関心領域上での指標関数を構築する。
• スパースな測定でも分解能と安定性を向上させるために、周波数帯域 $k \in (0, k_{\text{max}})$ の多周波数データを用いる。
• ヘルムホルツ方程式の漸近的挙動を活用し、指標関数を用いて高い値を示す領域を特定することで、源サポートを検出する。
• 本手法は、源と測定可能なデータを結ぶ遠方場パターン $u^\infty(\theta, k) = \int_D e^{-ik\theta \cdot y} F(y,k) dy$ に依存している。
• 離散的な $k$ と $\theta$ のサンプリングを用いた数値的実装により、前方問題の反復的解法を一切回避し、高速な計算を実現している。

実験結果

リサーチクエスチョン

• RQ1 反復処理を要しない手法を用いて、多周波数のスパースな遠方場または近方場データから音源のサポートを再構成できるか？
• RQ2 多周波数データが利用可能な場合、1つの観測方向からの情報が源サポートの特定にどの程度寄与できるか？
• RQ3 空間的に変化する分布や周波数依存の分布を示す源に対して、本手法はどのように動作するか？
• RQ4 源のサイズと波長の影響は、特に拡張されたまたは大きな物体の再構成精度にどのような影響を及えるか？
• RQ5 測定セットが正のレブゴフ測度を有するか、スパースである場合でも、本手法は源の位置とサイズを回復できるか？

主な発見

• 本手法は、1つの観測方向と多周波数データのみを用いても、音源のサポートを成功裏に再構成でき、指標関数上に明確に源の位置とサイズが特定できる。
• 2つの成分または複雑な形状の源に対しては、複数の観測方向からのストリップの重ね合わせが、正確に源サポートを囲む。
• 源関数 $F$ が波数 $k$ に依存する場合（例：$F_1(x,y;k) = k^2(x^2 - y^2 + 5)$ または $F_2(x,y;k) = e^{ik(x\cos 3\pi/2 + y\sin 3\pi/2)}(x^2 - y^2 + 5)$）であっても、再構成は正確に保たれる。
• 最短波長 ($\lambda_{\text{min}} \approx 0.628$) と同程度またはそれ以上のサイズの大きな物体は、複数の観測方向（$M=20$）を用いることで、形状・サイズ・位置を正確に再構成できる。
• 短い波長は再構成品質を向上させ、三角形や細いスラブ支持部の微細構造的詳細を正確に回復可能にする。
• 本手法は計算的に効率的であり、行列-ベクトル演算のみに依存しており、前方問題の高コストな反復的解法を回避するため、リアルタイムまたは大規模応用に適している。