Skip to content

定义Definition

Chebyshev 多项式与节点

Chebyshev polynomials and nodes · Chebyshev nodes

以余弦表示和端点聚集节点控制区间上一致逼近、插值条件性与高次振荡。

形式陈述 ​

第一类 Chebyshev 多项式是正交多项式族,其标准定义为

Tn(x)=cos⁡(narccos⁡x),−1≤x≤1.

其正交性使用的权是 w(x)=(1−x2)−1/2,不是区间上的常数权。用 x=cos⁡θ 换元,把收敛的反常积分变为余弦的正交积分:

∫−11Tm(x)Tn(x)1−x2dx=∫0πcos⁡(mθ)cos⁡(nθ)dθ={0,m≠n,π,m=n=0,π/2,m=n≥1.

例如不加权时 ∫−11T0T2dx=−2/3,不能沿用同一正交声明。

余弦倍角公式给出

T0(x)=1,T1(x)=x,Tn+1(x)=2xTn(x)−Tn−1(x)(n≥1).

当 n≥1 时,Tn 在 [−1,1] 上的极值点为 cos⁡(jπ/n)(0≤j≤n),函数值交替为 ±1;其零点为

cos⁡(2j+1)π2n,j=0,…,n−1.

次数 n 插值常用两套不同节点。Chebyshev roots 是 Tn+1 的零点

xj=cos⁡(2j+1)π2n+2,j=0,…,n;

Chebyshev–Lobatto 节点则在 n≥1 时包含两个端点,

xj=cos⁡jπn,j=0,…,n.

两者不能混称为同一公式。对 n≥1,Tn 的首项系数为 2n−1,首一多项式 Cn=21−nTn 的一致范数为 21−n,且在全部 n 次实首一多项式中最小。若另一首一多项式 q 的范数更小,则 q−Cn 在 n+1 个交替极值点上的符号也严格交替,每对相邻极值点之间至少有一个零点。这样得到 n 个不同零点,但首项相消使 deg⁡(q−Cn)≤n−1,矛盾。

节点的端点聚集由余弦映射自然产生,并使相应Lebesgue 常数只按对数级增长。

直觉

把角度 θ 上的均匀点通过 x=cos⁡θ 投到区间,会在 x=±1 附近自动挤密。高次多项式最容易在端点之间产生大摆动;多放一些端点附近的数据,相当于在最危险区域增加约束。

Chebyshev 多项式还可视为复单位圆上 zn 与 z−n 的对称组合:若 x=(z+z−1)/2,则 Tn(x)=(zn+z−n)/2。这个复数图像解释了余弦结构,也说明目标函数在复平面中离区间最近的奇异点会影响系数和逼近速度。

同一余弦倍角结构也出现在量子化中:把收缩 Hermitian 矩阵的特征值 x 变为二维旋转角 arccos⁡x,重复旋转的压缩块便给出 Tn(x)。量子信号处理在旋转之间加入可选相位来改变多项式响应;其复振幅还必须满足酉补全条件,不能仅从区间一致有界就断言可实现。

例子与边界

对 Runge 函数

f(x)=11+25x2,x∈[−1,1],

分别用 n+1 个等距节点 xj=−1+2j/n 与 Chebyshev–Lobatto 节点 xj=cos⁡(jπ/n) 做次数 n 插值。两组都用第二重心公式和各自的结构化权重在 binary64 中求值;在 [−1,1] 上取 20,001 个等距检验点,以

En=maxk|pn(tk)−f(tk)|

作为一致误差的高分辨率网格近似。按上述节点、结构化权重和检验网格逐一求值,独立运行脚本 在 NumPy 2.3.5 中得到下表:

次数 n 等距节点 En Chebyshev–Lobatto En 误差比
10 1.915659 1.321974×10−1 1.449089×101
20 5.982231×101 1.773782×10−2 3.372584×103
30 2.388281×103 2.425789×10−3 9.845376×105

表中数值是该次 binary64 实现的网格最大值,不是经过区间认证的连续上确界。尤其高次等距节点存在严重消去,求和顺序或数学库差异可改变末几位;比较节点趋势不应把全部打印数字都当作可移植精度。

两组使用同一目标、次数、求值公式和误差指标,因此差异主要来自节点。Runge 函数在复平面靠近实区间处有极点 ±i/5;等距插值的端点摆动和快速增长的 Lebesgue 常数放大这种困难,而 Lobatto 节点控制首一节点多项式并抑制端点放大。实验不是“换公式”的胜利,而是节点几何改变了插值问题。

Chebyshev 节点并非对所有目标、范数和约束绝对最优。周期函数可能更适合 Fourier 网格,局部尖峰可能需要分段或自适应节点,加权范数也会选择不同的正交族。即使节点良好,把高次多项式先展开成单项式系数仍可能引入严重舍入误差,应使用重心、Clenshaw 或离散余弦变换等结构化计算。

可以把“好节点”与最佳一致多项式直接比较。对 f=x3,两个 Chebyshev roots 为 ±1/2,一次插值是 x/2,端点误差达到 1/2;最佳一次式却是 3x/4,全区间误差只有 1/4,其交错证书可逐点验证。上文首一多项式的极小性优化的是节点乘积,不能据此宣布任何目标的节点插值都已最佳。已有 Chebyshev 系数的求值则交给Clenshaw 反向递推,它承担计算坐标的任务。

推论与应用

Chebyshev roots 与 Lobatto 节点都来自同一余弦结构,却分别适合不含端点和必须包含端点的任务。实现与文字必须明确采用哪一种,权重、端点条件和节点个数才能保持一致。

离散余弦变换可在节点值与 Chebyshev 系数之间快速转换,连接快速变换与谱方法;本页只建立节点和多项式结构,不把变换细节当作插值存在性的组成部分。

节点与多项式结构确定后,可沿三个计算方向继续:微分矩阵精确求插值导数,Clenshaw–Curtis 求积精确求插值积分,有理 Chebyshev 方法则通过映射把这套工具延伸到半无限区间。

参考资料
  • 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.
关系图谱18 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系