Traditionally the Volterra time and frequency domain analysis tools cannot be applied to severely non-linear systems. In this paper, a new method of building a time-domain NARX MISO model for a class of severely SISO non-linear systems that exhibit subharmonics is introduced and it is shown how this allows the Volterra time and frequency domain analysis to be extended to this class of non-linear systems. The new approach is based on decomposing the original single input based on a Fourier analysis to provide a set of modified input signals which have the same period as the output signal. A MISO NARX model can then be constructed from the decomposed multiple inputs and the single output signal. The resulting MISO model is shown to meet the basic requirement for the existence of a Volterra series representation from which important frequency domain properties can be derived, explained and discussed. This is done by first introducing the derivation of generalized frequency response functions (GFRFs) from time domain MISO NARX models. The steady state response synthesis problem using the input spectrum and the MISO GFRFs is then investigated in order to verify the effectiveness and accuracy of the MISO modelling approach for severely non-linear systems. Finally a new frequency domain analysis method is introduced for systems that exhibit subharmonic oscillations.