给定一个系统G
在LTI形式中,下面的命令将系统简化为您指定的任何所需顺序。判断是基于它的Hankel奇异值。
rng (1234“旋风”);4 G = rss(30日,3);%随机30状态模型%平衡截断模型与尺寸12:16[G1, info1] = balancmr (G, 12:16);%通过指定' 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”);标题(“误差范围和模型误差”)