We present a new model to predict diffraction patterns of femtosecond pulses through complex optical systems. The model is based on the extension of an ABCD matrix formalism combined with generalized Huygens-Fresnel transforms (already used in the CW regime) to the femtosecond regime. The model is tested to describe femtosecond digital in-line holography experiments realized in situ through a cylindrical Plexiglas pipe. The model allows us to establish analytical relations that link the holographic reconstruction process to the experimental parameters of the pipe and of the incident beam itself. Simulations and experimental results are in good concordance. Femtosecond digital in-line holography is shown to allow significant coherent noise reduction, and this model will be particularly efficient to describe a wide range of optical geometries. More generally, the model developed can be easily used in any experiment where the knowledge of the precise evolution of femtosecond transverse patterns is required.