【数学 函数空间】格拉姆-施密特正交化&切比雪夫多项式推导

  • 正交基切比雪夫多项式是线性代数和数值分析中非常重要的概念。接下来,我将解释如何构造一组正交基,以及如何使用切比雪夫多项式来构造正交基。
    切比雪夫多项式有两类:第一类第二类。以下为第一类切比雪夫多项式 T n ( x ) T_n(x) Tn(x)的推导。

一、正交基的定义与构造

在向量空间 V V V 中,如果一组向量 { v 1 , v 2 , . . . , v n } \{v_1, v_2, ..., v_n\} {v1,v2,...,vn} 彼此正交且非零,那么它们就构成 V V V 的一组正交基。若它们的范数均为 1,则称为标准正交基(即正交归一基)。

1. 通过格拉姆-施密特(Gram-Schmidt)正交化过程构造正交基

假设我们有一组线性无关的向量 { u 1 , u 2 , . . . , u n } \{u_1, u_2, ..., u_n\} {u1,u2,...,un},可以使用格拉姆-施密特方法将其变换为一组正交基 { v 1 , v 2 , . . . , v n } \{v_1, v_2, ..., v_n\} {v1,v2,...,vn}

步骤:

  1. 定义 v 1 = u 1 v_1 = u_1 v1=u1
  2. 对于 i > 1 i > 1 i>1,计算:
    v i = u i − ∑ j = 1 i − 1 ⟨ u i , v j ⟩ ⟨ v j , v j ⟩ v j v_i = u_i - \sum_{j=1}^{i-1} \frac{\langle u_i, v_j \rangle}{\langle v_j, v_j \rangle} v_j vi=uij=1i1vj,vjui,vjvj
  3. 归一化向量 v i v_i vi 以得到单位向量:
    e i = v i ∥ v i ∥ e_i = \frac{v_i}{\|v_i\|} ei=vivi

这将生成一组正交基 { e 1 , e 2 , . . . , e n } \{e_1, e_2, ..., e_n\} {e1,e2,...,en}

二、格拉姆-施密特正交化得到切比雪夫多项式

  • 可通过格拉姆-施密特正交化构造前两个(1, x x x)的正交基函数:

g 0 ( x ) = 1 2 , g 1 ( x ) = 3 2 x . g_0(x) = \frac{1}{\sqrt{2}}, \quad g_1(x) = \sqrt{\frac{3}{2}} x. g0(x)=2 1,g1(x)=23 x.

  • 现在,我们将继续推导出第三个正交基函数 g 2 ( x ) g_2(x) g2(x),从而构造出正交基的第三项。

  • 我们希望从初始函数集 { 1 , x , x 2 } \{1, x, x^2\} {1,x,x2} 中构造出正交的第三个函数 g 2 ( x ) g_2(x) g2(x),使得:
    ⟨ g 2 , g 0 ⟩ = 0 , ⟨ g 2 , g 1 ⟩ = 0. \langle g_2, g_0 \rangle = 0, \quad \langle g_2, g_1 \rangle = 0. g2,g0=0,g2,g1=0.

  • 我们需要将 f 2 ( x ) = x 2 f_2(x) = x^2 f2(x)=x2 正交化,以得到 g 2 ( x ) g_2(x) g2(x)

1. 从 f 2 ( x ) f_2(x) f2(x)中减去其在 g 0 ( x ) g_0(x) g0(x) g 1 ( x ) g_1(x) g1(x)上的投影

我们需要从 f 2 ( x ) f_2(x) f2(x) 中减去它在前两个正交基 g 0 ( x ) g_0(x) g0(x) g 1 ( x ) g_1(x) g1(x) 上的投影:
g 2 ( x ) = f 2 ( x ) − Proj g 0 ( f 2 ) − Proj g 1 ( f 2 ) . g_2(x) = f_2(x) - \text{Proj}_{g_0}(f_2) - \text{Proj}_{g_1}(f_2). g2(x)=f2(x)Projg0(f2)Projg1(f2).

首先,我们计算投影部分。

2. 计算投影到 g 0 ( x ) g_0(x) g0(x) 上的分量

Proj g 0 ( f 2 ) = ⟨ f 2 , g 0 ⟩ g 0 . \text{Proj}_{g_0}(f_2) = \langle f_2, g_0 \rangle g_0. Projg0(f2)=f2,g0g0.

内积 ⟨ f 2 , g 0 ⟩ \langle f_2, g_0 \rangle f2,g0 计算如下:
⟨ f 2 , g 0 ⟩ = ∫ − 1 1 x 2 ⋅ 1 2   d x = 1 2 ∫ − 1 1 x 2   d x . \langle f_2, g_0 \rangle = \int_{-1}^{1} x^2 \cdot \frac{1}{\sqrt{2}} \, dx = \frac{1}{\sqrt{2}} \int_{-1}^{1} x^2 \, dx. f2,g0=11x22 1dx=2 111x2dx.

计算积分:
∫ − 1 1 x 2   d x = [ x 3 3 ] − 1 1 = 2 3 . \int_{-1}^{1} x^2 \, dx = \left[\frac{x^3}{3}\right]_{-1}^{1} = \frac{2}{3}. 11x2dx=[3x3]11=32.
因此,
⟨ f 2 , g 0 ⟩ = 1 2 ⋅ 2 3 = 2 3 . \langle f_2, g_0 \rangle = \frac{1}{\sqrt{2}} \cdot \frac{2}{3} = \frac{\sqrt{2}}{3}. f2,g0=2 132=32 .
于是,
Proj g 0 ( f 2 ) = 2 3 ⋅ 1 2 = 1 3 . \text{Proj}_{g_0}(f_2) = \frac{\sqrt{2}}{3} \cdot \frac{1}{\sqrt{2}} = \frac{1}{3}. Projg0(f2)=32 2 1=31.

3. 计算投影到 g 1 ( x ) g_1(x) g1(x) 上的分量

Proj g 1 ( f 2 ) = ⟨ f 2 , g 1 ⟩ g 1 . \text{Proj}_{g_1}(f_2) = \langle f_2, g_1 \rangle g_1. Projg1(f2)=f2,g1g1.

内积 ⟨ f 2 , g 1 ⟩ \langle f_2, g_1 \rangle f2,g1 计算如下:
⟨ f 2 , g 1 ⟩ = ∫ − 1 1 x 2 ⋅ 3 2 x   d x = 3 2 ∫ − 1 1 x 3   d x . \langle f_2, g_1 \rangle = \int_{-1}^{1} x^2 \cdot \sqrt{\frac{3}{2}} x \, dx = \sqrt{\frac{3}{2}} \int_{-1}^{1} x^3 \, dx. f2,g1=11x223 xdx=23 11x3dx.

计算积分:
∫ − 1 1 x 3   d x = [ x 4 4 ] − 1 1 = 0. \int_{-1}^{1} x^3 \, dx = \left[\frac{x^4}{4}\right]_{-1}^{1} = 0. 11x3dx=[4x4]11=0.
因此,
⟨ f 2 , g 1 ⟩ = 3 2 ⋅ 0 = 0 , \langle f_2, g_1 \rangle = \sqrt{\frac{3}{2}} \cdot 0 = 0, f2,g1=23 0=0,

Proj g 1 ( f 2 ) = 0. \text{Proj}_{g_1}(f_2) = 0. Projg1(f2)=0.

4. 得到第三个正交基函数

将上面的结果代入,我们得到:
g 2 ( x ) = f 2 ( x ) − Proj g 0 ( f 2 ) − Proj g 1 ( f 2 ) = x 2 − 1 3 . g_2(x) = f_2(x) - \text{Proj}_{g_0}(f_2) - \text{Proj}_{g_1}(f_2) = x^2 - \frac{1}{3}. g2(x)=f2(x)Projg0(f2)Projg1(f2)=x231.

5. 标准化 g 2 ( x ) g_2(x) g2(x)

为了得到标准化的正交基,我们计算 g 2 ( x ) g_2(x) g2(x) 的范数:
∥ g 2 ( x ) ∥ = ∫ − 1 1 ( x 2 − 1 3 ) 2   d x . \|g_2(x)\| = \sqrt{\int_{-1}^{1} \left(x^2 - \frac{1}{3}\right)^2 \, dx}. g2(x)=11(x231)2dx .

展开并计算积分:
( x 2 − 1 3 ) 2 = x 4 − 2 3 x 2 + 1 9 , \left(x^2 - \frac{1}{3}\right)^2 = x^4 - \frac{2}{3}x^2 + \frac{1}{9}, (x231)2=x432x2+91,
因此,
∥ g 2 ( x ) ∥ 2 = ∫ − 1 1 ( x 4 − 2 3 x 2 + 1 9 )   d x . \|g_2(x)\|^2 = \int_{-1}^{1} \left(x^4 - \frac{2}{3}x^2 + \frac{1}{9}\right) \, dx. g2(x)2=11(x432x2+91)dx.

逐项计算积分:
∫ − 1 1 x 4   d x = 2 5 , ∫ − 1 1 x 2   d x = 2 3 , ∫ − 1 1 1 9   d x = 2 9 . \int_{-1}^{1} x^4 \, dx = \frac{2}{5}, \quad \int_{-1}^{1} x^2 \, dx = \frac{2}{3}, \quad \int_{-1}^{1} \frac{1}{9} \, dx = \frac{2}{9}. 11x4dx=52,11x2dx=32,1191dx=92.

因此,
∥ g 2 ( x ) ∥ 2 = 2 5 − 4 9 + 2 9 = 8 45 , \|g_2(x)\|^2 = \frac{2}{5} - \frac{4}{9} + \frac{2}{9} = \frac{8}{45}, g2(x)2=5294+92=458,
∥ g 2 ( x ) ∥ = 8 45 = 2 10 15 . \|g_2(x)\| = \sqrt{\frac{8}{45}} = \frac{2\sqrt{10}}{15}. g2(x)=458 =15210 .

最终得到标准化的 g 2 ( x ) g_2(x) g2(x)
g 2 ( x ) = x 2 − 1 3 2 10 15 = 15 2 10 ( x 2 − 1 3 ) = 3 10 4 ( x 2 − 1 3 ) . g_2(x) = \frac{x^2 - \frac{1}{3}}{\frac{2\sqrt{10}}{15}} = \frac{15}{2\sqrt{10}} \left(x^2 - \frac{1}{3}\right) = \frac{3\sqrt{10}}{4} \left(x^2 - \frac{1}{3}\right). g2(x)=15210 x231=210 15(x231)=4310 (x231).

三、结果

  • 通过格拉姆-施密特正交化,我们得到了第三个正交基函数:
    g 2 ( x ) = 3 10 4 ( x 2 − 1 3 ) . g_2(x) = \frac{3\sqrt{10}}{4} \left(x^2 - \frac{1}{3}\right). g2(x)=4310 (x231).

  • 至此,我们构造了一组在区间 [ − 1 , 1 ] [-1, 1] [1,1] 上的正交基:
    g 0 ( x ) = 1 2 , g 1 ( x ) = 3 2 x , g 2 ( x ) = 3 10 4 ( x 2 − 1 3 ) . g_0(x) = \frac{1}{\sqrt{2}}, \quad g_1(x) = \sqrt{\frac{3}{2}} x, \quad g_2(x) = \frac{3\sqrt{10}}{4} \left(x^2 - \frac{1}{3}\right). g0(x)=2 1,g1(x)=23 x,g2(x)=4310 (x231).

四、切比雪夫多项式性质

切比雪夫多项式是定义在区间 ([-1, 1]) 上的一组正交多项式。它们在数值分析、逼近理论和求解微分方程中有广泛应用。

1. 递推性

递归定义:
T 0 ( x ) = 1 , T 1 ( x ) = x T_0(x) = 1, \quad T_1(x) = x T0(x)=1,T1(x)=x
T n + 1 ( x ) = 2 x T n ( x ) − T n − 1 ( x ) for  n ≥ 1 T_{n+1}(x) = 2x T_n(x) - T_{n-1}(x) \quad \text{for } n \geq 1 Tn+1(x)=2xTn(x)Tn1(x)for n1

  • 利用递归关系来推导前几个切比雪夫多项式:

    • T 0 ( x ) = 1 T_0(x) = 1 T0(x)=1
    • T 1 ( x ) = x T_1(x) = x T1(x)=x
    • T 2 ( x ) = 2 x 2 − 1 T_2(x) = 2x^2 - 1 T2(x)=2x21
      T 2 ( x ) = 2 x T 1 ( x ) − T 0 ( x ) = 2 x 2 − 1 T_2(x) = 2x T_1(x) - T_0(x) = 2x^2 - 1 T2(x)=2xT1(x)T0(x)=2x21
    • T 3 ( x ) = 4 x 3 − 3 x T_3(x) = 4x^3 - 3x T3(x)=4x33x
      T 3 ( x ) = 2 x T 2 ( x ) − T 1 ( x ) = 4 x 3 − 3 x T_3(x) = 2x T_2(x) - T_1(x) = 4x^3 - 3x T3(x)=2xT2(x)T1(x)=4x33x
    • T 4 ( x ) = 8 x 4 − 8 x 2 + 1 T_4(x) = 8x^4 - 8x^2 + 1 T4(x)=8x48x2+1
      T 4 ( x ) = 2 x T 3 ( x ) − T 2 ( x ) = 8 x 4 − 8 x 2 + 1 T_4(x) = 2x T_3(x) - T_2(x) = 8x^4 - 8x^2 + 1 T4(x)=2xT3(x)T2(x)=8x48x2+1

显式公式:
T n ( ω ) = cos ⁡ ( n arccos ⁡ ( ω ) ) = cos ⁡ ( n cos ⁡ − 1 ( ω ) ) T_n(\omega) = \cos(n \arccos(\omega))= \cos(n \cos^{-1}(\omega)) Tn(ω)=cos(narccos(ω))=cos(ncos1(ω))

  • x = c o s − 1 ω x=cos^{-1}\omega x=cos1ω,则 T n ( ω ) = c o s ( n x ) T_n(\omega)=cos(nx) Tn(ω)=cos(nx)
  • 通过三角恒等式展开 c o s ( n x ) cos(nx) cos(nx),可以发现其与 T n ( ω ) T_n(\omega) Tn(ω)为恒等关系
  • 例如:
    • c o s ( 0 ∗ x ) = 1 cos(0*x)= 1 cos(0x)=1
    • c o s ( x ) = c o s ( c o s − 1 ω ) = ω cos(x)= cos(cos^{-1}\omega) = \omega cos(x)=cos(cos1ω)=ω
    • c o s ( 2 x ) = 2 c o s x − 1 = 2 ω 2 − 1 cos(2x)= 2cos^x-1 = 2\omega^2-1 cos(2x)=2cosx1=2ω21
2. 正交性

切比雪夫多项式在区间 ([-1, 1]) 上关于权函数 w ( x ) = 1 1 − x 2 w(x) = \frac{1}{\sqrt{1 - x^2}} w(x)=1x2 1 是正交的,即:
∫ − 1 1 T m ( x ) T n ( x ) w ( x )   d x = { 0 , m ≠ n π 2 , m = n ≠ 0 π , n = 0 \int_{-1}^{1} T_m(x) T_n(x) w(x) \, dx = \begin{cases} 0, & m \neq n \\ \frac{\pi}{2}, & m = n \neq 0 \\ \pi, & n = 0 \end{cases} 11Tm(x)Tn(x)w(x)dx= 0,2π,π,m=nm=n=0n=0

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值