Analytical solutions founded on three-dimensional theories play a crucial role in evaluating the credibility and precision of different plate theories and numerical methodologies. While Levy-type analytical solutions are widely recognized, they have been primarily confined to purely elastic plates. This study introduces a Levy-type analytical series solution for three-dimensional vibrations in a sandwich rectangular plate featuring a functionally graded material (FGM) core, along with piezoelectric material (PM) layers on the top and bottom surfaces. The behaviors of the FGM and PM layers were described using three-dimensional elasticity and piezoelasticity theories, respectively. In this study, the displacement functions and electric potential of each layer were expanded by Fourier series and polynomial auxiliary functions. An analytical series solution was then established by satisfying the governing equations of each layer, the mechanical and electric boundary conditions on the six faces of the plate, and the continuity conditions on the interfaces between the PM and FGM layers. To validate the proposed solutions, in-depth convergence studies were conducted for the vibration frequencies of the first six modes of sandwich square plates with various boundary conditions on the other pair of side faces. The well-converged results were then compared with published data based on various plate theories to verify the accuracy of these published data. Finally, accurate nondimensional frequencies were tabulated for the first six modes of sandwich rectangular plates with various aspect ratios, thickness-to-width ratios, PM-to-FGM layer thickness ratios, power law indices for the FGM layer, and six combinations of boundary conditions. These new numerical results when piezoelectric coupling is considered should be very useful to future analytical and numerical studies.