The occurrence and risk of recurrence of brain related injuries and diseases are difficult to characterize due to various factors including inter-individual variability. A useful approach is to analyze the brain electroencephalogram (EEG) for differences in brain frequency bands in the signals obtained from potentially injured and healthy normal subjects. However, significant shortcomings include: (1) contrary to empirical evidence, current spectral signal analysis based methods often assume that the EEG signal is linear and stationary; (2) nonlinear time series analysis methods are mostly numerical and do not possess any predictive features. In this work, we develop models based on stochastic differential equations that can output signals with similar frequency and magnitude characteristics of the brain EEG. Initially, a coupled linear oscillator model with a large number of degrees of freedom is developed and shown to capture the characteristics of the EEG signal in the major brain frequency bands. Then, a nonlinear stochastic model based on the Duffing oscillator with far fewer degrees of freedom is developed and shown to produce outputs that can closely match the EEG signal. It is shown that such a compact nonlinear model can provide better insight into EEG dynamics through only few parameters, which is a step towards developing a framework with predictive capabilities for addressing brain injuries.