Suppose $A \in \mathbb{C}^{m \times n}$. I know that if $m \ge n$ and $A$ has full rank, then $A^*A$ is invertible so $A^*A x = A^*b$ has a solution. However, looking at this question, I believe that this system has a solution even if $A$ is not of full rank, and even if $m \le n$.
I believe I have a solution using the SVD, but it is rather long and I'm sure that there is a better way. My question is: can you please give me a hint on how to prove this in an easier way?
Solution using SVD
Let $A^* = U \Sigma V^*$ be an SVD of $A^*$. That is $U, V$ are unitary and $$\Sigma = \begin{bmatrix} \sigma & 0\\ 0^T & 0 \end{bmatrix}$$ where $\sigma = diag(\sigma_1, ..., \sigma_r)$ where $r$ is the rank of $A^*$.
The equation $A^*A x = A^*b$ becomes $$U \Sigma \Sigma^TV^* x= U \Sigma V^* b$$ $$\Sigma(\Sigma^TV^* x- V^*b) = 0$$ $$\Sigma(\Sigma^Ty- V^*b) = 0$$ where $y = V^*x$. Since $V$ is invertible, we are free to select any $y$ that we want. Now, $$\Sigma^Ty = [\sigma_1 y_1 \cdots \sigma_r y_r \ \ 0 \cdots 0]^T$$ Select $y_1, ..., y_r$ so that the first $r$ components of $\Sigma^Ty$ match the first $r$ components of $V^*b$. Then the first $r$ components of $\Sigma^Ty- V^*b$ are $0$, and as a result, $$\Sigma(\Sigma^Ty- V^*b) = 0.$$