If we formally exponentiate the derivative operator $\frac{d}{dx}$ on $\mathbb{R}$, we get
$$e^\frac{d}{dx} = I+\frac{d}{dx}+\frac{1}{2!}\frac{d^2}{dx^2}+\frac{1}{3!}\frac{d^3}{dx^3}+ \cdots$$
Applying this operator to a real analytic function, we have
$$\begin{align*}e^\frac{d}{dx} f(x) &= f(x)+f'(x)+\frac{1}{2!}f''(x)+\cdots\\ &=f(x)+f'(x)((x+1)-x)+\frac{1}{2!}f''(x)((x+1)-x)^2+\cdots\\ &=f(x+1) \end{align*}$$
Does anyone have an explanation of why this should "morally" be true? I do not have a very good intuition for the matrix exponential which is probably holding me back here...
You can just see it as an identity: the shift operator can be expressed in terms of a Taylor series, and then we just compute its closed form.
There are other visualizations for this, though. You can think of $1 + \frac{d}{dx}$ as an infinitesimal shift operator, and exponentiation accumulates all of the infinitesimal shifts up into an actual shift.
In particular, if $E_k$ is the shift-by-$k$ operator, and $\Delta_k = E_k - 1$ then
$$ E_1 = (1 + \Delta_{1/n})^{n} $$
but we know that for small $k$, $\Delta_k f \approx k \frac{df}{dx}$. Thus,
$$ E_1 \approx \left(1 + \frac{1}{n} \frac{d}{dx} \right)^n $$
for large $n$. In fact, both of the following equalities turn out to be true:
$$ E_1 = \lim_{n \to \infty} \left(1 + \frac{1}{n} \frac{d}{dx} \right)^n = e^{d/dx}$$