In 21 cm experimental cosmology, accurate characterization of a radio telescope’s antenna beam response is essential to measure the 21 cm signal. Computational electromagnetic (CEM) simulations estimate the antenna beam pattern and frequency response by subjecting the EM model to different dependencies, or beam hyperparameters, such as soil dielectric constant or orientation with the environment. However, it is computationally expensive to search all possible parameter spaces to optimize the antenna design or accurately represent the beam to the level required for use as a systematic model in 21 cm cosmology. We therefore present the Model for Emulating Directivities and Electric fields of Antennas (MEDEA), an emulator that rapidly and accurately generates far-field radiation patterns over a large hyperparameter space. MEDEA takes a subset of beams simulated by CEM software, spatially decomposes them into coefficients on a complete, linear basis, and then interpolates them to form new beams at arbitrary hyperparameters. We test MEDEA on an analytical dipole and two numerical beams motivated by upcoming lunar lander missions, and then employ MEDEA as a model to fit mock radio spectrometer data to extract covariances on the input beam hyperparameters. We find that the interpolated beams have rms relative errors of at most 10−2 using 20 input beams or less, and that fits to mock data are able to recover the input beam hyperparameters when the model and mock are derived from the same set of beams. When a systematic bias is introduced into the mock data, extracted beam hyperparameters exhibit bias, as expected. We propose several extensions to MEDEA to potentially account for such bias.