B_MIN_SIMULATION estimates matrix B and the initial state x0 of the LDS by minimizing the simulation error. [B, x0, M1] = B_min_simulation(A,C,y,u,T); Input A : state matrix of the LDS C : matrix of the LDS y : output data organized as row vectors [ y(0) y(1) ... y(T-1)] u : input data organized as row vectors [ u(0) u(1) ... u(T-1)] T : length of the data Output B : matrix of the LDS system x0: initial Condition of the state of the LDS M1: Generalized Observability matrix of the LDS