Let H m×n be the set of all m × n matrices over the real quaternion algebra. We call that A ∈ H n×n is η-Hermitian if A = A η * , where A η * = −ηA * η, η ∈ {i, j, k}, i, j, k are the quaternion units. In this paper, we derive some solvability conditions and the general solution to a system of real quaternion matrix equations. As an application, we present some necessary and sufficient conditions for the existence of an η-Hermitian solution to some systems of real quaternion matrix equations. We also give the expressions of the general η-Hermitian solutions to these systems when they are solvable. Some numerical examples are given to illustrate the results of this paper