The aspirations of users to expand their data rate up to 1 Gbps in the fifth generation (5G) mobile networks open the horizon to exploit the row bandwidths in the millimeter wave (mmW) frequency bands. Several measurement campaigns have been launched to calculate the power delay profile and the large scale fading of the channels (path loss, shadowing) in mmW bands. In this paper, we design a 3D mmW multi-input multi-output channel model to extract the analytical narrowband impulse response. Based on a reference model, we establish a Sum-of-Sinusoids channel simulator in mmW frequency bands that is capable of generating the channel impulse response, which reflects the real conditions. The numerical results illustrate the outcomes of the simulator at different mmW frequency bands.