An Ornstein Uhlenbeck process $x_t$ satisfies the following stochastic differential equation:
$$dx_t = \theta (\mu - x_t ) dt + \sigma dW_t $$
where $\theta,\mu, \sigma > 0$, and $W_t$ denotes the Wiener process. When $\mu, \sigma$ and $\theta$ are constant, the analytic solution is:
$$x_t = x_0 e^{-\theta t} + \mu (1 - e^{-\theta t}) + \displaystyle\int_{0}^{t} \sigma e^{\theta(s-t)} dW_s$$
and the expectation is given by:
$$E(x_t) = x_o e^{-\theta t} + \mu (1 - e^{-\theta t})$$
However, I want to solve this SDE, and find the expectation, when $\mu = \mu(t)$, i.e. $\mu$ is a function of $t$.
Does anyone know of a solution or a reference for where a solution may be found? Thanks!