主要内容

佳能

规范状态空间实现

描述

例子

坐标系=佳能(sys类型)转换线性模型sys变成一个规范的状态空间模型坐标系类型指定是否坐标系为情态或伴语形式。

有关可控和可观察规范形式的信息,请参见规范的状态空间实现

例子

坐标系=佳能(sys“模态”,condt)指定上界condt关于块对角化变换的条件数。使用condt如果有接近的特征值坐标系

例子

坐标系T) =佳能(___)还返回状态坐标变换矩阵T它将状态空间模型的状态联系起来sys到美国坐标系

例子

全部折叠

aircraftPitchSSModel.mat包含飞机的状态空间矩阵,其中输入为电梯偏转角 δ 输出是飞机的俯仰角 θ

α ˙ ˙ θ ˙ - 0 3. 1 3. 5 6 7 0 - 0 0 1 3. 9 - 0 4 2 6 0 0 5 6 7 0 α θ + 0 2 3. 2 0 0 2 0 3. 0 δ y 0 0 1 α θ + 0 δ

将模型数据加载到工作区并创建状态空间模型sys

负载(“aircraftPitchSSModel.mat”);sys = ss(A,B,C,D)
sys = A = x1 x2 x3 x1 -0.313 56.7 0 x2 -0.0139 -0.426 0 x3 0 56.7 0 B = u1 x1 0.232 x2 0.0203 x3 0 C = x1 x2 x3 y1 0 01 D = u1 y1 0连续时间状态空间模型。

转换生成的状态空间模型sys与标准形式相伴。

Csys = canon(sys,“同伴”)
csys = A = x1 x2 x3 x1 0 0 -1.809e-16 x2 1 0 -0.9215 x3 0 1 -0.739 B = u1 x1 1 x2 0 x3 0 C = x1 x2 x3 y1 0 1.151 -0.6732 D = u1 y1 0连续时间状态空间模型。

坐标系的伴生规范形式是sys

pendulumCartSSModel.mat包含推车上倒立摆的状态空间模型,其中输出为推车位移x还有摆角 θ .控制输入u是推车上的水平力。

x ˙ x ¨ θ ˙ θ ¨ 0 1 0 0 0 - 0 1 3. 0 0 0 0 1 0 - 0 5 3. 0 0 x x ˙ θ θ ˙ + 0 2 0 5 u y 1 0 0 0 0 0 1 0 x x ˙ θ θ ˙ + 0 0 u

首先,加载状态空间模型sys到工作区。

负载(“pendulumCartSSModel.mat”“sys”);

转换sys模态标准形式,并提取变换矩阵。

[csys,T] = canon(sys,“模态”)
csys = A = x1 x2 x3 x4 x1 00 00 x2 0 -0.05 00 x3 00 -5.503 0 x4 00 0 5.453 B = u1 x1 1.875 x2 6.298 x3 12.8 x4 12.05 C = x1 x2 x3 x4 y1 16 -4.763 -0.003696 0.003652 y2 0 0.003969 -0.03663 0.03685 D = u1 y1 0 y2 0连续时间状态空间模型。
T =4×40.0625 1.2500 -0.0000 -0.1250 0 4.1986 0.0210 -0.4199 0 0.2285 -13.5873 2.4693 0 -0.2251 13.6287 2.4995

坐标系的模态规范形式是sys,而T的状态向量之间的变换sys而且坐标系

对于这个例子,考虑以下具有双极和近极集群的系统:

年代 y 年代 年代 ) 1 0 0 年代 - 1 ) 年代 + 1 ) 年代 年代 + 1 0 ) 年代 + 1 0 0 0 0 1 ) 年代 - 1 + ) ) 2 年代 - 1 - ) ) 2

创建一个zpk模型,并使用字符串将其转换为模态规范形式“模态”

sys = zpk([1],[0 -10 -10.0001 - 1 + 1我1-1i 1-1i 1 + 1], 100);Csys1 = canon(sys,“模态”);csys1。一个
ans =7×70000 000 1.0000 1.0000 0000 0 -1.0000 1.0000 3.2459 0000 000 1.0000 1.0000 0000 000 -1.0000 1.0000 0000 0000 -10.0000 4.0571 0000 0000 0 -10.0001
csys1。B
ans =7×10.1600 -0.0052 0.0427 -0.0975 0.5319 0 4.0095

sys有一对杆子在年代-10而且年代-10.0001,和两个复极的多重性2年代1 +我而且年代.因此,模态形式csys1状态空间模型的大小为2的两个极点是否接近年代-10,以及一个大小为4的复特征值块。

现在,把两根杆子分开年代-10通过增加块对角化变换的条件数。使用值为1 e10汽油对于这个例子。

Csys2 = canon(sys,“模态”1 e10);csys2。一个
ans =7×70000 0000 1.0000 1.0000 0000 0 -1.0000 1.0000 3.2459 0000 000 1.0000 1.0000 0000 000 -1.0000 1.0000 0000 0000 -10.0000 0000 0000 0 -10.0001
格式shortEcsys2。B
ans =7×10.000 -0.0000 0.0000 -0.0000 0.0000 1.6267 1.6267

一个矩阵的csys2包括相邻极点的独立对角线元素年代-10.中的一些非常大的值会导致条件数的增加B矩阵。

该文件icEngine.mat包含一个以0.04秒的采样率收集的1500个输入输出样本的数据集。输入u (t)是电压(V)控制旁通空气阀(BPAV),和输出y (t)为发动机转速(RPM/100)。

使用的数据icEngine.mat创建具有可识别参数的状态空间模型。

负载icEngine.matZ = iddata(y,u,0.04);Sys = n4sid(z,4,“InputDelay”2);

转换标识的状态空间模型sys与标准形式相伴。

Csys = canon(sys,“同伴”);

通过运行对模型参数的零迭代更新来获得结果表单的协方差。

opt = ssestOptions;opt.SearchOptions.MaxIterations = 0;Csys = ssest(z, Csys,opt);

比较的频率响应置信限sys坐标系

H = bodeploy (sys,csys,“r”。);showConfidence (h)

图中包含2个轴对象。从:u1到:y1的轴对象1包含2个line类型的对象。这些对象表示sys, csys。坐标轴对象2包含2个line类型的对象。这些对象表示sys, csys。

频率响应置信限相同。

输入参数

全部折叠

动态系统,指定为SISO或MIMO动态系统模型。您可以使用的动态系统包括:

  • 连续时间或离散时间数值LTI模型,例如特遣部队(控制系统工具箱)zpk(控制系统工具箱)党卫军(控制系统工具箱),或pid(控制系统工具箱)模型。

  • 广义的或不确定的LTI模型,如一族(控制系统工具箱)号航空母舰(鲁棒控制工具箱)模型。(使用不确定模型需要鲁棒控制工具箱™软件。)

    由此产生的规范状态空间模型假设

    • 可调控制设计块的可调组件的当前值。

    • 不确定控制设计块的标称模型值。

  • 确定的LTI模型,例如idtf中的难点idprocidpoly,idgrey模型。

不能使用频响数据模型,例如的朋友(控制系统工具箱)模型。

转换类型,指定为任意一种“模态”“同伴”.如果类型则未指定佳能默认情况下将指定的动态系统模型转换为模态规范形式。

伴生标准形式与可观测标准形式相同。有关可控和可观察规范形式的信息,请参见规范的状态空间实现

  • 模态形式

    在模态形式中,一个是一个块对角矩阵。对于实特征值,块大小通常是1乘1,对于复特征值,块大小通常是2乘2。但是,如果存在重复的特征值或附近特征值的集群,块大小可以更大。

    例如,对于一个具有特征值的系统 λ 1 σ ± j ω λ 2 ) ,模态一个矩阵是这样的

    λ 1 0 0 0 0 σ ω 0 0 ω σ 0 0 0 0 λ 2

  • 同伴的形式

    在相应的实现中,系统的特征多项式显式地出现在最右边的列中一个矩阵。对于具有特征多项式的系统

    P 年代 ) 年代 n + α 1 年代 n 1 + + α n 1 年代 + α n

    相应的同伴一个矩阵是

    一个 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 α n α n 1 α n 2 α n 3. α 1

    伴随变换要求系统从第一个输入开始是可控的。转换到伴生形式是基于可控性矩阵,几乎总是数值奇异的中等阶。因此,尽可能避免使用它。

    伴生标准形式与可观测标准形式相同。有关可观察和可控规范形式的更多信息,请参见规范的状态空间实现

块对角化变换的条件数的上界,指定为正标量。此参数仅当类型设置为“模态”

增加condt中特征值簇的大小一个矩阵的坐标系.设置condt =无穷对角化矩阵一个

输出参数

全部折叠

动态模型的规范状态空间形式,返回为党卫军(控制系统工具箱)模型对象。坐标系的状态空间实现sys以指定的标准形式类型

变换矩阵,返回为n × n矩阵,其中n是状态数。T状态向量之间的变换x状态空间模型sys状态向量xc坐标系

xcTx

此参数仅当sys是一个党卫军(控制系统工具箱)模型对象。

限制

  • 不能使用频率响应数据模型将其转换为规范状态空间形式。

  • 伴生形式对于大多数状态空间计算条件不佳,也就是说,到伴生形式的转换是基于可控性矩阵,对于中阶几乎总是数值奇异的。因此,尽可能避免使用它。

算法

佳能命令使用bdschur(控制系统工具箱)转换命令sys转换成模态形式并计算变换T.如果sys不是状态空间模型,佳能首先使用将其转换为状态空间党卫军

该方法采用基于可控性矩阵的状态相似度变换,简化为伴随形式[1]

参考文献

Kailath, T。线性系统, Prentice-Hall, 1980。

另请参阅

(控制系统工具箱)|(控制系统工具箱)|(控制系统工具箱)|(控制系统工具箱)|(控制系统工具箱)|(控制系统工具箱)|(控制系统工具箱)|(控制系统工具箱)|(鲁棒控制工具箱)|||||

R2006a之前介绍