本发明提出了一种基于切比雪夫多项式展开的时变神经动力学系统辨识方法。该方法首先用Volterra级数表征由仿真输入输出锋电位序列组成的时变神经系统，将前馈、反馈核函数用不同的Volterra核表示出来；接着，用Laguerre基函数对时变Volterra核进行展开，得到时变广义Laguerre–Volterra模型；然后，用切比雪夫多项式对时变广义Laguerre–Volterra模型的时变参数进行展开，将时变模型转化为时不变模型；最后，利用向前正交回归算法选出有意义的模型项，并利用广义线性拟合算法对时不变参数进行估计，进而通过反向求解得到时变参数及原始时变核函数。与现有的自适应滤波技术相比较，本方法对较强非平稳神经系统信号有更好的跟踪能力，可以实现对时变系统核函数的准确跟踪，为神经系统建模，特别是海量高维数据的系统建模提供新的研究方法，对揭示大脑完成信息处理的复杂神经动力学机制有重要的意义。