Propagation of optical field through an optical system characterized by an ABCD ray matrix has been investigated in the spatial frequency domain. The behavior of the angular spectrum of the optical field can be directly obtained by use of the diffraction integral in spatial frequency domain which has been derived by this paper. Some examples have shown that for some problems it is more convenient to study in frequency domain by using the frequency-domain diffraction integral formula than in spatial domain.