An analytical approach for the calculation of the dyadic Green's functions inside the rectangular cavity over a broad range of frequency is presented. Both vector potential and electric field dyadic Green's functions are considered. The method is based on the extraction of the Green's function at an imaginary wavenumber from itself to obtain a rapidly convergent eigenfunction expansion of the dyadic Green's function. The extracted term encompasses the singularity of the Green's function and is computed using spatial expansions. Results are illustrated for a rectangular cavity up to 5 wavelengths in size with thousands of cavity modes obtained by the 6th order convergent expansion. It is shown that for an accurate and broadband simulation, the proposed method is many times faster than the Ewald method.