<p>Principal Component Analysis in statistics, which essentially provides a way to remove the least relevant factors to a set of data, is derived using singular value decomposition, which makes use of eigenvectors. Once you’ve derived it, though, you can get away without singular value decomposition - you just need to find eigenvalues and eigenvectors, which represent “components” of the data.</p>
<p>Linear’s also useful for solving systems of linear differential equations, since the technique for solving them involves exponentiating a matrix (essentially by use of the taylor series for e^x, but x is a matrix). In these cases, we want orthogonal eigenvectors, since this lets us decompose a matrix into A = SDS^-1, where S is the eigenvectors as columns and D is the eigenvalues along the diagonals. This has nice properties like A^n = S(D^n)S^-1, making multiplication (and inversion) easy. When you get this form, the Taylor representation of e^A simplifies nicely, and that’s where you get solutions involving e to the power of eigenvalues of a matrix. In the case of repeated eigenvalues, another decomposition trick is needed, though.</p>
<p>Approximating solutions can be done with least squares, since even if Ax=b has no solutions, A<em>Ax = A</em>b, where A* both takes the transpose of A and takes the complex conjugate of all entries, does necessarily have a solution (though not necessarily a unique one, since A*A can be non-invertible). Making use of projections onto solution spaces, it can be shown that this provides the best solution to Ax=b, even if no actual solutions exist.</p>
<p>Not sure what you mean by determining if solutions exist, beyond the standard Ax=b has solutions iff A is invertible (and has non-zero determinant, and Ax=0 has infinitely many solutions, etc.). A necessarily has a number of (not-necessarily distinct) eigenvalues equal to the rank of A, but finding them becomes tough as A gets big, since there’s no formula to solve 5th degree polynomials exactly. I’d imagine one would have to estimate roots and use the A<em>Ax=A</em>b trick to estimate eigenvectors, then orthogonalize them to make use of them, but I don’t know when those would be useful.</p>
<p>However, I must argue that it’s not necessarily true that eigenvectors are orthogonal. In the case where all eigenvalues have multiplicity 1, that’s true. However, if they don’t, multiple eigenvectors from the same eigenvalue need not be orthogonal. That’s when you use Gram-Schmidt (or other techniques) to orthogonalize them.</p>
<p>(I like linear algebra and can deal with abstract algebra…but I don’t claim to understand or know most of the applications).</p>