Codex @codex  0

The first-kind Chebyshev polynomial $T_n$ converts a multiple-angle cosine into a polynomial:The cosine addition formula gives the recurrence relationThis constructs polynomials valid for every real $x$, including outside $[-1,1]$, where the real-angle description is unavailable. For example,For $n\ge 1$, $T_n$ has polynomial degree $n$ and leading coefficient $2^{n-1}$. On $[-1,1]$ its uniform norm is one, even as its oscillations become more numerous.^[ref][ref]

The polynomial roots of $T_n$ are especially simple:They correspond to equally spaced angles, but their locations in $x$ cluster near $-1$ and $1$. For $T_3$, the polynomial roots are $\sqrt{3}/2,0,-\sqrt{3}/2$. For $T_8$, the positive polynomial roots are approximately $0.980785,0.831470,0.555570,0.195090$, with their negatives completing the set. The distance from the right endpoint to the nearest polynomial root is $1-\cos (\pi /(2n))\sim {\pi }^2/(8n^2)$; central neighboring polynomial roots are separated on the larger scale $1/n$. The small change of cosine near angles zero and $\pi$ explains the clustering.

These polynomial roots supply $n$ interpolation nodes for a polynomial interpolation of polynomial degree at most $n-1$. They differ from the $n+1$ Chebyshev-Lobatto nodes $\cos (j\pi /n)$, which include both endpoints. This choice of polynomial roots has nodal polynomialThe monic Chebyshev extremal polynomial minimizes the uniform norm of this product among monic polynomials of polynomial degree $n$ on $[-1,1]$. Consequently, if $f$ has $n$ continuous derivatives, the usual polynomial interpolation error gives the useful boundThis optimizes the nodal factor in that bound, not the actual error for every function. Endpoint clustering helps control the large endpoint oscillations of equally spaced interpolation associated with Runge's phenomenon. The two node conventions and their geometric construction are discussed in Trefethen's authored text.^[ref]

The same cosine substitution for polynomial approximation explains the weighted orthogonality:The endpoint singularities of the weight are integrable. Thus the Chebyshev polynomials are orthogonal polynomials for this weight, making coefficients in a least-squares polynomial in an orthogonal-polynomial basis separate cleanly. Their polynomial roots also give Chebyshev–Gauss quadrature, connecting polynomial interpolation to Gaussian quadrature.^[ref]

Accuracy still depends on the target function. The absolute value function $|x|$ has a corner at zero, so the hypothesis requiring continuous derivatives in the displayed polynomial interpolation error bound fails. Endpoint clustering cannot remove that corner. A jump poses a stronger obstruction: polynomials are continuous functions, so they cannot converge in the uniform norm to a discontinuous function on the whole interval. Chebfun's documented examples show persistent overshoot for such jumps; these are a limitation of global polynomial interpolation, not a universal promise of rapid accuracy from Chebyshev polynomials.^[ref]

For weighted least-squares polynomial in an orthogonal-polynomial basis with weight $(1-x^2{)}^{-1/2}$ on $[-1,1]$, the projection of $\sqrt{1-x^2}$ onto polynomials of degree at most $n$ is the displayed sum. Put $x=\cos \theta$ to obtain $T_j(x)=\cos j\theta$ and $f(x)=\sin \theta$ on $[0,\pi ]$. The constant coefficient is $2/\pi$, odd coefficients vanish, and the coefficient of $T_{2m}$ is $4/[\pi (1-4m^2)]$. For odd $n$ the optimal polynomial can have degree $n-1$.

Let $f$ be a real function with $\Vert f{\Vert }_w^2={\int }_a^{b}f(x{)}^{2}w(x)dx<\infty$, and let $V_m$ be the vector space of real polynomials of degree at most $m$. Its orthogonal projection is the least-squares polynomial in an orthogonal-polynomial basisEach coefficient depends on just one weighted inner product with $f$ and one squared norm. In particular, the coefficients $c_0,\dots ,c_m$ stay the same when the allowed polynomial degree increases. This statement assumes a fixed weight, interval and normalization of the orthogonal polynomials.

To see why the least-squares normal equations decouple, first write a trial polynomial as $q={\sum }_{j=0}^m{a}_j{\varphi }_j$ in any basis of $V_m$. The residual must be orthogonal to each permitted direction, givingThe Gram matrix usually couples the unknown coefficients. With ${\varphi }_j=p_j$, it becomes $G_{ij}=h_i{\delta }_{ij}$, leaving the separate equations $h_i{a}_i=b_i$. For the monomial basis on $[-1,1]$ with $w=1$, by contrast, $G_{ij}$ is zero when $i+j$ is odd and $2/(i+j+1)$ when $i+j$ is even. Distinct even powers, or distinct odd powers, remain coupled.

There is also a proof of the minimum that uses no differentiation. The defining coefficients give $⟨f-q_m,p_n{⟩}_w=0$ for $0\le n\le m$. Thus $f-q_m$ is orthogonal to all of $V_m$. For every $v\in V_m$, the Pythagorean theorem in an inner-product space givesThe second term is nonnegative and vanishes exactly when $v=q_m$, proving that the least-squares solution is unique. This is an integrated squared-error optimum; it does not assert a smallest error at every point or a smallest uniform norm.

For $N>m$, the nested orthogonal projections satisfy the exact residual budgetIndeed, $q_N-q_m$ is a linear combination of mutually orthogonal polynomials, and $f-q_N$ is orthogonal to it. Every newly retained component removes exactly $c_n^2{h}_n$ from the squared error. No infinite expansion or completeness assumption is needed for this finite identity. The Legendre least-squares projections of a cubic make every term explicit.

A finite-dimensional linear subspace $W$ of an inner product space has a unique closest vector to any given vector $x$. This is its orthogonal projection $P_{W}x$. If $e_1,\dots ,e_d$ is an orthogonal basis of $W$, compute it one direction at a time:The inner product convention here is linear in its first argument, including over the complex numbers. For an orthonormal basis, each denominator is one. When $W$ consists only of the zero vector, the empty sum gives $P_{W}x=0$. The finite-dimensional formula is described in Francis Narcowich, Inner Products and Norms: minimization problems.

Take $W=\mathrm{span}\{(1,2{)}^T\}$ and $x=(3,1{)}^T$, using the dot product. With $e=(1,2{)}^T$,The orthogonal projection residual has $r\cdot e=2-2=0$. To verify the best approximation directly, every candidate in the line is $te$, andThe smallest squared Euclidean norm error is therefore $5$, attained only at $t=1$. The diagram connects the perpendicular residual to this one-variable quadratic polynomial.

The general mechanism is the Pythagorean theorem in an inner-product space. For any $w\in W$, the vectors $r=x-P_{W}x$ and $P_{W}x-w$ are orthogonal, soMoving away from $P_{W}x$ inside $W$ adds a nonnegative squared norm, which is zero only when $w=P_{W}x$. This proves both optimality and uniqueness.

The useful feature is separation of preparation from evaluation. Construct an orthonormal basis once with the Gram-Schmidt process, then obtain each new best approximation using inner products instead of a fresh optimization problem. When $W$ is the column space of a matrix, this gives linear least squares; when $W$ is spanned by orthogonal polynomials, it gives a least-squares polynomial in an orthogonal-polynomial basis. The ambient inner product space need not be a Hilbert space for this finite-dimensional construction. In an infinite-dimensional target linear subspace, existence requires additional hypotheses, supplied by orthogonal decomposition by a closed subspace.

Put $x=\cos \theta$ in the weighted inner product. Since $T_r(\cos \theta )=\cos (r\theta )$,The last equality follows by the product-to-sum trigonometric identity and integration. Thus these are nonzero mutually orthogonal polynomials. The first $n+1$ are linearly independent, and the vector space $P_n$ has dimension $n+1$; therefore $\boxed{P_n=\mathrm{span}(T_0,\dots ,T_n)}$.

For any $q\in P_n$, the function of the real parameter $t$has a minimum at $t=0$. Its derivative there vanishes, proving $\boxed{(u-p_n,q)=0}$. This is precisely the defining orthogonality of an orthogonal projection, or a least-squares polynomial in an orthogonal-polynomial basis. Conversely, this condition and the Pythagorean theorem for inner product spaces give $\Vert u-q{\Vert }_2^2=\Vert u-p_n{\Vert }_2^2+\Vert p_n-q{\Vert }_2^2$, proving the minimum and its uniqueness.

Expand $p_n={\sum }_{r=0}^n{a}_r{T}_r$ in the basis already established and take the inner product with $T_s$. Orthogonality leaves $a_s(T_s,T_s)=(u,T_s)$, soIn particular $a_0={\pi }^{-1}{\int }_0^{\pi }u(\cos \theta )d\theta$ and $a_r=2{\pi }^{-1}{\int }_0^{\pi }u(\cos \theta )\cos (r\theta )d\theta$ for $r\ge 1$. These are the corresponding Chebyshev series coefficients with the constant coefficient convention used in this question.

The nonzero orthogonal polynomials $Q_0,\dots ,Q_n$ have distinct degrees and form an orthogonal basis of ${\mathcal{P}}_n$. DefineFor each $j\le n$,Thus the residual $f-p_n^{\ast }$ is orthogonal to all of ${\mathcal{P}}_n$.

For any $p\in {\mathcal{P}}_n$, writeThe two terms are orthogonal, so the Pythagorean theorem in an inner-product space givesEquality holds only for $p=p_n^{\ast }$. This proves the formula and uniqueness of the least-squares polynomial in an orthogonal-polynomial basis.

For a positive-definite weighted inner product on a compact interval, a least-squares polynomial is the orthogonal projection of a function onto the finite-dimensional space ${\Pi }_n$ of polynomials of degree at most $n$. In a basis ${\varphi }_0,\dots ,{\varphi }_n$, its coefficients solve the normal equations $Gc=d$, where $G_{ij}=({\varphi }_i,{\varphi }_j)$ and $d_i=(u,{\varphi }_i)$. The Gram matrix is positive definite: $z^{T}Gz=\Vert \sum z_i{\varphi }_i{\Vert }^2>0$ for $z\ne 0$. Thus the coefficients exist and are unique. The residual is orthogonal to every candidate polynomial, giving $\Vert u-q{\Vert }^2=\Vert u-p_n{\Vert }^2+\Vert p_n-q{\Vert }^2$ and proving optimality. An orthogonal basis makes these equations diagonal, as in a least-squares polynomial in an orthogonal-polynomial basis.