The x-ray spectra of copper $K\ensuremath{\alpha}$ and $K\ensuremath{\beta}$ remain a topic of great interest due to the complex open-shell processes involved, with many discrepancies among theories and experiments. We present high-accuracy theory of copper $K\ensuremath{\alpha}$ and $K\ensuremath{\beta}$ diagram and satellite spectra in the multiconfiguration Dirac-Hartree-Fock method. Diagram spectra are expanded to $5s$ with simultaneous convergence of 28 000 configuration state functions (CSFs) $(K\ensuremath{\alpha})$ and to $6g$ with simultaneous convergence of 91 000 CSFs $(K\ensuremath{\beta})$, achieving eigenvalue convergence to $\ifmmode\pm\else\textpm\fi{}0.03$ eV or 0.000 25%, approximately a factor of 10 improvement over past work. It is necessary to invoke biorthogonalization, developments of the active space approach, analysis of markers for theoretical convergence of eigenvalues, and the question of self-consistency for both $K\ensuremath{\alpha}$ and $K\ensuremath{\beta}$ spectra. We make use of gauge convergence, eigenvalue convergence, and $A$-coefficient convergence. The satellite spectra are a major outcome---without these, it is not possible to make use of the increased accuracy of the diagram computations. The Cu $K\ensuremath{\alpha}\phantom{\rule{4pt}{0ex}}3{d}^{8}$ double shake satellite spectrum alone contains 1506 unique eigenvalues (transitions) and required simultaneous convergence of 593 000 CSFs. Ab initio shakeoff probabilities for $1s,\phantom{\rule{0.16em}{0ex}}2s,\phantom{\rule{0.16em}{0ex}}2p,\phantom{\rule{0.16em}{0ex}}3s,\phantom{\rule{0.16em}{0ex}}3p,\phantom{\rule{0.16em}{0ex}}3d$, and $4s$ subshells as a result of the $K$-shell photoionization process are presented. Portable spectral representations are provided in the supplemental material, discussed in the text. We investigate the meaning behind shake, shakeoff, and shakeup processes and the computational potential to investigate these at the current time. We separate the notions of total shake, single shake, and double shake, explain how these are observable in high-quality experimental spectra, and how these calculations can be used experimentally.