The free vibration of a circular plate with multiple perforations is analyzed by using the Rayleigh-Ritz method. Admissible functions are assumed to be separable functions of radial and tangential coordinates. Trigonometric functions are assumed in the circumferential direction. The radial shape functions are the boundary characteristic orthogonal polynomials generated following the Gram-Schmidt recurrence scheme. The assumed functions are used to estimate the kinetic and the potential energies of the plate depending on the number and the position of the perforations. The eigenvalues, representing the dimensionless natural frequencies, are compared with the results obtained using Bessel functions, where the exact solution is available. Moreover, the eigenvectors, which are the unknown coefficients of the Rayleigh-Ritz method, are used to present the mode shapes of the plate. To validate the analytical results of the plates with multiple perforations, experimental investigations are also performed. Two unique case studies that are not addressed in the existing literature are considered. The results of the Rayleigh-Ritz method are found to be in good agreement with those from the experiments. Although the method presented can be employed in the vibration analysis of plates with different boundary conditions and shapes of the perforations, circular perforations that are free on the edges are studied in this paper. The results are presented in terms of dimensionless frequencies and mode shapes.