拉盖尔多项式

来自testwiki
跳转到导航 跳转到搜索

Template:NoteTA

数学中,以法国数学家Template:Link-en命名的拉盖尔多项式定义为拉盖尔方程的标准解。

xy+(1x)y+ny=0

这是一个二阶线性微分方程

这个方程只有当n非负时,才有非奇异解。拉盖尔多项式可用在高斯积分法中,计算形如0f(x)dx的积分。

这些多项式(通常用L0L1等表示)构成一个多项式序列。这个多项式序列可以用罗德里格公式递推得到。

Ln(x)=exn!dndxn(exxn).

在按照下式定义的内积构成的内积空间中,拉盖尔多项式是正交多项式

f,g=0f(x)g(x)exdx.

拉盖尔多项式构成一个Template:Link-en

拉盖尔多项式在量子力学中有重要应用。氢原子薛定谔方程的解的径向部分,就是拉盖尔多项式。

物理学家通常采用另外一种拉盖尔多项式的定义形式,即在上面的形式的基础上乘上一个n!。

前几个拉盖尔多项式

前几个拉盖尔多项式的表达式与函数图像如下:

n Ln(x)
0 1
1 x+1
2 12(x24x+2)
3 16(x3+9x218x+6)
4 124(x416x3+72x296x+24)
5 1120(x5+25x4200x3+600x2600x+120)
6 1720(x636x5+450x42400x3+5400x24320x+720)
前六个拉盖尔多项式

递归定义

拉盖尔多项式也可以通过递归的方式进行定义。首先,规定前两个拉盖尔多项式为:

L0(x)=1
L1(x)=1x

然后运用下面的递推关系得到更高阶的多项式。

Lk+1(x)=1k+1((2k+1x)Lk(x)kLk1(x)).

广义拉盖尔多项式

上面提到的拉盖尔多项式的正交性,也可以用另外一种方式表达。即:如果X是一个服从指数分布随机变量(即,概率密度函数如下式):

f(x)={exif x>0,0if x<0,

那么:

E[Ln(X)Lm(X)]=0 whenever nm.

指数分布不是唯一的伽玛分布,对于任意的伽玛分布(概率密度函数如下,α > −1,参见Γ函数

f(x)={xαex/Γ(1+α)if x>0,0if x<0,

相应的正交多项式为形如下式的广义拉盖尔多项式(可以通过罗德里格公式得到):

Ln(α)(x)=xαexn!dndxn(exxn+α).

有时也将上面的多项式称为连带(联属,伴随)拉盖尔多项式。当取α = 0时,就回到拉盖尔多项式:

Ln(0)(x)=Ln(x).

广义拉盖尔多项式的性质与应用

  • 拉盖尔函数可以由合流超几何函数和Kummer变换得到: Ln(α)(x):=(n+αn)M(n,α+1,x)=(n+αn)i=0(1)i(ni)(α+ii)xi =ex(n+αn)M(α+n+1,α+1,x)=exsin(nπ)sin((n+α)π)Lαn1(α)(x)=exi=0(1)i(α+n+in)xii!.n为整数时,截断为n阶拉盖尔多项式。
  • n阶拉盖尔多项式可以通过将Template:Link-en应用在罗德里格公式上而得到,结果为Ln(α)(x)=i=0n(1)i(n+αni)xii!
  • n阶拉盖尔多项式的首项系数为(−1)n/n!;
  • 拉盖尔多项式在x=0的取值(常数项)为 Ln(α)(0)=(n+αn)nαΓ(α+1);
  • Ln(α)n的正(应该注意到 ((1)niLni(α))i=0n构成以施图姆序列),且这些根全部位于区间(0,n+α+(n1)n+α]中。
  • n很大,而α不变,x>0时,拉盖尔多项式的渐近行为如下:
Ln(α)(x)nα214πex2xα2+14cos(2x(n+α+12)π2(α+12)),以及
Ln(α)(x)nα2142πex2xα2+14exp(2x(n+α+12))[1]
  • 前几个广义拉盖尔多项式为:
L0(α)(x)=1
L1(α)(x)=x+α+1
L2(α)(x)=x22(α+2)x+(α+2)(α+1)2
L3(α)(x)=x36+(α+3)x22(α+2)(α+3)x2+(α+1)(α+2)(α+3)6
  • 根据拉盖尔多项式的定义,可以使用秦九韶算法计算拉盖尔多项式,程序代码如下:
 function LaguerreL(n, alpha, x) {
    LaguerreL:= 1; bin:= 1 
    for i:= n to 1 step -1 {
        bin:= bin* (alpha+ i)/ (n+ 1- i)
        LaguerreL:= bin- x* LaguerreL/ i
    }
    return LaguerreL;
 }

递推关系

拉盖尔多项式满足以下的递推关系:

Ln(α+β+1)(x+y)=i=0nLi(α)(x)Lni(β)(y),

特别地,有

Ln(α+1)(x)=i=0nLi(α)(x)以及Ln(α)(x)=i=0n(αβ+ni1ni)Li(β)(x),或Ln(α)(x)=i=0n(αβ+nni)Li(βi)(x);

还有

Ln(α)(x)j=0Δ1(n+αnj)(1)jxjj!=(1)ΔxΔ(Δ1)!i=0nΔ(n+αnΔi)(ni)(ni)Li(α+Δ)(x)=(1)ΔxΔ(Δ1)!i=0nΔ(n+αi1nΔi)(ni)(ni)Li(n+α+Δi)(x).

运用以上式子可以得到以下四条关系式:

  • Ln(α)(x)=Ln(α+1)(x)Ln1(α+1)(x)=j=0k(kj)Lnj(αk+j)(x),
  • nLn(α)(x)=(n+α)Ln1(α)(x)xLn1(α+1)(x), or xkk!Ln(α)(x)=i=0k(1)i(n+ii)(n+αki)Ln+i(αk)(x),
  • nLn(α+1)(x)=(nx)Ln1(α+1)(x)+(n+α)Ln1(α)(x)
  • xLn(α+1)=(n+α)Ln1α(x)(nx)Ln(α)(x);

将它们组合在一起,就得到了最常用的递推关系式:

Ln+1(α)(x)=1n+1((2n+1+αx)Ln(α)(x)(n+α)Ln1(α)(x)).

in均为整数时,拉盖尔多项式有以下的有趣性质:

(x)ii!Ln(in)(x)=(x)nn!Li(ni)(x);

进一步可以得到部分分式分解

Ln(α)(x)(n+αn)=1j=1nxjα+jLnj(j)(x)(j1)!=1xi=1nLni(α)(x)Li1(α+1)(x)α+i.

拉盖尔多项式的导函数

将拉盖尔多项式对自变量x求导k次,得到:

dkdxkLn(α)(x)=(1)kLnk(α+k)(x);

进一步有:

1k!dkdxkxαLn(α)(x)=(n+αk)xαkLn(αk)(x),

运用Template:Link-en可以得到:

Ln(α)(x)=(αα)(α+nαα)0xtα(xt)αα1xαLn(α)(t)dt.

将拉盖尔多项式对参变量α求导,得到下面的有意思的结果:

ddαLn(α)(x)=i=0n1Li(α)(x)ni.

广义拉盖尔多项式满足下面的微分方程:

xLn(α)(x)+(α+1x)Ln(α)(x)+nLn(α)(x)=0,

可以与拉盖尔多项式的k阶导数所满足的微分方程作一比较。

xLn(k)(x)+(k+1x)Ln(k)(x)+(nk)Ln(k)(x)=0,

仅在此式中,Ln(k)(x)dLn(x)dxk(后面这个符号又有了新的含义)。

于是,当α=0时,广义拉盖尔多项式可以用拉盖尔多项式的导数表示: Ln(k)(x)=(1)kdLn+k(x)dxk 式中的上标(k)容易与求导k次混淆。

正交性

伴随拉盖尔多项式在区间[0, ∞)上以权函数xα e −x正交:

0xαexLn(α)(x)Lm(α)(x)dx=Γ(n+α+1)n!δn,m,

这可由下式得到:

0xα1exLn(α)(x)dx=(αα+nn)Γ(α).

伴随对称核多项式可以用拉盖尔多项式表示为:

Kn(α)(x,y):=1Γ(α+1)i=0nLi(α)(x)Li(α)(y)(α+ii)=1Γ(α+1)Ln(α)(x)Ln+1(α)(y)Ln+1(α)(x)Ln(α)(y)xyn+1(n+αn)=1Γ(α+1)i=0nxii!Lni(α+i)(x)Lni(α+i+1)(y)(α+nn)(ni);

也有下面的递推关系:

Kn(α)(x,y)=yα+1Kn1(α+1)(x,y)+1Γ(α+1)Ln(α+1)(x)Ln(α)(y)(α+nn).

进一步地,在伴L2[0, ∞)空间上,有:

yαeyKn(α)(,y)δ(y),

在氢原子的量子力学处理中用到了下面的公式:

0xα+1ex[Ln(α)]2dx=(n+α)!n!(2n+α+1).

级数展开

设一个函数具有以下的级数展开形式:

f(x)=i=0fi(α)Li(α)(x).

则展开式的系数由下式给出

fi(α)=0Li(α)(x)(i+αi)xαexΓ(α+1)f(x)dx.

这个级数在Lp空间L2[0,)上收敛,当且仅当

fL22:=0xαexΓ(α+1)|f(x)|2dx=i=0(i+αi)|fi(α)|2<.

一个相关的展开式为:

f(x)=eγ1+γxi=0Li(α)(x1+γ)(1+γ)i+α+1n=0iγin(in)fn(α);

特别地

eγxLn(α)(x(1+γ))=i=nLi(α)(x)(1+γ)i+α+1γin(in),

这可由下式得到:

Ln(α)(x1+γ)=1(1+γ)ni=0nγni(n+αni)Li(α)(x).

还有,当Re(2αβ)>1时,

xαβf(x)Γ(αβ+1)=(αβ)i=0Li(β)(x)(β+ii)n=0i(1)in(αβin)(α+nn)fn(α),

这个结果可以由下式导出,

xαβLn(α)(x)Γ(αβ+1)=(αβ)(α+nn)i=n(1)in(αβin)Li(β)(x)(β+ii)

更多的例子

幂函数可以展开为:

xnn!=i=0n(1)i(n+αni)Li(α)(x)=(1)ni=0nLi(αi)(x)(αni),

二项式可以展开为:

(n+xn)=i=0nαii!Lni(x+i)(α).

进一步可以得到:

eγx=i=0γi(1+γ)i+α+1Li(α)(x)(当且仅当 Re(γ)>12时收敛)

更一般地

xβeγxΓ(β+1)=(α+βα)i=0Li(α)(x)(α+ii)j=0i(1)j(1+γ)α+β+j+1(α+β+jj)(α+iij).

对于非负的整数β,可以化简为:

xneγxn!=i=0γiLi(α)(x)(1+γ)i+n+α+1j=0n(1)njγj(n+αj)(inj),

γ=0时,可以化简为:

xβΓ(β+1)=(α+βα)i=0(1)i(βi)Li(α)(x)(α+ii),
xβLn(γ)(x)Γ(β+1)=(α+βα)i=0Li(α)(x)(α+ii)j=0n(1)ij(n+γnj)(β+ji)(α+β+jj).

雅可比Theta 函数有下面的表示:

kek2πx=i=0Li(α)(xt)k(k2πt)i(1+k2πt)i+α+1;

随意选定参量t,贝塞尔函数可以表示为: Jα(x)(x2)α=etΓ(α+1)i=0Li(α)(x24t)(i+αi)tii!; Γ函数可以展开为:

Γ(α)=xαi=0Li(α)(x)α+i((α)<12);

低阶不完全伽玛函数可展开为:

γ(s;z)tsΓ(s)=(zt)αΓ(α+1)i=0Li(α)(zt)(α+ii)j=0i(1)j(1+t)s+j(s1+jj)(α1+iij),
γ(s;z)tsΓ(s)=(α+sα+1)i=0(α+i+1i+1)Li+1(α)(zt)(α+i+1i)j=0i(1)j(1+t)α+1+s+j(α+s+jj)(α+i+1ij).

还有:

γ(s,z)=γsΓ(1s)i=0Li+1(s)(0)Li+1(s)(zγ)(1+γ)i+1n=0iγin(in)n+1s;

于是,高阶不完全伽玛函数就是:

Γ(s,z)zsez=k=0Lk(α)(z)(k+1)(k+1+αsk+1)((sα2)<14)=k=0Lk(α)(zt)2F1(1+α+k,1+k;2+α+ks;t1t)tk(k+1)(1+α+ks1+k)=tsk=0Lk(α)(zt)2F1(1s,1+αs;2+α+ks;t1t)(k+1)(1+α+ks1+k)=t1+αk=0Lk(α)(zt)2F1(1+α+k,1+αs;2+α+ks;1t)(k+1)(1+α+ks1+k),

2F1表示超几何函数

围道积分表示

拉盖尔多项式可以用围道积分表示,如下式所示:

Ln(α)(x)=12πiext1t(1t)α+1tn+1dt

积分方向逆时针绕原点一周。

与埃爾米特多項式的关系

广义拉盖尔多项式与埃爾米特多項式有下列关系:

H2n(x)=(1)n 22n n! Ln(1/2)(x2)

以及

H2n+1(x)=(1)n 22n+1 n! x Ln(1/2)(x2)

这里的Hn表示乘上了exp(−x2)的埃爾米特多項式(所谓的“物理学家形式”)。 正因为这样,广义拉盖尔多项式也在量子谐振子的量子力学处理中出现。

与超几何函数的关系

拉盖尔多项式可以用超几何函数来定义,具体地说,是用合流超几何函数定义:

Ln(α)(x)=(n+αn)M(n,α+1,x)=(α+1)nn!1F1(n,α+1,x)

(a)n阶乘幂,这里表示升阶乘

与贝塞尔函数的关系

拉盖尔多项式与变形贝塞尔函数之间有以下关系:

Ln(α)(x)=ex2(x4)n+122π(n+1)!(12n+1)k=0n(1)k+1(2n+1nk)(n+αn)(α+2n+1nk)(nk+αnk)(k+12)Kk+12(x2)=ex2(4x)n+α+12Γ(α+12)(α1n)(α12n)n!k=0n(2n12αkn)(2n1αkn)(α1kn)(α+12+k)Iα+12+k(x2),

进一步有:

Ln(α)(x)=24n(2n+1)(12n)k=0n(k+12)(2n+1nk)(nk)2(n+αk)(2n+α+1nk)xnk(nk)!Lk2k1(x).

外部链接

注释

参考文献

  • Abramowitz, Milton; Stegun, Irene A., eds. (1965), "Chapter 22 Template:Wayback", Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, New York: Dover, ISBN 0-486-61272-4 .
  • B Spain, M G Smith, Functions of mathematical physics, Van Nostrand Reinhold Company, London, 1970. Chapter 10 deals with Laguerre polynomials.
  • Eric W. Weisstein, "Laguerre Polynomial Template:Wayback", From MathWorld—A Wolfram Web Resource.
  • Template:Cite book
  • S. S. Bayin (2006), Mathematical Methods in Science and Engineering, Wiley, Chapter 3.