Recurrence Relation for Iterative Solution of Ax=b

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
1 reply · 2K views
Kidnapster
Messages
2
Reaction score
0
[itex]_nx[/itex] represents the nth iteration.

Homework Statement



Show that [itex]{_{n+2}}x = _nx + M^{-1}(b-A{_nx})[/itex]

Homework Equations



[itex]{_{n+1}}x = _nx + M_1^{-1}(b-A{_nx})[/itex]
[itex]{_{n+2}}x = {_{n+1}}x + M_2^{-1}(b-A{_{n+1}x})[/itex]
[itex]M = M_1(M_1+M_2-A)^{-1}M_2[/itex]

The Attempt at a Solution



I tried was adding the two equations, which gives [itex]_{n+2}x = M_2^{-1}(b-A{_{n+1}x}) + _nx + M_1^{-1}(b-A{_nx})[/itex]. But [itex]M_2^{-1}(b-A{_{n+1}}x)[/itex] devolves into a huge mess once I expand it. Subtraction is less helpful, as you are left without any [itex]_nx[/itex]. Directly substituting for [itex]_{n+1}x[/itex] in the second equation is even less helpful, as you are left with a huge mess of terms that can't really be factored. What am I missing?
 
Physics news on Phys.org
Never mind, figured it out. Here's my solution:

[tex] \begin{align*}<br /> {_{n+1}x} &= {_nx} + M_1^{-1}(b-A_nx)\\<br /> {_{n+2}x} &= {_{n+1}x} + M_1^{-1}(b-A_{n+1}x) \\<br /> {_{n+2}x} &= {_nx} + M_1^{-1}(b-A_nx) + M_2^{-1}(b-A[{_nx} + M_1^{-1}(b-A_nx)])\\<br /> {_{n+2}x} &= {_nx} + M_1^{-1}b - M_1^{-1}A_nx + M_2^{-1}b-M_2^{-1}A_nx - M_2^{-1}AM_1^{-1}b + M_2^{-1}AM_1^{-1}A_nx \\<br /> {_{n+2}x} &= {_nx} + M_1^{-1}(b-A_nx) + M_2^{-1}(b-A_nx) - M_2^{-1}AM_1^{-1}(b-A_nx) \\<br /> {_{n+2}x} &= {_nx} + M_1^{-1}(M_1+M_2-A)M_2^{-1}(b-A_nx)<br /> \end{align*}<br /> [/tex]