A technique for nonlinearization of the Lax pair for the scalar soliton equations in (1+1) dimensions is applied to the symmetric matrix KdV equation. As a result, a pair of finite-dimensional integrable Hamiltonian systems, which are of higher rank generalization of the classic Gaudin models, are obtained. The integrability of the systems are shown by the explicit Lax representations and r-matrix method.