What would you like to understand?

The solution to the normal equations is considered the least-squares solution because it satisfies the necessary and sufficient condition for minimizing the residual norm $\|\vec{b} - A\vec{x}\|$. The derivation shows that minimizing this norm is equivalent to requiring the residual $A\vec{x}^* - \vec{b}$ to be orthogonal to the column space $C(A)$.

Conditions: The original system $A\vec{x}=\vec{b}$ may be inconsistent; $A^T A \vec{x}^* = A^T \vec{b}$ has a solution

The solution to the normal equations is considered the least-squares solution because it satisfies the necessary and sufficient condition for minimizing the residual norm $\|\vec{b} - A\vec{x}\|$. The derivation shows that minimizing this norm is equivalent to requiring the residual $A\vec{x}^* - \vec{b}$ to be orthogonal to the column space $C(A)$.

Conditions: The original system $A\vec{x}=\vec{b}$ may be inconsistent; $A^T A \vec{x}^* = A^T \vec{b}$ has a solution