A numerically robust algorithm for computing ellipsoidal bounds on the state of a linear, discrete-time dynamic system was proposed.

When using the system identification methods to establish the dynamic thermal model of the building wall, the order of the discrete time model was usually set up to be constant value.

To start with, the system dynamic equation is transformed into discrete-time state equation in observable canonical form.