漫谈:Chebyshev多项式,Krylov子空间,Chebyshev迭代,共轭梯度方法
Posted zongy17
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了漫谈:Chebyshev多项式,Krylov子空间,Chebyshev迭代,共轭梯度方法相关的知识,希望对你有一定的参考价值。
漫谈:Chebyshev多项式,Krylov子空间,Chebyshev迭代,共轭梯度方法
Chebyshev迭代和共轭梯度方法的收敛速度(后者或称误差分析)都与Chebyshev多项式有着紧密联系,因此做一些整理,以期把其中的逻辑理清,推导理顺。只覆盖关键要点,不求面面俱到。
Chebyshev 多项式
标准Chebyshev多项式
与权函数
ρ
(
y
)
=
1
1
−
y
2
\\rho(y)=\\frac1\\sqrt1-y^2
ρ(y)=1−y21对应的正交多项式为Chebyshev多项式。标准的Chebyshev多项式分两段定义:
T
n
(
y
)
=
c
o
s
(
n
arccos
(
y
)
)
,
y
∈
[
−
1
,
1
]
1
2
[
(
y
+
y
2
−
1
)
n
+
(
y
−
y
2
−
1
)
n
]
,
y
∈
R
∖
[
−
1
,
1
]
T_n(y)=\\left\\ \\beginaligned &\\mathrmcos(n\\space\\arccos(y)) &\\space, y\\in [-1,1]\\\\ &\\frac12\\left[\\left(y+\\sqrty^2-1\\right)^n+\\left(y-\\sqrty^2-1\\right)^n\\right] &\\space, y \\in \\mathbfR \\setminus [-1,1] \\endaligned \\right.
Tn(y)=⎩⎨⎧cos(n arccos(y))21[(y+y2−1)n+(y−y2−1)n] ,y∈[−1,1] ,y∈R∖[−1,1]
概述一些基本性质。 n n n次的Chebyshev多项式在开区间 ( − 1 , 1 ) (-1,1) (−1,1)内有 n n n个零点,在 [ − 1 , 1 ] [-1,1] [−1,1]内有 n + 1 n+1 n+1个符号交错的极值点( − 1 -1 −1和 + 1 +1 +1交替)。Chebyshev多项式不是奇函数就是偶函数,且随 n n n交替变化。
对于任意位置区间
[
a
,
b
]
[a,b]
[a,b]的Chebyshev多项式,都可以通过平移和缩放来化为标准区间
[
−
1
,
1
]
[-1,1]
[−1,1]上的标准Chebyshev多项式。对于
x
∈
[
a
,
b
]
x\\in[a,b]
x∈[a,b],做变换
y
=
x
−
1
2
(
b
+
a
)
1
2
(
b
−
a
)
y=\\fracx-\\frac12(b+a)\\frac12(b-a)
y=21(b−a)x−21(b+a),即有
y
∈
[
−
1
,
1
]
y\\in[-1,1]
y∈[−1,1]:
T
n
~
(
x
)
=
T
n
(
y
=
x
−
1
2
(
b
+
a
)
1
2
(
b
−
a
)
)
\\widetildeT_n(x)=T_n\\left(y=\\fracx-\\frac12(b+a)\\frac12(b-a)\\right)
Tn
(x)=Tn(y=21(b−a)x−21(b+a))
相应地,权函数也做变换:
ρ
~
(
x
)
=
ρ
(
y
=
x
−
1
2
(
b
+
a
)
1
2
(
b
−
a
)
)
=
1
1
−
(
x
−
1
2
(
b
+
a
)
1
2
(
b
−
a
)
)
2
\\widetilde\\rho(x)=\\rho\\left(y=\\fracx-\\frac12(b+a)\\frac12(b-a)\\right)= \\frac1\\sqrt1-\\left( \\fracx-\\frac12(b+a)\\frac12(b-a) \\right)^2
ρ
(x)=ρ(y=21(b−a)x−21(b+a))=1−(21(b−a)x−21(b+a))21
三项递推关系
记
θ
=
a
r
c
c
o
s
y
\\theta=\\mathrmarccosy
θ=arccosy,利用和差化积公式有
cos
(
(
n
+
1
)
θ
)
+
cos
(
(
n
−
1
)
θ
)
=
2
cos
(
n
θ
)
cos
(
θ
)
\\cos((n+1)\\theta)+\\cos((n-1)\\theta)=2\\cos(n\\theta)\\cos(\\theta)
cos((n+1)θ)+cos((n−1)θ)=2cos(nθ)cos(θ)
故有
T
n
+
1
(
y
)
=
2
y
T
n
(
y
)
−
T
n
−
1
(
y
)
T_n+1(y)=2yT_n(y)-T_n-1(y)
Tn(Python)Markov,Chebyshev,Chernoff上界函数