System Resonse's Autocorrelation Function (using integral)

[RoshanBBQ]

Homework Statement

x(t) is wss.
$$R_x(\tau) = e^{-a|\tau|}$$
The process goes through a system with the impulse response
$$h(t) = u(t)e^{-bt}$$
What is the output's autocorrelation function?

The Attempt at a Solution

For starters, I already found a solution by finding the output power spectral density and finding its inverse Fourier transform. But I cannot make my integration method work. I want someone to tell me where my integration is going wrong:

$$R_y(\tau)=\int_{-\infty}^\infty \int_{-\infty}^\infty h(s)h(r)R_x(\tau +s -r)\,dr \,ds$$
$$R_y(\tau)=\int_{-\infty}^\infty \int_{-\infty}^\infty u(s)e^{-bs}u(r)e^{-br}e^{-a|\tau +s -r|}\,dr \,ds$$

So I separated the integral. Let us just ignore the second integral for now and focus on my method for the first half. If the first half is working, the second will work. So I write for one integral:
$$\tau +s -r > 0 \rightarrow \tau + s > r$$
So the first integral becomes, now eliminating the u(r) by incorporating it into the limits:
$$R_y(\tau)=\int_{-\infty}^\infty \int_{0}^{\tau + s} u(s)e^{-bs}e^{-br}e^{-a(\tau +s -r)}\,dr \,ds$$

But we then realize unless the below condition is true, the integration adds up to zero due to the unit step:
$$\tau + s > 0 \rightarrow s > -\tau$$

We can use this now as a limit of integration:

$$R_y(\tau)=\int_{-\tau}^\infty \int_{0}^{\tau + s} u(s)e^{-bs}e^{-br}e^{-a(\tau +s -r)}\,dr \,ds$$

But, due to u(s), we know the following condition must also be true or the answer is zero:
$$-\tau > 0$$
EDIT: I think the above is not true. If -tau < 0, the integration then goes from 0 to infinity. Is this the root of the error? If so, how do I handle this? uhh, I guess I'm just stuck.So I integragted (assuming this solution was for negative tau)
$$R_y(\tau)=\int_{-\tau}^\infty \int_{0}^{\tau + s} e^{-bs}e^{-br}e^{-a(\tau +s -r)}\,dr \,ds$$

and got (for tau < 0)
$$\frac{e^{b \tau}}{2(a+b)b}$$

So my total answer will be
$$\frac{e^{-b |\tau|}}{2(a+b)b} \ne \frac{1}{a^2-b^2} \left ( \frac{a}{b} e^{-b|\tau|} - e^{-a |\tau|} \right )$$Note, the answer it is not equal to is the one I got using power spectral density.

 

[RoshanBBQ]

OMG I GOT IT TO WORK USING THE REALIZATION I HAD IN MY EDIT!

So let's say I have the symbolic solution to that integral (after the second one) in terms of s. I will call it

x(s).

Let's say I set up a similar one for the other half and get that solution, calling it xx(r).

Then the solution for when tau < 0 is actually:
-x(tau) - xx(0)
(because x(inf) = 0 and xx(inf) = 0 -- I am just adding the two integrals, one from -tau to inf and the other from 0 to inf)

And this results in the same answer as from the power spectral density!