By using the method of matrix decomposition and expanding the aperture function into a finite sum of complex Gaussian functions, the closed-form propagation equ ations of Gaussian beams through a multi-apertured imaging system of B=0 are der ived and are illustrated with a few typical numerical examples. The advantages o f our treatment are pointed out.