Abstract:A standard way to solve linear algebraic systems $Au=f,\,\,(*)$ with ill-conditioned matrices $A$ is to use variational regularization. This leads to solving the equation $(A^*A+aI)u=A^*f_\d$, where $a$ is a regularization parameter, and $f_\d$ are noisy data, $||f-f_\d||\leq \d$. Numerically it requires to calculate products of matrices $A^*A$ and inversion of the matrix $A^*A+aI$ which is also ill-conditioned if $a>0$ is small. We propose a new method for solving (*) stably, given noisy data $f_\d$. This method, the DSM (Dynamical Systems Method) is developed in this paper for selfadjoint $A$. It consists in solving a Cauchy problem for systems of ordinary differential equations.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to solve the ill - conditioned system of linear algebraic equations \(Au = f\), where the matrix \(A\) is ill - conditioned, that is, the condition number \(k(A)\gg1\). The characteristic of this kind of problem is that a small input perturbation may lead to a huge change in the solution, making the problem ill - posed. The traditional solution method is to use variational regularization, but this involves calculating the product of the matrix \(A^*A\) and its inverse matrix \((A^*A + aI)^{-1}\), which is very difficult and time - consuming numerically, especially when \(a\) is very small.
To solve these problems, the author proposes a new method - the Dynamical Systems Method (DSM). DSM does not need to directly calculate the product and inverse of matrices, but indirectly solves the original problem by solving the Cauchy problem of ordinary differential equations (ODEs). Specifically, for a self - adjoint operator \(A\), the core idea of the DSM method is based on the following formula:
\[
B^{-1}(e^{Bt}-I)=\int_{0}^{t}e^{Bs}ds
\]
If
\[
\lim_{t\rightarrow\infty}\|e^{Bt}\| = 0
\]
then
\[
B^{-1}=-\lim_{t\rightarrow\infty}\int_{0}^{t}e^{Bs}ds
\]
By this method, solving \(A^{-1}f\) can be transformed into solving an initial - value problem of an ordinary differential equation:
\[
\dot{u}_{a}=i(A + ia)u_{a}+f,\quad u_{a}(0) = 0
\]
The main results of the paper include two theorems:
1. **Theorem 1**: When the noisy data \(f_{\delta}\) is close to the real data \(f\), the solution \(u_{a}(t)\) obtained by the DSM method converges to the minimum - norm solution \(y\) in the limits of \(a\rightarrow0\) and \(t\rightarrow\infty\).
2. **Theorem 2**: When appropriate parameters \(t_{\delta}\) and \(a(\delta)\) are selected, even in the presence of noisy data \(f_{\delta}\), the DSM method can still obtain a solution \(u_{\delta}\) close to the minimum - norm solution \(y\).
In conclusion, this paper proposes a new numerical method DSM to stably solve ill - conditioned systems of linear algebraic equations, avoiding the numerical difficulties in traditional methods, and has potential advantages in both theoretical and practical applications.