Based on Fraunhofer diffraction integral formula， using convolution theorem and the property of zero-order Hank transform， the analytical mathematical formula of light intensity distribution on the back focal plane of the lens after the fundamental transverse mode Gaussian beam perpendicularly enters the axicon-lens combination optical system is derived， and the mathematical relationship between the bottle beam intensity distribution and the main optical parameters in the axicon-lens combination is established， which provides theoretical support for the rapid simulation of bottle beam intensity distribution and the design of optical system for generating bottle beam. According to the derived formula， the numerical simulation and experiment are carried out， obtaining a uniformly distributed bottle beam with similar ring size， respectively， which verifies the correctness of the theoretical results.