In this paper, the moving least-squares differential quadrature (MLSDQ) method is employed for free vibration of thick antisymmetric laminates based on the first-order shear deformation theory. The generalized displacements of the laminates are independently approximated with the centered moving least-squares (MLS) technique within each domain of influence. The MLS nodal shape functions and their partial derivatives are computed quickly through back-substitutions after only one LU decomposition. Subsequently, the weighting coefficients in the MLSDQ discretization are determined with the nodal partial derivatives of the MLS shape functions. The MLSDQ method combines the merits of both the differential quadrature and meshless methods which can be conveniently applied to complex domains and irregular discretizations without loss of implementation efficiency and numerical accuracy. The natural frequencies of the laminates with various edge conditions, ply angles, and shapes are calculated and compared with the existing solutions to study the numerical accuracy and stability of the MLSDQ method. Effects of support size, order of completeness of basis functions, and node irregularity on the numerical accuracy are investigated in detail.