主要内容

渣油

计算和测试残差

描述

例子

残油(数据sys计算已识别模型的提前一步预测误差(残差),sys,并根据中的数据绘制剩余输入动态图如下所示数据

  • 对于时域数据,渣油绘制残差的自相关和残差与输入信号的相互关系。相关性是在滞后-25到25之间产生的。若要指定不同的最大滞后值,请使用residOptions.99%可信区域标记了统计上不显著的相关性,显示为x轴周围的阴影区域。

  • 对于频域数据,渣油绘制从输入信号到残差的频率响应的预兆图。99%置信区域标记了统计上不显著的响应,显示为x轴周围的区域。

要更改显示选项,右键单击图以访问上下文菜单。有关菜单的详细信息,请参见提示

例子

残油(数据sysLinespec设置线条样式、标记符号和颜色。

例子

残油(数据、sys1 sysn……)计算和绘制多个识别模型的残差sys1、……sysn

例子

残油(数据, Linespec1 sys1…,sysn Linespecn)为每个系统设置线条样式、标记符号和颜色。

例子

残油(___选项指定额外的剩余计算选项。使用选项使用任何前面的语法。

例子

残油(___类型指定绘图类型。使用类型使用任何前面的语法。

例子

ER) =渣油(数据sys返回计算的残差,E,和残差相关,R.不生成任何情节。

例子

全部折叠

时域数据加载。

负载iddata1数据= z1;

估计一个ARX模型。

Sys = arx(data,[1 1 0]);

绘制残差的自相关以及残差与输入之间的相互关系。

渣油(数据、系统)

图中包含2个轴。标题为AutoCorr的轴1包含两个类型为line的对象。该对象表示sys。标题为XCorr (u1)的轴2包含2个line类型的对象。该对象表示sys。

关联一直计算到默认的最大延迟25。99%可信区域标记了统计上不显著的相关性,显示为x轴周围的阴影区域。

转换数据到频域。

data2 = fft(数据);

计算已识别模型的残差,sys,以及频域数据。用红叉标出剩余响应。

渣油(data2 sys,“处方”

图中包含2个轴。标题为From: u1 To: e@y1的轴1包含一个line类型的对象。该对象表示sys。Axes 2包含一个类型为line的对象。该对象表示sys。

对于频域数据,渣油绘制波德图显示从输入到残差的频率响应。

时域数据加载。

负载iddata1

估计一个ARX模型。

syso2 = arx(z1,[1 1 0]);

估计传递函数模型。

sys2 =特遣部队(z1, 2);

绘制残差的相关性。

渣油(z1、sys1“b”sys2,“r”

图中包含2个轴。标题为AutoCorr的轴1包含4个类型为line的对象。这些对象表示sys1, sys2。标题为XCorr (u1)的轴2包含4个line类型的对象。这些对象表示sys1, sys2。

的残差之间的互相关sys2所有滞后的输入都在99%的置信区间。

时域数据加载。

负载iddata1

估计一个ARX模型。

Sys = arx(z1,[1 1 0]);

指定残差相关计算的最大延迟。

选择= residOptions (“MaxLag”35);

画出从输入到残差的脉冲响应。

渣油(z1、系统选择、“红外”

图中包含一个坐标轴。标题为From: u1 To: e@y1的轴包含两个line类型的对象。该对象表示sys。

时域数据加载。

负载iddata7

数据是一个双输入单输出的数据集。

估计一个ARX模型。

sys =特遣部队(z7 2);

计算残差及其与输入的自相关和相互关系。

[E R] =渣油(z7 sys);

R是一个26 × 3 × 3的相关矩阵。例如,

  • R (:, - 1, - 1)为滞后25前残差的自相关。

  • R (:, 1, 2)为残差与第一个输入的相互关系,直到滞后25。

E是一个iddata对象,其残差作为输出数据,验证数据(z7)的输入作为输入数据。您可以使用E识别误差模型并分析误差动力学。

画出错误。

情节(E)

图中包含3个轴。标题为e@y1的轴1包含一个类型为line的对象。标题为u1的坐标轴2包含一个类型为line的对象。这个对象代表E. Axes 3,标题u2包含一个类型为line的对象。这个对象表示E。

估计输入和残差之间的脉冲响应。用3个标准差置信区间画出来。

I =冲动(E);showConfidence (impulseplot(我,20),3)

图中包含2个轴。标题为From: u1的坐标轴1包含两个line类型的对象。这个对象代表I. Axes 2,标题为From: u2包含2个类型为line的对象。这个对象表示I。

输入参数

全部折叠

验证输入输出数据,指定为iddata对象。数据可以有多个输入输出通道。当sys是线性的,数据为时域或频域。当sys是非线性的,数据时域。

用于计算残差的系统,指定为已确定的线性或非线性模型。

例子:idpoly

作为字符向量指定的线样式、标记符号和颜色。有关更多信息,请参见情节.当类型被指定为“相关系数”,只使用线条风格。

例子:“Linespec”、“kx的

剩余分析选项,指定为residOptions选项设置。

Plot类型,指定为以下值之一:

  • “相关系数”-绘制残差的自相关图,e,残差与输入信号的互相关,u.相关性是在滞后-25到25之间产生的。使用residOptions指定一个不同的最大滞后值。99%可信区域表示统计上不显著的相关性,也显示为围绕x轴的阴影区域。置信区间的计算是在假设条件下完成的e白的,独立的u

    “相关系数”是时域数据的默认值。此绘图类型不适用于频域数据。

  • “红外”-绘制系统从输入到残差的滞后25的脉冲响应。的冲动命令首先估计脉冲响应模型e作为输出数据u作为输入。然后冲动计算估计模型的脉冲响应。99%置信区域标记了统计上不重要的响应显示为阴影区域。低震级表明模型可靠。

    此绘图类型不适用于频域数据。

  • “fr”-从输入到残差的频率响应(基于高阶FIR模型)显示为波德图。99%置信区域标记了统计上不重要的响应显示为阴影区域。在感兴趣的频率范围内的低幅度表示一个可靠的模型。

    “fr”是频域数据的默认值。

输出参数

全部折叠

模型残差,作为iddata对象。残差存储在E.OutputData,输入存储在E.InputData.使用E建立描述从输入到残差的动态的模型。如果sys是一个可靠的识别模型。

残差的相关性,返回为:

  • 双精度矩阵-用于时域数据

    R矩阵的大小是多少+ 1 - (纽约+ν)——- (纽约+ν).在那里,最大延迟是否在选项纽约是输出数,和ν为输入的个数。的默认值是25。

    在每一个滞后kk = 0: M),R (k, i, j)为期望值,< Z (t, i)还是Z (t + k - 1 j) >.在这里,Z = [E.OutputData E.InputData]

    例如,对于双输出单输入模型,Z = (e1, e2, u1).在那里,e1为第一次输出的残差,e2是剩余的第二个输出,和u1是输入。R是一个26——- - - - - -3.——- - - - - -3.矩阵,地点:

    • R(5、1、2)= < e1 (t) .e2 (t + 4) >相互关系是否滞后于4e1e2

    • R(5、1,3)= < e1 (t) .u1 (t + 4) >相互关系是否滞后于4e1u1

    • R(5、1、1)R (2, 2)R(5、3、3)自相关性是否处于滞后状态e1e2,u1,分别。

  • []—频域数据

提示

  • 右键单击情节打开上下文菜单,您可以访问以下选项:

    • 系统—选择系统查看残差相关或响应图。默认情况下,绘制所有系统。

    • 展示自信的地区-查看99%置信区域标记统计上不显著的相关性。仅适用于相关图。

    • 数据实验—仅支持多实验数据。切换不同实验的数据。

    • 特征—查看数据特征。不适用于相关图。

      • 峰值响应-查看数据的峰值响应。

      • 置信区域-查看99%置信区域标记统计上不重要的响应。

    • 显示-仅适用于频率响应图。

      • -查看频率响应的幅度。

      • 阶段-查看频率响应相位。

    • I / O分组—用于包含多个输入或输出通道的数据集。在图上选择输入和输出通道分组。不适用于相关图。

      • 没有一个-在各自独立的坐标轴上绘制输入输出通道。

      • 所有-将所有输入通道和所有输出通道组合在一起。

    • I / O选择器—用于包含多个输入或输出通道的数据集。选择要绘制的输入和输出通道的子集。默认情况下,绘制所有输出通道。

    • 网格-在情节中添加网格。

    • 正常化-将图中所有数据的y刻度归一化。不适用于频响数据。

    • 全视图-返回全视野。默认情况下,plot被缩放到全视图。

    • 初始条件—指定初始条件处理。

      指定为以下其中之一:

      • 估计—将初始条件作为估计参数。

      • —设置所有初始条件为零。

      • 吸收延迟和估计-将非零延迟吸收到模型系数中,并将初始条件作为估计参数。此选项仅适用于离散时间模型。

    • 属性-打开属性编辑器对话框自定义情节属性。

参考文献

[1] Ljung, L。系统识别:用户的理论.上鞍河,新泽西州:Prentice-Hall PTR, 1999,第16.6节。

之前介绍过的R2006a