对应教材 5.5.1–5.5.2、5.6.1–5.6.3、定理 5.6–5.7。老师课堂未讲,按教材补齐,标“教材有、课堂未讲”。
主线:用正交多项式作基,把“解正规方程”变成“算内积”;再讨论最佳一致/平方逼近。
设权函数 ρ(x)≥0,若
(φi,φj)=∫abρ(x)φi(x)φj(x)dx=0(i=j),
称 {φk} 在 [a,b] 上带权 ρ 正交。
由幂函数系 1,x,x2,… 经 Gram–Schmidt 正交化,可得正交多项式系。常用:
- Legendre 多项式:[−1,1],ρ=1;
- Chebyshev 多项式:[−1,1],ρ=1/1−x2;
- Laguerre 多项式:[0,∞),ρ=e−x;
- Hermite 多项式:(−∞,∞),ρ=e−x2。
定理 5.6(正交多项式基本性质):n 次正交多项式 pn 在 (a,b) 内有 n 个互异实根;正交多项式系线性无关;任意次数不超过 n 的多项式可用 p0,…,pn 线性表示。
三项递推:正交多项式满足
pk+1(x)=(x−αk)pk(x)−βkpk−1(x),
这是数值生成正交多项式、构造 Gauss 求积节点的基础。
在 C[a,b] 中,用 n 次多项式 pn 逼近 f,使
∥f−pn∥∞=a≤x≤bmax∣f(x)−pn(x)∣
最小。Chebyshev 交替定理给出最佳一致逼近多项式的特征:误差在 n+2 个点交替取 ±∥f−pn∥∞。实际计算常用 Chebyshev 展开/Remez 算法。
在 Lρ2(a,b) 中,使
∥f−p∥22=∫abρ(x)∣f(x)−p(x)∣2dx
最小。若 {φk} 正交,则
p(x)=k=0∑nckφk(x),ck=(φk,φk)(f,φk).
系数直接由内积给出,不用解线性方程组;这就是“换到正交基”的最大好处。Legendre/Chebyshev 展开是最常用的两种。
定理 5.7(最佳平方逼近的刻画):设 H=operatornameSpanarphi0,dots,arphin,p∗inH。则 p∗ 是 f 在 H 中的最佳平方逼近元,当且仅当
(f−p∗,arphik)=0,qquadk=0,dots,n.
等价地,系数 ck 满足法方程
sumj=0n(arphij,arphik)cj=(f,arphik),qquadk=0,dots,n.
若 arphik 线性无关,则解唯一。
证明要点:对任意 qinH,展开
lVertf−p∗−qVert22=lVertf−p∗Vert22−2(f−p∗,q)+lVertqVert22.
若 (f−p∗,arphik)=0,则 (f−p∗,q)=0,故 p∗ 最小;反之若某个内积不为零,取 q 沿该方向扰动可使误差更小。■
覆盖核对
本节覆盖教材 5.5.1–5.5.2、5.6.1–5.6.3、定理 5.6–5.7;老师课堂未讲。下一节:第 5 章习题。