Harmonic-Enriched Reproducing Kernel Approximation for Highly Oscillatory Differential Equations