In this paper, two constructions of mutually orthogonal zero correlation zone polyphase sequence sets are presented. The first one is based on DFT matrices and interleaving iteration. After each recursive step, the period of sequence and the length of zero-correlation zone are two times larger than that in the last step. The second method, based on DFT matrices and orthogonal matrices, can generate numbers of mutually orthogonal optimal ZCZ sequence sets whose parameters reach the theoretical bounds by using interleaving and shifting techniques. As a result, the algorithms proposed can provide more sequences for the QS-CDMA (quasi-synchronous CDMA) systems.