系统辨识课件

系统辨识1系统辨识1课程主要内容第一章概述

第二章过渡响应法和频率响应法第三章辨识线性系统脉冲响应函数的相关分析法结束第四章线性系统参数估计的最小二乘法第五章线性系统的状态估计法2课程主要内容第一章概述第二章过渡响应法和频率响第一章概述一、建模的必要性二、模型三、建模方法四、系统辨识的内容（或步骤）3第一章概述一、建模的必要性二、模型三、建模方法四、系一、建模的必要性课程的核心问题是建模，主要是辨识建模。系统辨识是研究辨识建模的理论和方法。数学模型的主要用途：控制理论与控制工程就一直围绕着

建立模型和控制器设计这两个主题来发展，它们相互依赖、相互渗透并相互发展。用来预报实际系统物理量研究实际系统往往需要事先知道一些物理量的数值，而其中有些量可能无法直接测量或测不准，所以需要建立数学模型来预报。第一章概述4一、建模的必要性课程的核心问题是建模，主要是辨识建模。系统辨为了设计控制系统目前,对被控系统的控制器的设计方法的选取,以及如何进行具体的控制结构和参数的设计都广泛依赖于对被控系统的理解及所建立的被控系统数学模型。用于分析实际系统工程上在分析一个新系统时，通常先进行数学仿真，仿真的前提必须有数学模型。第一章概述建模问题在控制器设计中起着非常重要的作用,是设计中首先需要解决的问题；是成功地进行控制器设计的关键之一。

5为了设计控制系统目前,对被控系统的控制器的设计方法的选取第一章概述系统的模型一般分物理模型与数学模型物理模型：指用物理、化学、生物等材料构成的用于描述系统中的关系和特征的实体模型。

模型：就是把系统实体的本质信息简缩成有用的描述形式，数学模型：描述系统中一些关系和特征的数据模型。控制领域的数学模型就是指能用来描述系统的动态或静态特性和行为的数学表达式或方程。是进行系统分析、预报、优化及控制系统设计的基础。

二、模型是一种简化描述。6第一章概述系统的模型一般分物理模型与数学模型物理模型三、建模方法1、理论建模法：通过对系统内在机理的分析，按照已知的一些物理定律导出各物理量关系来建立数学模型。理论建模法建立的模型称为机理模型。一般在理论建模中,根据模型应用的目的和精度要求,仅考虑系统中起主导作用的有限的几个因素即可。缺陷：当验前信息不足时，用理论建模法会遇到很大困难。对于比较复杂的过程，必须对机理模型简化，这就使得机理建模与实际过程间有一定的误差。第一章概述7三、建模方法1、理论建模法：通过对系统内在机理的分析，理论建模通常只能用以建立比较简单系统的模型（白箱问题）。由于许多系统的机理和所处的环境越来越复杂，因此，理论建模法的运用亦越来越困难，其局限性越来越大,需要建立新的建模方法。第一章概述在被建模的装置尚不存在（设计阶段）或虽存在但无法进行实验时，理论建模是取得模型的唯一途径，是验前问题中唯一可行的方法。理论建模的难点在于对有关学科知识及实际经验的掌握，故不属于课程的讨论范围。在理论建模方法难以进行或难以达到要求的情况下，系统辨识建模方法就幸运而生。8理论建模通常只能用以建立比较简单系统的模型（白箱问题）2、辨识建模法：

对被控系统进行测试，利用观测数据，通过辨识技术去构造系统模型的方法。系统辩识是研究怎样利用对未知系统的试验数据或在线运行数据(输入/输出数据)建立描述系统的数学模型的科学。系统辩识亦称为实验建模方法,它是“系统分析”和“控制系统设计”的逆问题。第一章概述是现代控制理论的一个分支。92、辨识建模法：对被控系统进行测试，利用观测数据，通过1）完全辨识问题：第一章概述完全不了解系统的任何基本特性（定常—时变；线性—非线性；确定—随机等）。这类问题称为黑箱问题。这是一个极难解决的问题，通常需要对系统作某些主观的先验假设。2）部分辨识问题：系统的某些基本特性假定是已知的，但不知动态模型的阶次或有关的系数。这类问题称为灰箱问题。显然比黑箱问题容易解决。根据对系统事先了解的程度（先验知识）可将辨识问题分成二类：完全辨识问题和部分辨识问题。101）完全辨识问题：第一章概述完全不了解系统的任大部分工程系统及工业过程都属于灰箱问题。通常对系统的结构会有很多了解，因此可推导得系统特定的数学模型。在这种情况下只要定阶和确定模型中的一组参数。从而模型化问题简化为参数估计。因此参数估计是一个最重要的问题。第一章概述有效的辨识策略：尽可能地掌握系统的先验知识，即尽可能地使系统“白化”；有效的辨识方法：“灰箱”方法。将两种方法结合起来，互为补充。对依然“黑”的部分，用理论建模方法不能确定的部分和参数，采用系统辨识方法。11大部分工程系统及工业过程都属于灰箱问题。通常对系统的第一章概述系统辨识的框图12第一章概述系统辨识的框图12模糊数学创始人L.A.Zadeh第一章概述1962年Zadeh从数学的角度定义:辨识就是在输入输出数据的基础上，从一组给定的模型类中，确定一个与所测系统等价的模型。1978年瑞典的李龙(Ljung)提出:系统辩识的三个要素——数据、模型类和准则。系统辩识是按照一个准则，在模型类中选择一个与数据拟合得最好的模型。拟合的好坏是一个不定的概念，所以要用准则来判别。

3、系统辨识的定义所谓辨识建模是从实验数据出发，根据辨识的目的以及对过程已有的验前知识，预先给出一个模型类（线性的、非线性的、定常的、时变的、连续的、离散的…）进行拟合。

13模糊数学创始人第一章概述1962年Zadeh从数学四、系统辨识的内容（或步骤）第一章概述它是一个迭代过程。大致包括：试验设计，模型结构确定，参数估计和模型验证。满意辨识目的及先验知识确定模型结构和准则模型的参数估计模型验证最终模型不满意数据预处理输入输出数据检测试验试验设计辨识的一般步骤14四、系统辨识的内容（或步骤）第一章概述它是一个迭代过大致包括：试验设计，模型结构确定，参数估计和模型验证。1、试验设计第一章概述

1）选择变量：以提取有效的信息（数据）为目的。首先根据试验对象，确定所要观测的变量。（u是人为给定的，y是观测的，y的选取不同会改变输出矩阵C的结构和数值。）通常为得到试验设计前的必要的知识，必须进行一些预备性试验（摸底）。四、系统辨识的内容（或步骤）15大致包括：试验设计，模型结构确定，参数估计和模型验证。第一章概述预备性试验：可用一些简单方法（阶跃响应，频率响应等）获得系统的如下信息：主要时间常数（系统频宽，与试验长度有关）允许的输入信号幅度（系统的线性范围）过程的非线性与时变性（有助于模型类的选择）噪声水平（以便用多大的输入，使得观测量有多大的信噪比）变量之间的延迟（滞后环节参数）2）输入信号的选择（阶跃、方波、脉冲、PRBS）。16第一章概述预备性试验：可用一些简单方法（阶跃响应，频第一章概述

3）采样速度的选择（要采集数据就有采样速度选择问题）。实际上先采用较短的采样间隔，在数据分析时，可根据需要隔几个取一个数据。

4）试验长度的确定（试验时间问题）。辨识精度与试验时间的长短有关。2、模型结构确定根据辨识的目的及对被辨识系统的先验知识，确定系统所属的模型类模型结构的选择主要取决于应用的目的及精度要求。通常模型精度与复杂性要折衷考虑。17第一章概述3）采样速度的选择（要采集数据就有采样第一章概述常用的模型类：

参数的或非参数的

线性的或非线性的

连续的或离散的

确定的或随机的

I/O的或状态的

时变的或定常(时不变)的集中参数的或分布参数的

频率域的或时间域的等等。18第一章概述常用的模型类：18第一章概述根据系统的空间、时间的离散化情况，模型可分为三类：1）集中参数的连续时间模型：空间变量是离散的，时间变量连续。如常微分方程，代数方程。2）集中参数的离散时间模型：时、空变量均离散。如差分方程，代数方程。3）分布参数模型：时、空变量均连续，如偏微分方程。它可以在空间上离散化，简化成分块集中参数，所以对它的辨识不介绍。a19第一章概述根据系统的空间、时间的离散化情况，模第一章概述3、参数估计模型结构确定后，其中未知部分就要通过观测数据进行估计。通常未知部分是以未知参数出现，故辨识工作就成了参数估计。4、模型验证一个模型辨出来后，是否可靠必须进行多次验证。参数估计的要求就是要辨识出来的模型与实际过程在某种意义下最“接近”。所以必须有个准则衡量。通常一个模型用一套数据进行辨识，然后用另一套数据来验证和修改。

20第一章概述3、参数估计模型结构确定后，其中未知第二章过渡响应法和频率响应法

§2

1过渡响应法（时域法）§2

2频率响应法（频域法）§2

3多输入多输出线性系统传函（矩阵）的辨识21第二章过渡响应法和频率响应法§21过渡响应法（时域模型可以有不同的形式，不同的模型适于不同的系统。古典辨识方法：采用时域法和频率法来辨识线性系统的传递函数。原则上只适用于SISO线性系统。SISO系统通常采用传递函数。MIMO系统通常采用状态空间表达式。由实验来建立数学模型——传递函数，可以为更复杂的系统辨识做预备性实验，它是现代系统辨识的基础，属于连续系统的数学模型的辨识领域。第二章过渡响应法和频率响应法22模型可以有不同的形式，不同的模型适于不同的系统。古典辨识方法第二章过渡响应法和频率响应法试验信号的选用：对系统模型的研究方法不同，输入试验信号也相应分成非周期的和周期的两种。用时域法建模：输入信号为非周期的。主要采用阶跃和方波（近似脉冲）函数。用频域法建模：输入信号用周期的。主要用正弦波，二进制周期函数。它们又分为单频和多频（组合正弦波及周期方波）23第二章过渡响应法和频率响应法试验信号的选用：对系统§2

1过渡响应法（时域法）采用非周期试验信号，通过系统的动态响应研究系统的模型。一、非参数模型的辨识在时域中建立线性系统非参数模型时，用很简便的方法就可得到脉冲响应曲线，阶跃响应曲线、方波响应曲线或它们的离散采样数据表。对于线性系统，脉冲响应，阶跃响应和方波响应之间是可以相互转换的。脉冲响应：可以采用幅值相当大，宽度很窄的方波来近似δ函数。第二章过渡响应法和频率响应法24§21过渡响应法（时域法）采用非周期试验信号，通过第二章过渡响应法和频率响应法二、由阶跃响应曲线辨识传函1、试探法

工业中常用的模型类：（即便是高阶系统也用低阶模型去逼近）由非参数模型转变成参数模型，包括确定传函的结构及参数。先观察试验所得响应曲线的形状特征，据此判断，从模型类中确定一种结构。然后进行参数估计，最后验证数据拟合程度，反复多次，直至误差e(t)最小（验证数据拟合可只取若干点）。25第二章过渡响应法和频率响应法二、由阶跃响应曲线辨识传函1第二章过渡响应法和频率响应法1）若阶跃响应曲线特征为：曲线逐渐上升到稳态值：可采用结构：待估参数为：K，T稳态增益：将试验曲线标么化，即,26第二章过渡响应法和频率响应法1）若阶跃响应曲线特征为：曲第二章过渡响应法和频率响应法要确定T，只要一对观测数据：y*(t1)，t1

则标么化后响应：可得：由若取y*(t1)=0.63，则T=t1验证数据拟合如何，可在t=T/2和t=2T二点进行：若拟合不好，则应另选模型结构类。127第二章过渡响应法和频率响应法要确定T，只要一对观测数第二章过渡响应法和频率响应法待估参数为：K，T，τ

稳态增益：将试验曲线标么化，即2）实验曲线是一条S形非周期曲线

ⅰ）可选用模型类：

则为了确定T和τ，必须将两个坐标值（观测值）代入，则28第二章过渡响应法和频率响应法待估参数为：K，T，τ稳态第二章过渡响应法和频率响应法两边同取对数得：根据两对观测值

y*(t1)和y*(t2)

，可求出T和τ。29第二章过渡响应法和频率响应法两边同取对数得：根据两对观测30第二章过渡响应法和频率响应法若选y*(t1)=0.39，y*(t2)=0.63，则模型验证：tT/20.8TT2Ty*(t)0.390.550.630.87由则3030第二章过渡响应法和频率响应法若选y*(t1)=0第二章过渡响应法和频率响应法待估参数为：K，T，

究竟选一阶惯性带延时的模型结构，还是选二阶模型，事先无法确定，完全看两种模型与试验曲线拟合程度，哪个精度高，选哪个。由于大多数工业过程的试验曲线是过阻尼的，即

1，∴只讨论此种情况，而

<1的传函辨识比较麻烦。ⅱ）也可选用模型类：

S形曲线本身就说明是过阻尼（

1）。若

<1,为振荡环节，响应曲线为振荡曲线，不是S形曲线。31第二章过渡响应法和频率响应法待估参数为：K，T，第二章过渡响应法和频率响应法稳态增益：将试验曲线标么化，即,由于直接求T、

不方便，所以转而取求传函的两个极点

1、2

，再由极点与参数T、

的关系求出T、

。从而转化成求中的T、

。分解因式：则∵

1，∴1、2>0且为实数。

32第二章过渡响应法和频率响应法稳态增益：将试验曲线标么化第二章过渡响应法和频率响应法代入可得：它的单位阶跃响应为：改写为：令

2=1(

>1)，代入上式得：两边同取对数得：33第二章过渡响应法和频率响应法代入可得：它的单位阶跃响第二章过渡响应法和频率响应法可见，当t

时，是一条直线。斜率：k=

1

，截距：在坐标纸上，根据数据y*(t)，画出t较大时的图形，作其渐进线，即可得斜率

k和截距

b

。则可得：

2

T，

当t

时,34第二章过渡响应法和频率响应法可见，当t时，是一第二章过渡响应法和频率响应法若用常用对数，则当t

时：则：缺点：计算G(s)时采用的点都是

t较大时的点，而当

t较大时，往往1

y*(t)

的值较小，这就会产生较大的误差。b35第二章过渡响应法和频率响应法若用常用对数，则当t第二章过渡响应法和频率响应法2、Laplace变换的极限定理法（终值定理法）利用

Laplace变换的极限定理，由非参数模型的单位阶跃响应，求参数模型——传递函数。它克服了试探法需选择模型类的不足，但它仅适用于下述一种模型类。设线性SISO定常系统的传函结构为：特点：系统只有极点、无零点。36第二章过渡响应法和频率响应法2、Laplace变换的极第二章过渡响应法和频率响应法当输入u(t)=1(t)时，输出y(t)为：终值定理为：对于阶跃响应:代入上式得：

K0

37第二章过渡响应法和频率响应法当输入u(t)=1(t)第二章过渡响应法和频率响应法在已存在的系统G(s)的基础上构造一个新系统G1(s)，当输入u(t)=1(t)时，其单位阶跃响应为：（y1(t)与y(t)的关系）求G1(s)的稳态增益K1

:

K1

当输入u(t)=1(t)时，输出y(t)为：38第二章过渡响应法和频率响应法在已存在的系统G(s第二章过渡响应法和频率响应法

G1(s)求拉氏变换：求得G1(s)与G(s)的关系:当输入u(t)=1(t)时，输出y(t)为：39第二章过渡响应法和频率响应法G1(s)求拉氏变第二章过渡响应法和频率响应法再利用终值定理可求得G1(s)的稳态增益K1：

a1

当输入u(t)=1(t)时，输出y(t)为：

把y1(t)定义成与y(t)有联系，当然G1(s)也与G(s)有联系，而输入均为u(t)=1(t)，再利用终值定理求G1(s)中的参数，从而也就求出G(s)中的参数。求得G1(s)的稳态增益K1与G(s)中的参数关系。40第二章过渡响应法和频率响应法再利用终值定理可求得G1(s第二章过渡响应法和频率响应法

K2

同理，在系统G1(s)基础上构造一个新系统G2(s)。G2(s)的单位阶跃响应为：求拉氏变换：当输入u(t)=1(t)时，输出y(t)为：41第二章过渡响应法和频率响应法K2同理，在系统G142第二章过渡响应法和频率响应法再利用终值定理可得：

a2

G2(s)求得G2(s)与G1(s)的关系:当输入u(t)=1(t)时，输出y(t)为：4242第二章过渡响应法和频率响应法再利用终值定理可得：第二章过渡响应法和频率响应法当输入u(t)=1(t)时，输出y(t)为：

Kr

同理，在系统Gr-1(s)基础上构造一个新系统Gr(s)。Gr(s)的单位阶跃响应为：再用终值定理，由数学归纳法可得：

ar

43第二章过渡响应法和频率响应法当输入u(t)=1(t)第二章过渡响应法和频率响应法特点：ⅰ）每求一次Ki，要计算一次面积，所以计算量大，而且误差随着积分次数增大而增大。故仅适用于低阶模型的辨识。ⅱ）使用过程受到一定的限制，仅适用于特定的模型结构（即传函G(s)只有极点，而没有零点的情况）。由上述（n+1）个方程可求出（n+1）个待估参数：

K，a1，…

，an44第二章过渡响应法和频率响应法特点：ⅱ）使用过程受到一定的第二章过渡响应法和频率响应法K1的物理意义：由可知，K1为阴影部分的面积。

（几何意义）将G(s)改写成极点形式：显然：（物理意义）145第二章过渡响应法和频率响应法K1的物理意义：由第二章过渡响应法和频率响应法三、由脉冲响应曲线辨识传函1、矩法脉冲响应g(t)可由单位阶跃响应微分后求得，也可用窄方波响应来近似。（方波宽度<<对象过渡过程时间）1）原理：令mi称为g(t)的i阶矩（时域的矩），可以由试验曲线g(t)的数值积分求得。实际上mi是G(s)泰勒展开式中的各个系数，于是可得G(s)的泰勒展开式：mi46第二章过渡响应法和频率响应法三、由脉冲响应曲线辨识传函第二章过渡响应法和频率响应法矩法就是找一个近似的传函，使得它的泰勒展开式式的系数尽可能与实际传函G(s)泰勒展开式的系数相等。使得2）定模型类：通常总是s的有理分式：47第二章过渡响应法和频率响应法矩法就是找一个近似的传第二章过渡响应法和频率响应法3）参数估计：待估参数：a1，…，an，b0，b1，…，bm，共

m+n+1个。求的泰勒级数展开式：令由（m+n+1）个方程可求得（m+n+1）个待估参数。可以由脉冲响应数据求得。48第二章过渡响应法和频率响应法3）参数估计：待估参数：a第二章过渡响应法和频率响应法若取则待估参数有4个：a1，a2，b0，b1。可令G(s)与的泰勒展开式的前四项系数相同。可得：解上述方程组可解出a1，a2，b0，b1。49第二章过渡响应法和频率响应法若取则待估参数有4个：a1第二章过渡响应法和频率响应法展开式的前n项系数也可采用升幂长除法求得。50第二章过渡响应法和频率响应法展开式的前n项系数也可采用升第二章过渡响应法和频率响应法特点：ⅰ）简单明了。（思路清晰）ⅱ）计算m0，m1，m2，…的计算量大，而且数值积分的误差越来越大。通常只适合于低阶系统。ⅲ）若取G(s)与的泰勒展开式的前n项的系数分别相等，则有:ⅳ）模型类为

s的有理分式，所以该法适应性广。∵频率特性∴系统稳态时，t

相当频率特性的低频段，即

0，故s=j

0。复域中s

0相当于时域中的t

，所以估计出的参数精度主要局限于t

的稳态情况。（稳态精度较好）51第二章过渡响应法和频率响应法特点：ⅱ）计算m0，m1，第二章过渡响应法和频率响应法2、Z变换法Z变换法是使得构造的模型传函

的脉冲响应与实验数据中的

g(t)在某些采样点上相等。Z变换的定义（脉冲传递函数G(z)是脉冲响应函数的采样序列g(kT)的Z变换）：1）原理：

Z变换法不是直接去求而是先求，即找一个近似的，使得。52第二章过渡响应法和频率响应法2、Z变换法Z变换法第二章过渡响应法和频率响应法2）定模型类：取为z

1的有理分式：3）参数估计：待估参数：a1，…，an，b0，b1，…，bm，共m+n+1个。求关于z

1的展开式：采用z

1的升幂长除法求得：显然关于z

k的系数就是在t=kT时的采样值。53第二章过渡响应法和频率响应法2）定模型类：取第二章过渡响应法和频率响应法令与g(t)在采样点t=kT（k=0，1，2，…，n+m）相等，即从而求得的待估参数：b0，b1，…，bm，a1，…，an

ⅲ）只能保证在某些采样点上与实验数据拟合，采样点之间无法保证相拟合。若用较简单的模型去逼近较复杂的实验曲线时。可能得到一个不稳定的传函。为了减小误差，必须缩短采样周期。4）由求：可查表求与的对应关系。特点：ⅰ）计算简单。（无需象矩法那样求数值积分）ⅱ）采用的模型类为z的有理分式，故适应性广。

a54第二章过渡响应法和频率响应法令与g(t)在采§2

2频率响应法（频域法）第二章过渡响应法和频率响应法采用周期试验信号。通过系统稳态响应研究系统模型。一、非参数模型的辨识1、正弦波法每次输入某一频率的正弦波，测量输出的稳态振荡波形，从而获得输出与输入的幅值比（幅频特性）和相位差（相频特性）。1）单频信号：每次试验只得频率特性上的一个点，费时间。55§22频率响应法（频域法）第二章过渡响应法和频率响第二章过渡响应法和频率响应法它是相对于某一单频信号而言，包含有不同频率的单频信号。然后对输出量进行谐波分析（富氏分析），从中分析出各个频率的输出波形。2）多频信号每次试验可得频率特性上的几个点。一般在系统的主要频区上测试15个点，就可以建立系统的非参数模型。3）特点该方法未考虑系统的干扰噪声的影响，所以所得的模型数据中都含有噪声的干扰信号的影响。56第二章过渡响应法和频率响应法它是相对于某一单频信号而实验数据（非参数模型）参数模型方法误差（计算误差、近似误差）引入噪声（过程噪声、量测噪声）第二章过渡响应法和频率响应法由于实际系统的测量都是载噪的，而且噪声对观测数据的影响常常达到不可忽略的地步，因此当噪声影响足以使得要求的精度达不到时，就必须考虑噪声的影响。2、相关法考虑噪声时，系统的方框图

n(t)：干扰噪声（过程噪声和量测噪声）。由于输入是人为的加入信号，所以可以认为是无噪声的。系统辨识的任务之一就是设法排除干扰的影响。采用相关方法可以提高对噪声的抗干扰性。57实验数据（非参数模型）参数模型方法误差（计算误差、近似误差）第二章过渡响应法和频率响应法随机过程1)平稳随机过程

X(t)：均值:均方值:自相关函数:是

的单值函数平稳随机过程的数字特征的特点：均值为常数，自相关函数为单变量（

=t2

t）的函数（即不随时间的平移而变化）。58第二章过渡响应法和频率响应法随机过程1)平稳随机过程第二章过渡响应法和频率响应法2)各态经历的平稳随机过程

X(t)：（集均值）

依概率1收敛于时间均值

依概率1成立

对于各态历经的平稳随机过程，它的集平均（均值，自相关函数等），可以用一个样本的时间平均来代替。因为它的任何一个样本函数都包含了随机过程的所有状态以及各状态的分布情况。样本函数：每进行一次观测得到一族时间t的函数x(t)就是随机过程X(t)的一个样本函数。59第二章过渡响应法和频率响应法2)各态经历的平稳随机过程第二章过渡响应法和频率响应法白噪声是一种理想化数学模型，实际上如果某种噪声在比有用带宽宽得多的范围内具有比较“平坦”的谱密度，就可以近似当作白噪声处理。

白噪声：零均值，谱密度为非零常数，自相关函数为

函数的平稳随机过程，且各态历经。维纳

辛钦公式：谱密度和自相关函数是一傅里叶变换对。60第二章过渡响应法和频率响应法白噪声是一种理想化数学模第二章过渡响应法和频率响应法频率特性测试仪BT6：61第二章过渡响应法和频率响应法频率特性测试仪BT6：61第二章过渡响应法和频率响应法同理可得：=0

正弦、余弦的正交性62第二章过渡响应法和频率响应法同理可得：=0正弦、余弦第二章过渡响应法和频率响应法由于n(t)具有零均值和各态历经性，所以在足够多的周期下测取数据能使载噪的正、余弦的所以可得系统频率特性的非参数模型：由于它具有一定的抗干扰性，所以基本上可得到系统真实传函的非参数模型，是获得高精度模型的有效方法。b63第二章过渡响应法和频率响应法由于n(t)具有零均值和第二章过渡响应法和频率响应法二、由实验频率特性辨识传函工业上常见的模型类的频率特性：由非参数模型辨识参数模型。64第二章过渡响应法和频率响应法二、由实验频率特性辨识传函第二章过渡响应法和频率响应法1、用一阶环节（模型）拟合1）当实验所得频率特性近似为一个半圆时，可用模型:待估参数为T和K。为了提高精度，可在实验曲线的主频区内多取n个

k和

(

k)，取它们的平均值，即65第二章过渡响应法和频率响应法1、用一阶环节（模型）拟合第二章过渡响应法和频率响应法2）带有纯滞后环节当实验曲线为：待估参数为：T、K和

。可考虑系统也许有纯滞后特性，即采用模型：∵纯滞后环节对幅频特性无影响，只影响相频特性，即产生一个附加的滞后角

。∴将带纯滞后环节与不带纯滞后环节的频率特性相比较，就可以求出相应的参数和纯滞后时间

。66第二章过渡响应法和频率响应法2）带有纯滞后环节当实验曲第二章过渡响应法和频率响应法2）带有纯滞后环节待估参数为：T、K和

。注意这里的

k，

k和

k均为负角度，67第二章过渡响应法和频率响应法2）带有纯滞后环节待估参数第二章过渡响应法和频率响应法2、用二阶模型拟合待估参数为K、T1和T2。1）当实验曲线分布在三、四两个象限。且，可采用模型:68第二章过渡响应法和频率响应法2、用二阶模型拟合待估参数第二章过渡响应法和频率响应法2、用二阶模型拟合待估参数为K、T1和T2。1）当实验曲线分布在三、四两个象限。且，可采用模型:取

=

k，并将：K=r(0)

G(j

)

=r(

k)

(

)=

(

k)代入可求得：K，T1和T2。理论上只要一对观测数据

k，r(

k)和

(

k)即可求得参数，实际上可多取几个数据，然后求参数的平均值。69第二章过渡响应法和频率响应法2、用二阶模型拟合待估参数第二章过渡响应法和频率响应法2）带有纯滞后环节可考虑采用模型：先作不带滞后环节的：可得G

(j

)的曲线。从而得：

同理作图，70第二章过渡响应法和频率响应法2）带有纯滞后环节可考虑采§2

3多输入多输出线性系统传函（矩阵）的辨识第二章过渡响应法和频率响应法设MIMO系统有：

r个输入，m

个输出，对应yi与uj之间的传函为：Gij(s)（i=1,2,…,m,j=1,2,…,r）系统的传递函数矩阵：71§23多输入多输出线性系统传函（矩阵）的辨识第二章第二章过渡响应法和频率响应法显然，输出yi与r个输入的变化有关。由于线性系统满足迭加原理，故可采用逐次地改变每一个输入量，同时保持其他各个输入量为恒值。仅改变uj值，同时采集y1、…、ym的值，采用SISO的建模方法，就可以求得m个传函：G1j(s)，G2j(s)，…，Gmj(s)。在实际中，研究MIMO系统，往往采用状态空间法(时域法)，而不用古典的频域法(传递函数法)。由现代控制理论可知，实际上两种方法有密切联系，即MIMO系统的状态空间描述和传递函数描述之间是可以相互转化的。然后，逐次改变uj，j由1到r，就可以求出全部Gij(s)。

272第二章过渡响应法和频率响应法显然，输出yi与r个输第三章辨识线性系统脉冲响应函数的相关分析法§3

1相关分析法（CorrelationFunctionMethods）§3

2伪随机二位式序列产生的方法及性质§3

3用M序列辨识线性系统的脉冲响应函数§3

4用逆重复M序列辨识线性系统的脉冲响应函数§3

5用伪随机序列辨识系统的步骤§3

6多变量系统的辨识73第三章辨识线性系统脉冲响应函数的相关分析法§31相关不能“在线”辨识，所处理的只是离线估计问题。时域法和频域法的缺点：

所得的结果没有递推公式，只适合于非时变系统，不宜用于实时估计和自适应控制。影响系统正常工作状态，甚至可能破坏设备。抗干扰性差，通常不考虑系统的噪声干扰。实际上系统总是载噪的（过程噪声，量测噪声），所以在辨识中存在输入信号的信噪比与实际所允许的输入信号幅值的矛盾。（即：增大幅值→信噪比↑，同时→影响实际工况。）采用伪随机二位式序列（PRBS）作为输入来辨识系统的脉冲响应函数，既可保证输入信号能激励系统的所有模态反应，又不会干扰系统的正常运行，而且可以进行在线辨识。第三章相关分析法采用PRBS信号加上相关分析法的辨识，可克服上述的缺点。74不能“在线”辨识，所处理的只是离线估计问题。时域法和频域§3

1相关分析法（CorrelationFunctionMethods）系统噪声实际上存在各种难以精确描述的因素：数学模型中未加以考虑的各种干扰作用，模型线性化及其他类似的假设所引起的误差等（系统误差）。第三章相关分析法这些难以用确定性模型来精确描述的随机因素对过程的影响不容忽视，但这种影响一般都不会处于主导地位，因此处理它们的方式应尽量简化。一般在确定性模型上以迭加输入的方式来考虑噪声的影响，并将多个噪声源综合成一个等效的噪声源来描述。输入输出的量测误差等（量测误差）。75§31相关分析法（CorrelationFunctiou(t)：系统输入变量z(t)：系统真正的输出y(t)：系统输出量的实际测量值n(t)：随机干扰（噪声）(量测误差和系统误差的总误差)第三章相关分析法设输入u(t)是一个平稳随机过程，则由相关理论可知z(t)也是一个平稳随机过程。由控制理论可知：SISO系统：76u(t)：系统输入变量第三章相关分

平稳随机过程的数字特征：均值=const，自相关函数为单变量函数（不随时间而变化。）第三章相关分析法维纳—霍甫方程维纳—霍甫方程是相关分析法的理论基础。g(t)是一个确定量可提到均值号前面通常求均值和积分允许交换顺序输入输出的互相关函数Ruz(

)等于输入自相关函数Ruu(

)与系统脉冲响应函数g(

)的卷积。77平稳随机过程的数字特征：均值=const，自相关函数为单设输入u(t)是各态历经的平稳随机过程，则输出z(t)也是各态历经的。第三章相关分析法则可得：自相关函数的集平均，可以用样本的时间平均来代替。78设输入u(t)是各态历经的平稳随机过程，则输出z(第三章相关分析法上式表明，通过相关分析得到实测的输出y(t)与输入u(t)之间的互相关函数Ruy(

)，在n(t)与u(t)不相关的条件下（此条件常常能满足），等价于真实输出z(t)与输入u(t)间的互相关函数Ruz(

)，因此,相关分析法有较强的抗干扰能力，这是相关法的一个突出优点。设n(t)与u(t)不相关，且n(t)为零均值，则：Run(

)=079第三章相关分析法上式表明，通过相关分析得到实测的输出相关法与古典法的差别之一：抗干扰强

第三章相关分析法①古典法：

直接用u(t)和y(t)通过求g(

)。②

相关法：

不直接用u(t)和y(t)，而要用u(t)，y(t)先算出Ruy(

)、Ruu(

)后，再解维纳—霍甫方程利用u(t)，y(t)计算Ruy(

)、Ruu(

)过程，即相关运算器具有滤除干扰n(t)的作用。

思路：维纳—霍甫方程

解决了抗干扰问题

80相关法与古典法的差别之一：抗干扰强第三章相关分析法解决办法：选择输入u(t)若采用白噪声作为输入，则维纳—霍甫方程的解就简单多了。第三章相关分析法思路：维纳—霍甫方程

解决了抗干扰问题

引出了解积分方程难的问题问题：维纳—霍甫方程是一个积分方程（算得Ruy(

)、Ruu(

)，求解），求解积分方程比较困难。81解决办法：选择输入u(t)第三章相关分析法思路：问题：8第三章相关分析法

白噪声是一零均值，谱密度为非零常数，自相关函数为δ函数的各态历经的平稳随机过程，即：思路：维纳—霍甫方程

解决了抗干扰问题

引出了解积分方程难的问题

采用白噪声作为u(t)解决之，使得Ruy(

)与成比例

显然，当采用白噪声输入时，Ruy(

)与成比例，这样估计的问题就简化成求

Ruy(

)了。82第三章相关分析法白噪声是一零均值，谱密度为非零常数，②相关法：为了使得，可采用白噪声输入，即：

u(t)={白噪声}，不干扰系统正常运行。

第三章相关分析法相关法与古典法的差别之二：不干扰系统正常工作

①古典法：

直接用u(t)和y(t)通过求g(

)。若采用脉冲输入：，则：，将上述积分方程简化，但干扰了系统正常运行。将直接输入

(t)转化成直接输入白噪声，通过求Ruu，Ruy，间接地将输入Ruu变成()。从而不干扰系统正常运行。a83②相关法：第三章相关分析法相关法与古典法的差别之二：第三章相关分析法问题：由于在理论上要求积分时间T为无穷大，所以要求观测时间足够长。解决办法：还是选择输入u(t)！采用具有周期性的，近似于白噪声的伪随机信号作为输入信号，以解决积分时间长的问题。思路：维纳—霍甫方程

解决了抗干扰问题

引出了解积分方程难的问题

采用白噪声作为u(t)解决之，使得Ruy(

)与成比例

尚存在求Ruy(

)积分时间长问题

采用周期性，近似白噪声伪随机信号解决之，并仍保持Ruy(

)与成比例84第三章相关分析法问题：解决办法：还是选择输入u(t)！第三章相关分析法1）周期性：u(t+T)=u(t)（解决用一个T周期的观测数据就可得到Ruy(

)的问题。）∴周期性平稳过程的自相关函数为周期函数。85第三章相关分析法1）周期性：u(t+T)=u(t)∴周第三章相关分析法∴Ruy(

)的积分范围只要在一个T周期内。2）具有白噪声的性质：（保持与Ruy(

)的简单比例关系）Ruu(

)为周期函数，Ruu(

)=Ruu(

+T)86第三章相关分析法∴Ruy()的积分范围只要在一个T第三章相关分析法若选择

T大于系统脉冲响应的过渡过程时间Ts，即：当

>T时，g(

)

0对于第一个周期的激励而言，

（时间

在0～T

区间）有：87第三章相关分析法若选择T大于系统脉冲响应的对第三章相关分析法采用具有上述二性质的输入信号后，即可保持采用白噪声信号所具有的优越性，又可以解决

Ruy(

)的积分时间太长的问题，理论上只要在一个周期

T内积分就可以了。

具有上述二性质的输入信号到底能找到吗？若有，又将如何产生呢？思路：维纳—霍甫方程

解决了抗干扰问题

引出了解积分方程难的问题

采用白噪声作为u(t)解决之，使得Ruy(

)与成比例

尚存在求Ruy(

)积分时间长问题

采用周期性，近似白噪声伪随机信号解决之，并仍保持Ruy(

)与成比例

88第三章相关分析法采用具有上述二性质的输入信号后，即可§3

2伪随机二位式序列产生的方法及性质一、M序列产生的方法及性质：第三章相关分析法随机地掷一枚硬币的随机试验，结果：正面：+1

；反面：

1

反复试验

得到以+1，

1两元素组成的随机序列{u(k)}。当实验次数N相当大时，该序列{u(k)}具有以下两性质：①序列中+1与

1出现的次数几乎相等；（E[u]=0）

②随机序列的自相关函数Ruu(0)=max，离开原点时，Ruu(

)=0。（Ruu

）显然该序列接近于白噪声，最好它还应该是一个周期序列，在一个周期内具有上述白噪声性质。它在一个周期内观测时是一个随机信号；若观测时间很长时是一个周期信号。由于序列只有+1，

1两元素，称为伪随机二位式序列

PRBS（PseudoRandomBinarySignal）序列(它有规律性，故称伪随机，且可以人为产生和复制。)89§32伪随机二位式序列产生的方法及性质一、M序列产生n=4，k=2

（初态为：1111）时，

码数NP=6n=4，k=3

（初态为：1111）时，M序列

最大长度的伪随机二位式序列。第三章相关分析法由

n个双稳态触发器顺序组成

n级移位寄存器，将第

k级与第

n级状态“异或”后，反馈到第一级输入端。究竟k选哪一级呢？它将影响输出的性质。若k选择合适，将得到一个M序列。（初态不能为全零，否则输出总是零）

码数NP=151、M序列的产生90n=4，k=2（初态为：1111）时，码数N一个n级移位寄存器的输出序列的最大长度=？除各级全0的状态外，共有（2n

1）种不同的组合状态。第三章相关分析法若NP=2n

1，则该序列为最大长度序列或M序列。

当

n>12时，大约半数序列要用二级反馈产生，其他的则要用4级反馈来产生M序列。2、M序列的性质1）是一个确定的周期性序列，它的周期长度NP=2n

12）一个周期内。“0”状态比“1”状态少1个。（∵避免出现全“0”状态）“1”状态：“0”状态：91一个n级移位寄存器的输出序列的最大长度=？除各级全第三章相关分析法3）若将序列中相邻状态不变的那一部分长度称“游程”（或“段”），则在一个周期内的游程总数为m。∵不允许

n个全零状态，∴只有一个

n个码为全“1”∴游程总数如：n=4

1111

0

1

0

11

00

1

000m=8根据概率论可知，若游程总数为m，则：游程长度为i的为：游程长度为n的为：92第三章相关分析法3）若将序列中相邻状态不变的那一部分长度第三章相关分析法4）移位相加性：若将一个M序列与将其延迟了r个码以后的序列，按模2加法原则相加，所得的新序列还是M序列，不过延迟了q个码，r、q均为整数，且1

r，q

Np

1。，1

r，q

Np

1例：5）M序列具有近似离散的白噪声性质。下面将详细讨论M序列的自相关函数和功率密度谱。b93第三章相关分析法4）移位相加性：，1r，q第三章相关分析法3、M序列的相关函数和功率密度谱：1）相关函数：定义:“1”状态的逻辑电平为“

a”，“0”状态的逻辑电平为“+a”（为负逻辑关系，反之结论一样）。由Ruu()定义可知：Tp=NpΔ

94第三章相关分析法3、M序列的相关函数和功率密度谱：1）相第三章相关分析法(1)离散情况(

=

/

):其中：

=

/

,

设

是

的整数倍，此时

Ruu(

)取值于

=0，1，2，…，Np―1Ruu()写成离散形式为：∵显然模2乘法的结果与模2加法的结果在逻辑上是完全一样的，都为异或关系。即：同号码+a2（0）——u(k)与u(k+

)码的电平符号相同。异号码

a2（1）——u(k)与u(k+

)码的电平符号相异。∴{（同号码个数）－（异号码个数）}95第三章相关分析法(1)离散情况(=/):

=1，2，…，Np―1第三章相关分析法①当

=1，2，…，Np―1时：u(k)u(k+)在逻辑状态上相当于原序列{u(k)}与另一延迟序列{u(k+)}按摸2加法原则相加。根据M序列的移位相加性质可知，所得的结果在逻辑状态上仍是一个M序列。∴同号码个数

=新序列“0”状态个数=异号码个数

=新序列“1”状态个数=当新的序列为“0”状态时，说明u(k)与u(k+

)是同号当新的序列为“1”状态时，说明u(k)与u(k+

)是异号96=1，2，…，Np―1第三章相关分析法①当第三章相关分析法②当

=0时：(同一个M序列自乘)

=1，2，…，Np―1∴∴97第三章相关分析法②当=0时：(同一第三章相关分析法（2）连续情况（

不是

的整数倍）：（为平均面积值）为一个周期内曲线

u(t)•u(t+τ)

所围成的面积。即：一个周期内曲线

u(t)•u(t+τ)所围成的：正面积－负面积98第三章相关分析法（2）连续情况（不是的整数倍第三章相关分析法①当

=0时：(0<

<

)M序列每出现一次状态转换，积分将出现一个负面积(τ•a2)由M序列的性质3可知，在一个周期内：状态转换的次数＝M序列游程总数＝负面积＝正面积＝99第三章相关分析法①当=0时：M序列每出现一次第三章相关分析法负面积

＝正面积＝100第三章相关分析法负面积＝正面积＝100第三章相关分析法它是

τ的线性函数，因此可确定其两点：当τ=0时，Ruu(τ)=a2;当τ=

时，Ruu(τ)=

∵Ruu(τ)=Ruu(-τ)

为偶函数，∴在－

<

<

内，Ruu(τ)

为一个

波。

101第三章相关分析法它是τ的线性函数，因此可确定其两点：第三章相关分析法②当

•

<

<(

+1)•

，且

=1，2，…，Np―1时：可像

=0时的做法一样，在的基础上加上或减去一些面积来得到的平均面积，可以发现从中加上的面积正好等于减去的面积。∴综上所述，可得：102第三章相关分析法②当•<<(+1)•第三章相关分析法103第三章相关分析法103第三章相关分析法

M序列的Ruu(

)是一串周期性的

波，它与白噪声的δ函数还有差别，但只要缩小

，三角波的宽度变窄，Ruu(

)就接近理想脉冲。

当Np→

时：即∴此时Ruu(μ)为离散白噪声的自相关函数。104第三章相关分析法M序列的Ruu()是一串周第三章相关分析法令：Ruu(

)=Ruu1(

)+Ruu2(

)

，则：3105第三章相关分析法令：Ruu()=Ruu1(第三章相关分析法2）功率密度谱：

Ruu(τ)是一个以TP为周期的周期函数，而周期函数的Fourier变换是一个离散的频谱，其基波频率fo=1/Tp

Ruu(τ)的周期TP=NpΔ，且有界，满足Dirichlet条件，故可以展开复数形式的Fourier级数：基波角频率：Fourier复系数：106第三章相关分析法2）功率密度谱：Ruu(τ)是一个第三章相关分析法根据维纳

辛钦公式（Fourier变换对）：

M序列的功率密度谱Φuu(ω)与自相关函数Ruu(τ)是一个Fourier变换对，即因为Ruu(τ)的Fourier级数一致收敛，所以积分与求和可以交换顺序。∴107第三章相关分析法根据维纳辛钦公式（Fourier变换对第三章相关分析法由δ函数的性质可知：f(τ)=1

与F(ω)=2πδ(ω)是一个Fourier变换对，即δ(ω–rωo)说明Φuu(ω)

是一个离散的频谱，ω只能在rωo处取值，所以：

108第三章相关分析法由δ函数的性质可知：f(τ)=第三章相关分析法∵Ruu(τ)=Ruu(-τ)

为τ的偶函数∴109第三章相关分析法∵Ruu(τ)=Ruu(-τ)第三章相关分析法三项之和为零110第三章相关分析法三项之和为零110第三章相关分析法Φuu(ω)

在离散频率上取值，并考虑到：(1)当r=0

时，ω=rωo=0，可得：∴111第三章相关分析法Φuu(ω)在离散频率上取值，并考虑第三章相关分析法(2)当r≠0

时，

∴最后得：Φuu(ω)是一个离散的频谱，基波角频率：112第三章相关分析法(2)当r≠0时，∴最后得：第三章相关分析法

M序列的功率密度谱

Φuu(ω)是一个离散的线条谱，且有一个

(Sinx/x)2形的包络线，它的第一次取零的频率就是时钟脉冲的频率1/Δ(ω=2π/Δ)。显然，Φuu(ω)=Φuu(-ω)是偶函数，是关于纵坐标对称的。a113第三章相关分析法M序列的功率密度谱Φuu(ω)是第三章相关分析法当ω满足时,功率密度谱Φuu(ω)下降3dB。

∵,∴则

∴3dB114第三章相关分析法当ω满足3dB第三章相关分析法性质：①M序列的有效频带是：～～②M序列的激励功率与试验信号的幅值平方

a2成正比，与序列长度

Np成反比。1153dB第三章相关分析法性质：①M序列的有效频带是：1、逆重复M序列的产生将2NP个码的M序列{u(k)}与2NP个码的方波信号{m(k)}，按模2加法规则相加，即可得到逆重复M序列{l(k)}。第三章相关分析法{l(k)}={u(k)}

{m(k)}1）逆重复M序列{l(k)}的周期

=2Tp（即为原来M序列{u(k)}周期的两倍，Tp=NPΔ）为偶数。二、逆重复M序列产生的方法及性质：2、M序列的性质M序列{u(k)}的NP=2n

1为奇数，而方波{m(k)}的周期为偶数（2个码），所以2NP内的逆重复M序列{l(k)}不会重复。{u(k)}{m(k)}{l(k)}1161、逆重复M序列的产生将2NP个码的M序列{u(k)}2）逆重复M序列的前、后半个周期是逆重复的，即第三章相关分析法

l(t)=-

l(t+Tp)3）一个周期内，“0”状态与“1”状态的个数相等，各为

NP。

若定义“1”状态的逻辑电平为“

a”，“0”状态的逻辑电平为“+a”，则在周期2Tp内为零均值（更接近于白噪声）。4）逆重复M序列{l(k)}与M序列{u(k)}不相关，即：Rul(μ)=0∵l(k)=-

l(k+Np)u(k)=u(k+Np)∴1172）逆重复M序列的前、后半个周期是逆重复的，即第三章相第三章相关分析法5）逆重复M序列自相关函数Rll(τ)与M序列自相关函数Ruu(τ)的关系为：Rll

(

)=令：Rll

(

)=Rll1(

)+Rll2(

)118第三章相关分析法5）逆重复M序列自相关函数Rll(τ)与第三章相关分析法Rll(

)=Rll1(

)+Rll2(

)119第三章相关分析法Rll()=Rll1()+Rl第三章相关分析法Rll(

)=Rll1(

)+Rll2(

)Ruu(

)=Ruu1(

)+Ruu2(

)120第三章相关分析法Rll()=Rll1()+Rl第三章相关分析法Rll1(

)=Ruu1()0τ

(Np-1)Δ121第三章相关分析法Rll1()=Ruu1()第三章相关分析法逆重复M序列的Rll

2(

)的均值为零，M序列的Ruu

2(

)的均值非零；在0τ

（Tp-Δ）内：逆重复M序列的Rll

1(

)等于M序列的Ruu

1(

)，即：Rll

1(

)=

Ruu

1()∴逆重复M序列更接近于白噪声。

b122第三章相关分析法逆重复M序列的Rll2()的均值为零§3

3用M序列辨识线性系统的脉冲响应函数第三章相关分析法一、脉冲响应函数的辨识在辨识线性系统的脉冲响应函数的试验中，需要将有关算法转变成离散形式。当采样周期T0与M序列的时钟脉冲同步时，即T0=，序列的长度Np满足（Np-1）>Ts（系统的过渡过程时间），即：

当t>(Np

-1)

时，

g(t)

0

维纳—霍甫方程的离散形式为：其中：123§33用M序列辨识线性系统的脉冲响应函数第三章相关分第三章相关分析法在离散情况下M序列的Ruu(

)为：∴——

它是由Ruu2(

)产生的令124第三章相关分析法在离散情况下M序列的Ruu()为：∴第三章相关分析法∵脉冲响应函数g(k)有界，∴C为有界常数，且C>0∴

1、作图法将

Ruy(

)上移–C就可以得到：-C一般可以通过对Ruy(

)的稳态值的目测得到。从而得到：125第三章相关分析法∵脉冲响应函数g(k)有界，∴第三章相关分析法2、解析法通过精确计算公式得到：∵方程两边同求和，可得：126第三章相关分析法2、解析法通过精确计算公式得到：∵方程两第三章相关分析法二、估计量的统计特征1、是无偏的，即方程两边同取均值，可得：127第三章相关分析法二、估计量的统计特征1、第三章相关分析法此时，又回到理论的维纳—霍甫方程，所以：=0此时g(μ)已经是一个确定的量了！128第三章相关分析法此时，又回到理论的维纳—霍甫方程，所以：第三章相关分析法2、是一致估计量，即方程两边同取二阶原点矩，可得：129第三章相关分析法2、是一致估计量，即方程两边同取第三章相关分析法同证明无偏性一样，可得：可得的方差为:自然说明了估计量的有效性。是一致估计量。130第三章相关分析法同证明无偏性一样，可得：可得的第三章相关分析法三、提高估计精度的方法的估计精度取决于Ruy()的精度。1、提高采样速率提高Ruy()的精度：取采样周期To=/

(=1～4),用更多的y(t)数据计算Ruy()。2、采用多个周期的M序列输入r+1个周期的M序列，测得r个周期的y(t)计算Ruy()。通常取

r=1～44131第三章相关分析法三、提高估计精度的方法的估计精度取决于R第三章相关分析法四、计算的方法（采用多个周期）1、一次完成法（离线计算法）定义：上式写出向量─矩阵形式：132第三章相关分析法四、计算的方法（采用多个周期）1第三章相关分析法UNp

x

r

NpYr

Np

x

1∴133第三章相关分析法UNpxrNpYrNpx1∴134第三章相关分析法2）需要输入r+1个周期的u(k)：u(-Np+1)～u(r

Np-1)

。特点：1）一次离线求出（μ=0,1，…，Np-1)。3）精度要求较高时，Ruy()的计算精度要高，r的数目要大，所以数据存储量大。4）不是递推公式，无法在线辨识。134134第三章相关分析法2）需要输入r+1个周期的u(k)第三章相关分析法2、递推算法设已获得M对I/O数据，且M≥μ（∵μ=0,1，…，Np-1，∴M≥Np-1），即必须先观测至少一个周期。

Ruy(μ)的递推公式，即全部的M对I/O数据的Ruy(μ,M)可以用过去的(M-1)对I/O数据算得的Ruy(μ,M-1)和第M次观测的最新数据y(M)和u(M)递推地计算出。

135第三章相关分析法2、递推算法设已获得M对I/O数据，第三章相关分析法可得向量─