On convergence rates of adaptive ensemble Kalman inversion for linear ill-posed problems
Fabian Parzer1, Otmar Scherzer1,2,3
1Faculty of Mathematics, University of Vienna, Oskar-Morgenstern-Platz 1, 1090 Vienna, Austria.
Ensemble Kalman inversion offers a novel regularization method for linear inverse problems. A new sampling scheme improves performance and creates an order-optimal method for accurate solutions.
Area of Science:
- Applied Mathematics
- Numerical Analysis
- Computational Science
Background:
- Linear inverse problems are fundamental in various scientific fields.
- Traditional regularization methods like Tikhonov regularization have limitations.
- Ensemble Kalman inversion (EKI) is an emerging technique for inverse problems.
Purpose of the Study:
- To present a deterministic Ensemble Kalman Inversion (EKI) as a regularization technique for linear inverse problems.
- To introduce a novel sampling scheme for EKI based on the Nyström method.
- To develop an adaptive EKI that couples sample size with the regularization parameter.
Main Methods:
- Interpreting EKI as a low-rank approximation of Tikhonov regularization.
- Implementing a Nyström-based sampling scheme for enhanced practical performance.
- Formulating an adaptive EKI with a coupled sample size and regularization parameter.
- Utilizing the discrepancy principle as a stopping criterion.
Main Results:
- The proposed Nyström-based sampling scheme improves the practical performance of EKI.
- The adaptive EKI formulation, coupled with the discrepancy principle, achieves order-optimal regularization.
- Numerical comparisons demonstrate the effectiveness of the proposed methods for Radon transform inverse problems.
Conclusions:
- Deterministic EKI provides an effective regularization strategy for linear inverse problems.
- The Nyström method enhances EKI's practical applicability.
- The adaptive EKI scheme offers theoretical guarantees for order-optimal solutions.
More Related Videos
13:07Optical Coherence Tomography Based Biomechanical Fluid-Structure Interaction Analysis of Coronary Atherosclerosis Progression
Published on: January 15, 2022
09:01Gain-compensation Methodology for a Sinusoidal Scan of a Galvanometer Mirror in Proportional-Integral-Differential Control Using Pre-emphasis Techniques
Published on: April 4, 2017
Related Concept Videos
Linear time-invariant Systems
The input-output behavior of an LTI system can be fully defined by its response to an impulsive excitation at its input. Once this impulse response is known, the system's reaction to any other input can be...
Linear Approximation in Time Domain
For a simple pendulum with a mass evenly distributed along its length and the center of mass located at half the pendulum's length,...
Mechanistic Models: Compartment Models in Algorithms for Numerical Problem Solving
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
One-Compartment Open Model: Wagner-Nelson and Loo Riegelman Method for ka Estimation
On...
Linear Approximation in Frequency Domain
In contrast, nonlinear systems do not inherently possess these properties. However, for small deviations around an operating point, a nonlinear system can often be approximated as linear....
Parameters Affecting Nonlinear Elimination: Zero-Order Input, First-Order Absorption and Two-Compartment Model
When a drug is administered through a constant intravenous infusion and eliminated via nonlinear pharmacokinetics, it follows zero-order input. For example, oral drugs undergo first-order absorption upon administration and are eliminated through nonlinear pharmacokinetics.
In the case of subcutaneously administered drugs,...
