Jacobi正交多项式

注：本文的内容主要根据文末中的参考文档[1,2,3]中的内容进行整理完成。

---

一些基本概念

设$I=[-1,1]$是实轴上的标准区间，定义在$I$上的正函数：
$${\omega }_{\alpha ,\beta }(x)=(1-x{)}^{\alpha }(1+x{)}^{\beta },\alpha >-1,\beta >-1$$
是权函数。
赋权的Sobolev空间记为$H_{{\omega }_{\alpha ,\beta }}^r(I)$。

• 当$r=0$时，$H_{{\omega }_{\alpha ,\beta }}^0(I)=L_{{\omega }_{\alpha ,\beta }}^2(I)$ 是平方可积的函数空间，其上的内积和范数定义为：
  $$\begin{gathered} ⟨u,v{⟩}_{{\omega }_{\alpha ,\beta }}={\int }_{I}u(x)v(x){\omega }_{\alpha ,\beta }(x)\mathrm{d}(x), \\ \Vert u{\Vert }_{{\omega }_{\alpha ,\beta }}=\sqrt{⟨u,u{⟩}_{{\omega }_{\alpha ,\beta }}}. \end{gathered}$$
• 当$r>0$且为整数时，
  $$H_{{\omega }_{\alpha ,\beta }}^r(I)=\left\{u\in L_{{\omega }_{\alpha ,\beta }}^2(I):\frac{{\mathrm{d}}^{k}u}{\mathrm{d}{x}^k}\in L_{{\omega }_{\alpha ,\beta }}^2(I),0\le k\le r\right\},$$
  它的内积，范数，半范数分别为：
  $$\begin{gathered} ⟨u,v{⟩}_{r,{\omega }_{\alpha ,\beta }}=\sum _{k=0}^r{\left(\frac{{\mathrm{d}}^{k}u}{\mathrm{d}{x}^k},\frac{{\mathrm{d}}^{k}v}{\mathrm{d}{x}^k}\right)}_{{\omega }_{\alpha ,\beta }}, \\ \Vert u{\Vert }_{r,{\omega }_{\alpha ,\beta }}=\sqrt{⟨u,u{⟩}_{r,{\omega }_{\alpha ,\beta }}}, \\ |u{|}_{r,{\omega }_{\alpha ,\beta }}=\sqrt{{\left(\frac{{\mathrm{d}}^{k}u}{\mathrm{d}{x}^k},\frac{{\mathrm{d}}^{k}v}{\mathrm{d}{x}^k}\right)}_{{\omega }_{\alpha ,\beta }}}. \end{gathered}$$
• 当$r>0$且为实数时，我们可以通过内插空间[4]来定义$H_{{\omega }_{\alpha ,\beta }}^r(I)$.

---

Jacobi 正交多项式的定义

本节给出Jacobi正交多项式的三种不同的定义[1]。

• 定义1： Jacobi正交多项式${\left\{J_n^{\alpha ,\beta }(x)\right\}}_{n=0}^{\infty }$可以通过正交化代数多项式基底$\left\{1,x,x^2,\cdots ,x^n,\cdots \right\}$得到，这里的正交化是在内积空间$L_{{\omega }_{\alpha ,\beta }}^2(I)$中进行的。我们称$J_n^{\alpha ,\beta }(x)$为$n$次Jacobi多项式。

  • 当$\alpha =\beta =-1/2$时，$J_n^{-\frac{1}{2},-\frac{1}{2}}$(x)是第一类Chebyshev多项式；
  • 当$\alpha =\beta =1/2$时，$J_n^{\frac{1}{2},\frac{1}{2}}(x)$是第二类Chebyshev多项式；
  • 当$\alpha =\beta =0$时，$J_n^{0,0}(x)$是Legendre多项式；
  • 当$\alpha =\beta$时，$J_n^{\alpha ,\alpha }(x)$也叫做Gegenbauer多项式或Ultraspherical多项式（特种球多项式或者超球多项式）。

• 定义2：Jacobi正交多项式$J_n^{\alpha ,\beta }(x)$是下面奇异Sturm-Liouville问题的特征函数
  $$\frac{\mathrm{d}}{\mathrm{d}x}\left({\omega }_{\alpha ,\beta }(x)\frac{\mathrm{d}}{\mathrm{d}x}u(x)\right)+{\lambda }_n^{\alpha ,\beta }{\omega }_{\alpha ,\beta }(x)u(x)=0,$$
  其中特征值为：
  $${\lambda }_n^{\alpha ,\beta }=n(n+\alpha +\beta +1),n\ge 0,\alpha ,\beta >-1.$$

• 定义3：Jacobi正交多项式$J_n^{\alpha ,\beta }(x)$可由下面的关系式生成：
  $${\omega }_{\alpha ,\beta }(x)J_n^{\alpha ,\beta }(x)=\frac{(-1{)}^n}{2^{n}n!}\frac{{\mathrm{d}}^n}{\mathrm{d}{x}^n}{\omega }_{n+\alpha ,n+\beta }(x),$$
  其中 $n\ge 0,\alpha ,\beta >-1$,
  或者下面更一般的关系式生成：对$n\ge m\ge 0,\alpha ,\beta >-1$,
  $${\omega }_{\alpha ,\beta }(x)J_n^{\alpha ,\beta }(x)=\frac{(-1{)}^m(n-m)!}{2^{n}n!}\frac{{\mathrm{d}}^m}{\mathrm{d}{x}^m}{\omega }_{m+\alpha ,m+\beta }(x)J_{n-m}^{m+\alpha ,m+\beta }(x),$$
  第一个零次Jacobi多项式定义为$J_0^{\alpha ,\beta }(x)=1$.

注：可以证明，关于Jacobi正交多项式的上面三个定义是相互等价的[5].

---

证明[-1,1]上的正交多项式是Jacobi正交多项式

（注：本小节主要说明[-1,1]上的正交多项式是Jacobi正交多项式，具体可以查看文献[3]中详细介绍。）
假设 { φ i ( t ) } \{\varphi_{i}(t)\} {φi​(t)}是$[-1,1]$上的一组关于权函数$\omega (t)$正交多项式。那么，${\phi }_n(t)$一定与不超过$n-1$次的多项式正交。这是由于${\phi }_n(t)$垂直于${\phi }_0(t),\cdots ,{\phi }_{n-1}(t)$所张成的空间，而任意不超过$n-1$次多项式都是这个空间的元素（可由${\phi }_0(t),\cdots ,{\phi }_{n-1}(t)$线性表出）。因此
$${\int }_{-1}^1\omega (t){\phi }_n(t)q_{n-1}(t)dt=0,\text{\hspace{1em}(1)}$$
令$r_n(t)$是$\omega (t){\phi }_n(t)$的第$n$重积分，也即
$$r_n^{(n)}(t)=\omega (t){\phi }_n(t).$$
利用分部积分法，对(1)进行求解
$${\left[r_n^{(n-1)}(t)q_{n-1}(t)-r_n^{(n-2)}(t)q_{n-1}^{\prime }(t)+\cdots +(-1{)}^{n-1}{r}_n(t)q_{n-1}^{n-1}(t)\right]}_{-1}^1+(-1{)}^n{\int }_{-1}^1{r}_n(t)q_{n-1}^{(n)}(t)\mathrm{d}t=0,$$
由于$q_{n-1}^{(n)}(t)=0$，于是有
$${\left[r_n^{(n-1)}(t)q_{n-1}(t)-r_n^{(n-2)}(t)q_{n-1}^{\prime }(t)+\cdots +(-1{)}^{n-1}{r}_n(t)q_{n-1}^{n-1}(t)\right]}_{-1}^1=0$$
由于$q_{n-1}(t)$是任意的多项式，因此
$$r_n^{(n-1)}(\pm 1)=r_n^{(n-2)}(\pm 1)=\cdots =r_n^{\prime }(\pm 1)=r_n(\pm 1)=0，\text{\hspace{1em}(2)}$$
这是由于如果多项式恒为0，那么多项式的系数一定为0。观察(2)式，对$r_n(t)$求导一直到第$n-1$次导数在$t$取$1$和$-1$时仍然都等于0，所以从第0次到第$n-1$次导函数中都含有$(1-t^2)$这个因子，更进一步地，$r_n(t)$含有一个$2n$次多项式的因子，使得它在$t$取$1$和$-1$时都等于零。设想：
$$r_n(t)=c_n(1-t^2{)}^{n+\alpha },(\alpha >-1)$$
可以验证$r_n(t)$满足(12)式。
由于$r_n(t)$的定义知道：$r_n^{(n)}(t)=\omega (t){\phi }_n(t)$, 也即
$${\phi }_n(t)=\omega (t{)}^{-1}{r}_n^{(n)}(t)={\phi }_n(t).$$
又由于${\phi }_n(t)$是$n$次多项式，因此为了满足$\omega (t{)}^{-1}{r}_n^{(n)}(t)$也是$n$次多项式，只需要选择
$$\omega (t)=(1-t^2{)}^{\alpha }.$$
于是，正交多项式${\phi }_n(t)$的表达式为
$${\phi }_n(t)=c_n\frac{1}{(1-t^2{)}^{\alpha }}\frac{{\mathrm{d}}^n}{\mathrm{d}{t}^n}(1-t^2{)}^{n+\alpha }.\text{\hspace{1em}(3)}$$
实际上，我们也可以取如下形式的$r_n(t)$:
$$r_n(t)=c_n(1-t{)}^{n+\alpha }(1+t{)}^{n+\beta },(\alpha >-1,\beta >-1).$$
对应的
$$\omega (t)=(1-t{)}^{\alpha }(1+t{)}^{\beta },$$
此时，正交多项式${\phi }_n(t)$的表达式为：
$${\phi }_n(t)=c_n\frac{1}{(1-t{)}^{\alpha }(1+t{)}^{\beta }}\frac{{\mathrm{d}}^n}{\mathrm{d}{t}^n}[(1-t{)}^{n+\alpha }(1+t{)}^{n+\beta }].\text{\hspace{1em}(4)}$$

---

参考文献：
[1] 赵廷刚，张建宏. Jacobi正交多项式的一些性质. 甘肃高师学报，14(5), 2009.
[2] Jacobi正交多项式种类. https://max.book118.com/html/2020/0815/5334111143002331.shtm
[3] 一般正交多项式性质. https://zhuanlan.zhihu.com/p/270620896.
[4] R. Askey. Orthogonal polynomials and special functions. Regional conderence seris in applied mathematics, SIAM, Phiadelphia, 1975.
[5] Guo Benyu and Wang Lilian. Jacobi interpolation approximations and their applications to singular differential equations. Advances in Computational Mathematics, 14: 227-276, 2000.