We mainly solve three problems. Firstly, by the decomposition of the (anti-)Hermitian generalized (anti-)Hamiltonian matrices, the necessary and sufficient conditions for the existence of and the expression for the (anti-)Hermitian generalized (anti-)Hamiltonian solutions to the system of matrix equationsAX=B,XC=Dare derived, respectively. Secondly, the optimal approximation solutionminX∈K∥X^-X∥is obtained, whereKis the (anti-)Hermitian generalized (anti-)Hamiltonian solution set of the above system andX^is the given matrix. Thirdly, the least squares (anti-)Hermitian generalized (anti-)Hamiltonian solutions are considered. In addition, algorithms about computing the least squares (anti-)Hermitian generalized (anti-)Hamiltonian solution and the corresponding numerical examples are presented.