Hyperbolic matrix functions are essential for solving hyperbolic coupled partial differential equations. In fact the best analytic-numerical approximations for resolving these equations come from the use of hyperbolic matrix functions. The hyperbolic matrix sine and cosine sh (A), ch (A)(A∈ Mr (C)) can be calculated using numerous different techniques. In this article we derive some explicit formulas of sh (tA) and ch (tA)(t∈ R) using the Fibonacci-Hörner and the polynomial decomposition, these decompositions are calculated using the generalized Fibonacci sequences combinatorial properties in the algebra of square matrices. Finally we introduce a third approach based on the homogeneous linear differential equations. And we provide some examples to illustrate your methods.