Let $X$ be a separable Hilbert space, and let $(e_i)_{i=1}^\infty$ be an orthonormal basis of $X$. For each $n\in \mathbb{N}$, let $X_n$ be the subspace spanned by $(e_i)_{i=1}^n$, and consider the operator $I_n:\mathbb{R}^n\to X$ be such that $I_n u=\sum_{i=1}^n u_i e_i$ for all $u=(u_i)_{i=1}^n$. Let $M$ be a symmetric positive semidefinite operator on $X$ (we can assume $M$ is of trace class if necessary).

Given $\lambda>0$ and $U\in X$. Is it true that $$ \lim_{n\to \infty} \langle I_n(I_n^* MI_n+\lambda \operatorname{id}_{\mathbb{R}^n})^{-1} I_n^* U, U\rangle_X = \langle ( M +\lambda \operatorname{id}_{X})^{-1} U,U\rangle_X. $$

The claim seems to be related to the convergence of projection operator in the weak operator topology. The matrix $I_n^* MI_n$ has a matrix representation $(\langle Me_i,e_j\rangle_X)_{ij}$. By using $I_n^*\operatorname{id}_{X} I_n=\operatorname{id}_{\mathbb{R}^n}$, we can rewrite $I_n^* MI_n+\lambda \operatorname{id}_{\mathbb{R}^n}$ as $I_n^* (M+\lambda\operatorname{id}_{X})I_n $. But I don't know how to handle the inverse operator, as $I_n:\mathbb{R}^n\to X$ is not invertible.