Gaussian beam decomposition of high frequency wave fields using expectation-maximization

Gil Ariel, Björn Engquist, Nicolay M. Tanushev, Richard Tsai

Research output: Contribution to journalArticlepeer-review

18 Scopus citations

Abstract

A new numerical method for approximating highly oscillatory wave fields as a superposition of Gaussian beams is presented. The method estimates the number of beams and their parameters automatically. This is achieved by an expectation-maximization algorithm that fits real, positive Gaussians to the energy of the highly oscillatory wave fields and its Fourier transform. Beam parameters are further refined by an optimization procedure that minimizes the difference between the Gaussian beam superposition and the highly oscillatory wave field in the energy norm.

Original languageEnglish
Pages (from-to)2303-2321
Number of pages19
JournalJournal of Computational Physics
Volume230
Issue number6
DOIs
StatePublished - 20 Mar 2011

Bibliographical note

Funding Information:
The authors were partially supported by the NSF under Grant No. DMS-0714612 . N.T. was also partially supported by the NSF under Grant No. DMS-0636586 (UT Austin RTG).

Funding

The authors were partially supported by the NSF under Grant No. DMS-0714612 . N.T. was also partially supported by the NSF under Grant No. DMS-0636586 (UT Austin RTG).

FundersFunder number
National Science FoundationDMS-0714612, DMS-0636586

    Keywords

    • Expectation-maximization
    • Gaussian beams
    • High frequency waves
    • Wave equation

    Fingerprint

    Dive into the research topics of 'Gaussian beam decomposition of high frequency wave fields using expectation-maximization'. Together they form a unique fingerprint.

    Cite this