We utilize multi-epoch MUSE spectroscopy to study binaries in the core of NGC 3201. Our sample consists of 3553 stars with 54883 spectra in total comprising 3200 main-sequence stars up to 4 magnitudes below the turn-off. Each star in our sample has between 3 and 63 (with a median of 14) reliable radial velocity (RV) measurements within five years of observations. We introduce a statistical method to determine the probability of a star showing RV variations based on the whole inhomogeneous RV sample. Using HST photometry and an advanced dynamical MOCCA simulation of this specific GC we overcome observational biases that previous spectroscopic studies had to deal with. This allows us to infer a binary frequency in the MUSE FoV and enables us to deduce the underlying true binary frequency of (6.75+/-0.72)% in NGC 3201. The comparison of the MUSE observations with the MOCCA simulation suggests a significant fraction of primordial binaries. We can also confirm a radial increase of the binary fraction towards the GC centre due to mass segregation. We discovered that in our sample at least (57.5+/-7.9)% of blue straggler stars (BSS) are in a binary system. For the first time in a study of GCs, we were able to fit Keplerian orbits to a significant sample of 95 binaries. We present the binary system properties of eleven BSS and show evidence that two BSS formation scenarios, the mass transfer in binary (or triple) star systems and the coalescence due to binary-binary interactions, are present in our data. We also describe the binary and spectroscopic properties of four sub-subgiant (or red straggler) stars. Furthermore, we discovered two new black hole (BH) candidates with minimum masses (Msini) of (7.68+/-0.50)M_{sun}_, (4.4+/-2.8)M_{sun}_, and refine the minimum mass estimate on the already published BH to (4.53+/-0.21)M_{sun}_. These BHs are consistent with an extensive BH subsystem hosted by NGC 3201.