The fractional Fourier transform (FRFT) is applied to treat the propagation of annular flat-topped beams. Based on the definition of FRFT in the cylindrical coordinate system, analytical formulae are derived for annular flat-topped beams through the FRFT optical systems. By using the formulae, the properties of annular flat-topped beams in the FRFT plane are illustrated numerically. The results show that the intensity distribution properties in the FRFT plane are closely related to the fractional order of the FRFT optical system and initial beam parameters.The fractional Fourier transform (FRFT) has become an international special research subject in optics since it was firstly introduced by Ozaktas, Mendlovic and Lohmann [1][2][3] in 1993. It has been extensively applied in signal and image processing, beam shaping and beam analysis, etc [4][5][6] . The optical signal processing and beam control by using FRFT are no longer limited by the space, compared with those by using the Fourier transform. The observing plane of an arbitrary fractional order can be chosen freely according to the practical need. The FRFT for various laser beams have been studied, such as flattened Gaussian beams [7] , elliptical Gaussian beams [8] , hollow Gaussian beams [9] , partially coherent Gaussian-Schell beams [10,11] and the beams generated by a Gaussian mirror resonator [12] . In the past, more attention was paid to the simple irradiance distributions such as those varying from Gaussian to flat-topped distributions. In reality, the unstable resonator configurations may produce some complex output intensity distributions in some high-average-power and high-beamquality laser sources. One of the most typical examples is that the copper-vapor lasers output intensity profile is a flat-topped distribution with an axial shadow, when the onaxis unstable resonators and injection-seeded copper-vapor laser oscillators are used [13] . To model the beam with a doughnut-shaped profile in the near field, Saghafi et al [14] proposed an analytical expression for it in the cylindrical coordinate system, which was expressed as a superposition of two flattopped Gaussian beams and was extremely in good agreement with the experimental result.In this paper, following Tovar's [15] research, where a multiGaussian beam shape was proposed as a model for laser beam profile that has a nearly flat top, we develop an analytical expression in the cylindrical coordinate system to describe the beam with a nearly flat-topped profile and an on-axial shadow simultaneously, namely, the annular flat-topped beam. The FRFT for an annular flat-topped beam is studied.As numerical examples, the intensity distribution of the annular flat-topped beam in the FRFT plane is calculated and a simple conclusion is obtained.The optical system for performing FRFT is depicted in Fig.1.