The generating functional, as always, is given by
$$ Z[J, \eta, \eta^{\dagger}] = \int \mathcal{D}\phi \mathcal{D} \psi \mathcal{D} \psi^{\dagger} \exp \left[ \frac{i}{\hbar} \int d^4 x \, \left( \mathcal{L}(\phi, \psi, \psi^{\dagger}) +J\phi + \eta^{\dagger}\psi + \psi^{\dagger}\eta\right)\right] = $$
$$ \exp \left[ -\frac{i g}{\hbar} \int d^4 x \, \frac{\delta}{\delta J} \frac{\delta}{\delta \eta} \frac{\delta}{\delta \eta^{\dagger}} \right] \; Z_0[J, \eta, \eta^{\dagger}], $$
where $Z_0$ is the generating functional of the free theory:
$$ Z_0[J, \eta, \eta^{\dagger}] = \hbar \int d^4 x \, \int d^4 y \, \left( \frac{1}{2} J(x) \Delta_M(x, y) J(y) + \eta^{\dagger}(x) \Delta_m(x, y) \eta(y) \right) $$
with $\Delta_m$ the Klein-Gordon propagator with mass $m$.
Your formula for correlations (aka Green's functions of the interacting theory) is correct. You have to Tailor-expand the exponential of the interacting term to the appropriate order in perturbation theory.
Just do the math and you will arrive at the correct expressions for correlation functions. Spoiler alert: the second one vanishes.
Btw you don't need to explicitly write the generating functional to calculate correlations. There's a simpler way: first, note that you can calculate expressions like
$$ \left< F[\phi, \psi, \psi^{\dagger}] \right>_0 = \mathcal{N}_0 \int\mathcal{D}\phi \mathcal{D}\psi \mathcal{D}\psi^{\dagger} \exp\left[ \frac{i}{\hbar} \int d^4 x \, \mathcal{L}_0(\phi, \psi, \psi^{\dagger}) \right] \cdot F[\phi, \psi, \psi^{\dagger}] $$
where $F$ is polynomial in fields with help of Wick's theorem. Also note the diagrammatic interpretation of the terms in the Wick expansion.
Then define
$$ \left< F[\phi, \psi, \psi^{\dagger}] \right> = \mathcal{N} \int\mathcal{D}\phi \mathcal{D}\psi \mathcal{D}\psi^{\dagger} \exp\left[ \frac{i}{\hbar} \int d^4 x \, \mathcal{L}(\phi, \psi, \psi^{\dagger}) \right] \cdot F[\phi, \psi, \psi^{\dagger}] = $$
$$ \frac{\mathcal{N}}{\mathcal{N}_0} \left< F \cdot \exp \left[ - \frac{i g}{\hbar} \int d^4 x \, \phi \psi^{\dagger} \psi \right] \right>_0. $$
and Tailor-expand the exponential of the interacting term to any order (given in advance). You will arrive at the polynomial expression (because both $F$ and the truncated Tailor series are polynomials in $\phi$, $\psi$ and $\psi^{\dagger}$). We already know how to compute these:
The expectation bracket, to each order in the perturbative series, is given by a sum of terms. Each term can be represented graphically as a Feynman diagram.
The expectation brackets are by definition normalized such that
$$ \left<1\right>_0 = \left<1\right> = 1, $$
which means that $\mathcal{N}_0$ and $\mathcal{N}$ are related to each other. There's a very generic result, which is: the ratio $\mathcal{N} /
\mathcal{N_0}$ corresponds to the product of all bubble graphs (the ones without external legs). These nicely factor out, giving a convenient way to calculate correlations:
The appropriate normalization can be accounted for by simply disregarding graphs with disconnected bubbles.