We investigate a class of exactly solvable quantum quench protocols with a finite quench rate in systems of one dimensional non-relativistic fermions in external harmonic oscillator or inverted harmonic oscillator potentials, with time dependent masses and frequencies. These hamiltonians arise, respectively, in harmonic traps, and the $c=1$ Matrix Model description of two dimensional string theory with time dependent string coupling. We show how the dynamics is determined by a single function of time which satisfies a generalized Ermakov-Pinney equation. The quench protocols we consider asymptote to constant masses and frequencies at early times, and cross or approach a gapless potential. In a right side up harmonic oscillator potential we determine the scaling behavior of the one point function and the entanglement entropy of a subregion by obtaining analytic approximations to the exact answers. The results are consistent with Kibble-Zurek scaling for slow quenches and with perturbation calculations for fast quenches. For cis-critical quench protocols the entanglement entropy oscillates at late times around its initial value. For end-critical protocols the entanglement entropy monotonically goes to zero inversely with time, reflecting the spread of fermions over the entire line. For the inverted harmonic oscillator potential, the dual collective field description is a scalar field in a time dependent metric and dilaton background.