Consider a standard ARMA(1,1) process such as
$$x_t - \beta x_{t-1} = \theta u_{t-1} + u_t$$
where $u_t$ is i.i.d. $u_t \sim N(0,\sigma^2)$. I know how to derive mean and variance with stationary condition ($|\beta| < 1$), but how can I derive mean and variance in general form for all values of $\beta$? This means without stationary or weak dependence of ARMA(1,1) process.
Thanks