We construct, for any given ℓ=12+N0, the second-order, linear partial differential equations (PDEs) which are invariant under the centrally extended conformal Galilei algebra. At the given ℓ, two invariant equations in one time and ℓ+12 space coordinates are obtained. The first equation possesses a continuum spectrum and generalizes the free Schrödinger equation (recovered for ℓ=12) in 1 + 1 dimension. The second equation (the “ℓ-oscillator”) possesses a discrete, positive spectrum. It generalizes the 1 + 1-dimensional harmonic oscillator (recovered for ℓ=12). The spectrum of the ℓ-oscillator, derived from a specific osp(1|2ℓ + 1) h.w.r., is explicitly presented. The two sets of invariant PDEs are determined by imposing (representation-dependent) on-shell invariant conditions both for degree 1 operators (those with continuum spectrum) and for degree 0 operators (those with discrete spectrum). The on-shell condition is better understood by enlarging the conformal Galilei algebras with the addition of certain second-order differential operators. Two compatible structures (the algebra/superalgebra duality) are defined for the enlarged set of operators.