I'd appreciate your advice on how to solve recurrence relations of the following kind using generating functions: general solution of:
\begin{cases}a_{n+1} &= 5a_n - 3b_n \\ b_{n+1} &= 4a_n - 2b_n\end{cases}
for $n>0$ and particular solution for $a_0 = 2, b_0 = 1$.
I know how to solve recursions where only one series, f.e. $(a_n)_n$ is involved, but not how to do this for $(a_n)_n$ and $(b_n)_n$ at once.
Akin to the matrix approach, simply eliminate one of the variables from the equation. For example,
$$3b_{n+1}=2a_{n+1}+2a_n$$
which can then be substituted into the definition of $a_{n+1}$ to get
$$a_{n+1}=3a_n-2a_{n-1}$$
which is now linear in one recurrence. Similarly, we have
$$4a_{n+1}=5b_{n+1}-2b_n$$
and hence,
$$b_{n+1}=3b_n-2b_{n-1}$$