  We present the gluon distribution functions of pion and kaon in small-$x$ and
large-$x$ regions, and compare them with the results obtained from lattice QCD
and continuum Schwinger function methods. Whether in the small-$x$ region or
the large-$x$ region, our gluon distribution of pion is consistent with the
results of lattice QCD and continuum Schwinger function method. In addition,
the first four moments of gluon distributions of pion and kaon at different
$Q^{2}$ scales are calculated. Furthermore, we present the mass decompositions
of pion and kaon with the dynamical parton distribution functions calculated by
the DGLAP equation with parton-parton recombination corrections. The mass
structures of pion and kaon are totally different from that of the proton.
