Based on Bohr model, we have presented a general formalism describing the collective motion for any deformed system, in which the collective Hamiltonian is expressed as vibrations in the body-fixed frame, rotation of whole system around the laboratory frame, and coupling between vibrations and rotation. Under the condition of decoupling approximation, we have derived the quantized Hamiltonian operator. Based on the operator, we have calculated the rotational spectra for some special octupole and hexadecapole deformed systems, and shown their dependencies on deformation. The result indicates that the contribution of octupole or hexadecapole deformations to the lowest band is regular, while that to higher bands is dramatic. These features reflecting octupole and hexadecapole deformations are helpful to recognize the properties of real nuclei with octupole and/or hexadecapole deformations coexisting with quadrupole deformations.