Abstract-In many detection and estimation problems associated with processing of second-order stationary random processes, the observation data are the sum of two zero-mean second-order stationary processes: the process of interest and the noise process. In particular, the main performance criterion is the signal-to-noise ratio (SNR). After linear filtering, the optimal SNR corresponds to the maximal value of a Rayleigh quotient which can be interpreted as the largest generalized eigenvalue of the covariance matrices associated with the signal and noise processes, which are block multilevel Toeplitz structured for -dimensional vector-valued second-order stationary -dimensional random processes x . In this paper, an extension of Szegö's theorem to the generalized eigenvalues of Hermitian block multilevel Toeplitz matrices is given, providing information about the asymptotic distribution of those generalized eigenvalues and in particular of the optimal SNR after linear filtering. A simple proof of this theorem, under the hypothesis of absolutely summable elements is given. The proof is based on the notion of multilevel asymptotic equivalence between block multilevel matrix sequences derived from the celebrated Gray approach. Finally, a short example in wideband space-time beamforming is given to illustrate this theorem.Index Terms-Asymptotic distribution, block multilevel Toeplitz matrix, generalized eigenvalues, multidimensional second-order vector valued stationary random process, Szegö's theorem.