Suppose $X$ and $Y$ are normally distributed random variables.
I know that if $X$ and $Y$ are independent, then $(X,Y)$ is jointly normal.
But in general case, $X$ and $Y$ are not independent, $(X,Y)$ is not necessary jointly normal.
I wonder if there is a vector of affine-transformed random vectors, for example
$$Z_1 = a_1X+b_1Y$$
$$Z_2 = a_2X +b_2Y$$
$Z_1$ and $Z_2$ are surely normal, and also dependent.
Is $(Z_1, Z_2)$ jointly normal? If yes, how to prove it?