一辩识离散线性系统脉冲响应特性的相关分析法.DOCVIP

实验一 辩识离散线性系统脉冲响应特性的相关分析法 实验目的 掌握利用相关分析法辩识离散线性系统冲激响应特性的基本步骤和使用要点； 了解递推和批量两种算法如何具体编程； 熟悉估计算法中各种参数选择对辩识效果的影响。 实验原理 待辩识系统的基本情况和辩识工作的要求 待辩识系统的运行情况如图1所示，设其在静态工作点（U0，Y0）附近作局部线性化所得动态模型为 y(k)+a1y(k-1)+a2y(k-2)=b1u(k-1-d)+b2u(k-2-d) (1) U0+ u U0+ u(k) Y(k) Y(k)=Y0+y(k) 待辩识系统U0 待辩识系统 U0 η(k)u(k) η(k) u(k) Y’ Y’(k)=Y0+y(k) +η(k) 图1 待辩识系统的基本运行情况 为了进行系统冲激响应特性的辩识，在对象输入端施加电平幅值V=1的M序列信号（其波形见图2）作为动态测试信号u(k)。本次实验的基本任务就是在输入、输出的实际观测值Y’(k)=Y0+y(k) +η(k) 和u(k)的基础上，用批量或递推算法估计出脉冲响应序列{h(k)}。 u(k) 1514131211100987654321 15 14 13 12 11 10 0 9 8 7 6 5 4 3 2 1 图2 M序列信号的波形 M序列测试信号的生成方法 M序列又称最大长度伪随机二进制序列，它实际上是一种特殊的周期序列，其电平值u(k)只取+V和-V (或者说其逻辑电平值只取1或0)。 M序列的当前值可以由过去的值进行异或运算得到。选取恰当的k0，便可以产生最长周期为Np=2n-1的M序列。很容易由计算机程序来产生M序列。如取n=4，k0=3，，便可按 依次推出，最后得出图2所示的 M序列。 M序列的循环长度为NP=2-1=15,由此得到N=4. M=[1 1 1 1 –1 –1 –1 1 –1 –1 1 1 –1 1 -1],寄存器初始状态为（1111）。 基本原理和基本关系式 在相关法中，采用伪随机信号M序列对系统模型进行在线辩识。由于M序列的自相关函数很接近于周期性的δ函数，可近似地当作白噪声信号叠加在系统输入信号上，按其与输出信号间的互相关函数，可得到系统的脉冲响应曲线。 设对象输入u为具有平稳遍历性的随机信号，其自相关函数为 (2) 当系统初始过渡过程结束后，输入和输出之间的互相关函数为 (3) 式中E[●]表示统计平均， ● 表示时间平均。由维纳 – 霍甫方程有 (4) 当u为伪随机信号M序列时，自相关函数的离散值为 (5) 式中,V、△和?Np分别为M序列的电平值、基本电平时间和离散周期长度，ι为时差的离散值。 相应的功率谱密度函数为 (6) 式中 (7) 若将M序列的自相关函数近似为 (8) 并取的周期大于对象的整定时间，即 (9) 则有 (10) 这样就可以得到脉冲响应曲线 (11) 批量算法和递推算法 批量算法（Batch algorithm）： (12) (13) 2) 递推算法（Recursive algorithm）： (14) 辩识对象运行情况的仿真 取系统对象模型中初值为。将M序列作为，即可求。设工作点（U0,Y0）分别为0.1和0.2，则对象的实际输入、输出U(k)和Y’(k)分别为 (18) 其中正态白噪声序列η(k)根据中心极限定理采用12个在[-0.5,0,5]中均匀分布的随机数叠加的方法生成，其方差可由遭受幅度因子调节。 在待辩识对象在大范围内是非线性的情况下，通过辩识得出的是对象在工作点附近的局部线性化动态模型，因此还应去处工作点Y0 ，以得出输出端动态分量的实测值。但工作点无法单独测量，只能根据最新得到的推算： (19) 或递推计算如下： (20) 式中Y0(k)的初值Y0(0)可设为0。于是，输出端动态分量的实测值可近似估算为 (21) 实验步骤 批量算法如图3： 初始化 初始化 生成数据 u’ 生成数据 u’(k-d) , Y’(k) , k=1-r-d,…,-1,0,…,Np-1 k k1 ← Np + r + d k k ← k +1 生成新的数据 k-k1 k-k1iTp? 打印结果 打印结果 图3 批量算法流程图 图3 批量算法流程图 第一步：产生M序列，如图所示 NP=15,N=4 第二步：求系统输出，编程