Our aim is the development and analysis of numerical schemes for approximately solving backward diffusion-wave problem, which involves a fractional derivative in time with order $\alpha\in(1,2)$. From terminal observations at two levels, i.e., $u(T_1)$ $u(T_2)$, we simultaneously recover initial data $u(0)$ $u_t(0)$ hence solution $u(t)$ all $t > 0$. First, existence, uniqueness, Lipschitz stab...