The proposed model generalizes the traditional linear dynamical system model. An efficient algorithm is designed for parameter estimation. The model works very well on high-dimensional neuro-imaging data. It is also a standard tool for big time-series data analysis in many domains.