主要内容

使用ARMAX估算模型

此示例显示了如何使用迭代估计方法估算具有用于三输入和单输出(MISO)系统的ARMAX结构的线性多项式模型armax.。有关工具箱中的所有可用估计命令的摘要,请参阅模型估计命令

加载示例数据集Z8有三个输入和一个输出,测量1-Second间隔并包含500个数据样本。

加载Iddata8.

采用armax.两个构造Idpoly.模型对象,估计参数:

一种 问: y T. = σ. 一世 = 1 N B. 一世 问: 一世 T. - N K. 一世 + C 问: E. T.

通常,您尝试不同的模型订单并比较结果,最终选择最能描述系统动态的最简单模型。以下命令指定估计数据集,Z8,以及订单一种B., 和C多项式NA.NB., 和NC., 分别。NK.[0 0 0]指定所有三个输入通道没有输入延迟。

选择= ARMAXOPTIONS;Opt.focus =.'模拟';opt.SearchOptions.maxItations = 50;opt.SearchOptions.Tolerance = 1E-5;na = 4;nb = [3 2 3];nc = 4;nk = [0 0 0];m_armax = ARMAX(Z8,[NA NB NC NK],OPT);

重点宽容, 和maxiter.是配置估计目标函数和搜索算法的属性的估计选项。这重点选项指定是否针对模拟或预测应用程序进行了优化。这宽容maxiter.搜索选项指定何时停止估计。有关这些属性的更多信息,请参阅armaxoptions.参考页面。

armax.是一个版本poly具有ARMAX模型结构的简化语法。这armax.方法都构造了Idpoly.模型对象并估计其参数。

查看结果模型对象的信息。

m_armax.
M_ARMAX =离散时间ARMAX模型:a(z)y(t)= b(z)u(z)+ c(z)e(z)a(z)= 1  -  1.284 z ^ -1 + 0.3048 z ^-2 + 0.2648 z^-3 - 0.05708 z^-4 B1(z) = -0.07547 + 1.087 z^-1 + 0.7166 z^-2 B2(z) = 1.019 + 0.1142 z^-1 B3(z) = -0.06739 + 0.06828 z^-1 + 0.5509 z^-2 C(z) = 1 - 0.06096 z^-1 - 0.1296 z^-2 + 0.02489 z^-3 - 0.04699 z^-4 Sample time: 1 seconds Parameterization: Polynomial orders: na=4 nb=[3 2 3] nc=4 nk=[0 0 0] Number of free coefficients: 16 Use "polydata", "getpvec", "getcov" for parameters and their uncertainties. Status: Estimated using ARMAX on time domain data "z8". Fit to estimation data: 80.86% (simulation focus) FPE: 2.888, MSE: 0.9868

m_armax.是一个Idpoly.模型对象。系数表示该多项式模型的估计参数。您可以使用礼物(m_armax)显示有关该模型的其他信息,包括参数不确定性。

查看此模型的所有属性值。

得到(m_armax)
答:[1 -1.2836 0.3048 0.2648 -0.0571] B:{[-0.0755 1.0870 0.1166] [1x3双]} C:[1-0.0610 -0.1296 0.0249 -0.0470] D:1 F:{[1][1] [1] [1]} IntegrateNoise:0变量:'z ^ -1'Iodelay:[0 0]结构:[1x1 pmodel.polynomial] noisevariance:2.7984 inputdelay:[3x1 double] outputdelay:0 ts:1 yountunit:“秒”InputName:{3x1 Cell} InputUnit:{3x1 Cell} InputGroup:[1x1 struct] OutputName:{'Y1'} OutputUnit:{''} OutputGroup:[1x1 struct]注意:[0x1字符串] UserData:[]名称:''samplinggrid:[1x1 struct]报告:[1x1 idesults.polyest]

报告Model属性包含有关估算结果的详细信息。查看内部的属性和值报告,使用点表示法。例如:

m_armax.report.
ANS =状态:'使用ARMAX进行仿真焦点的ARMAX'方法:'ARMAX'INITIONCONDITION:'ZERO'FIT:[1x1 struct]参数:[1x1 struct]选项已使用:[1x1 IDOPOPTIONS.polyest] RANDSTATE:[1x1 struct] datsumed:[1x1结构]终止:[1x1 struct]

此操作显示估计报告的内容,如模型质量措施(合身),搜索终止标准(终止)和估计数据的记录(dataused.和选项(选择)。

相关例子

更多关于