Se presenta un modelo para caracterizar sistemas planetarios. Para sistemas en tránsito, cuya caracterización es nuestra principal motivación científica, el modelo permite obtener estimaciones coherentes de algunos parámetros físicos (masas y radios) y orbitales (inclinación, excentricidad, longitud del semieje mayor, argumento del periastro) a partir de ajustes simultáneos de series temporales de datos, que son generalmente curvas de luz y curvas de velocidad radial. El modelo usado para realizar los ajustes se basa en un modelo de N cuerpos que se integra numéricamente, como alternativa a los modelos analíticos usados comúnmente. Su uso pretende alcanzar la máxima precisión posible en el análisis de sistemas planetarios múltiples en los que existe interacción dinámica apreciable entre dos o más planetas miembros. Hemos validado el modelo usándolo para extraer los parámetros de varios sistemas planetarios conocidos (alrededor de las estrellas huésped HD 209458, HAT-P-14 y WASP-44) y concluimos que el modelo se encuentra en un estado funcional. El modelo admite cierto margen de mejora, como la incorporación de algunos efectos físicos adicionales para un tratamiento más preciso de los datos y de mejoras en el método de integración con el objetivo de alcanzar una mayor eficiencia computacional. Este modelo es útil para el análisis de datos obtenidos con los espectrógrafos de nueva generación como, por ejemplo, CARMENES, cuyo objetivo más ambicioso es la búsqueda de planetas telúricos en la zona habitable de estrellas de baja masa. El modelo también servirá para el desarrollo del proyecto PHASES, que es un telescopio espacial de bajo costo para la obtención de espectros estelares de baja resolución calibrados de forma absoluta.