当前位置: 首页 > news >正文

切比雪夫多项式

切比雪夫多项式

约定

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

∥f∥∞=sup⁡x∈[a,b]∣f(x)∣ \|f\|_\infty = \sup_{x \in [a,b]} |f(x)| f=x[a,b]supf(x)

其中:

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

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

∥f∥∞=max⁡x∈[a,b]∣f(x)∣ \|f\|_\infty = \max_{x \in [a,b]} |f(x)| f=x[a,b]maxf(x)


一、定义

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

cos⁡((n+1)x)=2cos⁡x⋅cos⁡(nx)−cos⁡((n−1)x) \cos((n+1)x)=2\cos x\cdot\cos(nx)-\cos((n-1)x) cos((n+1)x)=2cosxcos(nx)cos((n1)x)

  • 递归定义
    T0(x)=1,T1(x)=x,Tn+1(x)=2xTn(x)−Tn−1(x) T_0(x) = 1,\quad T_1(x) = x,\quad T_{n+1}(x) = 2xT_n(x) - T_{n-1}(x) T0(x)=1,T1(x)=x,Tn+1(x)=2xTn(x)Tn1(x)

  • 三角函数定义(常用于区间 [−1,1][-1, 1][1,1]):
    Tn(x)=cos⁡(narccos⁡x) T_n(x) = \cos(n \arccos x) Tn(x)=cos(narccosx)


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

  • 递归定义
    U0(x)=1,U1(x)=2x,Un+1(x)=2xUn(x)−Un−1(x) U_0(x) = 1,\quad U_1(x) = 2x,\quad U_{n+1}(x) = 2xU_n(x) - U_{n-1}(x) U0(x)=1,U1(x)=2x,Un+1(x)=2xUn(x)Un1(x)

  • 三角函数定义
    Un(x)=sin⁡((n+1)arccos⁡x)1−x2 U_n(x) = \frac{\sin((n+1)\arccos x)}{\sqrt{1 - x^2}} Un(x)=1x2sin((n+1)arccosx)


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

定理内容

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

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

且其最大偏差为:

∥pn∥∞=12n−1 \left\| p_n \right\|_{\infty} = \frac{1}{2^{n-1}} pn=2n11


证明思路

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


步骤一:等振性质

第一类切比雪夫多项式:
Tn(x)=cos⁡(narccos⁡x) T_n(x) = \cos(n \arccos x) Tn(x)=cos(narccosx)

在区间 [−1,1][-1, 1][1,1] 上,有 n+1n+1n+1 个交错极值点 x0<x1<⋯<xnx_0 < x_1 < \dots < x_nx0<x1<<xn,满足:
Tn(xi)=(−1)i,i=0,1,…,n T_n(x_i) = (-1)^i,\quad i = 0,1,\dots,n Tn(xi)=(1)i,i=0,1,,n

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


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

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

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

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

∥pn∥∞=max⁡x∈[−1,1]∣12n−1Tn(x)∣=12n−1 \|p_n\|_{\infty} = \max_{x \in [-1,1]} \left| \frac{1}{2^{n-1}} T_n(x) \right| = \frac{1}{2^{n-1}} pn=x[1,1]max2n11Tn(x)=2n11


步骤三:反证法

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

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

定义误差函数:

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

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

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

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


结论

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


二、正交性质

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

∫−11Tm(x)Tn(x)11−x2dx={0,m≠nπ,n=0π2,n=m≠0 \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} 11Tm(x)Tn(x)1x21dx=0,π,2π,m=nn=0n=m=0

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

∫−11Um(x)Un(x)1−x2dx={0,m≠nπ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} 11Um(x)Un(x)1x2dx={0,2π,m=nm=n


三、性质与定理

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

[−1,1][-1, 1][1,1] 上,最小最大误差的插值点为:
xk=cos⁡(2k−12nπ),k=1,2,…,n x_k = \cos\left(\frac{2k-1}{2n} \pi\right),\quad k=1,2,\dots,n xk=cos(2n2k1π),k=1,2,,n

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

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

3. 零点分布

  • 第一类零点:
    xk=cos⁡(2k−12nπ),k=1,…,n x_k = \cos\left(\frac{2k-1}{2n} \pi\right), \quad k=1,\dots,n xk=cos(2n2k1π),k=1,,n
  • 第二类零点:
    xk=cos⁡(kn+1π),k=1,…,n x_k = \cos\left(\frac{k}{n+1} \pi\right), \quad k=1,\dots,n xk=cos(n+1kπ),k=1,,n

四、应用

1. 插值理论

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

2. 数值逼近

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

3. 谱方法(Spectral Methods)

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

4. 信号处理

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

5. 最优化与迭代方法

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

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

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

http://www.dtcms.com/a/279383.html

相关文章:

  • leetcode 1290. 二进制链表转整数 简单
  • C++类模版与友元
  • 进程、线程、协程
  • windows内核研究(进程与线程-进程结构体EPROCESS)
  • Django基础(一)———创建与启动
  • 【反转链表专题】【LeetCode206.反转链表】【LeetCode25.K个一组翻转链表】【LeetCode234.回文链表】
  • Spring Boot 自带的 JavaMail 集成
  • android Perfetto cpu分析教程及案例
  • 5G 到 6G通信技术的革新在哪里?
  • 腾讯云和火山云优劣势对比
  • 电力协议处理框架C++版(三)
  • CA-IS3082W 隔离485 收发器芯片可能存在硬件BUG
  • LTspic下载,帮助及演示电路
  • sfe_py的应力云图计算与显示step by step
  • 暑期自学嵌入式——Day02(C语言阶段)
  • 揭开图像的秘密:OpenCV直方图入门详解
  • 代数基本定理最简短的证明
  • 对于独热编码余弦相似度结果为0和词向量解决了词之间相似性问题的理解
  • ubuntu之坑(十五)——设备树
  • gRPC和http长轮询
  • 新手向:Python自动化办公批量重命名与整理文件系统
  • Dubbo 学习笔记
  • 谷歌收获成果:OpenAI收购Windsurf计划告吹
  • 工业软件加密锁复制:一场技术与安全的博弈
  • Mybatis05-参数和返回
  • 以太网供电(PoE)电源
  • 编程语言设计目的与侧重点全解析(主流语言深度总结)
  • vue中使用西瓜播放器xgplayer (封装)+xgplayer-hls 播放.m3u8格式视频
  • Spark 单机模式安装与测试全攻略​
  • STM32小实验1--点亮LED