[guided]The smoothness issue is the only part not settled by fibrewise linear algebra. A map between vector bundles can be checked locally in frames: if its coordinate matrix in smooth local frames has smooth entries, then the bundle map is smooth. Here the coordinate matrices will be identity matrices.
Choose an open set $U\subset N$ on which $E$ admits a smooth local frame $s_1,\dots,s_r$, with each $s_i:U\to E$ smooth and with $s_1(y),\dots,s_r(y)$ a basis of $E_y$ for every $y\in U$. Set $V=f^{-1}(U)$. Since $f:M\to N$ is smooth, $V$ is open in $M$. Define $\widetilde{s}_i:V\to f^*E$ by $\widetilde{s}_i(x)=(x,s_i(f(x)))$.
These sections are smooth because they are obtained by composing the smooth map $f$ with the smooth frame sections $s_i$, and they form a frame because every fibre $(f^*E)_x$ is $\{x\}\times E_{f(x)}$.
Now consider the dual bundles. Let $s_1^*,\dots,s_r^*$ be the dual frame of $E^*$ over $U$, characterized by
\begin{align*}
s_i^*(y)(s_j(y))=\delta_{ij}
\end{align*}
for all $y\in U$. The corresponding pulled-back frame of $f^*(E^*)$ consists of the sections $x\mapsto (x,s_i^*(f(x)))$. The dual frame of $(f^*E)^*$ is denoted by $\widetilde{s}_1^*,\dots,\widetilde{s}_r^*$ and satisfies
\begin{align*}
\widetilde{s}_i^*(x)(\widetilde{s}_j(x))=\delta_{ij}.
\end{align*}
Using the definition of $\Phi$, for every $x\in V$ and every $1\le j\le r$,
\begin{align*}
\Phi(x,s_i^*(f(x)))(\widetilde{s}_j(x))=\Phi(x,s_i^*(f(x)))(x,s_j(f(x)))=s_i^*(f(x))(s_j(f(x)))=\delta_{ij}.
\end{align*}
Thus $\Phi(x,s_i^*(f(x)))$ is exactly $\widetilde{s}_i^*(x)$. In the chosen frames, $\Phi$ sends each basis vector to the corresponding basis vector, so its coordinate matrix is the identity matrix. The identity matrix has smooth constant entries, and its inverse is again the identity matrix. Hence $\Phi$ is a smooth vector bundle isomorphism.
The same local-frame computation proves smoothness for exterior powers. For each strictly increasing $k$-tuple $I=(i_1,\dots,i_k)$, define $s_I:U\to \Lambda^k E$ by $s_I(y)=s_{i_1}(y)\wedge\cdots\wedge s_{i_k}(y)$. The sections $s_I$ form the standard exterior-power frame of $\Lambda^k E$. The corresponding frame of $\Lambda^k(f^*E)$ is given by the sections $\widetilde{s}_I:V\to \Lambda^k(f^*E)$ defined by $\widetilde{s}_I(x)=\widetilde{s}_{i_1}(x)\wedge\cdots\wedge \widetilde{s}_{i_k}(x)$.
By the defining formula for $\Psi_k$,
\begin{align*}
\Psi_k(x,s_I(f(x)))=\Psi_k(x,s_{i_1}(f(x))\wedge\cdots\wedge s_{i_k}(f(x)))=(x,s_{i_1}(f(x)))\wedge\cdots\wedge (x,s_{i_k}(f(x)))=\widetilde{s}_{i_1}(x)\wedge\cdots\wedge \widetilde{s}_{i_k}(x)=\widetilde{s}_I(x).
\end{align*}
Again the local coordinate matrix is the identity matrix. Therefore $\Psi_k$ is smooth, and the same identity-coordinate argument applies to its inverse. For $k=0$, both exterior-power bundles are the product line bundle in these local descriptions, and the map is the identity on each fibre.[/guided]