We consider non-linear sampled-data control systems consisting of a continuous linear part and a non-linear pulse modulator ( PM). The PM modulates a sequence of rectangular pulses by sign, frequency and duration, as a function of a linear combination of discrete values of the system coordinates. For the stated systems the concepts of the natural phase space E m with coordinates x n , x 1 n , …, x n ( m−1) and of the difference phase space D m with coordinates x n , x n+1 , …, x n+ m−1 are used. Non-linear difference equations of motion have been found in both the phase spaces, and the relation between them has been established. A sufficient condition for the asymptotic stability in-the-large of the systems being investigated is obtained on the basis of discrete analogs of the Lyapunov and LaSalle theorems. The condition has the form of a transcendental inequality, uniquely determined by the system parameters and by the Lyapunov function chosen. The problems of absolute stability and of the region of asymptotic stability also are solved with the aid of the inequality obtained. A procedure for determining the steady-state parameters of the system tracking a ramp input is given. For this steady-state mode sufficient conditions are obtained for stability in-the-large on the basis of Lyapunov's direct method. A discrete analog of the Yoshizawa theorem for the limit boundedness of non-linear sampled-data systems is formulated. On the basis of this theorem a procedure is proposed for determining the boundaries of the asymptotically stable set inside which all the phase trajectories of the system terminate.