DETERMINISTIC_SUBSPACE estimates the parameters of the Linear Dynamical System [A,B,C,H,Obs,x0,H1] = deterministic_subspace(y,u,i,n,large_min); Inputs y: matrix of measured outputs u: matrix of measured inputs i: number of block rows in Hankel matrices (i * #outputs) is the max. order that can be estimated Typically: i = 2 * (max order)/(#outputs) n: order of the system large_min: (true or false) is it a large scale minimization Outputs A,B,C : deterministic state space system x_{k+1) = A x_k + B u_k y_k = C x_k Example [A,B,C,] = det_stat(y,u,10,5, true); References Subspace Identification for Linear Systems Theory - Implementation - Applications Peter Van Overschee / Bart De Moor Kluwer Academic Publishers, 1996, Page 52 & Alessandro Chiuso , Giorgio Picci System Identification Notes 2005