In this work, a homotopy optimization method is proposed for reconstructing unknown inputs in non-linear dynamical systems. The unknown inputs are parametrized using a B-spline basis. This parameterization of inputs converts the unknown input identification problem into a parameter identification problem. The unknown parameters are identified through an optimization process. The proposed homotopy-based optimization method is designed to converge to the global optimal solution instead of a local minimum. The unknown parameters are obtained through a series of iterations guided by a homotopy parameter and an optimization algorithm. © 2019 R. Manikantan, et al.