Given commuting families of Hermitian matrices {A(1), ... , A(k)} and {B1, ... , B(k)}, conditions for the existence of a completely positive map Phi, such that Phi (Aj)B(j) for j = 1, ... , k, are studied. Additional properties such as unital or/and trace preserving on the map Phi are also considered. Connections of the study to dilation theory, matrix inequalities, unitary orbits and quantum information science are mentioned.