An approach is proposed to carry out multivariate calculations of the magnetic field distribution in the working gaps of a plate polygradient matrix of an electromagnetic separator, based on a combination of the advantages of two- and three-dimensional computer modeling. Two-dimensional geometric models of computational domains are developed, which differ in the geometric dimensions of the plate matrix elements and working air gaps. To determine the vector magnetic potential at the boundaries of two-dimensional computational domains, a computational 3D experiment is carried out. For this, three variants of the electromagnetic separator are selected, which differ in the size of the working air gaps of the polygradient matrices. For them, three-dimensional computer models are built, the spatial distribution of the magnetic field in the working intervals of the electromagnetic separator matrix and the obtained numerical values of the vector magnetic potential at the boundaries of the computational domains are investigated. The determination of the values of the vector magnetic potential for all other models is carried out by interpolation. The obtained values of the vector magnetic potential are used to set the boundary conditions in a computational 2D experiment. An approach to the choice of a rational version of a lamellar matrix is substantiated, which provides a solution to the problem according to the criterion of the effective area of the working area. Using the method of simple enumeration, a variant of the structure of a polygradient matrix with rational geometric parameters is selected. The productivity of the electromagnetic separator with rational geometric parameters of the matrix increased by 3–5 % with the same efficiency of extraction of ferromagnetic inclusions in comparison with the basic version of the device