Die vorliegende Arbeit untersucht die longitudinale Dynamik von Halbleiterlasern anhand eines Modells, in dem ein lineares hyperbolisches System partieller Differentialgleichungen mit gewöhnlichen Differentialgleichungen gekoppelt ist. Zunächst wird mit Hilfe der Theorie stark stetiger Halbgruppen die globale Existenz und Eindeutigkeit von Lösungen für das konkrete System gezeigt. Die anschließende Untersuchung des Langzeitverhaltens der Lösungen erfolgt in zwei Schritten. Zuerst wird ausgenutzt, dass Ladungsträger und optisches Feld sich auf unterschiedlichen Zeitskalen bewegen, um mit singulärer Störungstheorie invariante attrahierende Mannigfaltigkeiten niedriger Dimension zu finden. Der Fluss auf diesen Mannigfaltigkeiten kann näherungsweise durch Moden-Approximationen beschrieben werden. Deren Dimension und konkrete Gestalt ist von der Lage des Spektrums des linearen hyperbolischen Operators abhängig. Die zwei häufigsten Situationen werden dann einer ausführlichen numerischen und analytischen Bifurkationsanalyse unterzogen. Ausgehend von bekannten Resultaten für die Ein-Moden-Approximation, wird die Zwei-Moden-Approximation in dem speziellen Fall untersucht, dass die Phasendifferenz zwischen den beiden optischen Komponenten sehr schnell rotiert, so dass sie sich in erster Ordnung herausmittelt. Mit dem vereinfachten Modell können die Mechanismen verschiedener Phänomene, die bei der numerischen Simulation des kompletten Modells beobachtet wurden, erklärt werden. Darüber hinaus lässt sich die Existenz eines anderen stabilen Regimes voraussagen, das sich im gemittelten Modell als "bursting" darstellt. / We investigate the longitudinal dynamics of semiconductor lasers using a model which couples a linear hyperbolic system of partial differential equations with ordinary differential equations. We prove the global existence and uniqueness of solutions using the theory of strongly continuous semigroups. Subsequently, we analyse the long-time behavior of the solutions in two steps. First, we find attracting invariant manifolds of low dimension benefitting from the fact that the system is singularly perturbed, i. e., the optical and the electronic variables operate on different time-scales. The flow on these manifolds can be approximated by the so-called mode approximations. The dimension of these mode approximations depends on the number of critical eigenvalues of the linear hyperbolic operator. Next, we perform a detailed numerical and analytic bifurcation analysis for the two most common constellations. Starting from known results for the single-mode approximation, we investigate the two-mode approximation in the special case of a rapidly rotating phase difference between the two optical components. In this case, the first-order averaged model unveils the mechanisms for various phenomena observed in simulations of the complete system. Moreover, it predicts the existence of a more complex spatio-temporal behavior. In the scope of the averaged model, this is a bursting regime.