Cubic spline interpolation is able to recover temporally and spectrally resolved soft x-ray fluxes from an array of K-edge filtered x-ray diodes without the need for a priori assumptions about the spectrum or the geometry of the emitting volume. The mathematics of the cubic spline interpolation is discussed in detail. The analytic nature of the cubic spline solution allows for analytical error propagation, and the method of calculating the error for radiation temperature, spectral power, and confidence intervals of the unfolded spectrally resolved flux is explained. An unfold of a blackbody model demonstrates the accuracy of the cubic spline unfold. Tests of cubic spline performance using spectrally convolved detailed atomic model simulation results have been performed to measure the method's ability to conserve spectral power to within a factor of 2 or better in line-dominated regimes. The unfold is also demonstrated to work when information from the x-ray diode array is limited due to high signal-to-noise ratios or the lack of signal due to over-attenuation or over-filtration of the x-ray diode signal. The robustness of the unfold with respect to background subtraction and raw signal processing, signal alignment between diode traces, limited signal information, and initial conditions is discussed. Results from an example analysis of a halfraum drive are presented to demonstrate the capabilities of the unfold in comparison with previously established methods.