On the acceleration of Kaczmarz's method for inconsistent linear systems (Q920573): Difference between revisions

Es sei \(A\) eine rechteckige Matrix. Durch \(\| y\|_ w=(\sum^{m}_{i=1}w^ 2_ i| y_ i|^ 2)^{1/2}\) wird eine gewichtete Euklidische Norm definiert. Gesucht ist nun die gewichtete kleinste-Quadrate-Lösung \(x^+\) von \(Ax=b\) mit minimaler Euklidischer Norm. Dazu wird zunächst eine Verallgemeinerung des Iterationsverfahrens von Kaczmarz mit Iterationsmatrix \(T_{\omega}\) vorgeschlagen, das von einem Parameter \(\omega\) abhängt. Für alle \(\omega\) mit \(0<\omega <{\bar \omega}\) und z.B. Startvektor \(x_0=0\) konvergiert das Verfahren gegen einen Vektor \(x(\omega)=A^-b\), wobei \(A^-\) eine gewisse verallgemeinerte Inverse von \(A\) ist. Für \(\omega\to 0\) wird \(x(\omega)\to x^+\) gezeigt. Auf Grund einer theoretischen Aussage und numerischer Rechnungen kann man davon ausgehen, daß das Spektrum von \(T_{\omega}\) für kleine \(\omega >0\) in der komplexen Ebene sehr dicht am Intervall \([0,1]\) liegt. Dies legt es nahe, obiges Verfahren mit dem Tschebyscheffschen semiiterativen Verfahren (für nichtsymmetrische Iterationsmatrizen) zu beschleunigen. Eine Aussage über den Konvergenzfaktor des neuen Verfahrens sowie ein Vergleich mit einer variierten Beschleunigungsmethode wird gemacht.