I was wondering if anyone could explain the reasoning behind the $h$ normalization constant when calculating the partition function for a classical harmonic oscillator.
I know that the partition function is $Z=\sum e^{-\beta H(\vec{x},\vec{p})}$ where $H=\frac{p^2}{2m}+\frac{kx^2}{2}$. I also know that to compute this we need to do an integral instead of a sum (as $\vec{p}$ and $\vec{x}$ are (classically) continuous).
So I understand why we go from $Z=\sum e^{-\beta H(\vec{x},\vec{p})}$ to $$Z=\frac{1}{N^3}\int d^3p \int d^3x e^{-\beta H(\vec{x},\vec{p})}$$
where $[N] \sim \mbox{kg} ~ \mbox{m}^2 ~ \mbox{s}^{-1} $ for dimensions.
I just don't understand how we go from this to concluding that $N=h$.
I understand it has something to do with $dx ~ dp \sim h$, but I'm finding it difficult to motivate/justify this from a non-quantum perspective, seeing as this is a classical oscillator. Is there any way, or does this just have to be motivated from a quantum perspective?