The relation between chaotic dynamics of nonlinear Hamiltonian systems and equilibrium statistical mechanics in its canonical ensemble formulation has been investigated for two different nonlinear Hamiltonian systems. We have compared time averages obtained by means of numerical simulations of molecular dynamics type with analytically computed ensemble averages. The numerical simulation of the dynamic counterpart of the canonical ensemble is obtained by considering the behavior of a small part of a given system, described by a
microcanonical ensemble, in order to have fluctuations of the energy of the subsystem. The results for the Fermi-Pasta-Ulam model (i.e., a one-dimensional anharmonic solid) show a substantial agreement between time and ensemble averages independent of the degree of stochasticity of the dynamics. On the other hand, a very different behavior is observed for a chain of weakly coupled rotators, where linear exchange effects are absent. In the high-temperature limit (weak coupling) we have a strong disagreement between time and ensemble averages for the specific heat even if the dynamics is chaotic. This behavior is related to the presence of spatially localized chaos, which prevents the complete filling of the accessible phase space of the system. Localized chaos is detected by the distribution of all the characteristic Liapunov exponents.