Ordinary differential equations (ODEs) are widely used to model the dynamic behavior of a complex system. Parameter estimation and variable selection for a “Big System” with linear ODEs are very challenging due to the need of nonlinear optimization in an ultra-high dimensional parameter space. In this talk, we will introduce a parameter estimation and variable selection method based on a matrix-based approach, about its theoretical properties and dramatic numerical improvement compared to classical methods.