给予一个系统G
在LTI表单中,以下命令将系统减少到您指定的任何所需顺序。判断呼叫是基于其Hankel奇异值。
RNG(1234,'twister');g = RSS(30,4,3);%随机30状态模型尺寸平衡截断与尺寸12:16的型号[G1,INFO1] = BALANCMR(G,12:16);通过指定`maxerror'%schur平衡截断[G2,INFO2] = SCHURMR(G,'maxerror',[1,0.8,0.5,0.2]);西格玛(g,'b-',g1,'r--',g2,'G-。') 传奇('G'那'g1'那'g2'的)
该曲线比较原始模型G
随着较小的模型G1.
和G2.
。
为了确定是否满足理论误差绑定,计算原始系统增益与减少系统之间的频率的峰值差异。将其与存储在中的错误进行比较信息
结构体。
规范(G-G1(:,:,1),'inf'的)
ANS = 2.2965.
Info1.Errorbound(1)
ANS = 9.7120.
或者,通过以下命令绘制模型错误与错误绑定:
[sv,w] = sigma(g-g1(:,:,1));loglog(w,sv,w,Info1.Erlorbound(1)*那些(大小(w)))xlabel('rad / sec'); ylabel('sv');标题('错误绑定和模型错误'的)