The correspondence between the quantum spectra and the classical orbits of the two-dimensional elliptic billiards is investigated. In order to study systems without analytic solutions in the future, an expansion method for stationary states (EMSS) to get highly accurate numerical solutions is applied and tested. The closed-orbit theory was extended to the open orbit case. By comparing the Fourier-transformed spectra with the classical orbits, we find that the peak positions in the quantum spectra match with the lengths of the classical orbits with small numerical differences. This is another example showing that semiclassical method provides a bridge between quantum and classical mechanics.