I assume that you are workinig with the definition
\begin{align} L_X(\mu) = \left.\frac{d}{dt}\right\vert_{t=0}(\phi_t^X)^*\mu, \end{align}
where $\phi_t^X\in\mathrm{Diff}(M)$ is the flow generated by $X\in\mathcal T_M$, and $\mu\in\Gamma(M,E)$ is a section of a vector bundle $E\to M$ canonically associated with $M$ (e.g., the Berezinian bundle). Note here that $t$ should really be thought of as some parameter from the base ring $\mathcal O _S$ via the functor-of-points picture (i.e., we parametrise $M$ via a smooth submersion of supermanifolds $M\to S$).
The following two lemmas are useful for the problem at hand.
Lemma 1: Let $U\subset M$ such that there exists a coordinate chart $(x^\bullet):U\to U^{m|n}\subset\mathbb R^{m|n}$, and let $[-] : \left((\Omega^1)^{m|n}\right)^\times \to \mathcal B er(\Omega^1)$ be the Berezinian of a local basis of $\Omega^1$. It holds
$$L_{\partial/\partial x^i} [\mathrm{d} x^1 \cdots \mathrm{d} x^{m+n}] = 0.$$
Proof: Note that by definition of the Berezinian sheaf,
$$ \phi^* [\mathrm d x^\bullet] = ber(D\phi)[\mathrm d x^\bullet]$$
for all diffeomorphisms $\phi\in\mathrm{Diff}(M)$.
Now, to compute the Lie derivative, we need to compute the flow of the coordinate vector fields
$$
x^j \circ \phi^{\partial / \partial x^i}_t = x^j + \delta^{ij} t,
$$
or, in other words,
$$\phi^{\partial/\partial x^i}_t \begin{pmatrix}x^1\\ \vdots \\ x^{i-1} \\ x^i\\ x^{i+1}\\ \vdots \end{pmatrix} =
\begin{pmatrix}x^1\\ \vdots \\ x^{i-1} \\ x^i + t\\ x^{i+1}\\ \vdots \end{pmatrix}.$$
Here, $t$ should really be thought of as a parameter from the base via the funtor-of-points picture: In particular, whenever $x^i$ is an odd coordinate function, $t$ should also be an odd parameter.
In either case, for all $t\in\mathcal O _S$, the Jacobian is given by $D\phi^{\partial/\partial x^i}_t = \mathrm{id}$, such that
\begin{align}
L_{\partial/\partial x^i} [\mathrm d x^\bullet]
&= \left.\frac{d}{dt}\right\vert_{t=0}
\left(\phi^{\partial/\partial x^i}_t\right)^*[\mathrm d x^\bullet]\\
&= \left.\frac{d}{dt}\right\vert_{t=0}
ber(D\phi^{\partial/\partial x^i}_t ) [\mathrm d x^\bullet]\\
&= \left.\frac{d}{dt}\right\vert_{t=0}
ber(\mathrm{id}) [\mathrm d x^\bullet]\\
&=0,
\end{align}
which proves the claim.
Lemma 2: Let $E\to M$ be a vector bundle canonically associated with $TM$, and denote $\mathcal E = \Gamma(M,E)$. For all $\mu\in\mathcal E$, $X\in\mathcal T_M$ and $f\in\mathcal O_M$, it holds
$$ L_{fX}(\mu) = (-1)^{p(X)p(f)}L_X(f\mu) = (-1)^{p(X)p(f)}X(f)\mu + fL_X(\mu).$$
Proof: This is just a reincarnation of the Leibniz rule on $\mathbb R$ (or more precisely, the base ring $\mathcal O _S$), together with $L_Xf = X(f)$. Detailed proofs can be found in Deligne & Morgan's Notes on supersymmetry.
The solution to the problem at hand now simply follows from $$\frac{\partial \psi^2}{\partial x} = \frac{\partial f}{\partial \psi^1}=0$$
together with the two lemmas above.