Your version of $\psi$ is (I'm sure you know) derived from the Time Independent Schrodinger equation, $\hat{h}\psi=E\psi$. To find time-dependent solutions, we solve $$i\frac{\partial}{\partial t}\psi=\hat{H}\psi.$$ You were trying to solve for stationary states, and the entire point of those is that $|\psi(x)|^2$ does not change over time. Hence $$\frac{d}{dt}\left<x\right>=\left<p\right>=0; \frac{d}{dt}\left<p\right>=-\left<V^\prime (x)\right>=0,$$ which is in accordance with the Ehrenfest theorem (albeit uninformatively).
Be careful about the time dependence of your reported $\psi$: you're actually dealing with $\Psi(x, t)=\psi(x)e^{-itE/\hbar}$ for those stationary states. Of course, this doesn't relate to the Ehrnfest theorem, but it's something worth mentioning: the complex and real parts are oscillating, as shown by the pink and blue lines in this diagram (from the wikipedia page on QHO):
Referring to the same diagram, observe parts G and H: these represent coherent states, which can be understood using Ehrnfest's theorem.
See
Kanasugi, H., and H. Okada. “Systematic Treatment of General Time-Dependent Harmonic Oscillator in Classical and Quantum Mechanics.” Progress of Theoretical Physics, vol. 93, no. 5, 1995, pp. 949–960., doi:10.1143/ptp/93.5.949.