Se estudia la generación de múltiplos de marea y mareas compuestas en la laguna costera somero Ensenada de la Paz, Baja California Sur, México. Para ello, se hacen modelos numéricos unidimensionales y bidimensionales usando los métodos de diferencias finitas. Se usan mediciones de la elevación de la superficie del mar y de corrientes en distintos sitios de la laguna para validar y calibrar los modelos. Si se usa un canal rectangular, los campos de elevación y de velocidad discrepan notablemente de las mediciones. Si se incluye la variación de las áreas transversales en el término de la ecuación de continuidad, las discrepancias entre datos modelados y datos medidos disminuyen. La dinámica lineal resultó suficiente para simular el comportamiento de las constituyentes M₂ , S₂ , K₁ y O₁, las cuales son las más energéticas en la laguna. Con los estudios numéricos unidimensionales no-lineales resultó que en el interior de la laguna se generan ondas de marea tercidiurnas, cuartidiurnas y sextidiurnas, específicamente se generan SK₃, MK₃, 2MS₆, 2SM₆, MO₃ y M₄, en ese orden de importancia, lo cual concuerda con las mediciones. Se estudia la importancia relativa de los términos no-lineales en el proceso de generación. Continuidad no-lineal fue el mecanismo de generación más importante de las componentes tercidiurnas y cuartidiurnas y fricción cuadrática de las sextidiurnas. La magnitud de la velocidad de las componentes tercidiurnas en el canal resultaron del mismo orden de magnitud que las velocidades de las componentes semidiurnas en el interior de la laguna, ~2 cms⁻¹. El modelo numérico bidimensional elaborado fue validado con los estudios numéricos unidimensionales y con mediciones de campo. El modelo reproduce la dinámica no-lineal inducida por mareas en la laguna. El patrón horizontal de la componente de marea M₂ resultó similar al patrón observado en los estudios de campo. The generation of overtides and compound tides in the shallow coastal lagoon of Ensenada de la Paz, Baja California, Mexico is studied. One-dimension and two-dimension numerical models are implemented using finite differences. Sea surface elevation and horizontal currents observation are used to validate and to calibrate the models. If a rectangular channel is used, important discrepancies in the predicted values for elevation and velocity fields are found with respect to the measured ones. These discrepancies diminish by allowing into the continuity equation, variations of the transversal cross section. Linear dynamics proved to be enough to simulate the progress of the main constituents, M₂ , S₂ , K₁ and O₁, which are the more energetic in the lagoon. Using the one-dimensional nonlinear numerical model produced the generation of ter-diurnal, four-diurnal and sixth-diurnal harmonics SK₃, MK₃, 2MS₆, 2SM₆, MO₃ and M₄, in the interior of the lagoon, in agreement with the measurements. The relative importance of the non-linear terms in the generation process is studied. Non-linear continuity was the mechanism responsible for the generation of ter-diurnal and four-diurnal harmonics and cuadratic friction of sixth-diurnal harmonics. The velocity of the ter-diurnal constituent turned out to be of the same order of magnitude that the velocity of the semidiurnal constituent in the interior of the lagoon, ~2 cms⁻¹. The two-dimensional model was validated by using results from the one-dimensional model and field observations. it reproduced the lagoon's dynamics induced by tides. The horizontal pattern of the M₂ tide was in agreement with field observations.