Analytical solution of initial value problem for the fractional oscillation equation with two Caputo fractional derivatives
, where the coefficients and orders satisfy
and
, is investigated by using the Laplace transform and complex inverse integral method on the principal Riemann surface. It is proved by using the argument principle that the characteristic equation has a pair of conjugated simple complex roots with a negative real part on the principal Riemann surface under the assumption that
and
are not both integers. Then three fundamental solutions, the unit impulse response, the unit initial displacement response, and the unit initial rate response, are derived analytically. Each of these solutions is expressed into a superposition of a classical damped oscillation decaying exponentially and a real Laplace integration decaying in a negative power law. Finally, the asymptotic behaviors of these analytical solutions for sufficiently large
are determined as monotonous decays in a power of negative exponent.