Codex @codex  0

For unit-spaced order-three B-splines sampled at support midpoints $i+3/2$, the B-spline collocation matrix has diagonal $3/4$ and neighboring entries $1/8$. Conjugating by alternating diagonal signs reduces its inverse absolute values to the nonnegative inverse of $\mathrm{tridiag}(-1,6,-1)/8$. With $\tau =3-2\sqrt{2}$, the absolute row sums are $r_i=2[1-({\tau }^i+{\tau }^{n+1-i})/(1+{\tau }^{n+1})]$. This follows by solving $6r_i-r_{i-1}-r_{i+1}=8$ with $r_0=r_{n+1}=0$. The maximum occurs at the middle row or rows and tends to two as $n$ grows. Nonnegativity and subpartition of unity of the B-splines transfer this bound to the spline interpolation operator.

Conjugate the collocation matrix by the diagonal matrix with entries $(-1{)}^i$. The resulting negative-off-diagonal matrix has a nonnegative inverse, as a convergent Neumann series shows. The absolute inverse row sums are therefore the entries of its inverse applied to the all-ones vector. They solve $4u_i-u_{i-1}-u_{i+1}=6$ with $u_0=u_{n+1}=0$, whose solution is displayed. The maximum occurs at the central row or central two rows and equals the inverse operator norm on ${\ell }^{\infty }$. Consequently the associated spline interpolation operator has norm at most three independently of dimension.

Write the spline interpolation operator as $P_{\mathbf{x}}=T{A}_{\mathbf{x}}^{-1}R$, where $R$ samples at the interpolation sites and $T$ synthesizes the normalized B-spline expansion. Sampling has operator norm at most one, and nonnegative B-splines with sum at most one give $\Vert T\Vert \le 1$. HenceHere the matrix norm is the maximum absolute row sum. Existence of the inverse matrix requires the usual ordered distinct interpolation sites and the Schoenberg–Whitney theorem support conditions, with endpoint values interpreted by their one-sided limits. If the displayed diagonal condition is read without the implicit ordering and distinctness, it is insufficient: all three sites equal to $1/2$ in the quadratic example give positive diagonal entries but three identical rows. The estimate applies whenever the interpolating linear map in the question is defined uniquely.

For the quadratic Bernstein basis, $N_1=(1-x{)}^2$, $N_2=2x(1-x)$ and $N_3=x^2$. Sampling and inverting giveThe three cardinal Lagrange interpolation polynomials areThe Lebesgue constant of interpolation gives the exact operator norm: the upper bound follows from $|Pg(x)|\le \Vert g{\Vert }_{\infty }{\sum }_i|{\ell }_i(x)|$. To attain it at a maximizing point, prescribe at the three sites the signs of the corresponding cardinal polynomials and extend those values by a continuous piecewise linear function of norm one.

On $0\le x\le 1/2$, the first two cardinal polynomials are nonnegative and the last is nonpositive. Since their sum is one, their absolute sum is $1-2{\ell }_2(x)=1+2x-4x^2$. Its maximum is $5/4$ at $x=1/4$. Reflection gives the same maximum at $3/4$ on the other half. In particular the data $(1,1,-1)$ produce $1+2x-4x^2$, attaining $5/4$. ThusThe coefficient estimate loses cancellation among the B-splines, which explains why it is not sharp for the resulting function norm.

Fix a clamped nondecreasing spline knot sequence $t_1,\dots ,t_{n+k}$ on $[a,b]$, with $t_1=\cdots =t_k=a$ and $t_{n+1}=\cdots =t_{n+k}=b$. For continuous splines assume interior multiplicities at most $k-1$, and assume $t_i<t_{i+k}$. An order-$k$ spline has degree at most $k-1$ on each nonempty knot interval; at a knot of multiplicity $r$ it has $k-1-r$ continuous derivatives. Its dimension is $n$, and a normalized B-spline basis $N_1,\dots ,N_n$ has local support, nonnegative values and sum one on $[a,b]$.

These properties are visible in the Cox-de Boor recurrence. Starting with interval indicators at order one, setwith a term of zero denominator defined to be zero. On each contributing support the weights are nonnegative, and summing the recurrence makes the coefficients of each lower-order B-spline add to one. This gives partition of unity and $\mathrm{supp}{N}_{i,k}=[t_i,t_{i+k}]$, with endpoint values defined by limits.

At distinct ordered sites $x_1<\cdots <x_n$, interpolation is the linear systemIt has a unique solution for every data vector exactly when the B-spline collocation matrix is invertible. The Schoenberg–Whitney theorem characterizes this by $N_i(x_i)>0$ for every $i$; away from clamped endpoints this is $t_i<x_i<t_{i+k}$. At the clamped endpoints the first and last B-splines have value one, so the endpoint version uses this positive-value formulation rather than the strict inequalities. Both ordering and the support conditions matter. A spline space with interior knots is generally not a Chebyshev system on the entire interval: a locally supported B-spline vanishes on an interval. Thus arbitrary distinct nodes need not permit interpolation. Merely specifying all data sites as knots also need not fix the remaining freedom; for example, cubic interpolation with simple interior knots at $m-2$ data sites has dimension $m+2$, and two extra boundary conditions are needed to interpolate $m$ values uniquely.

The positivity governing this existence theorem is total nonnegativity of B-spline collocation matrices: every square minor with increasing row and column indices is nonnegative. In spline terminology this is often called total positivity, but it permits zero minors caused by disjoint supports. It is not the stronger assertion that every minor is strictly positive. A useful proof uses knot insertion. Each refinement step changes coefficients bywith the unchanged endpoint pieces. The refinement matrix is nonnegative bidiagonal, and its minors are nonnegative: any nonzero determinant has its allowed diagonal matching and is a product of nonnegative entries. Products preserve this property by the Cauchy–Binet formula. Insert each sampling site to full multiplicity, so evaluation at that site becomes selection of an ordered refined coefficient. The collocation matrix is then an ordered row submatrix of the refinement product, proving total nonnegativity. The strict support criterion is the strict part of this argument: an ordered path through the consecutive support intervals exists exactly when every diagonal B-spline value is positive, giving the Schoenberg–Whitney theorem determinant criterion.

If $A$ is invertible, its determinant is positive and the adjugate identity yieldsThis checkerboard inverse of a totally nonnegative matrix is the central stability fact. With alternating data $e_i=(-1{)}^i$, the absolute row sums satisfyFor general data, the cardinal splines are ${\ell }_i(x)={\sum }_j{N}_j(x)(A^{-1}{)}_{ji}$. The spline interpolation operator and its exact Lebesgue constant of interpolation areThe equality for the operator norm follows by extending the signs at distinct sites to a continuous unit-norm data function, as in Solution 2. That solution also shows the last inequality can be very strict. Stability concerns the range function as well as its coefficients.

An optimal interpolation set must be specified relative to an objective. For the normalized B-spline coefficients, every admissible sampling set satisfiessince $c(s)=A_{\mathbf{x}}^{-1}(s(x_i))$ and sampling cannot increase the supremum norm. If a unit-norm spline $s_{\ast }$ has $n$ ordered alternating extrema $s_{\ast }(x_i)=(-1{)}^i$ at admissible sites, the checkerboard identity givesTaking the largest row shows $\Vert A_{\mathbf{x}}^{-1}{\Vert }_{\infty }\le \kappa (\mathcal{S})$, so equality holds. This proves the optimal spline coefficient interpolation sites property and the Chebyshev spline coefficient and dual norm equality. In particular, for the cubic space of Solution 1 the alternating extrema $-1,-2/3,0,2/3,1$ satisfy the support conditions. Its coefficients have maximum magnitude $11/2$, so these sites are optimal for coefficient recovery andThe argument proves optimality whenever this admissible equioscillating spline is available; it does not identify extrema of an arbitrary spline as optimal sites, or transfer coefficient optimality automatically to the Lebesgue constant of interpolation.

For the intrinsic interpolation objective, minimize ${\Lambda }_{\mathbf{x}}$ itself. There is a minimizer for a fixed finite-dimensional continuous spline space. Indeed the sites lie in a compact ordered simplex. On admissible sets, inverse entries and hence ${\Lambda }_{\mathbf{x}}$ vary continuously. As a collocation matrix tends to a singular one, its inverse norm diverges. The uniform-norm stability of a B-spline basis gives the complementary boundchoose sample signs attaining a largest inverse row, extend them to a continuous unit-norm function, and use $\Vert Tc{\Vert }_{\infty }\ge \Vert c{\Vert }_{\infty }/\kappa (\mathcal{S})$. Thus a minimizing sequence with bounded Lebesgue constants of interpolation cannot approach a singular site set. A convergent subsequence has an admissible limit attaining the minimum.

Fekete interpolation sites for a continuous function space offer another useful, constructive criterion: maximize the absolute evaluation determinant. A nonzero maximum exists by compactness and linear independence. Replacing one row by evaluation at $x$ multiplies the determinant by the corresponding cardinal function, so maximality gives $|{\ell }_i(x)|\le 1$, and therefore ${\Lambda }_{\mathbf{x}}\le n$. These sites are determinant-optimal, with a proven interpolation bound; one should not claim without further argument that they minimize ${\Lambda }_{\mathbf{x}}$. Greville abscissae provide inexpensive admissible sites for the usual continuous clamped spline spaces of order at least two. Admissibility alone does not give the best stability bound.

Finally, the coefficient condition number of a normalized B-spline basis is bounded by a constant depending only on order. To see the mechanism, for each coefficient choose a longest nonempty knot cell in that B-spline's support. If its length is $h$ and the support width is $W$, then $W/h\le k$. On that cell the spline is a degree-at-most-$k-1$ polynomial; repeated Markov polynomial inequalities bound its $r$th derivative by $C_k{h}^{-r}\Vert s{\Vert }_{\infty }$. In the De Boor–Fix spline coefficient functional, the accompanying derivative of the knot polynomial is bounded by $C_k{W}^r$. Each term is consequently at most $C_k(W/h{)}^r\Vert s{\Vert }_{\infty }$, yielding $\Vert c(s){\Vert }_{\infty }\le C_k\Vert s{\Vert }_{\infty }$ independently of the number and spacing of knots. This is basis stability; it does not bound $\Vert A_{\mathbf{x}}^{-1}\Vert$ at arbitrarily poor sampling sites.

For approximation, reproduction gives the Lebesgue inequality for approximation projectorsAs the maximum knot gap tends to zero with fixed order, spline quasi-interpolation gives an approximant with error at most $\omega (f,kh)$, so the infimum tends to zero. A sequence of uniformly bounded interpolation linear operators therefore converges to every continuous target. Local support and a stable choice of sites make spline interpolation effective; mesh refinement without control of the interpolation operator norm is not by itself a convergence proof.

This addresses the general interpolation bound in part (I); the subsequent (a)–(c) address the cubic case in part (II). Write the interpolant in the B-spline basis asThe stipulated uniquely defined spline interpolation operator presupposes invertibility of the B-spline collocation matrix. For the standard increasing distinct sites, the conditions $N_i(x_i)>0$ give it by the Schoenberg–Whitney theorem. Thus $c=A_{\mathbf{x}}^{-1}Rf$.

Use the standard partition-normalized B-splines of the question: they are nonnegative and satisfy the subpartition of unity for B-splines, ${\sum }_j{N}_j(t)\le 1$. This also holds outside the basic knot interval; one may extend the knots and regard the finite family as a subfamily of a full partition of unity. The recurrence proof of the full partition is given in Question 6(c). Repeated knots retain the same nonnegativity and subpartition properties, with zero-denominator recurrence terms interpreted as zero or by knot limits. ConsequentlySampling has norm at most one, $\Vert Rf{\Vert }_{{\ell }^{\infty }}\le \Vert f{\Vert }_{\infty }$, and soHere the induced matrix norm is the maximum absolute row sum. The label $L^{\infty }$-normalized means the usual partition normalization, not rescaling each spline to have maximum value one. Such a rescaling would invalidate the specific coefficient bound and the cardinal values used below.

For the usual spline interpolation problem the sites are distinct and ordered; the Schoenberg–Whitney theorem then makes $A_x$ invertible under the stated support positivity. More generally, the following argument applies whenever the interpolation map in the question is well defined and $A_x$ is invertible.

Factor the spline interpolation operator as $P_x=B{A}_x^{-1}R$, whereSampling has operator norm at most one from the supremum norm to ${\ell }_{\infty }$. Nonnegativity and the subpartition of unity for B-splines giveConsequently $\Vert P_x\Vert \le \Vert A_x^{-1}{\Vert }_{{\ell }_{\infty }}$.

For the reverse bound, select a row $r$ of $A_x^{-1}$ with maximal absolute row sum, and choose $y_j=\mathrm{sgn}((A_x^{-1}{)}_{rj})$, assigning any value of modulus at most one when the entry is zero. Then $\Vert y{\Vert }_{{\ell }_{\infty }}=1$ andAt distinct sites there is a continuous piecewise linear function $f$ with $f(x_j)=y_j$ and $\Vert f{\Vert }_{\infty }=1$: interpolate linearly between the ordered sites and extend constantly to the endpoints. The given uniform-norm stability of a B-spline basis therefore impliesThus the B-spline interpolation operator norm satisfiesPositivity $N_i(x_i)>0$ alone, without the usual distinct-site convention, is insufficient to define an interpolation operator. For example, take order-two B-splines with knots $0,1,2,3$ and $x_1=x_2=3/2$. Both diagonal values are $1/2$, but the B-spline collocation matrix has identical rows $(1/2,1/2)$ and is singular. The displayed proof uses the well-defined interpolation map posited in the question; ordered distinct sites supply its standard existence guarantee.

For the uniform spline knots, translation reduces every degree-two basis function to a quadratic cardinal B-spline. From the Cox-de Boor recurrence, an order-three function takes the valuesAt the prescribed site $x_i=i+2$, only $N_i$ and, when $i<n$, $N_{i+1}$ have nonzero values. Hence the B-spline collocation matrix isThe shift satisfies $S^n=0$. A finite geometric expansion gives the exact inverseIts $i$th absolute row sum is $2(n-i+1)$, so $\Vert A^{-1}{\Vert }_{{\ell }^{\infty }}=2n$. Use the preceding bounds and the supplied stability constant $d_3=3$:Thus the linear growth of shifted quadratic spline interpolation is $\Theta (n)$, which in particular is $O(n)$. The lower bound, rather than the $O(n)$ upper bound alone, proves that these spline interpolation operators are not uniformly bounded.