数学物理方程第六章 勒让德多项式
勒让德多项式
y2su2inz2u2u0
1
r 2 sin2
2u
2
0
( R(r)( )( ) 代人(6.1)
1 r 2
d dr
r 2
dR dr
1
R r 2 sin
d
d
sin
(1
x2
)
d2P dx2
2x
dP dx
n(n
1)
1
m2 x2
P
0
(6.7)
如 果u(r, , )和无 关,则( ) 常数,
根
据
d 2
d 2
m2
0,
得
到 m
0
(1
x2)
d2P dx2
2x
dP dx
n(n
1)P
0
(6.8)
(6.8)称为勒让德方程
11 11
sisni2n2
dd2dd222
n(ln
1)
(6.2) (6.3)
为了以后需要,将常数l记为n(n+1), 其中n可以为实数
也可以为复数。
把 (6.2)化简
R1r 2drd22rR2dd2rR22r d2dRrr ddRrn(nn1(n)R1)0
1
sin
d
d
sin
d
d
n(n 1)sin2
1
d 2
d 2
左端只和 有关,右端只和 有关,因此两端只能是常数
1 sin
数学物理方程课件第六章勒让德多项式
2 (2n)!
2n n!
2n n! 2n n!2n 1 2n 153
2 (2n)!
2n 1!
2 2n 1
数学物理方程与特殊函数
第6章勒让德多项式
性质2 递推公式
(n 1)Pn1 (x) (2n 1)xPn (x) nPn1 (x) 0
Pn1 (x) Pn1 (x) 2n 1Pn (x)
n0
Cn
2n 1 2
1 1
x Pn (x)dx
C0
1 2
1
1 x P0 (x)dx
1 2
1
x dx
1
1 2
C2n1 0
C2n
4n 1 2
1 1
x
P2n
(x)dx
4n
1
1 0
xP2n
( x)dx
4n 1
22n 2n!
1 d2n 0 x dx2n
(x2 1)2n dx
4n 1 22n 2n !
数学物理方程与特殊函数
第6章勒让德多项式
三 勒让德多项式
y APn (x) BQn (x)
Pn
(x)
M
(1)m
m0
2n 2m!
2n m!(n m)!(n
2m)!
xn2m
Pn
1 2n n!
dn dx n
(x2
1)n
当n为偶数时M
n 2
当n为奇数时 M
n 1 2
P0 (x) 1
P1(x) x
2)(n 1)(n 4!
3)
x4
]
c 1 c0
y2
a1[ x
(n
1)(n 3!
2)
勒让德多项式递推公式证明
勒让德多项式递推公式证明以勒让德多项式是数学中一类重要的特殊函数,其递推公式是证明其性质的关键。
本文将通过介绍以勒让德多项式的定义、性质和递推公式的证明,来解释这一标题。
以勒让德多项式是数学中的一类正交多项式,它们是解决物理和工程问题中的常微分方程的重要工具。
以勒让德多项式的定义如下:$$P_n(x) = \frac{1}{2^n n!} \frac{d^n}{dx^n} \left[(x^2 - 1)^n\right]$$其中,$n$为非负整数,$P_n(x)$表示以勒让德多项式的第$n$阶,$x$为自变量。
以勒让德多项式具有一系列重要的性质,如正交性、归一性等,这些性质使其在数学和物理学中得到广泛应用。
以勒让德多项式的递推公式是证明其性质的关键。
递推公式的形式如下:$$(n+1)P_{n+1}(x) = (2n+1)xP_n(x) - nP_{n-1}(x)$$下面我们来证明这个递推公式。
我们将以勒让德多项式的定义代入递推公式中,得到:$$(n+1)\left(\frac{1}{2^{n+1} (n+1)!} \frac{d^{n+1}}{dx^{n+1}} \left[(x^2 - 1)^{n+1}\right]\right) = (2n+1)x\left(\frac{1}{2^n n!} \frac{d^n}{dx^n} \left[(x^2 - 1)^n\right]\right) - n\left(\frac{1}{2^{n-1} (n-1)!} \frac{d^{n-1}}{dx^{n-1}} \left[(x^2 - 1)^{n-1}\right]\right) $$化简上式,可以得到:$$\frac{1}{2^{n+1} (n+1)!} \frac{d^{n+1}}{dx^{n+1}} \left[(x^2 - 1)^{n+1}\right] = \frac{2n+1}{2^n n!}x\frac{d^n}{dx^n} \left[(x^2 - 1)^n\right] - \frac{n}{2^{n-1} (n-1)!}\frac{d^{n-1}}{dx^{n-1}} \left[(x^2 - 1)^{n-1}\right]$$我们将上式中的$n+1$分布到第一项中,并利用导数的链式法则进行化简,得到:$$\frac{1}{2^{n+1} (n+1)!} \frac{d}{dx}\left[(2n+1)x(x^2-1)^n\right] = \frac{(2n+1)}{2^n n!}x\frac{d^n}{dx^n} \left[(x^2 - 1)^n\right] - \frac{n}{2^{n-1} (n-1)!}\frac{d^{n-1}}{dx^{n-1}} \left[(x^2 - 1)^{n-1}\right] $$再次化简上式,得到:$$\frac{1}{2^{n+1} (n+1)!} (2n+1)\left[x\frac{d}{dx}\left[(x^2-1)^n\right] + (x^2-1)^n\right] = \frac{(2n+1)}{2^n n!}x\frac{d^n}{dx^n} \left[(x^2 - 1)^n\right] - \frac{n}{2^{n-1} (n-1)!}\frac{d^{n-1}}{dx^{n-1}} \left[(x^2 - 1)^{n-1}\right] $$继续化简上式,可以得到:$$\frac{1}{2^{n+1} (n+1)!} (2n+1)\left[x\frac{d}{dx}\left[(x^2-1)^n\right] + (x^2-1)^n\right] = \frac{(2n+1)}{2^n n!}x\frac{d^n}{dx^n} \left[(x^2 - 1)^n\right] - \frac{n}{2^{n-1} (n-1)!}\frac{d^{n-1}}{dx^{n-1}} \left[(x^2 - 1)^{n-1}\right]$$再次化简上式,得到:$$\frac{1}{2^{n+1} (n+1)!} (2n+1)\left[x\frac{d}{dx}\left[(x^2-1)^n\right] + (x^2-1)^n\right] = \frac{(2n+1)}{2^n n!}x\frac{d^n}{dx^n} \left[(x^2 - 1)^n\right] - \frac{n}{2^{n-1} (n-1)!}\frac{d^{n-1}}{dx^{n-1}} \left[(x^2 - 1)^{n-1}\right] $$继续化简上式,可以得到:$$\frac{1}{2^{n+1} (n+1)!} (2n+1)\left[x\frac{d}{dx}\left[(x^2-1)^n\right] + (x^2-1)^n\right] = \frac{(2n+1)}{2^n n!}x\frac{d^n}{dx^n} \left[(x^2 - 1)^n\right] - \frac{n}{2^{n-1} (n-1)!}\frac{d^{n-1}}{dx^{n-1}} \left[(x^2 - 1)^{n-1}\right] $$再次化简上式,得到:$$\frac{1}{2^{n+1} (n+1)!} (2n+1)\left[x\frac{d}{dx}\left[(x^2-1)^n\right] + (x^2-1)^n\right] = \frac{(2n+1)}{2^n n!}x\frac{d^n}{dx^n} \left[(x^2 - 1)^n\right] - \frac{n}{2^{n-1} (n-1)!}\frac{d^{n-1}}{dx^{n-1}} \left[(x^2 - 1)^{n-1}\right] $$通过以上推导,我们证明了以勒让德多项式的递推公式。
勒让德多项式
例1:将 x 2 在[-1,1]内展成勒让德多项式的级数形式
x 2 Cn Pn (x) n0
Cn
2n 1 2
1 1
x
2
Pn
(
x)dx
1 1
xk
Pn
( x)dx
0
n2
4 1
C2 2
1 x2 1 (3x2 -1)dx 5
1 2
4
1 3x4 x2
1
dx
5 6 2 2 45 3 3
第6章勒让德多项式
例2:将Pl(x) 在[-1,1]内展成勒让德多项式的级数形式
解:方法一
l 1
(l 1) / 2
Pl(x) CnPn (x) CnPn (x)
Cl2n1Pl2n1 ( x)
n0
n0
n0
2l 4n 1
Cl2n1
2
1
1 Pl(x)Pl2n1(x)dx
2l 4n 1 2
1 0
xd
d 2n1 dx 2 n 1
(x2
1)2n
4n 22n
1 2n
!
x
d 2 n 1 dx 2 n 1
(x2
1)2n|10源自1 0d 2 n1 dx 2 n 1
(x2
1)2n
dx
4n 22n
1 2n
!
d 1 2n1 0 dx2n1
(x2
1)2n dx
4n 22n
1 2n
!
d2n2 dx 2 n 2
0
0
0
/ 2 sin 2n1 d 2n / 2 sin 2n1 d
0
2n 1 0
1 P2n (x)dx 1
第六章 勒让德多项式
y1 ( x ) = ∑ m = 0 a2 m x 2 m ,
∞
y2 ( x ) = a1 x + a3 x 3 + a5 x 5
西安理工大学应用数学系
不妨取n为非负整数,那么对应多项式结构如何? 不妨取 为非负整数,那么对应多项式结构如何?这时 为非负整数
an+2 = an+4 =⋯= 0 ak ≠ 0, k ≤ n
( n − 1)( n + 2) a3 = − a1 3⋅ 2 ( n − 3)( n + 4) ( n − 1)( n − 3)( n + 2)( n + 4) a5 = − a3 = a1 5⋅4 5!
西安理工大学应用数学系
( n − 3)( n + 4) ( n − 1)( n − 3)( n + 2)( n + 4) a5 = − a3 = a1 5⋅4 5!
y2 ( x ) 中有
西安理工大学应用数学系
( k − n)( k + n + 1) ak + 2 = ak k = 0,1, 2,⋯ ( k + 1)( k + 2) m n( n − 2)⋯ ( n − 2 m + 2)( n + 1)( n + 3)⋯ ( n + 2 m − 1) a2 m = ( −1) a0 (2m )! m ( n − 1)( n − 3)⋯( n − 2m + 1)( n + 2)( n + 4)⋯( n + 2m ) a2 m +1 = ( −1) a1 (2m + 1)!
( k − n)( k + n + 1) ak + 2 = ak ( k + 1)( k + 2) ( k + 1)( k + 2) ak = ak + 2 k ≤ n−2 ( k − n)( k + n + 1) n( n − 1) an − 2 = − an 2(2n − 1) ( n − 2)( n − 3) n( n − 1)( n − 2)( n − 3) an − 4 = − an − 2 = an 4(2n − 3) 2 ⋅ 4(2n − 1)(2n − 3) n( n − 1)( n − 2)⋯( n − 2m + 1) m an − 2 m = ( −1) an 2 ⋅ 4⋯ ⋅ 2m (2n − 1)⋯ (2n − 2m + 1)
勒让德多项式
(1) k (2l 2k )! l 2 k Pl ( x) l x k 0 2 k!(l k )!(l 2k )!
n
一. 特殊值、奇偶性和图形
l 2 l 1 n 当l为奇数时 2
当l为偶数时 n
Pl (1) 1,
P2 n (0) c0 (1) n
六. 勒让德多项式的正交性、完备性与模
0, lk 2 1 Pl ( x)Pk ( x)dx Nl2 , l k 2l 1
1
勒让德多项式完备性 若f(x)是定义在[-1,1]区间上任意一个平方可积的函数,
那么
f ( x) cl Pl ( x)
l 0
(l 1) P l 1 ( x) lP l 1 ( x) (2l 1) xP l ( x)
2. P l ( x) P l 1 ( x) 2 xP l ( x) P l 1 ( x)
3. 4.
P l 1 ( x) xP l ( x) (l 1) P( x) Pl 1 ( x) P l 1 ( x) 2l 1P l ( x)
1 1 2rx r xr
2
r Pl ( x)
l l 0 2
(l 1) P l 1 ( x) lP l 1 ( x) (2l 1) xP l ( x)
1 2rx r 2
(1 2rx r ) lr l 1Pl ( x)
l 0
( x r ) r l Pl ( x) (1 2rx r 2 ) lr l 1Pl ( x)
证
2 l
1 dl 2 l Pl ( x) l ( x 1 ) 2 l! dx l
勒让德多项式
从而得到
1
ห้องสมุดไป่ตู้
Θ
sin θ
d dΘ (sin θ ) + n( n + 1) sin 2 θ = m 2 dθ dθ
( 6. 4 )
( 6. 5 )
1 d 2Φ + m2 = 0 2 Φ dϕ
(2 勒让德多项式的一些性 质; )
有关的定解问题。 (3 会用勒让德多项式求解 有关的定解问题。 )
§6.1
勒让德方程的引出
u xx + u yy + uzz = 0
在第四章中, 域内的迪利克雷问题: 在第四章中,我们用格 林函数法解决了球形区 域内的迪利克雷问题:
{
球函数
z
θ
●
拉普拉斯方程 第一类边界条件
数学物理方法
第六章 勒让德多项式 ( Legendre polynomials )
勒让德( 勒让德(1752~1833) ~ ) Legendre . Adrien-Marie 阿德利昂·玛利 埃 勒让德 公元1752─公元1833 为法国数学家, 勒让德( 1752─公元1833) 阿德利昂 玛利·埃·勒让德(公元1752─公元1833)为法国数学家,生于 玛利 巴黎,卒于巴黎。 1770年毕业于马扎兰学院 1775年任巴黎军事学院数学 年毕业于马扎兰学院。 巴黎,卒于巴黎。约1770年毕业于马扎兰学院。1775年任巴黎军事学院数学 教授。1782年以 关於阻尼介质中的弹道研究》获柏林科学院奖金, 年以《 教授。1782年以《关於阻尼介质中的弹道研究》获柏林科学院奖金,次年当 选为巴黎科学院院士。1787年成为伦敦皇家学会会员 年成为伦敦皇家学会会员。 选为巴黎科学院院士。1787年成为伦敦皇家学会会员。 曾与拉格朗日( )、拉普拉斯 拉普拉斯( 勒让德 (Legendre) 曾与拉格朗日(Lagrange)、拉普拉斯(Laplace) 并列为法国数学界的“ 世纪末19世纪初法国数学的复兴, 并列为法国数学界的“三 L ”,为18世纪末19世纪初法国数学的复兴,做出了 , 18世纪末19世纪初法国数学的复兴 卓越的贡献。 卓越的贡献。
第六章---数理方程勒让德多项式
y2
x
(n
1)(n 3!
2)
x3
(n
1)(n
3)(n 5!
2)(n
4)
x5
(2k 1 n)(2k 3 n) (1 n)(n 2) (n 2k) x2k1 (2k 1)!
6. 3 勒让德多项式
6. 3 勒让德多项式
将6.2中的递推公式写成
ak
(k 2) (k 1) (n k)(k n 1)
2)!
6. 3 勒让德多项式
an4
(2n 4)! 2!2n (n 2)!(n
4)!
一般地当 n 2k 0 时,有
6.1 勒让德方程的引出
第二个方程为
d 2
d 2
cot
d
d
n
n
1
m2
sin2
0
令 x cos ,并记 P( x) (cos )
1 x2
d2P dx 2
2x
dP dx
n
n
1
m2 1 x2
P
0
k0
(k c 2)(k c 1)ak2 [(k c)(k c 1) n(n 1)ak 0
a k0 k2
(k
(k+c)(ck)(k
cc1)a1k)xk
c2
n(n
(k+c 1)(k c 2)
0 1)
ak
大学物理-勒让德多项式
(三) 勒让德多项式的正交性与正交归一关系式 (I) 勒让德多项式的正交性 勒让德方程与有限性条件构成施—刘型本征值问题:
勒让德方程可化为施—刘型方程,其中 (x) = 1,a =
–1,b = 1,而本征函数 yl (x) = Pl (x) 是实函数,故有如下正 交关系
另一种形式:
(II) 勒让德多项式的完备性 根据施—刘型本征值问题的完备性定理,对于在区间
示——微分表示
证明:由二项式展开定理得到
所以
注意到:凡是指数 (2l – 2s) < l 的项经 l 次求导后为 0,故 只剩下 (2l – 2s) ≥ l 的项,即 s ≤ l / 2,于是得到
罗德里格斯公式的应用之一 因为
所以
对于 Pl (1),利用上式,n > 0 各项的值均为 0,只有当 n = 0 的项,其值不为 0,因此
之间存在着一定的关系。具体如下:
(1) (2) (3) (4)
证明:(I) 由生成函数关系式 两边对 r 求导,有
改写为 两边乘以 (1–2 x r + r2),再将生成函数关系式代入,得到
比较两边 rl 的系数,有 整理上式,得到 当 l = 0 时,P0(x) = 1,P1(x) = x, 所以 x P0(x) = P1(x) (II) 将生成函数关系式两边对 x 求导,得到
例3 将函数 f (x) = xn 在区间 [–1,1] 上按完备正交函数系
Pl (x) 展开为广义傅里叶级数。
解:所给函数满足斯–刘型本征函数完备性定理的条件,故
xn Cl Pl (x) l 0
广义傅里叶系数为
Cl
2l 1 2
1 1
x
n
Pl
第六章 勒让德多项式1
( 2 1)l ( x)
l 1
C
dx
(6.1.11)
z 平面上围绕 z x 点的任一闭合回路,
并取正方向.这叫作勒让德多项式的施列夫利积分表示式.
式(6.1.11)还可以进一步表为下述拉普拉斯积分.
Pl ( x) 1 π 0
π
( x i 1 x 2 cos )l d
l
阶勒让德多项式.勒让德多项式也称为第一类勒让德函数.
式(6.1.7)即为勒让德多项式的级数表示. 注意到 x cos , 故可方便地得出前几个勒让德多项式:
P0 ( x) 1
P1 ( x) x cos
P2 ( x) 1 2 (3x 2 1) 1 4 (3cos 2 1)
(6.1.12)
x2 1
【证明】 取 C 为圆周,圆心在 z x ,半径为 .在
C
上有:
x x 1e
2
i
d i x 2 1ei d i( x)d
并注意到
2 1 ( x x 2 1ei )2 1 ( x 2 1)(1 ei2 ) 2 x x 2 1ei
P3 ( x)
1
1 (5 x3 3x) (5cos 3 3cos ) 2 8
1 1 P4 ( x) (35 x 4 30 x 2 3) (35cos 4 20 cos 2 9) 8 64 1 1 P5 ( x) (63x5 70 x3 15 x) (63cos 5 35cos 3 30 cos ) 8 128
dl dx
l
x
2l 2 k
(2l 2k )(2l 2k 1) [2l 2k (l 1)]x
《数学物理方法》第六章勒让德函数
①为了讨论系数的解析性质,以判定z0=0是方程的 常点、正则奇点还是非正则奇点,必须将p(x)及q(x) 分别延拓为
但为叙述与书写方便,仍采用x⇔z的记号
12
2. 系数递推公式 由此得系数递推公式
13
3. 由递推公式求系数,得通解
14
勒让德多项式 微分表达式-罗德里格斯(Rodrigues)公式; 母函数; 积分表达式—施列夫利公式和拉普拉斯
积分 递推公式.
6.2.1 勒让德多项式的微分表达式—罗德里格 斯公式 证明 从罗德里格斯公式右边出发来证明.
二项式展开定理为
32
对(x2-1)l求l阶导数后除以(2ll!)得到
为何求和指标的最大值为[l/2],因为对于指 数(2l-2s)<l的项,在求l 阶导数后均为零,故: 只含(2l-2s)≥l的项,即:s≤ l/2的项.这样当 l 为偶数时,l/2为最大值; l为奇数时,(l-1)/2 为最大值。用简写符号表示就是 [l/2]
证明 (1)在|t|<1内,将w(x,t)展开为泰勒级数
其中al为泰勒系数, C为在|t|<1内包围t=0点的回路
①奇点
的|t12|<1
36
(2)为证明al =Pl(x),作变换(u为复变数)
37
代入al ,便有
其中u平面的曲线Cʹ是在式(6.2.5)的变换下t平面曲线 C的像.当t=0时,由式 (6.2.6)得到u=x.既然t=0在 曲线C的内部,因此u=x在曲线Cʹ的内部.
式(6.1.17)乘以任意常数仍为勒让德方程的解 历史上为了让这个多项式与函数(1-2xt+t2)-1/2
的展开系数一致,选择最高次幂项的系数Cl 为
勒让德多项式pn(0)项取值
勒让德多项式是数学中的一个重要概念,它在物理学和工程学中都有广泛的应用。
勒让德多项式的第n个零阶项pn(0)取值对于理解和应用勒让德多项式至关重要。
本文将从勒让德多项式的定义、性质和应用角度对pn(0)项取值进行深入探讨。
一、勒让德多项式的定义1. 勒让德多项式是一个经典的正交多项式族,它们是勒让德方程的解。
勒让德多项式Pn(x)定义为以下多项式序列的一部分:P0(x) = 1P1(x) = xP2(x) = (1/2)(3x^2-1)P3(x) = (1/2)(5x^3-3x)…Pn(x) = (1/2^n)(d^n/dx^n)((x^2-1)^n)2. 勒让德多项式具有许多重要的性质,如正交性、递推关系和母函数等,这些性质使得勒让德多项式在实际问题中具有广泛的应用价值。
二、勒让德多项式的性质及pn(0)项取值1. 正交性:勒让德多项式满足以下正交性质:∫[-1,1] Pm(x)Pn(x)dx = (2/(2n+1))δmn其中δmn是克罗内克δ 符号,当m=n时取值为1,否则为0。
利用这一性质,可以得到pn(0)项取值的一般表达式。
2. 递推关系:勒让德多项式之间存在递推关系,即Pn+1(x)和Pn(x)之间存在一定的关系,这使得我们可以通过递推的方法求得pn(0)的取值。
3. pn(0)的表达式:根据正交性质和递推关系,可以得到勒让德多项式的pn(0)的具体表达式:pn(0) = (-1)^n(2n)!!其中(2n)!!是双阶乘的符号,表示连乘2n, 2n-2, 2n-4, …, 2或1。
三、勒让德多项式的应用1. 在物理学中,勒让德多项式被广泛应用于描述球对称的问题,如原子轨道、分子振动等。
pn(0)项的取值对于求解这些问题具有重要意义。
2. 在工程学中,勒让德多项式被用于信号处理、控制系统等领域。
了解pn(0)项的取值可以帮助工程师更好地理解和应用勒让德多项式。
勒让德多项式pn(0)项的取值对于理解和应用勒让德多项式具有重要意义。
第六章 勒让德多项式
第六章 勒让德多项式在这一章,我们将通过在球坐标系中对Laplace 方程进行分离变量,引出§2.6中曾指出过的勒让德方程,并讨论这个方程的解法及解的有关性质。
勒让德方程在区间[1,1]-上的有界解构成了另一类正交函数系-勒让德多项式。
§6.1 勒让德方程的引出现在对球坐标系中的Laplace 方程进行分离变量,在球坐标系中Laplace 方程为2222222111()(sin )0sin sin u u ur r r r r r θθθθθϕ∂∂∂∂∂++=∂∂∂∂∂ (6.1) 令 (,,)()()()u r R r θϕθϕ=ΘΦ, 代入(6.1)得2222222111()(sin )0sin sin d dR d d d r R R r dr dr r d d r d θθθθθϕΘΦΘΦ+Φ+Θ= 以2r R ΦΘ乘上式各项得 2222111()(sin )0sin sin d dR d d d r R dr dr d d d θθθθθϕΘΦ++=ΘΦ 或2222111()(sin )sin sin d dR d d d r R dr dr d d d θθθθθϕΘΦ=--ΘΦ 上式左端只与r 有关,右端只与θ,ϕ有关,要它们相等只有当它们都是常数时才有可能。
为了以后的需要,我们把这个常数写成(1)n n +的形式(这是可以做到的,因为任何一个实数总可以写成这种形式,这里的n 可能为实数,也有可能为复数),则得21()(1)d dRr n n R dr dr=+ (6.2)22211(sin )(1)sin sin d d d n n d d d θθθθθϕΘΦ+=-+ΘΦ (6.3) 将方程(6.2)左端的导数计算出来,即有2222(1)0d R dRr r n n R dr dr+-+= 这是一个欧拉方程,它的通解为(1)12()n n R r A r A r -+=+其中12,A A 为任意常数。
勒让德多项式
课程设计报告n 阶勒让德(Legendre)多项式一、设计任务与目标n 阶勒让德(Legendre)多项式可以递归定义如下:⎪⎩⎪⎨⎧>---===--1/))(*)1()(**)12((101)(21n nx p n x p x n n xn x p n n n(1) 输入n 和x 的数值,输出此时勒让德多项式的数值。
例如输入2,1,应输出1/2。
(2)输入n 的数值, 输出此时的勒让德多项式。
例如输入2,应输出3/2 x 2 - 1/2。
本次上机实践所使用的平台和相关软件。
平台:Windows xp 相关软件:VC6.0二、方案设计与论证对于这个题目,我分析了一下,第一问是要求我要用递归方法去求最终的值,所以我在程序中编写了子函数treat ,并在主函数main 中调用,在子函数中不断调用自己本身。
第二问,由于不能按照常规来做,只能够想一些特别的方法,例如:利用字符串输出,但这种方法不行。
经老师提醒,先做好这个表达式的每一项的情况,然后再将他们整合输出,于是我选择了这个方法并向着这个方向去做,后来在网上找了相关的资料,我发现了这么一条公式:∑⎥⎦⎤⎢⎣⎡=-----=202)!2()!(!2)!22()1()(n m mn nmn x m n m n m m n x p ,这一条公式可以求出表达式的每一项,我利用四个数组,第一个数组是记录m)1(-的结果;第二个数组是记录)!2()!22(m n m n --约简后的结果;第三个数组是记录)!(!2m n m n -约简后的结果;第四个数组是记录m n 2-的结果。
最后输出每一项并整合最终的结果。
在计算)!2()!22(m n m n --之前,我采用了没有约简的方法去做,结果数值超出了我设定的int 型数据的范围,导致我只能够输出n=6的情况,n=7输出错误。
后来利用约简的方法,于是结果达到了n=8的情况。
接着我采用了double 型,结果能输出n=10的情况,但是在运行的过程中发现,输出很慢。
勒让德多项式
六、结论与心得
在调试的过程中,一开始我发现原来的程序进行到n=6时 就发生了溢出,经过老师的指点后,我调整了程序后,将一部 分数据进行约分,减少大数据的出现,有了进一步的发展,能 够输到n=8,然而我没有放弃,利用double型和int型的区别去 想,继续调整,终于达到n=10的阶段。但是后来的调整时发现 数据很大,难以快速得到结果,因此我还没有时间去想如何调 整到一个高效的程序。在这个过程中,我知道了有些事情,只 有经历过了,才会知道自己的不足,我决定以后学习好专业语 言,努力改进自己编写过的程序,将程序写得更好!
况,n=7输出错误。后来利用约简的方法,于是结果达到了n=8
的情况。接着我采用了double型,结果能输出n=10的情况,但
是在运行的过程中发现,输出很慢。
三、程序框图或流程图,程序清单与调用关系
第一问如图一所示:
图1 — 递归调用流程图 Y Y N N
开始 输入n x n=0 ?
n=1? 返回 1 返回 1 递归调用treat() 输出结果 结束
第二问如图二所示:
输出整合的结果
Y N 开始 输入 n m=0
m﹤n/2 ?
结束
m=m+1
图2 — 求表达式的流程图
四、全部源程序清单
#include<stdio.h>
#define MAX 1000
#include<math.h>
float treat(int n,float x)
{//计算n阶勒让德多项式的值
{ if(D[i]!=0&&B[i]!=1)//输出首项的表达式 printf("%.0f/%.0fx^%.0f",A[i]*B[i],C[i],D[i]); else { if(D[i]==0&&B[i]==1)//输出n=0的表达式 printf("%.0f",A[i]); else if(D[i]!=0&&B[i]==1) printf("x");//输出n=1的表达式 else printf("%.0f/%.0f",A[i]*B[i],C[i]);//输出
勒让德多项式
都是原方程的解。 都是原方程的解。当 a0 ≠ 0, a1 ≠ 0 ,则
y1 , y2 线性无关。 线性无关。
得到原方程
d2y dy (1 − x ) 2 − 2 x + n( n + 1) y = 0 dx dx
2
的两个线性无关的解, 的两个线性无关的解,则原方程的通解可表示成如下形式
n( n + 1) 2 n( n − 2)( n + 1)( n + 3) 4 y = a0 1 − x + x + ... 2! 4! ( n − 1)( n + 2) 3 ( n − 1)( n − 3)( n + 2)( n + 4) 5 x + x + ... + a1 x − 3! 5!
3o + k + 2)(c + k + 1)a k + 2 (c − [(c + k )(c + k + 1) − n( n + 1)]a k = 0 k = 0,1,2... ( )
1o c(c − 1)a0 = 0 a0 ≠ 0
c = 0,1
2o c(c + 1)a1 = 0 设a1 ≠ 0
0
2o c(c + 1)a1 = 0
3o + k + 2)(c + k + 1)a k + 2 (c − [(c + k )(c + k + 1) − n( n + 1)]a k = 0 k = 0,1,2,3,... ( )
(6.13) (6.14)
1o c(c − 1)a0 = 0
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
(
)
n
n n! 1 1 n 2 − = x x2 ) ( ) ( 1 ∑ n n 2 n! 2 n! m =0 (n − m )!m!
a n −6 = −
2
n
一般说来,当 n − 2m ≥ 0 时,有
M
a n − 2 m = (− 1)
m
2
n
(2n − 2m )! m!(n − m )!(n − 2m )!
(2n − 2m )! x n−2m m!(n − m )!(n − 2m )! (2n − 2m )! x n−2m m!(n − m )!(n − 2m )!
2 2 2
(6.2.1)
的解为
y = ∑ ak x k
k =0
∞
(6.2.2)
,整理得 对上式求导,得出 y ′, y ′′ 的级数表达式,连同式(6.2.2)一齐代入式(6.2.1)
∑ {(k + 1)(k + 2)a
k =0
∞
k +2
+ [n(n + 1) − k (k + 1)]a k }x k = 0
(3x (5x
2
−1
) ) ) )
3
− 3x
4
(35x (63x
− 30 x 2 + 3
5
− 70 x 3 + 15 x
它们的图形如图 6-1 所示。
为了应用上的方便,我们将 Pn ( x ) 表示为
Pn ( x ) =
n 1 dn 2 ( x − 1) n n 2 n! dx
(6.3.2)
的形式。称式(6.3.2)为勒让德多项式的罗德里格斯(Rodrigues)表达式。该公式的证明如下。 证明:用二项式定理把 x − 1 展开,有
式中, C1 , C 2 为任意常数。 对式(6.1.5)进行整理,有
d 2Φ dΦ ⎡ m2 ⎤ ( ) n n + cot + + 1 − ϕ ⎢ ⎥Φ = 0 dϕ ⎣ dϕ 2 sin 2 ϕ ⎦
这个方程称为连带的勒让德方程。
(6.1.6)
为了表达上的方便,我们引入新的变量 x = cos ϕ 。由于 0 ≤ ϕ ≤ π ,所以 − 1 ≤ x ≤ 1 ,并记
n(n − 1)(n − 2)(n − 3)(n − 4)(n − 5) an 2 ⋅ 4 ⋅ 6(2n − 1)(2n − 3)(2n − 5)
M
为了形式上好看,并使所得的多项式在 x = 1 处值为 1,我们将任意常数 a n 的值取为
an =
(2n )! 2 2 n (n!)
=
1 ⋅ 3 ⋅ 5L (2n − 1) n!
将变量 R (r ), Y (θ , ϕ ) 分离,得
]R = 0
∂ ⎛ ∂Y ⎞ ∂ 2Y 1 d ⎛ 2 dR ⎞ 1 1 1 ⎜ ⎟ + 2 sin ϕ ⎜r ⎟=− [ 2 2 2 R dr ⎝ dr ⎠ Y r sin ϕ ∂ϕ ⎜ ∂ϕ ⎟ ⎝ ⎠ r sin ϕ ∂θ
]
上式左端只于 r 有关,右端只与 θ , ϕ 有关,所以二者都是常数时才能恒等。为了方便后续的讨论, 我们把这个常数写成 n(n + 1) 的形式(这里的 n 可以是实数,也可以是复数) ,于是有
(n = 1,2,L)
从而有
an−2 = −
(2n − 2)! n(n − 1) (2n )! =− n 2 n 2(2n − 1) 2 (n!) 2 (n − 1)!(n − 2 )!
2
n
an−4 = −
(2n − 4)! 2!(n − 2 )!(n − 4 )! (2n − 6)! 3!(n − 3)!(n − 6 )!
(6.2.5) (6.2.6)
显然,级数 y1 , y 2 都是勒让德方程式(6.2.1)的解,且 y1 与 y 2 是线性无关的。由达朗贝尔(比值)判 别法和式(6.2.6)的收敛区间为 x < 1 ,因此在 x < 1 内,勒让德方程式(6.2.1)的通解为
y = C1 y1 + C 2 y 2
dy y (1 − x ) d − 2 x + n(n + 1) y = 0 dx dx
2 2 2
(6.1.8)
式(6.1.8)称为勒让德方程。一部分定解问题的求解,最后都归纳为勒让德方程的求解。 6.2 勒让德方程的求解
和求贝赛尔方程一样,我们设勒让德方程
y dy (1 − x ) d − 2 x + n(n + 1) y = 0 dx dx
y = Φ (ϕ ) ,于是有
dΦ dy dx dy dy = = − sin ϕ = − 1− x2 dϕ dx dϕ dx dx
dΦ dy dx d ( − 1 − x 2 dy = = dx dϕ dx dϕ dx
将 x, y, 及 y 的导数代入式(6.1.6)中,整理得
) dx
dϕ
= −x
式中, a 0 , a1 为任意常数。若用 y1 , y 2 分别表示式(6.2.4)中的两个级数,即
(6.2.4)
n(n + 1) 2 n(n − 2 )(n + 1)(n + 3) 4 x + x −L 2! 4! (n − 1)(n + 2) x 3 + (n − 1)(n − 3)(n + 2)(n + 4) x 5 − L y2 = x − 3! 5! y1 = 1 −
∂ 2u ∂ 2u ∂ 2u + + =0 ∂x 2 ∂y 2 ∂z 2
的求解问题。本节,在球坐标系下,我们应用分离变量的方法来处理拉普拉斯方程。在球坐标系中, 拉普拉斯方程为
1 ∂ ⎛ 2 ∂u ⎞ 1 ⎛ ∂u ⎞ 1 ∂ 2u ⎜ sin ϕ ⎟+ 2 =0 ⎜r ⎟+ 2 2 ∂ϕ ⎟ r 2 ∂r ⎝ ∂r ⎠ r 2 sin ϕ ⎜ ⎝ ⎠ r sin ϕ ∂θ
式中, y1 , y 2 分别表示级数式(6.2.5)和式(6.2.6) ; C1与C 2 为任意常数。 6.3 勒让德多项式
上面我们解出了勒让德方程式(6.2.1)的解,由级数式(6.2.5)和式(6.2.6)可知,当 n 不是整 数时, y1与y 2 都是无穷级数,它们在 x < 1 内绝对收敛,在 x = + 1 时趋向于无穷。
d2y dy + 1− x2 dx dx 2
(
)
m dy ⎡ y (1 − x ) d − 2 x + ⎢n(n + 1) − dx 1− x dx
2 2 2 2
⎣
2
⎤ ⎥y = 0 ⎦
(6.1.7)
若 u (r , θ , ϕ ) 与 θ 无关,则由式(6.1.4)可知, Θ(θ ) 是常数,则 m = 0 。这时,式(6.1.7)简化为
(k ≤ n − 2)
这样,系数 a n − 2 , a n − 4 , L 都可以通过最高次项的系数 a n 来表示,即
an−2 = −
n(n − 1) an 2(2n − 1) n(n − 1)(n − 2 )(n − 3) an 2 ⋅ 4(2n − 1)(2n − 3)
an−4 = −
a n −6 = −
r2
d 2R dR + 2r − n(n + 1)R = 0 2 dr dr
(6.1.2)
∂Y ⎞ 1 ∂ 2Y 1 ∂ ⎛ ⎜ ⎟ sin ϕ + + n(n + 1)Y = 0 2 2 sin ϕ ∂ϕ ⎜ ∂ϕ ⎟ ⎝ ⎠ sin ϕ ∂θ
式(6.1.3)的解 Y (θ , ϕ ) 与半径 r 无关,故称之为球面函数,或简称为球函数。 式(6.1.2)是欧拉方程,其通解为
式中,0 ≤ ϕ ≤ π , 0 ≤ θ ≤ 2π 。 令式(6.1.1)的解为 u (r , θ , ϕ ) = R (r )Y (θ , ϕ ), 代入式(6.1.1) ,整理得
(6.1.1)
1 d ⎛ 2 dR ⎞ 1 ∂ ⎛ ∂Y ⎞ 1 ∂ 2Y ⎜ ⎟ ϕ r Y + [ sin + ⎜ ⎟ 2 2 2 ∂ϕ ⎟ r 2 dr ⎝ dr ⎠ r 2 sin ϕ ∂ϕ ⎜ ⎝ ⎠ r sin ϕ ∂θ
3!
1
a5 = −
a7 = −
(n − 1)(n − 3)(n + 2)(n + 4) a
5!
7!
1
(n − 1)(n − 3)(n − 5)(n + 2)(n + 4)(n + 6) a
1
将这些值代入式(6.2.2) ,则得到了勒让德方程式(6.2.1)的解
⎡ n(n + 1) 2 n(n − 2)(n + 1)(n + 3) 4 ⎤ y = a 0 ⎢1 − x + x − L⎥ 2! 4! ⎣ ⎦ (n − 1)(n + 2) x 3 + (n − 1)(n − 3)(n + 2)(n + 4) x 5 − L⎤ ⎡ + a1 ⎢ x − ⎥ 3! 5! ⎣ ⎦
(6.1.3)
R(r ) = Ar n + Br − (n +1)
式中,A,B 为任意常数。 式 (6.1.3) 中喊有两个自变量 θ , ϕ , 再次应用分离变量的方法, 令 Y (θ , ϕ ) = Θ(θ )Φ (ϕ ) , 代入 (6.1.3) 式中,整理得
dΦ ⎞ 1 d ⎛ 1 d 2Θ ⎜ sin ϕ ⎟Θ + Φ + n(n + 1)ΘΦ = 0 dϕ ⎟ sin ϕ dϕ ⎜ sin 2 ϕ dθ 2 ⎝ ⎠