主要内容

渣滓

计算和测试残差

描述

例子

渣油(数据sys计算识别的模型的一步预测误差(残差),sys,并根据数据中的数据,将剩余输入动态图绘制为以下内容之一数据

  • 对于时域数据,渣滓绘制残差的自相关和与输入信号的残差的互相关。为LAG -25至25生成相关性。要指定不同的最大滞后值,请使用Resentions..标记统计上微不足道的相关性的99%置信区显示为X轴周围的阴影区域。

  • 对于频域数据,渣滓绘制从输入信号到残差的频率响应的波德图。标记统计上不显著响应的99%置信区域显示为一个围绕x轴的区域。

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

例子

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

例子

渣油(数据、sys1 sysn……)计算并绘制多个识别模型的残余SYS1.,......,SYSN.

例子

渣油(数据,sys1,linepec1,...,sysn,linepecn)为每个系统设置线型,标记符号和颜色。

例子

渣油(___选项指定其他剩余计算选项。用选项使用以前的任何语法。

例子

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

例子

[E.R.] =渣油(数据sys返回计算的残差,E.,剩余相关,R..没有情节产生。

例子

全部收缩

加载时间域数据。

加载Iddata1.data = z1;

估计ARX模型。

系统= arx(数据,[1 1 0]);

绘制残差和残差与输入之间的交叉相关性。

渣油(数据、系统)

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

相关性被计算直到默认的最大延迟,25。标记统计上微不足道的相关性的99%置信区显示为X轴周围的阴影区域。

将数据转换为频域。

data2 = fft(数据);

计算已识别模型的残差,sys和频域数据。用红十字绘制残余响应。

resid(Data2,Sys,'rx'

图中包含2个轴。具有标题的轴1来自:U1至:E @ Y1包含类型线的对象。该对象表示sys。轴2包含类型线的对象。此对象表示SYS。

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

加载时间域数据。

加载Iddata1.

估计ARX模型。

drawtext (z1,[1]) = 0, '买进'),colorred;

估计传递函数模型。

sys2 = tfest(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]);

指定剩余相关计算的最大滞后。

opt = Resentions(“MaxLag”,35);

从输入到残差的脉冲响应绘制。

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

图包含轴。具有标题的轴:U1至:E @ Y1包含2个类型的对象。该对象表示sys。

加载时间域数据。

加载Iddata7.

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

估计ARX模型。

sys = tfest(z7,2);

计算残差及其自相关和与输入的交叉相关性。

[e,r] =渣油(z7,sys);

R.是一个26×3-3×3的相关性。例如,

  • R (:, - 1, - 1)是残差的自相关,直到滞后25。

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

E.是IDData对象,具有残差作为输出数据以及验证数据(Z7)的输入作为输入数据。您可以使用E.要识别错误模型并分析错误动态。

绘制错误。

情节(e)

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

估计输入和残差之间的脉冲响应。用3个标准偏差置信区绘制它们。

我=梗死(e);ShowConfidence(Impulseplot(I,20),3)

图中包含2个轴。具有标题的轴1来自:U1包含2个类型的型号。此对象表示I.带有标题的轴2:U2包含2个类型线的2个对象。这个对象代表I.

输入参数

全部收缩

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

用于计算剩余量的系统,指定为一个确定的线性或非线性模型。

例子:Idpoly.

线条样式,标记符号和颜色,指定为字符向量。有关更多信息,请参阅阴谋.什么时候类型被指定为“相关系数”,仅使用线条样式。

例子:'linespec','kx'

残差分析选项,指定为Resentions.选项集。

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

  • “相关系数”-绘制残差的自相关图,E.,以及残差与输入信号的互相关,.相关性是在滞后-25到25时产生的。用Resentions.指定不同的最大滞后值。标记统计上微不足道的99%置信区也显示为X轴周围的阴影区域。假设置信区域的计算E.成为独立的白人

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

  • “红外”- 将系统的脉冲响应绘制到从输入到残差的系统的LAG 25。这冲动命令首先估计脉冲响应模型E.作为输出数据和作为输入。然后冲动计算估计模型的脉冲响应。标记统计上微不足道的响应的99%置信区显示为阴影区域。低幅度表示可靠的模型。

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

  • 'fr'- 从输入到残差的频率响应(基于高阶FIR模型)被示为Bode图。标记统计上微不足道的响应的99%置信区显示为阴影区域。频率范围内的低幅度表示可靠的模型。

    'fr'默认为频域数据。

输出参数

全部收缩

模型残差,返回为iddata.对象。残留物存储在e.outputdata.,并且输入存储在E.Inputdata..用E.构建描述从输入到残差的动态的模型。动态可以忽略不计sys是一个可靠的识别模型。

残差的相关性,作为以下内容之一返回:

  • 双打矩阵 - 时间域数据

    R.矩阵的大小m+ 1-by-(纽约+ν)-经过-(纽约+ν)。在那里,m是指定的最大滞后选项纽约是输出的数量,和ν是输入的数量。默认值m是25。

    每个滞后K.k = 0:m),r(k,i,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) >是滞后4之间的互相关E1.E2.

    • R(5、1,3)= < e1 (t) .u1 (t + 4) >是滞后4之间的互相关E1.u1

    • r(5,1,1)r(5,2,2)r(5,3,3)自相关是滞后的吗E1.E2., 和u1, 分别。

  • []- 用于频域数据

尖端

  • 右键单击绘图打开上下文菜单,您可以在其中访问以下选项:

    • 系统- 选择系统以查看剩余相关或响应图。默认情况下,绘制所有系统。

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

    • 数据实验—仅用于多实验数据。切换来自不同实验的数据。

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

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

      • 信心地区- 查看标志着统计上微不足道的99%的置信区。

    • 表演- 仅适用于频率响应图。

      • - 查看频率响应的大小。

      • 阶段- 查看频率响应的阶段。

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

      • 没有任何- 在自己的单独轴上绘制输入 - 输出通道。

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

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

    • 网格- 将网格添加到图中。

    • 正常化- 标准化图中所有数据的y级。不适用于频率响应数据。

    • 全视图- 返回完整视图。默认情况下,绘图缩放到完整视图。

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

      指定为下列之一:

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

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

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

    • 属性- 打开“属性编辑器”对话框以自定义绘图属性。

参考

[1] Ljung,L.系统识别:用户理论.上鞍河,新罕布什尔州:Prentice-Hall Ptr,1999,第16.6节。

之前介绍过的R2006a