This paper studies the correlation properties of Walsh functions. It is found that Walsh functions generated by a Hadamard matrix could be divided into several cyclically orthogonal subsets. An approach to obtain such subsets is proposed as a theorem and is proved through the mathematical induction. The application of the cyclically orthogonal subsets of Walsh functions to DS-CDMA communications is then discussed.