Tengo un método iterativo \begin{eqnarray} X_{k+1}=(1+\beta)X_k-\beta X_k A X_k~~~~~ k = 0,1,\ldots \end{eqnarray} con aproximación inicial X0=βA∗ (β es un escalar que va de 0 a 1) que converge a la pseudo-inversa de Moore-Penrose A+. {Xk} son una secuencia de aproximaciones. Debo demostrar que {Xk} definido por las ecuaciones anteriores satisface R(Xk)=R(A∗) y N(Xk)=N(A∗) para k≥0. Donde R(Xk), N(Xk) denotan el rango y el espacio nulo de {Xk}. Estoy usando inducción matemática para probar el resultado anterior. Esto es lo que hice.
Trivialmente se cumple para k=0 (ya que para k=0,X0=βA∗). Supongamos que el resultado es cierto para algún k. Sea y∈N(Xk) un vector arbitrario. De la ecuación anterior tenemos Xk+1y=(1+β)Xky−βXkAXky=0 (ya que Xky=0). Esto nos da y∈ N(Xk+1) por lo tanto N(Xk)⊆N(Xk+1) es verdadero. Ahora mi pregunta es cómo se puede demostrar de manera análoga la afirmación R(Xk)⊇R(Xk+1). He intentado mucho pero no logro descubrir cómo probar R(Xk)⊇R(Xk+1). Un pequeño indicio funcionará para mí. Estaré muy agradecido.