In this study, we obtained results for the computation of eigen-pairs, singular value decomposition, pseudoinverse, and the least squares problem for elliptic quaternion matrices. Moreover, we established algorithms based on these results and provided illustrative numerical experiments to substantiate the accuracy of our conclusions. In the experiments, it was observed that the p-value in the algebra of elliptic quaternions directly affects the performance of the problem under consideration. Selecting the optimal p-value for problem-solving and the elliptic behavior of many physical systems make this number system advantageous in applied sciences.