We can, of course, always write a system of equations as a matrix problem of the form Ax= b. In particular, if A is a square matrix, mapping, say, Rn to Rn, there are three possibilities for the equation Ax= b.
1) the equation has a unique solution.
2) the equation has NO solution.
3) the equation has an infinite number of solutions.
If A is invertible, we can multiply both sides of Ax= b by its inverse to get A^{-1}Ax= x= A^{-1}b. That clearly is a solution to the equation and, because x must be equal to that, the solution is unique.
If A is not invertible, then it has rank, k, less than n. That means that A maps all of Rn into a subspace of Rn of dimension k< n. It also means that A has null space (space of vectors, v, such that Av= 0) of dimension n- k.
In that case, if b does not happen to lie in that subspace, the equation, Ax= b, has no solution. For any x, Ax is in that subspace and b is not. If b does happen to lie in that subspace, there is a solution, a vector x such that Ax= b, but it is also true that for any v in the null space of A, A(x+ v)= Ax+ Av= b+ 0= b so there are an infinite number of solutions.
Since that is true for square matrices, what does it have to do with this problem? Well, if we have m equations, in n unknowns, with m< n, we can represent it has the matrix equation Ax= b where A is Not square. It has m rows and n columns. But we can make it square by adding n- m rows of all zeros and adding n- m 0s to the vector b. Since that matrix obviously has determinant 0, it is not invertible, which means it has either no solution or an infinite number of solutions. Because the system is homogeneous, Ax= 0, it has the obvious solution, x= 0, which means "no solution" cannot apply. We are left with an infinite number of solutions.