By the quadratic form theory, we succeed in diagonalizing the Hamiltonian for three-dimensional anisotropic coupled harmonic oscillators, through three symplectic transformations keeping the commutation relations unchanged. The quantized energy spectrum and the exact wave function of the system are given.