We have modeled some magnetic systems, which consists of a number of sublattices, in the mean field approximation. This is possible in crystalline systems formed by two or more magnetic ions coupled by specific interactions such as the crystal field, exchange, among others. The main idea is to solve the microscopic Hamiltonian that models a given magnetic system in order to obtain their magnetic equation of state: M (H, T). For this, we use the appropriate arrangements to different magnetic (ferro-, ferri- and antiferro- magnetic) scheme sublattices. From the solutions of the Hamiltonian (eigenvalues and eigenvectors), physical quantities of interest were determined. In principle we consider systems with localized magnetism due to 3d and 4f electrons with the participation of non-magnetic ligands including 3d-4f systems with the presence of crystal field. In this dissertation we use the model of two-and three sublattices in order to obtain the equation of state for the following systems: RKKY exchange in RNi2B2C, superexchange in (Y3-zRz)(T1xFe1-x)(T2yFe3-y)O12, LixFe3- xO4 and (NixMn1-x)1.5[Cr(CN)6]. In these formulas, R represents a rare earth ion, T1 and T2 represent non-magnetic ions. Some representative cases are presented to illustrate the different equations of state and behavior of sublattices, metamagnetism, temperature compensation, etc. The extension to other similar systems can be direct or need to incorporate additional phenomenological parameters. / Consideramos neste estudo sistemas magnéticos modeláveis mediante sub- redes na aproximação do campo médio. Isto é possível em sistemas cristalinos formados por dois ou mais íons magnéticos acoplados por interações específicas como as do campo cristalino, troca, entre outros. A ideia central é resolver o hamiltoniano microscópico que modela um determinado sistema magnético de maneira a se obter sua equação de estado magnética: M(H,T). Para isto usamos o esquema de sub- redes adequado a diferentes arranjos magnéticos (ferro-, ferri- e antiferro- magnéticos). Com as soluções do hamiltoniano (autovalores e autovetores), grandezas físicas de interesse foram determinadas. Em princípio, consideramos sistemas com magnetismo localizado devido a elétrons 3d e 4f com participação de ligandos não magnéticos incluindo sistemas 3d-4f com presença de campo cristalino. Neste trabalho de dissertação empregamos o modelo de duas e três sub- redes para obter as equações de estado nos seguintes sistemas: troca tipo RKKY em RNi2B2C, supertroca em (Y3-zRz)(T1xFe1-x)(T2yFe3-y)O12, LixFe3-xO4 e (NixMn1- x)1.5[Cr(CN)6]. Nestas fórmulas, R representa um íon de terra rara, T1 e T2 representam íons não magnéticos. Alguns casos representativos são apresentados para exemplificar as diferentes equações de estado e comportamento das sub- redes, metamagnetismo, temperatura de compensação etc. Em princípio, a extensão para outros sistemas semelhantes, usando os modelos apresentados aqui, pode ser direta ou precisar incorporar parâmetros fenomenológicos adicionais.