形式陈述
第一类 Chebyshev 多项式是正交多项式理路正交多项式Orthogonal polynomials在带正权的函数内积下构造逐次升阶的正交多项式,并由三项递推与交错零点组织计算。族,其标准定义为
其正交性使用的权是 ,不是区间上的常数权。用 换元,把收敛的反常积分变为余弦的正交积分:
例如不加权时 ,不能沿用同一正交声明。
余弦倍角公式给出
当 时, 在 上的极值点为 (),函数值交替为 ;其零点为
次数 插值常用两套不同节点。Chebyshev roots 是 的零点
Chebyshev–Lobatto 节点则在 时包含两个端点,
两者不能混称为同一公式。对 , 的首项系数为 ,首一多项式 的一致范数为 ,且在全部 次实首一多项式中最小。若另一首一多项式 的范数更小,则 在 个交替极值点上的符号也严格交替,每对相邻极值点之间至少有一个零点。这样得到 个不同零点,但首项相消使 ,矛盾。
节点的端点聚集由余弦映射自然产生,并使相应Lebesgue 常数理路Lebesgue 常数与插值条件性Lebesgue constant · Interpolation conditioning用插值算子的无穷范数刻画节点集合对采样扰动和最佳逼近误差的放大。只按对数级增长。
直觉
把角度 上的均匀点通过 投到区间,会在 附近自动挤密。高次多项式最容易在端点之间产生大摆动;多放一些端点附近的数据,相当于在最危险区域增加约束。
Chebyshev 多项式还可视为复单位圆上 与 的对称组合:若 ,则 。这个复数理路复数Complex number形如 a+bi 的数,按坐标规则构成实数域的二次扩张。图像解释了余弦结构,也说明目标函数在复平面中离区间最近的奇异点会影响系数和逼近速度。
同一余弦倍角结构也出现在量子化理路Qubitization 与二维信号子空间Qubitization从Hermitian信号块构造两反射量子行走,逐本征值推导二维旋转、端点退化和Chebyshev响应,并说明一般块编码所需的额外受控访问。中:把收缩 Hermitian 矩阵的特征值 变为二维旋转角 ,重复旋转的压缩块便给出 。量子信号处理理路量子信号处理Quantum signal processing · QSP给出固定单qubit信号约定下完整的多项式对可实现条件,解释降阶合成,并通过一阶、Chebyshev三阶和x立方相位序列区分完整矩阵元与实部响应。在旋转之间加入可选相位来改变多项式响应;其复振幅还必须满足酉补全条件,不能仅从区间一致有界就断言可实现。
例子与边界
对 Runge 函数
分别用 个等距节点 与 Chebyshev–Lobatto 节点 做次数 插值。两组都用第二重心公式理路重心 Lagrange 插值Barycentric Lagrange interpolation · Barycentric interpolation预计算重心权后以线性成本求值 Lagrange 插值多项式,并说明节点命中、权重尺度与数值稳定性的条件。和各自的结构化权重在 binary64 中求值;在 上取 20,001 个等距检验点,以
作为一致误差的高分辨率网格近似。按上述节点、结构化权重和检验网格逐一求值,独立运行脚本 在 NumPy 2.3.5 中得到下表:
| 次数 |
等距节点 |
Chebyshev–Lobatto |
误差比 |
| 10 |
|
|
|
| 20 |
|
|
|
| 30 |
|
|
|
表中数值是该次 binary64 实现的网格最大值,不是经过区间认证的连续上确界。尤其高次等距节点存在严重消去,求和顺序或数学库差异可改变末几位;比较节点趋势不应把全部打印数字都当作可移植精度。
两组使用同一目标、次数、求值公式和误差指标,因此差异主要来自节点。Runge 函数在复平面靠近实区间处有极点 ;等距插值的端点摆动和快速增长的 Lebesgue 常数放大这种困难,而 Lobatto 节点控制首一节点多项式并抑制端点放大。实验不是“换公式”的胜利,而是节点几何改变了插值问题。
Chebyshev 节点并非对所有目标、范数和约束绝对最优。周期函数可能更适合 Fourier 网格,局部尖峰可能需要分段或自适应节点,加权范数也会选择不同的正交族。即使节点良好,把高次多项式先展开成单项式系数仍可能引入严重舍入误差,应使用重心、Clenshaw 或离散余弦变换等结构化计算。
可以把“好节点”与最佳一致多项式理路最佳一致多项式逼近Best uniform polynomial approximation · Minimax polynomial approximation · 极小极大多项式逼近在固定次数内最小化整区间最大误差,证明最优解存在,并区分插值、积分投影与一致最佳解。直接比较。对 ,两个 Chebyshev roots 为 ,一次插值是 ,端点误差达到 ;最佳一次式却是 ,全区间误差只有 ,其交错证书理路Chebyshev 交错定理Chebyshev alternation theorem · Equioscillation theorem · 等振荡定理用交替达到最大误差的点认证最佳多项式,完整证明必要、充分、唯一性及近最佳下界。可逐点验证。上文首一多项式的极小性优化的是节点乘积,不能据此宣布任何目标的节点插值都已最佳。已有 Chebyshev 系数的求值则交给Clenshaw 反向递推理路Clenshaw 三项递推求和Clenshaw algorithm · Clenshaw summation反向递推计算 Chebyshev 展开,证明收尾公式,并区分系数约定、算法误差与逼近误差。,它承担计算坐标的任务。
推论与应用
Chebyshev roots 与 Lobatto 节点都来自同一余弦结构,却分别适合不含端点和必须包含端点的任务。实现与文字必须明确采用哪一种,权重、端点条件和节点个数才能保持一致。
离散余弦变换可在节点值与 Chebyshev 系数之间快速转换,连接快速变换与谱方法;本页只建立节点和多项式结构,不把变换细节当作插值存在性的组成部分。
节点与多项式结构确定后,可沿三个计算方向继续:微分矩阵理路Chebyshev 微分矩阵Chebyshev differentiation matrix从 Lobatto 插值基函数导出稠密微分矩阵,并用多项式精确性与混叠反例核验算子。精确求插值导数,Clenshaw–Curtis 求积理路Clenshaw–Curtis 求积Clenshaw-Curtis quadrature精确积分 Chebyshev–Lobatto 插值多项式,得到可嵌套求积并区分加权正交与普通积分。精确求插值积分,有理 Chebyshev 方法理路无界区域的有理 Chebyshev 方法Rational Chebyshev method以有理映射把半轴变成有限区间,显式变换导数并检查尺度、无穷远条件和解析性变化。则通过映射把这套工具延伸到半无限区间。
参考资料
- Lloyd N. Trefethen, “Six Myths of Polynomial Interpolation and Quadrature,” Mathematics Today 47(4), 2011.
- Lloyd N. Trefethen, Approximation Theory and Approximation Practice, extended ed., SIAM, 2019.
- Jean-Paul Berrut and Lloyd N. Trefethen, “Barycentric Lagrange Interpolation,” SIAM Review 46(3), 2004.
- NIST Digital Library of Mathematical Functions, §18.3 Definitions, accessed 2026.