We derive the partition function of the one-body and two-body systems of classical noncommutative
harmonic oscillator in two dimensions. Then, we employ the path integral approach to the quantum noncommutative
harmonic oscillator and derive the partition function of the both systems at finite temperature.