加性误差法近似植物模型
给定一个系统G
在LTI形式中,以下命令将系统降低到您指定的任意所需顺序。判断是基于它的汉克尔奇异值。
rng (1234“旋风”);G = rss(30,4,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] = sigma(G-G1(:,:,1));重对数(w, sv, w, info1.ErrorBound(1) *(大小(w)))的包含(“rad /秒”); ylabel (“SV”);标题(“误差范围和模型误差”)