Analytical solution to a Linear advection-reaction  PDE I am looking for an analytical solution for the linear PDE
$(1)\qquad\qquad \qquad f_t+ A f_x + B f = 0,  $
Where $A$ and $B$ are constant matrices and $f=f(x,t)$ is a vector. 
Clearly each one of $f_t + A f_x = 0,$ and $f_t+ B f = 0,$ has an analytical solution involving eigenvectors and eigenvalues of $A$ and $B$ separately, but I'm not sure about equation (1).
My question is:

What are the minimal conditions on $A$ and $B$ such that there is an analytical solution to the problem in (1), (clearly if they share eigenvectors, I can write a solution...) and what is the solution in that case?

 A: By analytical, I presume that you mean explicit, or in close form. The known case so far is when $A$ and $B$ can be diagonalized in the same basis. Notice that the case where some eigenvalues come as complex conjugate pairs is subtle.
To see the importance of the condition $[A,B]=0_n$, let us make the Fourier transform of your equation:
$$\partial_t\hat f=(B-i\xi A)\hat f(t,\xi).$$
Even if you can compute the exponential of $t(B-i\xi A)$ (unlikely), it will be a complicated function of $\xi$, whose backward Fourier transform is not explicit at all. When $A$ and $B$ commutte, we have
$$\exp t(B-i\xi A)=e^{tB}e^{it\xi A},$$
and this explains that everything works well in this case. A less obvious good situation is when the commutator $[A,B]$ is a linear combination of $A$ and $B$. This argument was used (in an infinite dimensional setting) to compute explicitly the solution of the Cauchy problem for the harmonic oscillator:
$$\partial_tu-\Delta u+|x|^2u=0.$$
If you accept a less explicit formula, and if you assume that $A$ is diagonalisable with real eigenvalues (we say that the operator $\partial_t+A\partial_x$ is hyperbolic), then the Cauchy problem is well posed. You can construct the solution by using the Laplace transorm in time and Fourier transform in space. You obtain an integral formula involving the inverses $(B-i\xi A)^{-1}$.
A: Use the ansatz $f(x,t) = exp(a x) exp(c t) v$ where v is a constant vector.  You have a solution whenever $(c + a A + B) v = 0$, i.e. a nontrivial solution if v is an eigenvector 
of $ a A + B $ for eigenvalue -c.   
A: We have been working on a special case of your system (1) recently in an economic problem. Special case means that we have a two-dimensional system and A has a special structure. We obtain examples for closed form solutions whose structure resembles gamma-distributions, i.e. the structure is $f(x,t) = (c_0 + c_1x+c_2t)exp(c_3x+c_4t)$. For more details, see our pdf-file, also available at www.waelde.com/pub. We would expect that there are many generalizations of this solution. We also believe that it would work for more general structures of A than we used so far. (In fact, any hint to the mathematical literature on generalizations of this solution would be very welcome.)
