This paper presents a method of time-domain synthesis of dynamical systems treating the process signals as distributions or generalized functions, in the manner originally established by L. Schwartz. The technique of synthesizing the transfer function or the state-space model of a system employs exponentially weighted series of the generalized time derivatives of the impulse distribution, known also as the Poisson moment functional expansion. General algorithms are developed and their application is illustrated in typical eases.