[guided]Because $\mu_Q$ and $\mu_P$ are real scalars, the scalar multiples $\mu_Q I$ and $\mu_P I$ commute with $Q$ and $P$. Therefore the centered commutator is the same as the original commutator:
\begin{align*}
[(Q-\mu_Q I),(P-\mu_P I)]=[Q,P].
\end{align*}
We now compute $[Q,P]$ on the Schwartz function $\psi$. The position-after-momentum expression is
\begin{align*}
(QP\psi)(x)=-i\hbar x\psi'(x).
\end{align*}
The momentum-after-position expression uses the product rule:
\begin{align*}
(PQ\psi)(x)=-i\hbar(\psi(x)+x\psi'(x)).
\end{align*}
Subtracting these expressions gives
\begin{align*}
([Q,P]\psi)(x)=i\hbar\psi(x).
\end{align*}
Thus $[Q,P]\psi=i\hbar\psi$. Since $\|\psi\|_{L^2}=1$, this implies
\begin{align*}
([Q,P]\psi,\psi)_{L^2}=i\hbar.
\end{align*}
The space $\mathcal S(\mathbb R)$ is invariant under multiplication by $x$, differentiation, and scalar shifts. Therefore it is invariant under $Q$, $P$, $Q-\mu_Q I$, and $P-\mu_P I$, and all compositions used in the commutator are defined on $\psi$. For symmetric operators $S$ and $T$ on a common invariant domain, with the $L^2$ inner product linear in the first variable, the identity
\begin{align*}
2i\,\operatorname{Im}(S\psi,T\psi)_{L^2}=-(([S,T]\psi,\psi)_{L^2})
\end{align*}
follows by expanding the commutator and using symmetry. With $S=Q-\mu_Q I$ and $T=P-\mu_P I$, we obtain
\begin{align*}
|\operatorname{Im}(u,v)_{L^2}|=\frac{\hbar}{2}.
\end{align*}
The previous step already proved
\begin{align*}
\Delta_\psi Q\,\Delta_\psi P\geq |\operatorname{Im}(u,v)_{L^2}|.
\end{align*}
Combining the two displays gives
\begin{align*}
\Delta_\psi Q\,\Delta_\psi P\geq \frac{\hbar}{2}.
\end{align*}[/guided]