This paper uses a new framework for solving a class of linear matrix differential equations. For doing so, the operational matrix of the derivative based on the shifted Bernstein polynomials together with the collocation method are exploited to decrease the principal problem to system of linear matrix equations. An error estimation of this method is provided. Numerical experiments are reported to show the applicably and efficiency of the propounded method. Manuscript profile