In this paper, by using the hybrid algorithm merging hill-climbing with simulated annealing, the precise design of the diffractive optical element (DOE) for beam smoothing in fractional Fourier transform domain is completed. True beam smoothing has been obtained, the intensity distribution of any point, besides the used sampling point in the optimization, on the output plane is consistent with the demand for beam smoothing.