With the development of laser technology in the field of optics, ultra-fast optics has become an important research field. Compared with the traditional technology, ultrafast optics can be realized not only under shorter pulse function, but also on a smaller scale, which can more quickly reflect the dynamic process. We present an analytical calculation of the full three-dimensional (3D) coherent spectrum with a finite duration two-dimensional (2D) Gaussian pulse envelope. Our starting point is the solution of the optical Bloch equations for three-level potassium atomic gas in the 3D time domain by using the projection-slice theorem, error function and Fourier-shift theorem of 3D Fourier transform. These principles are used to calculate and simplify the third-order polarization equation generated by the device, and the analytical calculation of three-dimensional Fourier transform frequency spectrum at T = 0 is obtained. We simulate the analytic solution by using mathematics software. By comparing the simulations with the experimental results, with the homogeneous line-width fixed, we can obtain the relationship among the in-homogeneous broadening, the correlation diagonal coefficients and the three-dimensional spectrum characteristics, which can be identified quantitatively by fitting the slices of three-dimensional Fourier transform spectrum peaks in an appropriate direction. The results show that the three-dimensional Fourier transform spectrum will extend along the diagonal direction with the increasing of the in-homogeneous broadening, and the spectrogram progressively becomes a circle with the increasing of the diagonal correlation coefficient, and the amplitude also gradually turns smaller. According to the analytical solution, we give a complete two-dimensional spectrum of the T = 0 interface. The results can be fit to the experimental 3D coherent spectrum for arbitrary inhomogeneity.