A Note on Randomized Kaczmarz Algorithm for Solving Doubly-Noisy Linear Systems

08/31/2023
by   El Houcine Bergou, et al.
0

Large-scale linear systems, Ax=b, frequently arise in practice and demand effective iterative solvers. Often, these systems are noisy due to operational errors or faulty data-collection processes. In the past decade, the randomized Kaczmarz (RK) algorithm has been studied extensively as an efficient iterative solver for such systems. However, the convergence study of RK in the noisy regime is limited and considers measurement noise in the right-hand side vector, b. Unfortunately, in practice, that is not always the case; the coefficient matrix A can also be noisy. In this paper, we analyze the convergence of RK for noisy linear systems when the coefficient matrix, A, is corrupted with both additive and multiplicative noise, along with the noisy vector, b. In our analyses, the quantity R̃=Ã^†_2^2 Ã_F^2 influences the convergence of RK, where à represents a noisy version of A. We claim that our analysis is robust and realistically applicable, as we do not require information about the noiseless coefficient matrix, A, and considering different conditions on noise, we can control the convergence of RK. We substantiate our theoretical findings by performing comprehensive numerical experiments.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
08/15/2023

On Subsampled Quantile Randomized Kaczmarz

When solving noisy linear systems Ax = b + c, the theoretical and empiri...
research
10/12/2019

Structured condition number for multiple right-hand side linear systems with parameterized quasiseparable coefficient matrix

In this paper, we consider the structured perturbation analysis for mult...
research
04/08/2023

A comparison of Krylov methods for Shifted Skew-Symmetric Systems

It is well known that for general linear systems, only optimal Krylov me...
research
08/04/2021

QuantileRK: Solving Large-Scale Linear Systems with Corrupted, Noisy Data

Measurement data in linear systems arising from real-world applications ...
research
04/09/2014

Noisy Optimization: Convergence with a Fixed Number of Resamplings

It is known that evolution strategies in continuous domains might not co...
research
03/16/2022

Convergence Acceleration of Preconditioned CG Solver Based on Error Vector Sampling for a Sequence of Linear Systems

In this paper, we focus on solving a sequence of linear systems with an ...
research
04/22/2021

Operator Augmentation for General Noisy Matrix Systems

In the computational sciences, one must often estimate model parameters ...

Please sign up or login with your details

Forgot password? Click here to reset