Can you express the doublet as a sum of unit step functions, eahc of which are appropriately scaled and delayed?
If so, and you can find the output of the sytstem in response to a unit step input, then you can use the linearity and time invariance properties of the system to construct the output in response to the doublet.