Пример построения и анализа алгоритма оценивания в пакете Octave
sys=ss(A,B,C,D);
sys.InputName='Actuator';
sys.OutputName={'Pendulum angle','Flywheel angle'};
sys.StateName={'Pendulum angle','Pendulum rate','Flywheel angle','Flywheel rate'};
%est=estim(sys2,L')
An=A-L'*C;
Cn=eye(length(A));
est=ss(An,L',Cn,zeros(length(A),2))
syscl=series(sys,est);