主要内容

LagOp

创建滞后算子多项式

描述

创建一个p度,-维滞后算子多项式一个l) =一个0+一个1l1+一个2l2+……+一个plp通过指定系数矩阵一个0、……一个p以及相应的滞后。l滞后(或退档)操作符使ljytyt- - - - - -j

LagOp对象的功能使您能够处理指定的多项式。例如,您可以通过一个多项式过滤时间序列数据,确定一个是否稳定,或者通过执行多项式代数(包括加、减、乘和除)来组合多个多项式。

为了将一个包含滞后算子多项式的动态模型拟合到数据上,需要创建合适的模型对象,然后将其拟合到数据上。关于单变量模型,请参见华宇电脑估计;有关多变量模型,请参见varm估计.为了进一步分析,您可以创建LagOp对象的估计系数。

创建

描述

例子

一个= LagOp (系数创建一个滞后算子多项式一个与系数系数,并设置系数财产。

例子

一个= LagOp (系数名称,值使用一个或多个名称-值对参数指定其他选项。例如,LagOp(系数,“滞后”,[0 4 8],“宽容”,1平台以及)将系数与滞后相关联04,8,并设置滞后包容公差为1平台以及

输入参数

全部展开

滞后算子多项式系数,指定为数值矢量,-借-的平方数值矩阵,或单元向量-借-平方数字矩阵。

价值 返回的多项式LagOp
数值长度向量p+ 1 一个p-次,一维滞后算子多项式,其中系数(j是滞后系数吗j– 1
-借-广场数字矩阵 0度,-D滞后算子多项式,其中系数一个0,滞后系数为0
长度p+ 1细胞向量-借-广场数字矩阵 一个p-D滞后算子多项式,其中系数(j滞后系数矩阵是多少j– 1

例子:LagOp (1:3)创建一个多项式一个l) = 1 + 2l1+ 3l2

名称-值对的观点

指定可选的逗号分隔的对名称,值论据。的名字参数名和价值为对应值。的名字必须出现在引号内。可以以任意顺序指定多个名称和值对参数名称1,值1,…,名称,值

例子:8“滞后”,[4],“宽容”,1平台以及将系数与滞后相关联48,并设置系数大小公差为1平台以及

与多项式系数相关的滞后,指定为逗号分隔对,由“滞后”和一个唯一的非负整数向量。滞后一定有努姆科夫元素,努姆科夫是以下内容之一:

  • 如果系数是数字或单元格向量,努姆科夫元素个数(系数).滞后系数滞后(j系数(j

  • 如果系数是一个矩阵,努姆科夫1

例子:‘滞后’,[0 4 8 12]

滞后包含公差,指定为逗号分隔对组成“宽容”和一个非负的数值标量。

c为系数元素j以最大的震级。如果c宽容LagOp去除系数j从多项式。因此,MATLAB®执行以下操作:

  • 从向量中移除相应的滞后滞后财产。

  • 替换数组中对应的系数系数财产零(

  • 如果一个被移除的滞后是多项式的最终滞后,MATLAB将多项式的程度降低到多项式中出现的第二大滞后的程度。例如,如果MATLAB从4次多项式中去除滞后3和4,因为它们的系数大小低于宽容,得到的多项式的次数为2。

例子:“宽容”,1平台以及

属性

全部展开

滞后运算符多项式系数,指定为数值标量的滞后索引单元阵列或-借-矩阵。系数{j是滞后系数吗j,其中整数j0.例如,A.Coefficients {0}存储的滞后系数0一个

滞后属性存储的非零系数的指标系数

此属性是只读的。

多项式的学位p,指定为非负数值标量。最大值(滞后),与非零系数相关的最大滞后。

数据类型:双重的

此属性是只读的。

多项式的维度,指定为正数值标量。你可以应用多项式一个只有一个-D时间序列变量。

数据类型:双重的

此属性是只读的。

与非零系数相关的多项式滞后,指定为非负整数的滞后索引向量。

工作的价值滞后,输入以下代码将其转换为标准MATLAB向量。

滞后= A.Lags;

数据类型:双重的

对象的功能

滤器 应用滞后算子多项式对时间序列进行滤波
isEqLagOp 确定两个LagOp对象是相同的数学多项式
isNonZero 求与的非零系数相关的滞后LagOp对象
不稳定 确定滞后算子多项式的稳定性
- 滞后算子多项式减法
mldivide 滞后算子多项式左除法
mrdivide 滞后算子多项式右除法
mtimes 滞后算子多项式乘法
+ 滞后算子多项式加法
反映 反映滞后零点附近的滞后算子多项式系数
托克拉里 将滞后算子多项式对象转换为单元阵列

例子

全部折叠

创建一个LagOp对象,该对象表示滞后算子多项式

一个 l 1 - 0 6 l + 0 0 8 l 2

A=LagOp([1-0.60.08])
一个=一维滞后算子多项式 : ----------------------------- 系数:-0.6 - 0.08[1]滞后:[0 1 2)学位:2维:1

显示滞后系数0。

a0 = A.Coefficients {0}
a0 = 1

给第三个滞后赋一个非零系数:

A.系数{3}=0.5
一个=一维滞后算子多项式 : ----------------------------- 系数:[1 -0.6 0.08 0.5]滞后:[0 1 2 3]学位:3维:1

多项式次数增加到3.

创建一个LagOp对象,该对象表示滞后算子多项式

一个 l 1 + 0 2 5 l 4 + 0 1 l 8 + 0 0 5 l 1 2

nonzeroCoeffs = [1 0.25 0.1 0.05];时滞= [0 4 8 12];= LagOp (nonzeroCoeffs“滞后”,滞后)
A=1-D滞后算子多项式:------------------------------系数:[1 0.25 0.1 0.05]滞后:[0 4 8 12]度:12维:1

从滞后算子多项式中提取系数,并显示从滞后0到滞后12的所有系数。

allCoeffs = toCellArray(一个);%提取系数allCoeffs = cell2mat (allCoeffs ');allLags = 0: A.Degree;%准备显示延迟表(allCoeffs,“RowNames”“滞后”+字符串(allLags))
ans=13×1表allCoeffs _________滞后0 1滞后10滞后2 0滞后3 0滞后4 0.25滞后5 0滞后6 0滞后7 0滞后8 0.1滞后9 0滞后10 0滞后11 0滞后12 0.05

创建一个LagOp对象,该对象表示滞后算子多项式

一个 l 0 5 0 0 0 1 0 0 0 - 0 5 + 1 0 2 5 0 1 - 0 5 1 - 0 5 0 1 5 - 0 2 1 l 4

Phi0 = [0.5 0 0;...0 1 0;...0 0 -0.5];Phi4 = [1 0.25 0.1;...-0.5 1 -0.5;...0.15 - -0.2 1);= {Phi0 Phi4};时滞= [0 4];一个= LagOp(φ,“滞后”,滞后)
A=3-D滞后算子多项式:------------------------------系数:[具有2个非零系数的滞后索引单元数组]滞后:[0 4]度:4维:3

创建多元滞后算子多项式

创建2-D 3次滞后算子多项式

一个 l 1 0 0 1 + 0 5 0 2 5 - 0 1 0 4 l + 0 0 5 0 0 2 5 - 0 0 1 0 0 4 l 3.

m = 2;A0 =眼(m);A1 = [0.5 0.25;-0.1 - 0.4);A3 = 0.1 * A1;Coeffs = {A0 A1 A3};时滞= [0 1 3];一个= LagOp(多项式系数,“滞后”,滞后)
A =二维滞后算子多项式:-----------------------------系数:[带有3个非零系数的滞后索引单元阵列]滞后:[0 1 3]度:3维:2

确定多项式稳定性

一个时滞算子多项式是稳定的,如果它的特征多项式的所有特征值的大小小于1。

确定多项式是否稳定,并返回其特征多项式的特征值。

(tf,测评)=趋于稳定(A)
tf =逻辑1
测评=6×1复杂-0.5820 - 0.1330i 0.0821 + 0.2824i 0.0821 - 0.2824i 0.0499 + 0.2655i

逆多项式

计算的倒数 一个 l 通过将2×2的单位矩阵除以 一个 l

Ainv = mrdivide(眼睛(A.Dimension))
Ainv = 2- d滞后算子多项式:-----------------------------系数:[带有10个非零系数的滞后索引单元阵列

Ainv是一个LagOp表示函数的倒数的对象 一个 l ,一个9次滞后算子多项式。理论上,滞后算子多项式的逆是无限次的,但是mrdivide系数幅值公差截断多项式。

一个 l 和它的逆矩阵。

checkinv=Ainv*A
checkinv=2-D滞后算子多项式:-------------------------------系数:[具有4个非零系数的滞后索引单元数组]滞后:[0 10 11 12]度:12维:2

因为逆计算返回截断的理论逆,即乘积checkinv包含代表剩余部分的滞后。你可以减少mrdivide系数大小公差,以获得一个更精确的反多项式。

过滤时间序列

生成二维时间序列 y t 一个 l e t 一个 0 e t + 一个 1 e t - 1 + 一个 2 e t - 2 + 一个 3. e t - 3. 通过过滤二维序列的100随机标准高斯偏差 e t 通过多项式。

T = 100;e = randn (T, m);y =过滤器(A, e);情节(A。度+ 1):T,y) title(“过滤系列”

图中包含一个坐标轴。标题为“过滤系列”的轴包含两个类型为line的对象。

y一个97乘2的矩阵代表什么 y t yp观察结果少于e因为滤器需要第一个p的观察e生成时初始化动态序列的步骤: y (1:4)

介绍了R2010a