[guided]Fix one index $j\in\{1,\dots,m\}$. The point is to reduce the statement for the function $f_j$ to the ordinary one-dimensional Strong Law of Large Numbers. For every $i\in\mathbb N$, define the real-valued random variable
\begin{align*}
Y_{i,j}:\Omega\to\mathbb R,\qquad Y_{i,j}(\omega):=f_j(X_i(\omega)).
\end{align*}
This is measurable because it is the composition of the measurable map $X_i:(\Omega,\mathcal A)\to(S,\mathcal S)$ with the measurable map $f_j:(S,\mathcal S)\to(\mathbb R,\mathcal B(\mathbb R))$.
We must also verify integrability before applying the Strong Law of Large Numbers. Since $X_i$ has distribution $P$, the change-of-law identity for expectations gives
\begin{align*}
\mathbb E[|Y_{i,j}|]=\mathbb E[|f_j(X_i)|]=\int_S |f_j(x)|\,dP(x)=P|f_j|<\infty.
\end{align*}
The sequence $(Y_{i,j})_{i\ge 1}$ is i.i.d.: independence follows because the variables $X_i$ are independent and each $Y_{i,j}$ is a [measurable function](/page/Measurable%20Function) of the corresponding $X_i$, and identical distribution follows because all $X_i$ have distribution $P$.
Therefore the Strong Law of Large Numbers (citing a result not yet in the wiki: Strong Law of Large Numbers) applies to the integrable i.i.d. sequence $(Y_{i,j})_{i\ge 1}$. Hence there exists an event $A_j\in\mathcal A$ with $\mathbb P(A_j)=1$ such that, for every $\omega\in A_j$,
\begin{align*}
\frac{1}{n}\sum_{i=1}^{n}Y_{i,j}(\omega)\to \mathbb E[Y_{1,j}].
\end{align*}
Substituting the definition of $Y_{i,j}$ gives
\begin{align*}
\frac{1}{n}\sum_{i=1}^{n}Y_{i,j}(\omega)=\frac{1}{n}\sum_{i=1}^{n}f_j(X_i(\omega))=P_n f_j(\omega).
\end{align*}
Since $X_1$ has distribution $P$,
\begin{align*}
\mathbb E[Y_{1,j}]=\mathbb E[f_j(X_1)]=\int_S f_j(x)\,dP(x)=P f_j.
\end{align*}
Thus, on the event $A_j$, the empirical average $P_n f_j$ converges to the population average $P f_j$.[/guided]