主要内容

robstab

不确定系统的鲁棒稳定性

描述

例子

(stabmarg,wcu)= robstab ()计算不确定系统的鲁棒稳定性裕度。这种稳定裕度相对于不确定性水平中指定。鲁棒稳定性的优势大于1意味着系统是稳定的所有值的建模不确定性。鲁棒稳定性的优势小于1意味着一些值的系统变得不稳定不确定因素在其指定的范围。例如,一个边缘0.5意味着以下:

  • 保持稳定,只要不确定元素值保持在0.5归一化单位的名义值。

  • 有不稳定扰动大小0.5归一化单位。

结构stabmarg包含上、下界限实际稳定裕度和稳定裕度的临界频率最小。结构wcu包含不稳定值的确定因素。

例子

(stabmarg,wcu)= robstab (,w)限制了鲁棒稳定性裕度计算指定的频率w

  • 如果w是一个单元阵列的形式{wmin, wmax},然后robstab限制了稳定裕度计算之间的时间间隔wminwmax

  • 如果w是一个向量的频率呢robstab计算鲁棒稳定性裕度只在指定的频率。

例子

(stabmarg,wcu)= robstab (___,选择)计算指定附加选项。使用robOptions创建选择。您可以使用该语法与任何以前的输入参数组合。

例子

(stabmarg,wcu,信息)= robstab (___)返回一个结构与额外的稳定利润的信息和不稳定扰动。看到信息对这种结构的详细信息。您可以使用该语法与任何以前的输入参数组合。

例子

全部折叠

考虑控制系统的植物含有参数不确定性和动态不确定性。创建一个模型的使用不确定的元素。

k =尿素的(“k”10“百分比”,40);δ= ultidyn (“δ”[1]);G =特遣部队(18日1.8 k[1]) *(1 + 0.5 *δ);

创建一个模型的控制器,构建闭环传递函数。

C = pid (2.3、3、0.38、0.001);CL =反馈(G * C, 1);

闭环系统阶跃响应图表明,名义上是稳定的。

步骤(CL.NominalValue)

图包含一个坐标轴对象。坐标轴对象包含一个类型的对象。该对象代表untitled1。

检查闭环系统的鲁棒稳定性。

[stabmarg, wcu] = robstab (CL);stabmarg
stabmarg =结构体字段:下界:1.5960 UpperBound: 1.5993 CriticalFrequency: 4.8627

下界UpperBound领域的stabmarg表明闭环系统的鲁棒稳定性裕度大约是1.6。这个结果意味着系统可以承受更多的不确定性比约60%不确定元素中指定的不稳定。

您可以使用uscale规模系统稳定裕度的不确定性,为全方位的安全检查系统响应的不确定性。规模的不确定性在CL创建一个系统的鲁棒稳定性裕度的最大容许量的不确定性。

CLmaxunc = uscale (CL, stabmarg.UpperBound);CLmaxunc.Uncertainty.delta
不确定线性时不变动态“δ”1输出,1输入,获得小于1.6。
CLmaxunc.Uncertainty.k
不确定的实际参数“k”名义价值10和可变性(-64、64)%。

不确定的元素CLmaxunc射程范围约为1.6倍的原始建模不确定性CL

输出wcu是一个结构,包含最小的微扰kδ使系统不稳定。确认,用这些值不稳定的闭环模型和检查杆的位置。

CLunst = usubs (CL, wcu);极(CLunst)
ans =8×1复杂102我0.0000×-9.9314 + -0.1027 + 0.1009 -0.1027 - 0.1009我0.0000 + 0.0486 0.0000 - 0.0486 -0.0115 -0.0403我-0.0216 + 0.0000 + 0.0000 + 0.0000

结果系统的无阻尼对复杂波兰人与固有频率4.89,这使它不稳定。的CriticalFrequency领域的stabmarg包含相同的值,也就是频率的CL是最近不稳定。

检查的相对灵敏度的鲁棒稳定性裕度不确定系统的元素。考虑一个包含不确定因素的控制系统模型。

k =尿素的(“k”10“百分比”,40);δ= ultidyn (“δ”[1]);G =特遣部队(18日1.8 k[1]) *(1 + 0.25 *δ);C = pid (2.3、3、0.38、0.001);CL =反馈(G * C, 1);

创建一个选项设置robstab,使灵敏度的计算。

选择= robOptions (“敏感”,“上”);

计算鲁棒稳定性,指定信息输出访问关于计算的附加信息。

[stabmarg、wcu信息]= robstab (CL,选择);

检查灵敏度领域的信息

info.Sensitivity
ans =结构体字段:三角洲:80 k: 20

这个字段中的值显示多少改变了每个元素上的归一化摄动影响稳定裕度。例如,灵敏度k是21。这个值意味着给定变化dk在归一化范围的不确定性k导致约21%的变化,或0.21 * dk在稳定裕度。保证金在这种情况下要敏感得多δ,变化幅度约81%的规范化不确定性的变化范围。

考虑一个包含不确定因素的控制系统模型。

k =尿素的(“k”10“百分比”,40);δ= ultidyn (“δ”[1]);G =特遣部队(18日1.8 k[1]) *(1 + 0.5 *δ);C = pid (2.3、3、0.38、0.001);CL =反馈(G * C, 1);

默认情况下,robstab所有频率的稳定裕度计算只有最弱。如何稳定裕度随频率、使用“VaryFrequency”选择robOptions。例如,计算系统的稳定裕度之间的频率点0.1和10 rad / s。

选择= robOptions (“VaryFrequency”,“上”);[stabmarg、wcu信息]= robstab (CL,{0.1, 10},选择);信息
信息=结构体字段:模型:1频率:[19 x1双]范围:[19 x2双]WorstPerturbation: [19 x1 struct]敏感性:[1 x1 struct]

robstab返回向量的频率信息的输出,频率字段。info.Bounds包含稳定裕度的上界和下界在每个频率。使用这些值绘制频率依赖的稳定裕度。

semilogx (info.Frequency info.Bounds)标题(稳定裕度与频率的)ylabel (“保证金”)包含(“频率”)传说(“下限”,“上界”)

图包含一个坐标轴对象。坐标轴对象与标题稳定裕度与频率,包含频率、ylabel边缘包含2线类型的对象。这些对象代表下界,上界。

当你使用“VaryFrequency”选项,robstab选择自动频率点。它选择的频率保证包括稳定裕度的频率最低(在指定的范围内)。显示返回的频率值确认包括临界频率。

info.Frequency
ans =19×10.1000 0.1061 0.1425 0.1914 0.2572 0.3455 0.4642 0.6236 0.8377 1.1253⋮
stabmarg.CriticalFrequency
ans = 4.8269

或者,而不是使用“VaryFrequency”,您可以指定特定的频率计算的鲁棒稳定性的利润率。info.Bounds包含利润在所有指定的频率。然而,这些结果并不保证包括最弱的边缘,这可能指定的频点之间。

w = logspace (1、1,25);[stabmarg、wcu信息]= robstab (CL, w);semilogx (w, info.Bounds)标题(稳定裕度与频率的)ylabel (“保证金”)包含(“频率”)传说(“下限”,“上界”)

图包含一个坐标轴对象。坐标轴对象与标题稳定裕度与频率,包含频率、ylabel边缘包含2线类型的对象。这些对象代表下界,上界。

输入参数

全部折叠

动态系统与不确定性,指定为一个号航空母舰,ufrd,一族,或genfrd模型包含不确定元素。为一族genfrd模型,robstab使用任何可调模块的当前值和折叠成已知的(不确定)模型的一部分。

对于频率特性模型,ufrdgenfrd,robstab假设系统是名义上的稳定。

也可以不确定模型的一个数组。在这种情况下,robstab返回数组中最小的利润在所有模型,和信息输出包含相应的指数模型。

频率计算鲁棒稳定性的利润率,指定为单元阵列{wmin, wmax}或一个向量的频率值。

  • 如果w是一个单元阵列的形式{wmin, wmax},那么函数计算利润之间的频率范围wminwmax

  • 如果w的频率是一个矢量,然后在每个指定频率函数计算利润。例如,使用logspace生成一个行向量与对数间隔的频率值。

    号航空母舰一族模型,当w是一个矢量,robstab(忙,w)相当于robstab (ufrd(忙,w))。因此,名义上必须稳定。

在rad /单位指定的频率TimeUnit,在那里TimeUnitTimeUnit模型的属性。

选择鲁棒稳定性计算的利润率,指定为你创建一个对象robOptions。可用的选项包括设置让你:

  • 提取频率相关稳定利润。

  • 检查的敏感性,每个不确定元素的利润率。

  • 提高稳定裕度计算的结果为底层通过设置特定的选项mussv计算。特别是,设置选项“MussvOptions”“锰”可以减少之间的差距下界和上界。N是重启的数量。

所有可用选项的更多信息,请参阅robOptions

例子:robOptions(“敏感性”,“对”,“MussvOptions”,“m3”)

输出参数

全部折叠

鲁棒稳定性裕度和临界频率作为结构返回包含以下字段:

描述

下界

下界的实际模型的鲁棒稳定性裕度,作为一个标量值返回。稳定裕度的保证不小于下界。换句话说,所有建模不确定性的规范化程度下界,该系统是保证稳定。

UpperBound

上限的实际模型的鲁棒稳定性裕度,作为一个标量值返回。稳定裕度的保证不超过UpperBound。换句话说,存在一些uncertain-element值与此相关大小导致不稳定。robstab返回一个这样的实例wcu

CriticalFrequency

频率的稳定幅度是最小的,在rad /TimeUnit,在那里TimeUnitTimeUnit的属性

鲁棒稳定性的优势大于1说明是稳定的所有值的建模不确定性。一个健壮的利润率小于1意味着稳定变得不稳定的一些值不确定元素在其指定的范围。例如,一个边缘0.5意味着以下:

  • 保持稳定,只要不确定元素值保持在0.5归一化单位的名义值。

  • 有不稳定扰动大小0.5归一化单位。

使用uscale规模的不确定性稳定裕度的检查实际的容许范围的不确定性。

如果的名义价值不稳定,稳定裕度是0。如果是一个ufrdgenfrd模型中,robstab假设它名义上是稳定的。

如果是一个数组不确定的模型,stabmarg包含数组中最小的利润在所有模型。在这种情况下,信息输出包含在其相应的索引模型模型字段。

最小扰动的不确定导致的不稳定的元素,作为结构返回的字段的名称不确定的元素。每个字段包含实际的不稳定值为每一个不确定的因素。例如,如果包括一个不确定的矩阵和输出不确定动力学δ,然后wcu.M是一个数字矩阵和wcu.delta是一个状态空间模型的输出。

使用wcu usubs(忙)这些不确定元素的值替代获取不稳定的动态系统,偏离标称系统的最少。使用actual2normalized将这些实际不确定性值归一化单位稳定裕度的表达。

尿素的参数的范围不是围绕标称值,robstab使调整后的目的分析:

  • 当最严重的扰动(最小的不稳定扰动)不在有效范围内的actual-to-normalized转换(见getLimits),然后robstab设置相应的条目wcu到最近的有效值。换句话说,如果actpert实际单位是最严重的扰动,robgainwcu到最近的价值区间内ActLims返回的getLimits

  • 当没有不稳定扰动robstab设置相应的条目wcu的名义价值尿素的参数。

关于鲁棒稳定性利润率的附加信息,作为结构返回以下字段:

描述

模型

指数模型的稳定裕度最弱,什么时候是一个数组的模型。

频率

频率点robstab返回的鲁棒稳定性,作为一个向量返回。

  • 如果“VaryFrequency”选择robOptions“关闭”,然后info.Frequency临界频率,最小的边缘发生的频率。如果最小下界和上界最小的稳定裕度发生在不同的频率,然后info.Frequency包含这两个频率是一个矢量。

  • 如果“VaryFrequency”选择robOptions“上”,然后info.Frequency包含了频率选择robstab。这些频率保证包括稳定裕度的频率最小。

  • 如果你指定一个向量的频率w在计算稳定的利润info.Frequency= w。当你指定一个频率向量,这些频率不能保证包括稳定裕度的频率最小。

“VaryFrequency”只选择是有意义的号航空母舰一族模型。robstab忽略了选择ufrdgenfrd模型。

界限

上下边界的实际模型的鲁棒稳定性裕度,作为一个数组返回。info.Bounds (: 1)包含在每个相应的频率下限info.Frequency,info.Bounds (: 2)包含相应的上界。

WorstPerturbation

最小的不稳定扰动在每个频率点info.Frequency,作为一个结构数组返回。等领域的info.WorstPerturbation不确定元素的名称吗,每个字段包含的不稳定值在每个频率对应的元素。例如,如果包括一个不确定的参数p和输出不确定动力学δ,然后info.WorstPerturbation.p是一家集数值和info.WorstPerturbation.delta是一家集的输出状态空间模型。

灵敏度

每个不确定元素的稳定利润的敏感性,作为结构时返回“敏感”选择robOptions“上”。等领域的info.Sensitivity不确定元素的名称吗。每个字段包含一个百分比,措施多少对应元素的不确定性影响了稳定裕度。例如,如果info.Sensitivity.p是50,然后给定分数变化的不确定性范围p原因的一半分数稳定裕度的变化。

如果“敏感”选择robOptions(默认设置),然后呢info.Sensitivity

算法

计算鲁棒性在特定频率相当于保证金计算结构奇异值,μ,对于一些适当的块结构(μ分析)。

号航空母舰一族模型,robstab(忙)robstab(忙,{wmin, wmax})使用一个算法,发现最小的保证金在频率。该算法不依赖于频率网格和不间断的不利影响μ结构奇异值。看到获得可靠的估计的鲁棒性为更多的信息。

ufrdgenfrd模型,robstab计算μ在每个频率点的上下界限。这个计算频率点之间并没有提供担保,可以准确如果有不连续或急剧的山峰μ。的语法robstab (uss, w),在那里w是一个向量的频率点,是一样的吗robstab (ufrd (uss, w))也依赖于频率的网格计算。

一般来说,状态空间模型的算法比frequency-gridding更快和更安全的方法。然而,在某些情况下,状态方程算法需要许多μ计算。在这些情况下,指定频率的网格作为一个向量w可以更快,健壮性,提供保证金顺利,频率不同。这种平滑变化与动态不确定性是典型的系统。

版本历史

介绍了R2016b