Parseval's theorem says that:
If $f(x)\in\mathcal{L}_2(-\infty,\infty)$
$\int_{-\infty}^{\infty} dx |f(x)|^2=\frac{1}{2\pi}\int_{-\infty}^{\infty} dk |\tilde f(k)|^2$
where $\tilde f(k)$ is the Fourier transform of $f(x)$.
I looked up Parseval's theorem on Wikipedia, but there the proof is not given. Some articles just say that it proves that total energy is same in position space and momentum space (i.e. the physics explanation).
Could some provide any sources or suggestions about the proof of the theorem?
I would say it usually comes with the 'package' that shows that Fourier transforms extend to $L^2$ functions.
One starting point is to consider smooth functions and knowing one particular Fourier transform. For $\sigma>0$, let $$ \phi_\sigma (t) = \exp( -\sigma^2 t^2 /2) $$ Then somehow (you would typically need some complex analysis to do so) show that: $$ \widehat{\phi}_\sigma(x) = \frac{1}{\sqrt{2\pi} \sigma} \exp(-x^2/ (2\sigma^2))$$ When $f,g$ are Schwarz functions (smooth, rapid decays), the triple integral: $$ \int\int\int e^{it(x-y)} \phi_\sigma(t) f(x)g(y) \,dt \,dx\, dy $$ converges absolutely and Fubini applies. Carrying out the $t$ integration you get $$ \int \int \widehat{\phi}_\sigma(x-y) f(x)g(y)\,dx \,dy$$ and you show that carrying out the $x$ integral and letting $\sigma\rightarrow 0$ this converges to $\int f(y) g(y) \, dy$. While if you first carry out the $x$ and $y$ integrals you get: $$ \int \int \phi_\sigma(t) \hat{f}(t)\hat{g}(-t) \, dt$$ which converges to $\int \hat{f}(t) \hat{g}(-t) \, dt$. This looks nicer if you introduce the complex $L^2$ scalar product $(f,g) = \int f(x)\overline{g(x)} dx$ as it takes the simple form of an isometry. Writing ${\cal F}$ for the Fourier transform (including a factor $1/\sqrt{2\pi}$): $$ ({\cal F} f, {\cal F} g) = (f,g) $$ This shows that ${\cal F}$ is uniformly continuous on Schwarz functions and as they are dense in $L^2$, the identity extends by continuity to all $L^2$ functions (with a more abstract interpretation).
From the above identity all the standard properties of ${\cal F}$ follows: Parseval/Plancherel (let $f=g$) as well as the Fourier inversion formulas, noting e.g. that $ ({\cal F} f, {\cal F} g) = ({\cal F}^* {\cal F}f, g) = (f,g) $