第6章脉冲响应函数.docVIP

脉冲响应函数的辨识 6.1辨识问题的提法 下图所示，、将作用在系统上的一切随机干扰和噪声，用一个作用于系统输出的等效随机干扰源来代替。其中，输入信号是过程的运行操作信号，是可以直接观测的确定性变量；是过程的实际输出，是不能被观测到的；是过程的观测输出，混有随机噪声。 由此可以提出辨识问题： 在已知输入、输出的观测量、以及（可以根据脉冲响应过渡历程时间的先验知识作粗略估计）的情况下，要求估计出脉冲响应函数。 下面介绍两种辨识脉冲响应函数的常用方法：相关分析法和最小二乘法。 6.2用相关分析法辨识脉冲响应函数 相关函数是基于一种统计的描述，是由输出信号同其余变量之间的关系确定脉冲响应函数。假定噪声是一个零均值平稳随机过程，并与不相关，且过程是线性时不变的、因果性的系统，过程的未知脉冲响应函数为，则过程的输入、输出和脉冲响应函数之间的基本关系如下： （6.1） (6.2) 把变量用代换，得 (6.3) 由于已经假设与输入信号不相关，因此对应的相关系数，是可得维纳-霍夫方程。 (6.4) 若将（6.4） (6.5) 式中为的采样周期，； (6.6) (6.7) M为足够大的整数，为计算起点。 6.3 用最小二乘法辨识脉冲响应函数 假定对连续信号和以为周期进行采样，当足够小时，和在采样期间可看成是常数。 在，并令时，有 (6.8) 将观测数据代入，可形成m+1线性方程组，写成矩阵方程形式： (6.9) 定义误差矢量为： (6.10) 要求指标函数J相对于h达到最小： (6.11) 其定义为一个矢量函数，使得由于干扰噪声引起的误差的平方和最小。可直接得到最小二乘估计： (6.12) 6.4 最小二乘辨识与相关分析法辨识的关系 已知脉冲响应的最小二乘估计为： (6.13) 假设输入输出序列，是平稳的和遍历的随机序列，由m个数据构成输入自相关函数，输入输出的互相关函数关系为，则 (6.14) (6.15) 其中为任意正整数，表示在时间序列中相关函数计算的起点。相关分析法辨识具有最小二乘估计的性质，最小二乘的解也可以建议起离散型Wiener-Holf方程。 6.5激励信号的选择 选择激励信号是为了提高辨识精度和缩短辨识时间。 1 随机白噪声作激励信号 白噪声的相关函数为： (6.16) 其中，为Kronecke函数，为白噪声的方差。 采用白噪声作为激励信号的相关分析方法辨识过程的脉冲响应函数的结构图如下： 图中为采样周期，为方差。 2随机信号作为激励信号 用随机信号作为激励信号计算互相关函数，理论上要用无限长时间的观测数据。为了减少时间，可用“周期性的随机信号”，即伪随机信号作为激励信号。 在（0，T）时间内为白噪声，在此时间以外是周期函数，相关函数的计算如下： (6.17) 伪随机信号的自相关函数不需要在无穷大时间内计算，只需在一个周期为T内计算。类似可得： (6.18) 用这样的信号作为激励信号，可得 (6.19) 若T≥，当时，有，则 (6.20) 3随机二位式序列（PRBS）作激励信号 目前实际使用最多的是所谓伪随机二位式序列(PRBS)。其中，最大长度二位式序列，简称“M序列”。由自相关函数定义，可以计算出“M”序列的自相关函数： (6.21) 可以看出，如果，且时，“M序列”的自相关函数就接近于伪随机白噪声信号的自相关系数。 6.6 用伪随机二位式序列（PRBS）辨识脉冲响应函数 假设输入、输出信号的采样周期与“M序列”的时钟脉冲同步，即相同，且令 。取“M序列”的长度为 (6.22) 根据周期函数的相关函数计算公式和离散化相关函数计算公式，并取起点时刻为，即可得到 (6.23) (6.24) 计算“序列”的自相关函数 (6.25) 式中a为“M序列”的幅值。 由于， (6.26) 并由，可得到脉冲响应序列的计算公式： (6.27) 令，则可使互相关函数的计算更为简便，得到： (6.28) 关于“M序列”参数的确定可参考如下三个因素： 时钟脉冲周期；序列的长度；幅值。 - 5 - ＋ ＋ 过程g(t