We embed Feynman integrals in the subvarieties of Grassmannians through homogenization of the integrands in projective space, then obtain Gel'fand-Kapranov-Zelevinsky systems satisfied by those scalar integrals. The Feynman integral can be written as linear combinations of the hypergeometric functions of a fundamental solution system in neighborhoods of regular singularities of the Gel'fand-Kapranov-Zelevinsky system, whose linear combination coefficients are determined by the integral on an ordinary point or some regular singularities. Taking some Feynman diagrams as examples, we elucidate in detail how to obtain the fundamental solution systems of Feynman integrals in neighborhoods of regular singularities. Furthermore we also present the parametric representations of Feynman integrals of the two-loop self-energy diagrams that are convenient to embed in the subvarieties of Grassmannians.