切比雪夫多项式

78 阅读1分钟

切比雪夫多项式

约定

设函数 f(x)f(x) 定义在区间 [a,b][a, b] 上,则其无穷范数定义为:

f=supx[a,b]f(x)\|f\|_\infty = \sup_{x \in [a,b]} |f(x)|

其中:

  • sup\sup 表示上确界(least upper bound)
  • f(x)f(x)[a,b][a,b] 上连续时,这个上确界就是最大值。

因此若 f(x)f(x) 连续,有:

f=maxx[a,b]f(x)\|f\|_\infty = \max_{x \in [a,b]} |f(x)|

一、定义

1. 第一类切比雪夫多项式 Tn(x)T_n(x)

cos((n+1)x)=2cosxcos(nx)cos((n1)x)\cos((n+1)x)=2\cos x\cdot\cos(nx)-\cos((n-1)x)
  • 递归定义

    T0(x)=1,T1(x)=x,Tn+1(x)=2xTn(x)Tn1(x)T_0(x) = 1,\quad T_1(x) = x,\quad T_{n+1}(x) = 2xT_n(x) - T_{n-1}(x)
  • 三角函数定义(常用于区间 [1,1][-1, 1]):

    Tn(x)=cos(narccosx)T_n(x) = \cos(n \arccos x)

2. 第二类切比雪夫多项式 Un(x)U_n(x)

  • 递归定义

    U0(x)=1,U1(x)=2x,Un+1(x)=2xUn(x)Un1(x)U_0(x) = 1,\quad U_1(x) = 2x,\quad U_{n+1}(x) = 2xU_n(x) - U_{n-1}(x)
  • 三角函数定义

    Un(x)=sin((n+1)arccosx)1x2U_n(x) = \frac{\sin((n+1)\arccos x)}{\sqrt{1 - x^2}}

切比雪夫极值定理证明(最小最大误差定理)

定理内容

在所有最高项系数为 1 的 nn 次多项式中,
在区间 [1,1][-1,1] 上与 0 的最大偏差最小的多项式是:

pn(x)=12n1Tn(x)p_n(x) = \frac{1}{2^{n-1}} T_n(x)

且其最大偏差为:

pn=12n1\left\| p_n \right\|_{\infty} = \frac{1}{2^{n-1}}

证明思路

采用反证法,结合切比雪夫多项式的等振性质(Equioscillation)来证明其最小最大逼近性。


步骤一:等振性质

第一类切比雪夫多项式:

Tn(x)=cos(narccosx)T_n(x) = \cos(n \arccos x)

在区间 [1,1][-1, 1] 上,有 n+1n+1 个交错极值点 x0<x1<<xnx_0 < x_1 < \dots < x_n,满足:

Tn(xi)=(1)i,i=0,1,,nT_n(x_i) = (-1)^i,\quad i = 0,1,\dots,n

这是“等振性”(Equioscillation)的体现。


步骤二:构造归一化多项式

因为 Tn(x)T_n(x) 的最高项系数是 2n12^{n-1},所以定义单位最高项多项式:

pn(x)=12n1Tn(x)p_n(x) = \frac{1}{2^{n-1}} T_n(x)

它的最高次项系数为 1,且最大模为:

pn=maxx[1,1]12n1Tn(x)=12n1\|p_n\|_{\infty} = \max_{x \in [-1,1]} \left| \frac{1}{2^{n-1}} T_n(x) \right| = \frac{1}{2^{n-1}}

步骤三:反证法

假设存在另一个多项式 qn(x)q_n(x),也是次数为 nn,最高项系数为 1,且满足:

qn(x)<pn(x)=12n1\|q_n(x)\|_\infty < \|p_n(x)\|_\infty = \frac{1}{2^{n-1}}

定义误差函数:

E(x)=pn(x)qn(x)E(x) = p_n(x) - q_n(x)

E(x)E(x) 也是次数不超过 nn 的多项式,且:

  • x0,x1,,xnx_0, x_1, \dots, x_n 上:
    E(xi)=12n1(1)iqn(xi)E(x_i) = \frac{1}{2^{n-1}} (-1)^i - q_n(x_i)
  • 因为 qn(x)<12n1|q_n(x)| < \frac{1}{2^{n-1}},所以 E(xi)E(x_i) 的符号与 (1)i(-1)^i 相同;
  • E(x)E(x) 在区间上至少有 nn 个变号点 → 至少有 nn 个零点;
  • deg(E)n\deg(E) \le n,所以这 nn 个零点是全部根;
  • 所以 E(x)=ci=1n(xξi)E(x) = c \prod_{i=1}^n (x - \xi_i),是一个 nn 次多项式。

但注意:由于 pn(x)p_n(x)qn(x)q_n(x) 的最高项系数都是 1, 则 E(x)E(x) 的最高项系数为 0 —— 这与它是 nn 次多项式矛盾!


结论

矛盾说明假设不成立,
所以 12n1Tn(x)\boxed{\frac{1}{2^{n-1}} T_n(x)} 是所有单位最高项多项式中与 0 偏差最小者,
且最小偏差正好为 12n1\boxed{\frac{1}{2^{n-1}}}


二、正交性质

1. 第一类 Tn(x)T_n(x) 正交性(权函数 w(x)=11x2w(x) = \frac{1}{\sqrt{1 - x^2}}):

11Tm(x)Tn(x)11x2dx={0,mnπ,n=0π2,n=m0\int_{-1}^{1} T_m(x) T_n(x) \frac{1}{\sqrt{1 - x^2}} dx = \begin{cases} 0, & m \ne n \\ \pi, & n = 0 \\ \frac{\pi}{2}, & n = m \ne 0 \end{cases}

2. 第二类 Un(x)U_n(x) 正交性(权函数 w(x)=1x2w(x) = \sqrt{1 - x^2}):

11Um(x)Un(x)1x2dx={0,mnπ2,m=n\int_{-1}^{1} U_m(x) U_n(x) \sqrt{1 - x^2} dx = \begin{cases} 0, & m \ne n \\ \frac{\pi}{2}, & m = n \end{cases}

三、性质与定理

1. 切比雪夫极值定理(最小最大误差)

[1,1][-1, 1] 上,最小最大误差的插值点为:

xk=cos(2k12nπ),k=1,2,,nx_k = \cos\left(\frac{2k-1}{2n} \pi\right),\quad k=1,2,\dots,n

2. 最佳一致逼近定理(Chebyshev Approximation)

用切比雪夫多项式进行逼近能得到一致逼近误差最小的结果。

3. 零点分布

  • 第一类零点:
    xk=cos(2k12nπ),k=1,,nx_k = \cos\left(\frac{2k-1}{2n} \pi\right), \quad k=1,\dots,n
  • 第二类零点:
    xk=cos(kn+1π),k=1,,nx_k = \cos\left(\frac{k}{n+1} \pi\right), \quad k=1,\dots,n

四、应用

1. 插值理论

  • 切比雪夫节点用于插值可显著减小误差,避免龙格现象。
  • 比等距插值稳定得多。

2. 数值逼近

  • 切比雪夫多项式用于构造最佳一致逼近。
  • 广泛用于最小最大逼近问题。

3. 谱方法(Spectral Methods)

  • 切比雪夫多项式用于偏微分方程的数值解(如Navier-Stokes)。
  • 常配合 FFT 使用。

4. 信号处理

  • 用作滤波器设计中的窗函数。
  • 应用于 FIR/IIR 滤波器逼近。

5. 最优化与迭代方法

  • 常出现在误差函数或最优逼近算法中。

五、与其他正交多项式比较

多项式类型权函数定义域应用方向
切比雪夫一类 TnT_n11x2\frac{1}{\sqrt{1 - x^2}}[1,1][-1, 1]最小最大逼近
切比雪夫二类 UnU_n1x2\sqrt{1 - x^2}[1,1][-1, 1]谱方法、物理建模
勒让德多项式11[1,1][-1, 1]通用正交逼近
拉盖尔多项式exe^{-x}[0,+)[0, +\infty)量子力学、谱分析
埃尔米特多项式ex2e^{-x^2}(,+)(-\infty, +\infty)概率统计、量子力学