Simulating wave propagation is one of the fundamental problems in scientific computing. In this paper, we consider one-dimensional two-way wave equations, and investigate a family of L stable high order discontinuous Galerkin methods defined through a general form of numerical fluxes. For these L stable methods, we systematically establish stability (hence energy conservation), error estimates (in both L and negative-order norms), and dispersion analysis. One novelty of this work is to identify a sub-family of the numerical fluxes, termed αβ-fluxes. Discontinuous Galerkin methods with αβ-fluxes are proven to have optimal L error estimates and superconvergence properties. Moreover, both the upwind and alternating fluxes belong to this sub-family. Dispersion analysis, which examines both the physical and spurious modes, provides insights into the sub-optimal accuracy of the methods using the central flux and the odd degree polynomials, and demonstrates the importance of numerical initialization for the proposed non-dissipative schemes. Numerical examples are presented to illustrate the accuracy and the long-term behavior of the methods under consideration.