A complex system contains generally many elements that are networked by their couplings. The time series of output records of the system’s dynamical process is subsequently a cooperative result of the couplings. Discovering the coupling structure stored in the time series is an essential task in time series analysis. However, in the currently used methods for time series analysis the structural information is merged completely by the procedure of statistical average. We propose a concept called mode network to preserve the structural information. Firstly, a time series is decomposed into intrinsic mode functions and residue by means of the empirical mode decomposition solution. The mode functions are employed to represent the contributions from different elements of the system. Each mode function is regarded as a mono-variate time series. All the mode functions form a multivariate time series. Secondly, the co-occurrences between all the mode functions are then used to construct a threshold network (mode network) to display the coupling structure. This method is illustrated by investigating gait time series. It is found that a walk trial can be separated into three stages. In the beginning stage, the residue component dominates the series, which is replaced by the mode function numbered M14 with peaks covering ～680 strides (～12 min) in the second stage. In the final stage more and more mode functions join into the backbone. The changes of coupling structure are mainly induced by the co-occurrent strengths of the mode functions numbered as M11, M12, M13, and M14, with peaks covering 200–700 strides. Hence, the mode network can display the rich and dynamical patterns of the coupling structure. This approach can be extended to investigate other complex systems such as the oil price and the stock market price series.