Template:NoteTA
在数学 中,以法国 数学家 Template:Link-en 命名的拉盖尔多项式 定义为拉盖尔方程 的标准解。
x y ″ + ( 1 − x ) y ′ + n y = 0
这是一个二阶线性微分方程 。
这个方程只有当n 非负时,才有非奇异解。拉盖尔多项式可用在高斯积分法 中,计算形如∫ 0 ∞ f ( x ) d x 的积分。
这些多项式(通常用L 0 , L 1 等表示)构成一个多项式序列 。这个多项式序列可以用罗德里格公式 递推得到。
L n ( x ) = e x n ! d n d x n ( e − x x n ) .
在按照下式定义的内积构成的内积空间 中,拉盖尔多项式是正交多项式 。
⟨ f , g ⟩ = ∫ 0 ∞ f ( x ) g ( x ) e − x d x .
拉盖尔多项式构成一个Template:Link-en 。
拉盖尔多项式在量子力学中有重要应用。氢原子薛定谔方程 的解的径向部分,就是拉盖尔多项式。
物理学家通常采用另外一种拉盖尔多项式的定义形式,即在上面的形式的基础上乘上一个n !。
前几个拉盖尔多项式
前几个拉盖尔多项式的表达式与函数图像如下:
n
L n ( x )
0
1
1
− x + 1
2
1 2 ( x 2 − 4 x + 2 )
3
1 6 ( − x 3 + 9 x 2 − 1 8 x + 6 )
4
1 2 4 ( x 4 − 1 6 x 3 + 7 2 x 2 − 9 6 x + 2 4 )
5
1 1 2 0 ( − x 5 + 2 5 x 4 − 2 0 0 x 3 + 6 0 0 x 2 − 6 0 0 x + 1 2 0 )
6
1 7 2 0 ( x 6 − 3 6 x 5 + 4 5 0 x 4 − 2 4 0 0 x 3 + 5 4 0 0 x 2 − 4 3 2 0 x + 7 2 0 )
前六个拉盖尔多项式
递归定义
拉盖尔多项式也可以通过递归的方式进行定义。首先,规定前两个拉盖尔多项式为:
L 0 ( x ) = 1
L 1 ( x ) = 1 − x
然后运用下面的递推关系 得到更高阶的多项式。
L k + 1 ( x ) = 1 k + 1 ( ( 2 k + 1 − x ) L k ( x ) − k L k − 1 ( x ) ) .
广义拉盖尔多项式
上面提到的拉盖尔多项式的正交性,也可以用另外一种方式表达。即:如果X 是一个服从指数分布 的随机变量 (即,概率密度函数 如下式):
f ( x ) = { e − x if x > 0 , 0 if x < 0 ,
那么:
E [ L n ( X ) L m ( X ) ] = 0 whenever n ≠ m .
指数分布不是唯一的伽玛分布 ,对于任意的伽玛分布(概率密度函数如下,α > −1,参见Γ函数 )
f ( x ) = { x α e − x / Γ ( 1 + α ) if x > 0 , 0 if x < 0 ,
相应的正交多项式为形如下式的广义拉盖尔多项式 (可以通过罗德里格公式 得到):
L n ( α ) ( x ) = x − α e x n ! d n d x n ( e − x x n + α ) .
有时也将上面的多项式称为连带(联属,伴随)拉盖尔多项式 。当取α = 0时,就回到拉盖尔多项式:
L n ( 0 ) ( x ) = L n ( x ) .
广义拉盖尔多项式的性质与应用
拉盖尔函数可以由合流超几何函数 和Kummer变换得到: L n ( α ) ( x ) : = ( n + α n ) M ( − n , α + 1 , x ) = ( n + α n ) ∑ i = 0 ( − 1 ) i ( n i ) ( α + i i ) x i = e x ⋅ ( n + α n ) M ( α + n + 1 , α + 1 , − x ) = e x sin ( n π ) sin ( ( n + α ) π ) L − α − n − 1 ( α ) ( − x ) = e x ⋅ ∑ i = 0 ( − 1 ) i ( α + n + i n ) x i i ! . 当n 为整数时,截断为n 阶拉盖尔多项式。
n 阶拉盖尔多项式可以通过将Template:Link-en 应用在罗德里格公式上而得到,结果为L n ( α ) ( x ) = ∑ i = 0 n ( − 1 ) i ( n + α n − i ) x i i ! 。
n 阶拉盖尔多项式的首项系数为(−1)n /n !;
拉盖尔多项式在x=0的取值(常数项 )为 L n ( α ) ( 0 ) = ( n + α n ) ≈ n α Γ ( α + 1 ) ;
L n (α ) 有n 个实 的正根 (应该注意到 ( ( − 1 ) n − i L n − i ( α ) ) i = 0 n 构成以施图姆序列 ),且这些根全部位于区间 ( 0 , n + α + ( n − 1 ) n + α ] 中。
当n 很大,而α 不变,x > 0 时,拉盖尔多项式的渐近行为如下:
L n ( α ) ( x ) ≈ n α 2 − 1 4 π e x 2 x α 2 + 1 4 cos ( 2 x ( n + α + 1 2 ) − π 2 ( α + 1 2 ) ) ,以及
L n ( α ) ( − x ) ≈ n α 2 − 1 4 2 π e − x 2 x α 2 + 1 4 exp ( 2 x ( n + α + 1 2 ) ) 。[ 1]
L 0 ( α ) ( x ) = 1
L 1 ( α ) ( x ) = − x + α + 1
L 2 ( α ) ( x ) = x 2 2 − ( α + 2 ) x + ( α + 2 ) ( α + 1 ) 2
L 3 ( α ) ( x ) = − x 3 6 + ( α + 3 ) x 2 2 − ( α + 2 ) ( α + 3 ) x 2 + ( α + 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;
}
递推关系
拉盖尔多项式满足以下的递推关系:
L n ( α + β + 1 ) ( x + y ) = ∑ i = 0 n L i ( α ) ( x ) L n − i ( β ) ( y ) ,
特别地,有
L n ( α + 1 ) ( x ) = ∑ i = 0 n L i ( α ) ( x ) 以及L n ( α ) ( x ) = ∑ i = 0 n ( α − β + n − i − 1 n − i ) L i ( β ) ( x ) ,或L n ( α ) ( x ) = ∑ i = 0 n ( α − β + n n − i ) L i ( β − i ) ( x ) ;
还有
L n ( α ) ( x ) − ∑ j = 0 Δ − 1 ( n + α n − j ) ( − 1 ) j x j j ! = ( − 1 ) Δ x Δ ( Δ − 1 ) ! ∑ i = 0 n − Δ ( n + α n − Δ − i ) ( n − i ) ( n i ) L i ( α + Δ ) ( x ) = ( − 1 ) Δ x Δ ( Δ − 1 ) ! ∑ i = 0 n − Δ ( n + α − i − 1 n − Δ − i ) ( n − i ) ( n i ) L i ( n + α + Δ − i ) ( x ) .
运用以上式子可以得到以下四条关系式:
L n ( α ) ( x ) = L n ( α + 1 ) ( x ) − L n − 1 ( α + 1 ) ( x ) = ∑ j = 0 k ( k j ) L n − j ( α − k + j ) ( x ) ,
n L n ( α ) ( x ) = ( n + α ) L n − 1 ( α ) ( x ) − x L n − 1 ( α + 1 ) ( x ) , or x k k ! L n ( α ) ( x ) = ∑ i = 0 k ( − 1 ) i ( n + i i ) ( n + α k − i ) L n + i ( α − k ) ( x ) ,
n L n ( α + 1 ) ( x ) = ( n − x ) L n − 1 ( α + 1 ) ( x ) + ( n + α ) L n − 1 ( α ) ( x )
x L n ( α + 1 ) = ( n + α ) L n − 1 α ( x ) − ( n − x ) L n ( α ) ( x ) ;
将它们组合在一起,就得到了最常用的递推关系式:
L n + 1 ( α ) ( x ) = 1 n + 1 ( ( 2 n + 1 + α − x ) L n ( α ) ( x ) − ( n + α ) L n − 1 ( α ) ( x ) ) .
当i 与n 均为整数时,拉盖尔多项式有以下的有趣性质:
( − x ) i i ! L n ( i − n ) ( x ) = ( − x ) n n ! L i ( n − i ) ( x ) ;
进一步可以得到部分分式分解 :
L n ( α ) ( x ) ( n + α n ) = 1 − ∑ j = 1 n x j α + j L n − j ( j ) ( x ) ( j − 1 ) ! = 1 − x ∑ i = 1 n L n − i ( − α ) ( x ) L i − 1 ( α + 1 ) ( − x ) α + i .
拉盖尔多项式的导函数
将拉盖尔多项式对自变量x 求导k 次,得到:
d k d x k L n ( α ) ( x ) = ( − 1 ) k L n − k ( α + k ) ( x ) ;
进一步有:
1 k ! d k d x k x α L n ( α ) ( x ) = ( n + α k ) x α − k L n ( α − k ) ( x ) ,
运用Template:Link-en 可以得到:
L n ( α ′ ) ( x ) = ( α ′ − α ) ( α ′ + n α ′ − α ) ∫ 0 x t α ( x − t ) α ′ − α − 1 x α ′ L n ( α ) ( t ) d t .
将拉盖尔多项式对参变量α 求导,得到下面的有意思的结果:
d d α L n ( α ) ( x ) = ∑ i = 0 n − 1 L i ( α ) ( x ) n − i .
广义拉盖尔多项式满足下面的微分方程:
x L n ( α ) ′ ′ ( x ) + ( α + 1 − x ) L n ( α ) ′ ( x ) + n L n ( α ) ( x ) = 0 ,
可以与拉盖尔多项式的k阶导数所满足的微分方程作一比较。
x L n ( k ) ′ ′ ( x ) + ( k + 1 − x ) L n ( k ) ′ ( x ) + ( n − k ) L n ( k ) ( x ) = 0 ,
仅在此式中,L n ( k ) ( x ) ≡ d L n ( x ) d x k (后面这个符号又有了新的含义)。
于是,当α = 0 时,广义拉盖尔多项式可以用拉盖尔多项式的导数表示:
L n ( k ) ( x ) = ( − 1 ) k d L n + k ( x ) d x k
式中的上标(k)容易与求导k次混淆。
正交性
伴随拉盖尔多项式在区间[0, ∞)上以权函数x α e −x 正交:
∫ 0 ∞ x α e − x L n ( α ) ( x ) L m ( α ) ( x ) d x = Γ ( n + α + 1 ) n ! δ n , m ,
这可由下式得到:
∫ 0 ∞ x α ′ − 1 e − x L n ( α ) ( x ) d x = ( α − α ′ + n n ) Γ ( α ′ ) .
伴随对称核多项式可以用拉盖尔多项式表示为:
K n ( α ) ( x , y ) : = 1 Γ ( α + 1 ) ∑ i = 0 n L i ( α ) ( x ) L i ( α ) ( y ) ( α + i i ) = 1 Γ ( α + 1 ) L n ( α ) ( x ) L n + 1 ( α ) ( y ) − L n + 1 ( α ) ( x ) L n ( α ) ( y ) x − y n + 1 ( n + α n ) = 1 Γ ( α + 1 ) ∑ i = 0 n x i i ! L n − i ( α + i ) ( x ) L n − i ( α + i + 1 ) ( y ) ( α + n n ) ( n i ) ;
也有下面的递推关系:
K n ( α ) ( x , y ) = y α + 1 K n − 1 ( α + 1 ) ( x , y ) + 1 Γ ( α + 1 ) L n ( α + 1 ) ( x ) L n ( α ) ( y ) ( α + n n ) .
进一步地,在伴L 2 [0, ∞)空间上,有:
y α e − y K n ( α ) ( ⋅ , y ) → δ ( y − ⋅ ) ,
在氢原子的量子力学处理中用到了下面的公式:
∫ 0 ∞ x α + 1 e − x [ L n ( α ) ] 2 d x = ( n + α ) ! n ! ( 2 n + α + 1 ) .
级数展开
设一个函数具有以下的级数展开形式:
f ( x ) = ∑ i = 0 f i ( α ) L i ( α ) ( x ) .
则展开式的系数由下式给出
f i ( α ) = ∫ 0 ∞ L i ( α ) ( x ) ( i + α i ) ⋅ x α e − x Γ ( α + 1 ) ⋅ f ( x ) d x .
这个级数在Lp空间 L 2 [ 0 , ∞ ) 上收敛,当且仅当
‖ f ‖ L 2 2 : = ∫ 0 ∞ x α e − x Γ ( α + 1 ) | f ( x ) | 2 d x = ∑ i = 0 ( i + α i ) | f i ( α ) | 2 < ∞ .
一个相关的展开式为:
f ( x ) = e γ 1 + γ x ⋅ ∑ i = 0 L i ( α ) ( x 1 + γ ) ( 1 + γ ) i + α + 1 ∑ n = 0 i γ i − n ( i n ) f n ( α ) ;
特别地
e − γ x ⋅ L n ( α ) ( x ( 1 + γ ) ) = ∑ i = n L i ( α ) ( x ) ( 1 + γ ) i + α + 1 γ i − n ( i n ) ,
这可由下式得到:
L n ( α ) ( x 1 + γ ) = 1 ( 1 + γ ) n ∑ i = 0 n γ n − i ( n + α n − i ) L i ( α ) ( x ) .
还有,当Re ( 2 α − β ) > − 1 时,
x α − β f ( x ) Γ ( α − β + 1 ) = ( α β ) ∑ i = 0 L i ( β ) ( x ) ( β + i i ) ∑ n = 0 i ( − 1 ) i − n ( α − β i − n ) ( α + n n ) f n ( α ) ,
这个结果可以由下式导出,
x α − β L n ( α ) ( x ) Γ ( α − β + 1 ) = ( α β ) ( α + n n ) ∑ i = n ( − 1 ) i − n ( α − β i − n ) L i ( β ) ( x ) ( β + i i )
更多的例子
幂函数 可以展开为:
x n n ! = ∑ i = 0 n ( − 1 ) i ( n + α n − i ) L i ( α ) ( x ) = ( − 1 ) n ∑ i = 0 n L i ( α − i ) ( x ) ( − α n − i ) ,
二项式 可以展开为:
( n + x n ) = ∑ i = 0 n α i i ! L n − i ( x + i ) ( α ) .
进一步可以得到:
e − γ x = ∑ i = 0 γ i ( 1 + γ ) i + α + 1 L i ( α ) ( x ) (当且仅当 Re ( γ ) > − 1 2 时收敛)
更一般地
x β e − γ x Γ ( β + 1 ) = ( α + β α ) ∑ i = 0 L i ( α ) ( x ) ( α + i i ) ∑ j = 0 i ( − 1 ) j ( 1 + γ ) α + β + j + 1 ( α + β + j j ) ( α + i i − j ) .
对于非负的整数β ,可以化简为:
x n e − γ x n ! = ∑ i = 0 γ i L i ( α ) ( x ) ( 1 + γ ) i + n + α + 1 ∑ j = 0 n ( − 1 ) n − j γ j ( n + α j ) ( i n − j ) ,
当γ = 0 时,可以化简为:
x β Γ ( β + 1 ) = ( α + β α ) ∑ i = 0 ( − 1 ) i ( β i ) L i ( α ) ( x ) ( α + i i ) , 或
x β L n ( γ ) ( x ) Γ ( β + 1 ) = ( α + β α ) ∑ i = 0 L i ( α ) ( x ) ( α + i i ) ∑ j = 0 n ( − 1 ) i − j ( n + γ n − j ) ( β + j i ) ( α + β + j j ) .
雅可比Theta 函数 有下面的表示:
∑ k ∈ ℤ e − k 2 π x = ∑ i = 0 L i ( α ) ( x t ) ∑ k ∈ ℤ ( k 2 π t ) i ( 1 + k 2 π t ) i + α + 1 ;
随意选定参量t,贝塞尔函数 可以表示为:
J α ( x ) ( x 2 ) α = e − t Γ ( α + 1 ) ∑ i = 0 L i ( α ) ( x 2 4 t ) ( i + α i ) t i i ! ;
Γ函数 可以展开为:
Γ ( α ) = x α ∑ i = 0 L i ( α ) ( x ) α + i ( ℜ ( α ) < 1 2 ) ;
低阶不完全伽玛函数 可展开为:
γ ( s ; z ) t s Γ ( s ) = ( z t ) α Γ ( α + 1 ) ∑ i = 0 L i ( α ) ( z t ) ( α + i i ) ∑ j = 0 i ( − 1 ) j ( 1 + t ) s + j ( s − 1 + j j ) ( α − 1 + i i − j ) ,
γ ( s ; z ) t s Γ ( s ) = ( α + s α + 1 ) ∑ i = 0 ( α + i + 1 i + 1 ) − L i + 1 ( α ) ( z t ) ( α + i + 1 i ) ∑ j = 0 i ( − 1 ) j ( 1 + t ) α + 1 + s + j ( α + s + j j ) ( α + i + 1 i − j ) .
还有:
γ ( s , z ) = γ s Γ ( 1 − s ) ∑ i = 0 L i + 1 ( − s ) ( 0 ) − L i + 1 ( − s ) ( z γ ) ( 1 + γ ) i + 1 ∑ n = 0 i γ i − n ( i n ) n + 1 − s ;
于是,高阶不完全伽玛函数就是:
Γ ( s , z ) z s e − z = ∑ k = 0 L k ( α ) ( z ) ( k + 1 ) ( k + 1 + α − s k + 1 ) ( ℜ ( s − α 2 ) < 1 4 ) = ∑ k = 0 L k ( α ) ( z t ) ⋅ 2 F 1 ( 1 + α + k , 1 + k ; 2 + α + k − s ; t − 1 t ) t k ( k + 1 ) ( 1 + α + k − s 1 + k ) = t s ∑ k = 0 L k ( α ) ( z t ) ⋅ 2 F 1 ( 1 − s , 1 + α − s ; 2 + α + k − s ; t − 1 t ) ( k + 1 ) ( 1 + α + k − s 1 + k ) = t 1 + α ∑ k = 0 L k ( α ) ( z t ) ⋅ 2 F 1 ( 1 + α + k , 1 + α − s ; 2 + α + k − s ; 1 − t ) ( k + 1 ) ( 1 + α + k − s 1 + k ) ,
2 F 1 表示超几何函数 。
围道积分表示
拉盖尔多项式可以用围道积分 表示,如下式所示:
L n ( α ) ( x ) = 1 2 π i ∮ e − x t 1 − t ( 1 − t ) α + 1 t n + 1 d t
积分方向逆时针绕原点一周。
与埃爾米特多項式的关系
广义拉盖尔多项式与埃爾米特多項式有下列关系:
H 2 n ( x ) = ( − 1 ) n 2 2 n n ! L n ( − 1 / 2 ) ( x 2 )
以及
H 2 n + 1 ( x ) = ( − 1 ) n 2 2 n + 1 n ! x L n ( 1 / 2 ) ( x 2 )
这里的H n 表示乘上了exp(−x 2 )的埃爾米特多項式 (所谓的“物理学家形式”)。
正因为这样,广义拉盖尔多项式也在量子谐振子 的量子力学处理中出现。
与超几何函数的关系
拉盖尔多项式可以用超几何函数 来定义,具体地说,是用合流超几何函数 定义:
L n ( α ) ( x ) = ( n + α n ) M ( − n , α + 1 , x ) = ( α + 1 ) n n ! 1 F 1 ( − n , α + 1 , x )
( a ) n 是阶乘幂 ,这里表示升阶乘 。
与贝塞尔函数的关系
拉盖尔多项式与变形贝塞尔函数 之间有以下关系:
L n ( α ) ( x ) = e x 2 ( x 4 ) n + 1 2 2 π ( n + 1 ) ! ( − 1 2 n + 1 ) ⋅ ⋅ ∑ k = 0 n ( − 1 ) k + 1 ( 2 n + 1 n − k ) ( n + α n ) ( α + 2 n + 1 n − k ) ( n − k + α n − k ) ( k + 1 2 ) K k + 1 2 ( x 2 ) = e x 2 ( 4 x ) n + α + 1 2 Γ ( α + 1 2 ) ( − α − 1 n ) ( − α − 1 2 n ) ⋅ ⋅ n ! ∑ k = 0 n ( − 2 n − 1 − 2 α k − n ) ( − 2 n − 1 − α k − n ) ( − α − 1 k − n ) ( α + 1 2 + k ) I α + 1 2 + k ( x 2 ) ,
进一步有:
L n ( α ) ( x ) = 2 4 n ( 2 n + 1 ) ( − 1 2 n ) ∑ k = 0 n ( k + 1 2 ) ( 2 n + 1 n − k ) ( n k ) 2 ( n + α k ) ( 2 n + α + 1 n − k ) x n − k ( n − k ) ! L k − 2 k − 1 ( 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.