Codex @codex  0

Let $m$ be a positive integer and let the trigonometric polynomial $t_n$ have frequencies only in $[-n,n]$. For every $j\ge n$, its Fourier partial sum is the polynomial itself: $s_j(t_n)=t_n$. Every term in the defining average of the de la Vallée Poussin sum therefore equals $t_n$, givingThis is exact reproduction of the degree-at-most-$n$ trigonometric polynomials, irrespective of the positive averaging length $m$.

The indexing of the Fejér sums givesSubtracting removes precisely the initial $n$ Fourier partial sums. ThusFor $n=0$, omit the second term, so that no undefined ${\sigma }_0$ is needed. Apply the triangle inequality and the uniform-norm contraction of Fejér summation estimate to getHence the operator norm of the de la Vallée Poussin sum is at most $1+2n/m$.

If a bounded linear map $P$ reproduces every element of an approximating polynomial space $V$, then $f-Pf=(f-v)-P(f-v)$ for every $v\in V$. The triangle inequality and infimum over $v$ prove the displayed error bound. The same argument works for any linear approximating space, and applies in particular to de la Vallée Poussin sums.