Two classes of multidimensional phase control systems with differentiable vector periodic functions are considered, the class of continuous systems described by ordinary differential equations and the class of discrete systems described by difference equations. The number of cycle slips for angular coordinates in phase systems with differentiable nonlinearities is studied. The study is based on the direct Lyapunov method and uses periodic Lyapunov functions, extensions of the phase space of the system, and the Yakubovich-Kalman lemma. This lemma provides necessary and sufficient conditions for the existence of Lyapunov functions by using the transfer matrix of the linear part of the system. As a result, for phase systems possessing global asymptotics, frequency criteria making it possible to sharpen estimates for the deviation of angular coordinates from their initial values are obtained. These criteria contain multiparameter frequency inequalities with variable parameters satisfying certain algebraic inequalities.