International audienceLaser beam can carry multi-scales properties in space and time that impact the beam quality. The study of their evolution along complex optical sequences is of crucial interest, particularly in high intensity laser chains. For such analysis, results obtained with standard numerical calculus are strongly dependent of the sampling. In this paper, we develop an analytic model for a sinusoidal phase modulation inside a sequence of first order optics elements based on the Wigner matrix formalism. A Bessel decomposition of the Wigner function gives pseudo-Wigner functions that obey to the general ABCD matrix law transformation without approximations and sampling considerations. Applied to a Gaussian beam, explicit expression...