Matrix Li-Yau-Hamilton estimates for the heat equation on Kähler manifolds