跳转至内容

三角函数/爱好者/切比雪夫多项式

来自维基教科书,开放的书籍,开放的世界


切比雪夫多项式,以 帕夫努提·切比雪夫 命名,[1] 是与三角多角公式相关的多项式序列。

我们通常区分

  • 第一类切比雪夫多项式,记为Tn,与 密切相关,以及
  • 第二类切比雪夫多项式,记为Un,与 密切相关。

字母T 的使用是因为切比雪夫 的另一种转写方式是Tchebycheff(法语)或Tschebyschow(德语)。

切比雪夫多项式TnUnn 次多项式,并且任何一种切比雪夫多项式序列都构成一个“多项式序列”。

在域 −1 < x < 1 中,前几个第一类切比雪夫多项式:平坦T0, T1, T2, T3, T4T5.

前几个第一类切比雪夫多项式是

在域 −1 < x < 1 中,第二类切比雪夫多项式的最初几个:图中的扁平U0, U1, U2, U3, U4U5. 虽然图像中没有显示出来,但 Un(1) = n + 1Un(−1) = (n + 1)(−1)n.

第二类切比雪夫多项式的最初几个是

的定义

[edit | edit source]

定义 有许多不同的方法,它们都得出相同的多项式。我们用来定义 的定义是

换句话说, 是用 表示 的多项式。

例如

直接来自


使用此定义,切比雪夫多项式的递推关系立即得出


递推来自关系

这是从余弦加法公式推导的关系的重新排列

的特殊情况,其中

我们在将两个波加在一起时看到的拍频。

用途

[edit | edit source]

如果在一个区间 上逼近某个函数,可以使用多项式逼近,例如

为了近似它的值。因为 在 0 到 1 之间, 项从左到右越来越小。当项数足够多时,通常比这里多,我们通常可以对行为良好的函数得到非常好的近似值。在我们的例子中,我们将级数截断在第四项,并且忽略了 及更高阶的项。

在上面的例子多项式中,实际的数字只是虚构的数字,作为示例,没有特别的意义——如果你想知道为什么使用这些特定数字。

事实证明,从某种意义上说,截断级数不是用三次多项式近似原始函数的实际值的最佳方法。一种更复杂但更好的方法使用切比雪夫多项式。


如果相反,我们可以将函数表示为

其中 是常数,那么如果我们将这个级数截断在第四项,即 T3,我们又得到一个没有 或更高阶项的多项式。在我们展开并收集系数后,这些系数通常不会与仅仅截断级数时的系数相同。这听起来可能很疯狂,所以一个实际例子可以更清楚地说明这一点。

实际例子:展开 (2-x)-1

快速检查一下...这个公式看起来合理吗?

代入 ,我们得到

代入 ,我们得到

.

这个公式看起来合理。

如果我们只取四项,那么在 处的误差是

现在用切比雪夫多项式来表示这个公式

展开后得到

待办事项:添加图表显示最大误差更小。不幸的是,当前计算结果没有显示这一点。

嗯……看起来我需要重新计算切比雪夫系数。



切比雪夫多项式在逼近理论中很重要,因为切比雪夫多项式 *T**n* 的根被用作多项式插值的节点。得到的插值多项式最小化了龙格现象的问题,并在最大范数下提供了一个接近连续函数最佳逼近多项式的逼近。



显式公式

[edit | edit source]

定义切比雪夫多项式的不同方法会导致不同的显式公式,例如



其中 是一个超几何函数.

切比雪夫多项式在逼近理论中很重要,因为切比雪夫多项式 *T**n* 的根被用作多项式插值的节点。得到的插值多项式最小化了龙格现象的问题,并在最大范数下提供了一个接近连续函数最佳逼近多项式的逼近。

微分方程研究中,它们是切比雪夫微分方程的解。

分别对应第一类和第二类多项式。这些方程是 Sturm–Liouville 微分方程 的特例。

定义

[edit | edit source]

第一类切比雪夫多项式 由递归关系定义

Tn 的传统 生成函数

指数生成函数为

第二类切比雪夫多项式递归关系 定义

Un 的一个例子是 生成函数

三角定义

[edit | edit source]

第一类切比雪夫多项式可以通过 三角恒等式 定义

因此

对于 n = 0, 1, 2, 3, ..., 而第二类多项式满足

这在结构上与狄利克雷核 非常相似。

cos(nx) 是 cos(x) 的 n 次多项式,这可以通过观察 cos(nx) 是 棣莫弗公式 一侧的实部,而另一侧的实部是 cos(x) 和 sin(x) 的多项式,其中 sin(x) 的所有幂都是偶数,因此可以通过恒等式 cos2(x) + sin2(x) = 1 替换。

这个恒等式与递归生成公式结合起来非常有用,因为它可以让人只用基本角的余弦来计算任何整数倍角的余弦。 评估前两个切比雪夫多项式

可以很容易地确定

等等。为了简单地检查结果是否合理,将等号两边的系数相加(即设置 等于零,此时余弦为 1),可以看到前一个表达式中 1 = 2 − 1,后一个表达式中 1 = 4 − 3

两个直接推论是合成恒等式(或“嵌套属性”)

以及复指数形式的切比雪夫多项式表达式:给定 z = a + bi

佩尔方程定义

[编辑 | 编辑源代码]

切比雪夫多项式也可以定义为佩尔方程的解

在环 R[x] 中。[2] 因此,它们可以通过佩尔方程的标准技术生成,该技术采用基本解的幂

第一类和第二类切比雪夫多项式之间的关系

[编辑 | 编辑源代码]

第一类和第二类切比雪夫多项式之间有着密切的关系,它们可以通过以下方程来表示

切比雪夫多项式导数的递归关系可以从这些关系中推导出

此关系式在求解微分方程的切比雪夫谱方法中使用。

等效地,这两个序列也可以从一对互递归方程定义:

这些可以从三角公式推导出;例如,如果,那么

请注意,如果我们像某些著作一样,遵循将我们的Un(n次多项式)用Un+1表示的替代约定,则这些方程和三角方程都将采用更简单的形式。


根和极值

[编辑 | 编辑来源]

任何一种 n 次切比雪夫多项式在区间 [−1,1] 内有 n 个不同的简单根,称为切比雪夫根。这些根有时被称为切比雪夫节点,因为它们用作多项式插值的节点。使用三角定义和以下事实

可以很容易地证明Tn的根是

类似地,Un的根是

第一类切比雪夫多项式的一个独特的性质是在区间 −1 ≤ x ≤ 1 上,所有 极值 的值都为 −1 或 1。因此,这些多项式只有两个有限的 临界值,这是 Shabat 多项式 的定义属性。第一类和第二类切比雪夫多项式都在端点处有极值,分别为

微分与积分

[edit | edit source]

多项式的导数可能不像直观那样简单。通过对多项式在三角形式上的微分,很容易证明

由于在 x = 1x = −1 处会出现零除(0/0 不确定形式,具体来说),最后两个公式在数值上可能会有问题。可以证明

证明

第一类 切比雪夫多项式 的二阶导数为

如果如上所示进行评估,会产生问题,因为函数在x = ±1处不确定。由于该函数是多项式,(所有)导数必须对所有实数存在,因此对上述表达式求极限应得到期望的值。

其中目前只考虑。对分母进行因式分解

由于整体极限必须存在,因此分子和分母的极限必须独立存在,并且

分母(仍然)趋于零,这意味着分子必须趋于零,即 ,这在后面将很有用。由于分子和分母都趋于零,因此洛必达法则适用。

证明 类似,需要注意

实际上,以下更一般的公式成立:

该结果在特征值问题的数值解法中非常有用。

关于积分,Tn的一阶导数意味着

第一类多项式的递推关系涉及导数,因此可以建立

正交性

[编辑 | 编辑源代码]

共同构成了一系列的正交多项式。第一类多项式关于权重

在区间 上是正交的,即有

这可以通过令 并使用恒等式 来证明。

类似地,第二类多项式关于权重

在区间 上是正交的,即有

(注意,权重 在归一化常数范围内是维格纳半圆分布 的密度)。

也满足离散正交条件

其中 高斯勒让德 零点

最小 ∞-范数

[编辑 | 编辑源代码]

对于任何给定的 n ≥ 1,在所有首项系数为 1 的 n 次多项式中,

是那个在区间 [−1, 1] 上的最大绝对值最小的多项式。

这个最大绝对值为

并且|ƒ(x)| 恰好在 n + 1 个地方达到这个最大值:在

假设 是一个首项系数为 1 的 n 次多项式,它在区间 [−1, 1] 上的最大绝对值小于 .

我们定义

因为在 的极值点处,我们有

是一个 n - 1 次多项式,根据 介值定理,它至少有 n 个根,但这对于一个 n - 1 次多项式来说是不可能的。

其他性质

[编辑 | 编辑源代码]

切比雪夫多项式是超球面多项式或 盖根鲍尔多项式 的特例,而盖根鲍尔多项式本身是 雅可比多项式 的特例。

对于每一个非负整数 nTn(x) 和 Un(x) 都是 n 次多项式。它们是 x偶函数或奇函数,因为 n 是偶数或奇数,所以当它们被写成 x 的多项式时,它们分别只有偶数次项或奇数次项。事实上,

Tn 的最高次项系数为 2n − 1,如果 1 ≤ n,但如果 0 = n,则为 1。

Tn李萨如曲线 的特例,频率比为 n

一些多项式序列,如 卢卡斯多项式 (Ln)、迪克森多项式(Dn)、斐波那契多项式(Fn),都与切比雪夫多项式 TnUn 有关。

第一类切比雪夫多项式满足以下关系

这很容易从余弦的积化和公式推导出来。第二类多项式满足类似的关系

.

类似于公式

我们有类似的公式

.


作为基集

[编辑 | 编辑源代码]
非光滑函数(顶部)y = −x3H(−x),其中 HHeaviside 阶跃函数,以及(底部)其 Chebyshev 展开的第 5 个偏和。第 7 个和在图形分辨率下与原始函数无法区分。

在适当的Sobolev 空间中,Chebyshev 多项式集形成一个完整的基集,因此相同空间中的函数可以在 −1 ≤ x ≤ 1 上通过以下展开式表示:[3]

此外,如前所述,Chebyshev 多项式形成了一个正交的基,这意味着(除其他外)系数 an 可以通过应用内积轻松确定。这个和被称为Chebyshev 级数Chebyshev 展开

由于 Chebyshev 级数与傅里叶余弦级数通过变量变化相关联,因此所有适用于傅里叶级数的定理、恒等式等都有一个 Chebyshev 对应物。[3] 这些属性包括

  • Chebyshev 多项式构成一个完整的正交系统。
  • 如果函数是分段光滑连续,则 Chebyshev 级数收敛于 ƒ(x)。在大多数情况下,光滑度要求可以放宽 - 只要 ƒ(x) 及其导数的间断点有限。
  • 在间断点处,级数将收敛于左右极限的平均值。

傅里叶级数继承的大量定理和恒等式使 Chebyshev 多项式成为数值分析中的重要工具;例如,它们是谱方法中最流行的通用基函数[3],通常优于三角级数,因为对于连续函数来说收敛速度通常更快(吉布斯现象仍然是一个问题)。

考虑 的 Chebyshev 展开。我们可以表达

可以通过应用内积或离散正交性条件来找到系数。对于内积,

得到

或者,当您无法评估您尝试逼近的函数的内积时,离散正交性条件给出

其中克罗内克德尔塔 函数,而 高斯-勒让德 零点。

这使得我们可以通过离散余弦变换非常有效地计算系数

示例 2

[edit | edit source]

为了提供另一个示例


部分和

[编辑 | 编辑源代码]

的部分和

在各种函数的近似和微分方程的求解中非常有用(参见谱方法)。确定系数 *a**n* 的两种常用方法是通过使用内积,如在伽辽金方法中,以及通过使用与插值相关的配置方法。

作为插值函数,(*N* − 1)th 部分和的 *N* 个系数通常是在切比雪夫-高斯-勒让德点(或勒让德网格)上获得的,这会导致最小误差并避免与均匀网格相关的龙格现象。

切比雪夫形式的多项式

[编辑 | 编辑源代码]

任意 *N* 次多项式可以用第一类切比雪夫多项式表示。这样的多项式 *p*(*x*) 具有以下形式

可以用克伦肖算法对切比雪夫形式的多项式进行求值。

扩展多项式

[编辑 | 编辑源代码]

扩展多项式在某种意义上等效于第一类切比雪夫多项式,但它们使人们能够在某些情况下避免平方根和常规三角函数,尤其是在有理三角学中。


  1. 切比雪夫多项式首次出现在:P. L. Chebyshev (1854) "Théorie des mécanismes connus sous le nom de parallélogrammes," Mémoires des Savants étrangers présentés à l’Académie de Saint-Pétersbourg, vol. 7, pages 539-586.
  2. Jeroen Demeyer Diophantine Sets over Polynomial Rings and Hilbert's Tenth Problem for Function Fields, Ph.D. theses (2007), p.70.
  3. a b c Boyd, John P. (2001). 切比雪夫与傅立叶谱方法 (PDF) (第二版). 多佛. ISBN 0486411834.
  4. 切比雪夫插值:互动之旅
华夏公益教科书