minimal balanced realization
slb=minreal(sl [,tol])
:sl,slb syslin lists : :tol real (threshold) :
[ae,be,ce]=minreal(a,b,c,domain [,tol]) returns the balanced realization of linear system sl ( syslin list).
sl is assumed stable.
tol threshold used in equil1.
A=[-`eye`_(2,2),`rand`_(2,2);`zeros`_(2,2),-2*`eye`_(2,2)];
B=[`rand`_(2,2);`zeros`_(2,2)];C=`rand`_(2,4);
sl=`syslin`_('c',A,B,C);
slb=minreal(sl);
`ss2tf`_(sl)
`ss2tf`_(slb)
`ctr_gram`_(sl)
`clean`_(`ctr_gram`_(slb))
`clean`_(`obs_gram`_(slb))