给定一个系统G
在LTI形式中,以下命令将系统缩减为您指定的任何所需顺序。判断调用是基于其汉克尔奇异值。
rng (1234“旋风”);4 G = rss(30日,3);%随机30状态模型%平衡截断模型尺寸12:16[G1, info1] = balancmr (G, 12:16);% Schur通过指定“MaxError”来平衡截断[G2, info2] = schurmr (G,“MaxError”, 0.8, 0.5, 0.2);σ(G,“b -”G1,“r——”G2,“g -”。)传说(‘G’,G1的,“G2”)
图中比较了原始模型G
通过简化模型G1
和G2
.
为了确定理论误差界限是否满足,计算原始系统增益与简化系统增益之间的峰值频率差。将其与存储在信息
结构。
规范(G-G1 (:,: 1),“正”)
ans = 2.2965
info1.ErrorBound (1)
ans = 9.7120
或者,通过以下命令绘制模型错误和错误界限:
[sv, w] =σ(G-G1 (:,: 1));重对数(w, sv, w, info1.ErrorBound(1) *(大小(w)))的包含(“rad /秒”); ylabel (“SV”);标题(误差界与模型误差)