A method for setting up a time-periodic Lyapunov function for a linear system with periodic coefficients based on an auxiliary matrix-valued function is proposed and developed. New sufficient conditions for the stability of large-scale periodic systems decomposed into an even number of subsystems are formulated Keywords: large-scale periodic system, linear system with periodic coefficients, block-diagonal matrix-valued function, time-periodic Lyapunov function, asymptotic stability conditions Introduction. In natural sciences and engineering, motions repeating in equal and approximately equal time intervals are attributed to oscillatory phenomena [10,14,21]. Only feasible motions of this kind are of practical interest in mechanics. The principle of choice (from the set of all theoretically possible solutions of the corresponding equations of motion) tells us that such motions correspond to stable solutions.The positive experience of using the direct Lyapunov method in developing reliable approaches to the stability analysis of large-scale systems and systems with an infinite number of degrees of freedom and the justified universality of the Lyapunov function method [1,6,8,12,22] with its generalizations and modifications [3,5,7,9,11], which are currently used to develop the method of multicomponent functions [4,19,20] and establish new stability conditions for some classes of systems [13,15,16,18], give importance to the problem of setting up an approximating V-function for systems of linear differential equations with periodic coefficients [2].We will employ the result obtained in [2,17] as to a method of setting up a Lyapunov function based on the general theorem of asymptotic stability of periodic system (based on a matrix-valued auxiliary function). Generalizing this result, we will formulate sufficient conditions for the uniform asymptotic stability of large-scale periodic systems. We will also analyze the general case of a system of arbitrary order that has matrix trigonometric polynomials as coefficients and consists of an even number of subsystems. The idea of forming the off-diagonal elements of a block-diagonal matrix-valued function by solving Fredholm equations of the second kind will be applied to the system under consideration. The derivation of diagonal elements is reduced to quadratures.1. Lyapunov function. Stability Conditions. Consider a system