The aim of this work is to derive an accurate model of two-dimensional switched control heating system from data generated by a Finite Element solver. The nonintrusive approach should be able to capture both temperature fields, dynamics and the underlying switching control rule. To achieve this goal, the algorithm proposed in this paper will make use of three main ingredients: proper orthogonal decomposition (POD), dynamic mode decomposition (DMD) and artificial neural networks (ANN). Some numerical results will be presented and compared to the high-fidelity numerical solutions to demonstrate the capability of the method to reproduce the dynamics.