The Schrödinger solutions for a three-dimensional central potential system whose Hamiltonian is composed of a time-dependent harmonic plus an inverse harmonic potential are investigated. Because of the time-dependence of parameters, we cannot solve the Schrödinger solutions relying only on the conventional method of separation of variables. To overcome this difficulty, special mathematical methods, which are the invariant operator method, the unitary transformation method, and the Nikiforov-Uvarov method, are used when we derive solutions of the Schrödinger equation for the system. In particular, the Nikiforov-Uvarov method with an appropriate coordinate transformation enabled us to reduce the eigenvalue equation of the invariant operator, which is a second-order differential equation, to a hypergeometric-type equation that is convenient to treat. Through this procedure, we derived exact Schrödinger solutions (wave functions) of the system. It is confirmed that the wave functions are represented in terms of time-dependent radial functions, spherical harmonics, and general time-varying global phases. Such wave functions are useful for studying various quantum properties of the system. As an example, the uncertainty relations for position and momentum are derived by taking advantage of the wave functions.