A Novel Algorithm to Solve for an Underwater Line Source Sound Field Based on the Theory of Global Matrix Coupled Modes and the Chebyshev–Tau Spectral Method