说明：  对复杂系统进行辨识，寻找降阶实现。 M序列-伪随机序列（相比白噪声，工程上更容易实现）输入作系统辨识， 这里有几个处理： 1、 原始数据中的输出要减去y(1)，消除零漂。 2、 M序列取几个周期，因为刚开始不稳定，多以从第二个周期开始，显然，输出也是N周期的。 3、 算自相关和互相关的时候，最好自己写代码，因为是周期信号做自/互相关，用指令只能得到有限信号的自/互相关，需要作拓展，但是其实也是近似的。最好还是用代码自己写，清晰一点。 4、 脉冲响应得到Hankel阵时，g(0)是不取的（课本P201：It is important to mention that even if , does not appear in the Hankel matrix.）。 5、 取Hankel阵46*46（足够大的维数），求基于Hankel阵的平衡实现，然后直接截断成5阶。（用Hankel matrix说明原因。） 6、 得到是离散的模型，进行连续化，同时修改增益令辨识结果的低频增益与原始数据的低频增益相同。
(Identify the complex system and find the realization of order reduction)