Using Richards-Wolf vectorial diffraction integral,three-dimensional light field complex amplitude function is presented when linearly polarized Bessel-Gaussian beam passes through a high numerical-aperture system with primary coma aberration.In the case of linearly polarized Bessel-Gaussian beam,the intensity distribution of polarization components at the focal plane are simulated under different comatic aberration coefficients.It is shown that the primary coma aberration makes the intensity distribution transform into uniaxial symmetry from biaxial symmetry,and the position of all high intensity lobes is taken a shift.For large primary coma aberration,interference fringe presents.