The randomized Kaczmarz (RK) method is an efficient iterative projection algorithm with low computational complexity for solving consistent linear systems. However, noise is inevitable in real-world applications, and both the coefficient matrix and the right-hand side vector may be contaminated by noise. The convergence analysis of RK-type methods for doubly noisy linear systems, where both the system matrix and the measurement vector are perturbed, remains relatively limited. In this paper, we investigate the limiting behavior of the RK algorithm for solving doubly noisy inconsistent linear systems without imposing any additional initial assumptions. Furthermore, to the best of our knowledge, this work provides the first convergence analysis of the randomized extended Kaczmarz (REK), randomized block Kaczmarz (RBK), and randomized double block Kaczmarz (RDBK) algorithms for doubly noisy linear systems. We prove that these algorithms converge to a neighborhood of the least-squares solution of the underlying noiseless system. Compared with existing theoretical estimates, the proposed bounds effectively characterize the convergence behavior of these algorithms when applied to doubly noisy linear systems. Finally, numerical experiments are conducted to validate the theoretical results.
The classical method of Kaczmarz, introduced in 1937, is a textbook iterative method for solving linear systems $A x = b.$ Despite its widespread use, particularly in the context of solving inverse problems where it is commonly included in software packages within Matlab, Python, and Julia, the precise characterization...
Stochastic gradient methods have gained increasing attention for solving large-scale inverse problems due to their computational efficiency. However, their theoretical justification in the context of ill-posed problems remains underdeveloped, particularly regarding convergence rate analysis, where existing results...
Qi-Nian Jin, Xi-Liang Lu, Lin Tian· Numerische Mathematik· 0 citations
Abstract.
We propose an adaptive randomized truncation estimator for Krylov subspace methods that optimizes the trade-off between the solution variance and the computational cost while remaining unbiased. The estimator solves a constrained optimization problem to compute the truncation probabilities on the fly with mi...
Qi Luo, Florian T. Schäfer· SIAM Journal on Scientific C...· 0 citations
A globally convergent regularized Newton method with positive definite regularization for solving nonsmooth optimization problems that replaces the identity matrix in traditional algorithms with a general positive-definite symmetric matrix to regularize the generalized Hessian.
In order to solve large sparse linear systems, we introduce a relaxation factor into the iterative formula based on the randomized Kaczmarz method. According to this idea, a parameterized randomized Kaczmarz method is constructed. Then, the convergence of this method is proved when the linear systems are consistent. Nu...
Chen-Xiao Gao, Fang Chen· Journal of Mathematical Stud...· 0 citations
We use cookies to run the site and, with your consent, for analytics and to show ads.
See our Cookie Policy.