The research generalizes the traditional asymptotic waveform evaluation method only with first order poles to the method with higher order poles case. We discuss the invertibility of the induced generalized Vandermonde matrix with the determinant formula of the matrix. Via establishing the equivalence of the multiplicity of the poles and the coefficients of the introduced matrix, the generalized asymptotic waveform evaluation algorithm is built. Then the real system realization of the lower order poles-residues form is asserted. At last two experiments are provided to check the validity of our algorithm.