It's possible, because that's partly how VSOP was made, but the amount of effort involved is immense. What you are describing is at the very least a doctoral thesis, if not several of them.
Why not just use the SPICE kernels directly?
Analytic ephemerides are generally power series of trig functions of polynomials in orbital resonances, as described for the VSOP family of models in Simon, Francou, Fienga, and Manche (2013), "New analytical planetary theories VSOP2013 and TOP2013", Astronomy & Astrophysics 557 A49. Understanding what's going on in that paper, and what the words they use actually mean, requires good working knowledge of celestial mechanics at the level of Richard Battin's An Introduction to the Mathematics and Methods of Astrodynamics (1999), which begins with continued fraction expansions for elliptic integrals, and gets harder from there.
All the VSOP theories are equations for the six equinoctial elements (a, $ \lambda$, k, h, q, p) for each planet (originally eight, but expanded to include Pluto in 2010). The theory starts out with lots of parameters, but in the end each equation uses only time as a variable, because all of the other things had their values adjusted to fit measured data. Actually, in VSOP's case, they fit not to the measurements directly, but instead to the output of a numerical integration, which was fit to actual measurements. VSOP2000 fit to DE403, VSOP2010 fit to DE405, and VSOP2013 fit to the Paris Observatory's INPOP10. For more information, consult Fienga (2009), "Evolution of INPOP planetary ephemerides", and Hilton and Hohenkerk (2011), "A comparison of the high accuracy planetary ephemerides DE421, EPM2008, and INPOP08".
The full 54 equations of VSP2013 can be found at https://ftp.imcce.fr/pub/ephem/planets/vsop2013/solution/ . VSOP2013-secular.dat is the easiest part to understand. The full theory, using all of the files in that directory, involves more than 2.5 million coefficients. Choosing good criteria for discarding many of those terms is the topic of the last page of https://ftp.imcce.fr/pub/ephem/planets/vsop2013/solution/README.pdf , the effects of which on the number of coefficients is summarized in Table 7 of Simon, et al. (the first link above). There is also a Chebyshev polynomial version of the VSOP2013 model, in the neighboring https://ftp.imcce.fr/pub/ephem/planets/vsop2013/ephemerides/ .