数学物理方法_第5章 Legendre多项式(上)
legendre多项式推导

legendre多项式推导
摘要:
1.Legendre 多项式的定义和基本概念
2.Legendre 多项式的性质和应用
3.Legendre 多项式的推导过程
正文:
Legendre 多项式是数学中的一种多项式,它的定义和基本概念如下:Legendre 多项式是指形如P_n(x) = (x - x_1)(x - x_2)...(x - x_n) 的多项式,其中x_1, x_2,..., x_n 是多项式的根。
它也可以写成如下形式:P_n(x) = A_n x^n + A_{n-1} x^{n-1} +...+ A_1 x + A_0
其中A_n, A_{n-1},..., A_1, A_0 是多项式的系数。
Legendre 多项式具有很多重要的性质和应用,比如:
1.Legendre 多项式的根是均匀分布在实数轴上的。
2.Legendre 多项式的系数A_n, A_{n-1},..., A_1, A_0 都可以通过Vieta 定理求出。
3.Legendre 多项式在数值分析和计算机图形学中有广泛应用,比如在求解数值积分和计算曲线积分时,常常使用Legendre 多项式作为基函数。
数学物理方程-第五章格林函数法

第五章 格林函数法在第二章中利用分离变量法求出了矩形区域和圆域上位势方程Dirichlet 问题的解.本章利用Green 函数法求解一些平面或空间区域上位势方程Dirichlet 问题. 另外,也简单介绍利用Green 函数法求解一维热传导方程和波动方程半无界问题. 应指出的是:Green 函数法不仅可用于求解一些偏微分方程边值问题或初边值问题,特别重要的是,它在偏微分方程理论研究中起着非常重要的作用.§5⋅1 格林公式在研究Laplace 方程或Poisson 方程边值问题时,要经常利用格林(Green )公式,它是高等数学中高斯(Gauss )公式的直接推广.设Ω为3R 中的区域,∂Ω充分光滑. 设k 为非负整数,以下用()k C Ω表示在Ω上具有k 阶连续偏导的实函数全体,()k C Ω表示在Ω上具有k 阶连续偏导的实函数全体. 如()10()()()()u C C C C ∈Ω⋂ΩΩ=Ω,表示(,,)u x y z 在Ω具有一阶连续偏导数而在Ω上连续. 另外,为书写简单起见,下面有时将函数的变量略去.如将(,,)P x y z 简记为P ,(,,)P x y z x ∂∂简记为Px∂∂或x P 等等.设(,,)P x y z ,(,,)Q x y z 和(,,)R x y z 1()C ∈Ω,则成立如下的Gauss 公式()P Q R dV Pdydz Qdydx Rdxdy x y z Ω∂Ω∂∂∂++=++∂∂∂⎰⎰⎰⎰⎰ (1.1) 或者()(cos cos cos )P Q R dV P Q R ds x y z αβγΩ∂Ω∂∂∂++=++∂∂∂⎰⎰⎰⎰⎰ (1.2) 如果引入哈米尔顿(Hamilton )算子: (,,)x y z∂∂∂∇=∂∂∂,并记(,,)F P Q R =,则Gauss 公式具有如下简洁形式⎰⎰⎰⎰⎰∂⋅=⋅∇ΩΩds n F dv F(1.3)其中(cos ,cos ,cos )n αβγ=为∂Ω的单位外法向量.注1 Hamilton 算子是一个向量性算子,它作用于向量函数(,,)F P Q R =时,其运算定义为(,,)(,,) ,F P Q R x y zP Q Rx y z∂∂∂∇⋅=⋅∂∂∂∂∂∂=++∂∂∂形式上相当于两个向量作点乘运算,此即向量F 的散度div F . 而作用于数量函数(,,)f x y z 时,其运算定义为(,,)(,,)f f ff f x y z x y z∂∂∂∂∂∂∇==∂∂∂∂∂∂, 形式上相当于向量的数乘运算,此即数量函数f 的梯度grad f .设(,,)u x y z ,2(,,)()v x y z C ∈Ω,在(1.3)中取F u v =∇得 ()u v dV u v nds Ω∂Ω∇⋅∇=∇⋅⎰⎰⎰⎰⎰ (1.4)直接计算可得v u v u v u ∇∇+=∇⋅∇∆)( (1.5)其中xx yy zz v v v v ∆=++. 将(1.5)代入到(1.4)中并整理得vu vdV uds u vdV n Ω∂ΩΩ∂∆=-∇⋅∇∂⎰⎰⎰⎰⎰⎰⎰⎰ (1.6) (1.6)称为Green 第一公式.在(1.6)中将函数u ,v 的位置互换得uv udV vds v udV n Ω∂ΩΩ∂∆=-∇⋅∇∂⎰⎰⎰⎰⎰⎰⎰⎰ (1.7) 自(1.6)减去(1.7)得()()v uu v v u dV uv ds n nΩ∂Ω∂∂∆-∆=-∂∂⎰⎰⎰⎰⎰ (1.8) (1.8)称为Green 第二公式.设点0(,,)P ξηζ∈Ω,点3(,,)P x y z R ∈,||00P P r P P -==引入函数 001(,)4P PP P r πΓ=,注意0(,)P P Γ是关于六个变元(,,)x y z 和(,,)ξης的函数且00(,)(,)P P P P Γ=Γ. 如无特别说明, 对b 求导均指关于变量(,,)x y z 的偏导数. 直接计算可得00(,)0, P P P P ∆Γ=≠即0(,)P P Γ在3R 中除点0P 外处处满足Laplace 方程.设0ε>充分小使得00(,){(,,) ||}B B P P x y z P P εε==-≤⊂Ω. 记\G B =Ω,则G B ∂=∂Ω⋃∂. 在Green 第二公式中取0(,)v P P =Γ,G Ω=. 由于在区域G 内有0∆Γ=,故有()GGu udV uds n n∂∂Γ∂-Γ∆=-Γ∂∂⎰⎰⎰⎰⎰ 或者()()GBu u udV uds u ds n n n n ∂Ω∂∂Γ∂∂Γ∂-Γ∆=-Γ+-Γ∂∂∂∂⎰⎰⎰⎰⎰⎰⎰ (1.9) 在球面B ∂上,021()414P P r n rr r ππ∂∂Γ∂Γ=-=-=∂∂∂, 因此21(,,)4BBuuds ds u x y z n πε∂∂∂Γ==∂⎰⎰⎰⎰ (1.10)其中(,,)P x y z B ∈∂.同理可得 14BBu u ds ds n n πε∂∂∂∂Γ=∂∂⎰⎰⎰⎰(,,)ux y z n ε∂'''=∂ (1.11) 其中(,,)P x y z B '''∈∂.将(1.10)和 (1.11)代入到(1.9)中并令0ε+→,此时有0(,,)(,,)P x y z P ξηζ→,(,,)0ux y z n ε∂'''→∂,并且区域G 趋向于区域Ω,因此可得()(,,)uudV uds u n nξηζΩ∂Ω∂Γ∂-Γ∆=-Γ+∂∂⎰⎰⎰⎰⎰, 即(,,)()u u u ds udV n n ξηζ∂ΩΩ∂∂Γ=Γ--Γ∆∂∂⎰⎰⎰⎰⎰ (1.12) (1.12)称为Green 第三公式. 它表明函数u 在Ω内的值可用Ω内的u ∆值与边界∂Ω上u 及nu∂∂的值表示.注2 在二维情形,Green 第一公式和Green 第二公式也成立. 而对于Green第三公式, 需要取011(,)ln 2P P rπΓ=,其中0(,)P ξη∈Ω,2(,)P x y R ∈,r =0P P r =0||P P -=此时Green 第三公式也成立.§5⋅2 Laplace 方程基本解和Green 函数基本解在研究偏微分方程时起着重要的作用. 本节介绍Laplace 方程的基本解,并在一些特殊区域上由基本解生成Green 函数,由此给出相应区域上Laplace 方程或Poisson 方程边值问题解的表达式. 下面以Dirichlet 问题为例介绍Laplace 方程的基本解和Green 函数方法的基本思想.5.2.1 基本解设30(,,)P R ξηζ∈,若在点0P 放置一单位正电荷,则该电荷在空间产生的电位分布为(舍去常数0ε)001(,,)(,)4P Pu x y z P P r π=Γ=(2.1) 易证: 0(,)P P Γ在30\{}R P 满足0 .u -∆= 进一步还可以证明[1],在广义函数的意义下0(,)P P Γ满足方程0(,)u P P δ-∆= (2.2)其中0(,)()()()P P x y z δδξδηδζ=---. 0(,)P P Γ称为三维Laplace 方程的基本解.当n =2时,二维Laplace 方程的基本解为0011(,)ln2P PP P r πΓ=(2.3)其中0(,)P ξη,2(,)P x y R ∈,0P P r =同理可证,0(,)P P Γ在平面上除点0(,)P ξη外满足方程0 u -∆=,而在广义函数意义下0(,)P P Γ满足方程0(,)u P P δ-∆= (2.4)其中0(,)()()P P x y δδξδη=--.注1 根据Laplace 方程的基本解的物理意义可以由方程(2.2)和(2.4)直接求出(2.1)和(2.3),作为练习将这些内容放在本章习题中. 另外,也可以利用Fourier 变换求解方程(2.2)和(2.4)而得到Laplace 方程的基本解.5.2.2 Green 函数 考虑如下定解问题(,,), (,,) (2.5)(,,)(,,), (,,)(2.6)u f x y z x y z u x y z x y z x y z ϕ-∆=∈Ω⎧⎨=∈∂Ω⎩设0(,,)P ξηζ∈Ω,21(,,)()()u x y z C C ∈Ω⋂Ω是(2.5)— (2.6)的解,则由Green 第三公式可得(,,)()u u u ds udV n n ξηζ∂ΩΩ∂∂Γ=Γ--Γ∆∂∂⎰⎰⎰⎰⎰ (2.7) 在公式(2.7)的右端,其中有两项可由定解问题(2.5)—(2.6)的边值和自由项求出,即有uds ds n n ϕ∂Ω∂Ω∂Γ∂Γ=∂∂⎰⎰⎰⎰udV f dV ΩΩΓ∆=-Γ⎰⎰⎰⎰⎰⎰.而在u ds n ∂Ω∂Γ∂⎰⎰中,un∂∂在边界∂Ω上的值是未知的. 因此须做进一步处理. 注2 若要求解Neumann 问题,即将(2.6)中边界条件换为(,,)ux y z nϕ∂=∂.此时,在方程(2.7)右端第二项uds n∂Ω∂Γ∂⎰⎰中,u 在边界∂Ω上的值是未知的,而其余两项可由相应定解问题的边值和自由项求出.如何由(2.7)得到定解问题(2.5)-(2.6)的解?Green 的想法就是要消去(2.7) 右端第一项uds n∂Ω∂Γ∂⎰⎰. 为此,要用下面的Green 函数取代(2.7)中的基本解. 设h 为如下定解问题的解0,(,,)(2.8),(,,)(2.9)h x y z h x y z -∆=∈Ω⎧⎨=-Γ∈∂Ω⎩在Green 第二公式中取v h =得()h u h udV uh ds n nΩ∂Ω∂∂-∆=-∂∂⎰⎰⎰⎰⎰ 或者0()u hhu ds h udV n n ∂ΩΩ∂∂=--∆∂∂⎰⎰⎰⎰⎰ (2.10) 将(2.7)和(2.10)相加得(,,)()u G u Gu ds G udV n n ξηζ∂ΩΩ∂∂=--∆∂∂⎰⎰⎰⎰⎰ (2.11) 其中0(,)G P P h =Γ+.由(2.2)和(2.8)—(2.9)可得,0(,)G P P 是如下定解问题的解00(,), (,,)(2.12)(,)0, (,,)(2.13)G P P P x y z G P P P x y z δ-∆=∈Ω⎧⎨=∈∂Ω⎩0(,)G P P 称为Laplace 方程在区域Ω的Green 函数.由于G 在∂Ω上恒为零,由(2.11)可得(,,)Gu uds G udV n ξηζ∂ΩΩ∂=--∆∂⎰⎰⎰⎰⎰ Gds GfdV n ϕ∂ΩΩ∂=-+∂⎰⎰⎰⎰⎰. (2.14) 因此,若求出了区域Ω的Green 函数0(,)G P P ,则(2.14)便是定解问题(2.5)— (2.6)的解.§5⋅3 半空间及圆域上的Dirichlet 问题由第二节讨论可知,只要求出了给定区域Ω上的Green 函数,就可以得到该区域Poisson 方程Dirichlet 问题的解. 对一般区域,求Green 函数并非易事. 但对于某些特殊区域,Green 函数可借助于基本解的物理意义利用对称法而得出. 下面以半空间和圆域为例介绍此方法.5.3.1 半空间上Dirichlet 问题设{(,,)|0},{(,,)|0}x y z z x y z z Ω=>∂Ω==. 考虑定解问题2(,,),(,,) (3.1)(,,0)(,),(,) (3.2)u f x y z x y z u x y x y x y Rϕ-∆=∈Ω⎧⎨=∈⎩设0(,,),P ξηζ∈Ω则1(,,)P ξηζ-为0P 关于∂Ω的对称点. 若在0P ,1P 两点各放置一个单位正电荷,则由三维Laplace 方程的基本解知,它们在空间产生的电位分别为00111(,)41(,)4P P r P P r ππΓ=Γ=其中0011||,||r P P r P P =-=-. 由于0P 和1P 关于∂Ω对称,且1P ∉Ω,故有01001[(,)(,)](,), (,)(,)0,.P P P P P P P P P P PP δ-∆Γ-Γ=∈Ω⎧⎨Γ-Γ=∈∂Ω⎩ 即001(,)(,)(,)G P P P P P P =Γ-Γ为上半空间的Green 函数,且有001(,)(,)(,)G P P P P P P =Γ-Γ011114r r π⎛⎫=- ⎪⎝⎭14π⎡⎤= (3.3)直接计算可得3/2222012()()z G G n zx y ζπξηζ∂Ω=∂∂=-=-∂∂⎡⎤-+-+⎣⎦(3.4)将(3.3)—(3.4)代入到公式(2.14)得(,,)Gu ds Gfd n ξηζϕν∂ΩΩ∂=-+∂⎰⎰⎰⎰⎰ 3/2222001(,)2()() (,)(,,)x y dxdyx y G P P f x y z dxdydzϕζπξηζ∞∞-∞-∞∞∞∞-∞-∞=⎡⎤-+-+⎣⎦+⎰⎰⎰⎰⎰上式便是定解问题(3.1)— (3.2)的解.5.3.2 圆域上Dirichlet 问题设222{(,)|}x y x y R Ω=+<,则222{(,)|}x y x y R ∂Ω=+=. 考虑圆域Ω上的Dirichlet 问题(,), (,) (3.5)(,)(,), (,)(3.6)u f x y x y u x y g x y x y -∆=∈Ω⎧⎨=∈∂Ω⎩设0(,)P ξη∈Ω,1(,)P ξη为0(,)P ξη关于圆周∂Ω的对称点,即201,OP OP R =如图3-1所示 . 由于201OP OP R =,因此对任意M ∈∂Ω有 01~OP M OMP ∆∆ROP r r MP M P ||010=1P01011||P MPMR r OP r =图3.1因此有0101111ln ln 022||P M PMR r OP r ππ-= (3.7) 上式说明函数01001111(,)ln ln22||P P P PR G P P r OP r ππ=- (3.8) 在∂Ω上恒为零. 又由于1P ∉Ω,故有000(,)(,),(,)0,.G P P P P P G P P P δ-∆=∈Ω⎧⎨=∈∂Ω⎩即0(;)G P P 是圆域上的Green 函数.引入极坐标(,)P ρθ,设0000(,)(,)P P ξηρθ=,则21100(,)(,)R P P ξηθρ=. 用α表示0OP 与OP 的夹角,则有000cos cos cos sin sin cos()αθθθθθθ=+=-利用余弦定理可得0P P r =(3.9)1P P r =(3.10)将(3.9)和(3.10)代入到(3.8)中并整理得22222000042220002cos()1(,)ln 42cos()R R R G P P R R ρρρρθθπρρρρθθ+--=-+-- (3.11) 直接计算可得RG Gn ρρ∂Ω=∂∂=∂∂2222000122cos()R R R R ρπρρθθ-=-+-- . (3.12) 记()(cos ,sin )g R R ϕθθθ=,则有00(,)Gu ds Gfd n ρθϕσ∂ΩΩ∂=-+∂⎰⎰⎰ 222022000()()122cos()R d R R πρϕθθπρρθθ-=+--⎰- 222222000422200002cos()1(cos ,sin )ln 42cos()R R R R f d d R R πρρρρθθρθρθρρθπρρρρθθ+--+--⎰⎰(3.13)(3.13)便是定解问题(3.5)—(3.6)的解.注1 当0f =时(3.13)称为圆域上调和函数的Poisson 公式.注2 利用复变函数的保角映射,可以将许多平面区域变换为圆域或半平面.因此,与保角映射结合使用,可以扩大对称法以及Green 函数法的应用范围. 在本章习题中有一些这类题目,Green 函数法更多的应用可查阅参考文献[13].§5⋅4* 一维热传导方程和波动方程半无界问题5.4.1 一维热传导方程半无界问题为简单起见,仅考虑以下齐次方程定解问题20 , 0 , 0 (4.1)(0,)0 , 0 (4.2)(,0)() , 0 t xx u a u x t u t t u x x x ϕ-=<<∞>=≥=<<∞ (4.3)⎧⎪⎨⎪⎩该定解问题称为半无界问题, 这是一个混合问题,边界条件为(4.2). 类似于上节Poisson 方程在半空间和圆域上Dirichlet 问题的求解思想,也要以热方程的基本解为基础,使用对称法求出问题(4.1)—(4.3)的Green 函数,并利用所得到的Green 函数给出该问题的解.一维热传导方程的基本解为224(,)() .x a tx t H t -Γ=(,)x t Γ是如下问题的解20, , 0 (4.4)(,0)(), . (4.5)t xx u a u x t u x x x δ⎧-=-∞<<∞>⎨=-∞<<∞⎩相当于在初始时刻0t =,在0x =点处置放一单位点热源所产生的温度分布.若将上面定解问题中的初始条件换为(,0)()u x x δξ=-,只要利用平移变换'x x ξ=-易得此时(4.4)—(4.5)的解为(,)x t ξΓ-.为求解定解问题(4.1)—(4.3),先考虑()()x x ϕδξ=-,其中ξ为x 轴正半轴上的任意一点. 此时,相当于在x ξ=点处置放一单位点热源. 则此单位点热源在x 轴正半轴上产生的温度分布,如果满足边界条件(4.2),它便是(4.1)—(4.3)的解,即为该问题的Green 函数. 为此,设想再在x ξ=-点,此点为x ξ=关于坐标原点的对称点,处置放一单位单位负热源,这时在x ξ=点处置放的单位点热源产生的温度分布(,)x t ξΓ-和在x ξ=-处置放的单位负热源产生的温度分布(,)x t ξ-Γ+在0x =处相互抵消,从而在0x =处的温度恒为零. 因此,问题(4.1)—(4.3)的Green 函数为(,)(,)(,) G x t x t x t ξξξ-=Γ--Γ+ (4.6) 利用叠加原理可得原问题的解为(,)() (,)u x t G x t d ϕξξξ∞=-⎰ . (4.7)若将(4.2)中的边界条件换为(0,)()u t g t =或(0,)0x u t =,请同学们考虑如何求解相应的定解问题.5.4.2 一维波动方程半无界问题 考虑以下齐次方程定解问题20, 0, 0 (4.8)(0,)0, 0 (4.9)(,0)0, (,0)(), 0 tt xx t u a u x t u t t u x u x x x ψ-=<<∞>=≥==<<∞ (4.10)⎧⎪⎨⎪⎩一维波动方程的基本解(,)x t Γ为1, 2(;) 0, .x ata x t x at ⎧<⎪Γ=⎨⎪≥⎩完全类似于上小节的分析,可得该问题的Green 函数为(,)(,)(,)G x t x t x t ξξξ-=Γ--Γ+, (4.11) 其中0ξ>. 因此,该定解问题的解便可表示为(,)() (,)u x t G x t d ψξξξ∞=-⎰. (4.12)注意到(,)x t ξΓ-的具体表示式为1, 2(;) 0, x at ax t x at ξξξ⎧-<⎪Γ-=⎨⎪-≥⎩类似地有1, 2(;) 0, x at ax t x at ξξξ⎧+<⎪Γ+=⎨⎪+≥⎩将上面两式代入到(4.12)中并整理可得1(), 0 2(,)1(), 0.2x atx at x at at xd x at a u x t d x at a ψξξψξξ+-+-⎧-≥⎪⎪=⎨⎪-<⎪⎩⎰⎰ 若将(4.9)中的边界条件换为(0,)0x u t =,请同学们考虑如何求解相应的定解问题.注1 对一维波动方程半无界问题,除上面使用的Green 函数法以外,也可以用延拓法或特征线法求解[1]. 相比之下,Green 函数法最简单.注2 类似于本章前两节,对一维热传导方程和波动方程初边值问题,也可以建立起解的Green 公式表达式,相当于本章第二节中的(2.14), 并以此为基础而给出上面(4.7)和(4.12)两式的严格证明[2]. 由于本章主要是通过对一些比较简单的偏微分方程定解问题的求解,重点介绍Green 函数法的基本思想和一些特殊区域Green 函数的具体求法,故略去了(4.7)和(4.12)两式的推导过程.习 题 五1.设3R Ω⊂为有界区域,∂Ω充分光滑,21()()u C C ∈Ω⋂Ω. 证明(1)u udV ds n Ω∂Ω∂∆=∂⎰⎰⎰⎰⎰. (2)2u u udV u ds u dV n Ω∂ΩΩ∂∆=-∇∂⎰⎰⎰⎰⎰⎰⎰⎰. 2. 设3R Ω⊂为有界区域,∂Ω充分光滑,21()()u C C ∈Ω⋂Ω满足下面问题0, (,,)(,,)0, (,,).xx yy zz u u u u x y z u x y z x y z ∆=++=∈Ω⎧⎨=∈∂Ω⎩ 证明 (,,)0u x y z ≡,并由此推出Poisson 方程Dirichlet 问题解的唯一性.若将定解问题中的边界条件换为0, (,,),u x y z n∂=∈∂Ω∂问(,,)u x y z 在Ω中等于什么? Poisson 方程Neumann 问题的解是否具有唯一性? 3*设3R Ω⊂为有界区域,∂Ω充分光滑,21()()u C C ∈Ω⋂Ω满足下面问题(,,)(,,), (,,)(,,)(,,), (,,).u c x y z u f x y z x y z u x y z x y z x y z ϕ-∆+=∈Ω⎧⎨=∈∂Ω⎩其中 (,,)c x y z 在闭域Ω非负有界且不恒为零. 证明或求解以下各题(1) 如果0,(,,), 0,(,,),f x y z x y z ϕ=∈Ω=∈∂Ω证明(,,)0u x y z ≡.(2)如果0,(,,),f x y z =∈Ω而边界条件换为0, (,,),u x y z n∂=∈∂Ω∂问(,,)u x y z 在区域Ω中等于什么?4.(1) 验证0∆Γ=,0P P ≠,其中0(,) 3P P n Γ==01(,)22P P n πΓ== (2)设()u u r =, 22y x r +=, 求0,0xx yy u u r +=≠,并且满足(1)0, u = (0,) 1B u n ds δ∂∇⋅=-⎰的解, 其中(0,)B δ是以原点为圆心δ为半径的圆形域,n 为(0,)B δ∂的单位外法向量.(3) 设()u u r =, 222z y x r ++=, 求0=++zz yy xx u u u ,0≠r ,并且满足B(0,)lim ()0, 1r u r u nds δ→∞∂=∇⋅=-⎰⎰的解, 其中(0,)B δ是以原点为球心δ为半径的球形域,n 为(0,)B δ∂的单位外法向量.5. 设2R Ω⊂有界区域,∂Ω充分光滑,21()()u C C ∈Ω⋂Ω. 证明(,)()u u u ds ud n n ξησ∂ΩΩ∂∂Γ=Γ--Γ∆∂∂⎰⎰⎰ 其中0(,)P ξη∈Ω,0(,)P P Γ如第4题所示.6. 设2R Ω⊂有界区域,∂Ω充分光滑,0(,)P ξη∈Ω,2(,)P x y R ∈,0(,)P P Γ为二维Laplace 方程的基本解. 考虑定解问题(,), (,)(,)(,), (,)u f x y x y u x y x y x y ϕ-∆=∈Ω⎧⎨=∈∂Ω⎩ 若(,)h x y 是如下定解问题的解00, (,)(,)(,),(,)h x y h x y P P x y ∆=∈Ω⎧⎨=-Γ∈∂Ω⎩证明 若21(,)()()u x y C C ∈Ω⋂Ω,则有(,)G u ds Gfd n ξηϕσ∂ΩΩ∂=-+∂⎰⎰⎰, 其中G h =Γ+.7. 设3R Ω⊂有界区域,∂Ω充分光滑, 考虑定解问题(,,), (,,)(,,), (,,).u f x y z x y z u x y z x y z nϕ-∆=∈Ω⎧⎪∂⎨=∈∂Ω⎪∂⎩ 证明该问题可解的必要条件为0f dV ds ϕΩ∂Ω+=⎰⎰⎰⎰⎰. 8* 证明上半空间Laplace 方程Dirichlet 问题的Green 函数0(,)G P P 满足 020010(,), (,),0, .4P PG P P x y R z P P r π<<∈>≠ 对平面上圆域Laplace 方程Dirichlet 问题的Green 函数0(,)G P P ,给出类似结果.9. 利用对称法求二维Laplace 方程Dirichlet 问题在上半平面的Green 函数, 并由此求解下面定解问题0, (,),0(,0)(), (,).u x y u x x x ϕ-∆=∈-∞∞>⎧⎨=∈-∞∞⎩10. 求二维Laplace 方程在下列区域上 Dirichlet 问题的Green 函数.(1) {(,)|}x y x y Ω=>. (2) {(,)|0,0}x y x y Ω=>>.11. 设222{(,)|,0}x y x y R y Ω=+<>. 考虑半圆域Dirichlet 问题0,(,)(,)(,), (,).u x y u x y x y x y ϕ-∆=∈Ω⎧⎨=∈∂Ω⎩ 应用对称法求区域Ω上的Green 函数.12*求解定解问题 0,(,,)(,,)(,,),(,,).u x y z u x y z g x y z x y z -∆=∈Ω⎧⎨=∈∂Ω⎩ 其中32222,(0,){(,,)|}xx yy zz u u u u B R x y z R x y z R ∆=++Ω==∈++<.13.[解对边值的连续依赖性]设Ω为半径等于R 的圆域,考虑如下问题(,), (,)(,)(,),(,) 1,2.k k k u f x y x y u x y g x y x y k -∆=∈Ω⎧⎨=∈∂Ω=⎩ 利用Poisson 公式证明2121(,)(,)max{(,)(,)(,)}u x y u x y g x y g x y x y -≤-∈∂Ω14*证明在广义函数的意义下,11(,0)ln 2P rπΓ=满足 ()()u x y δδ-∆=, 其中xx yy r u u u =∆=+.15*设Ω为半径等于R 的圆域,考虑如下问题 0, (,)(,)(,),(,) .u x y u x y g x y x y -∆=∈Ω⎧⎨=∈∂Ω⎩ 如果(,)g x y 在∂Ω连续,证明由Poisson 公式给出的解是该问题的古典解(真解).16*设(,)u x y 为平面上区域Ω上的调和函数,000(,)P x y ∈Ω且0(,)B P R ⊂Ω. 证明调和函数的平均值公式00002(,)(,)11(,)(,)(,)2B P R B P R u x y u x y ds u x y dxdy R R ππ∂==⎰⎰⎰ 17*[极值原理]设2R Ω⊂有界区域,边界充分光滑,2()()u C C ∈Ω⋂Ω为Ω 内的调和函数,并且在某点000(,)P x y ∈Ω达到u 在闭域Ω上的最大(小)值,利用平均值公式证明u 为常数.18*[极值原理]设2R Ω⊂有界区域,边界∂Ω充分光滑, 2()()u C C ∈Ω⋂Ω. 如果u 在区域Ω内调和且不等于常数,则u 在闭域Ω上的最大值和最小值只能在区域的边界∂Ω上达到.19*利用第12题的结果,建立在3R Ω⊂内调和函数的平均值公式,并证明和第16题类似的结果.20*设2R Ω⊂有界区域,2()(), (),1,2,k k u C C g C k ∈Ω⋂Ω∈∂Ω=满足 (,), (,)(,)(,),(,) k kk u f x y x y u x y g x y x y -∆=∈Ω⎧⎨=∈∂Ω⎩ 证明 2121(,)(,)max{(,)(,)(,)}u x y u x y g x y g x y x y -≤-∈∂Ω.21.设D 和Ω为平面上的两个区域,()(,)(,)f z x y i x y ϕψ=+在区域D 内解析且不等于常数,()f D =Ω,即f 将区域D 保形映射到区域Ω.证明 如果(,)u x y 在区域Ω内调和,则((,),(,))u x y x y ϕψ在区域D 内调和.22.(1)找一个在上半平面解析的函数()f z ,在边界{(,),0}x y x R y ∈=上满足00(),, (),,f x A x x f x B x x =>=<其中A 和B 为实常数.(2)求下面定解问题的一个解0, 0,0(,0)0,0, (0,)10,0.xx yy u u x y u x x u y y +=>>⎧⎨=>=>⎩23*求下面定解问题的一个解22220, 1(,)0,0, (,)1,0, 1.xx yy u u x y u x y y u x y y x y ⎧+=+<⎪⎨=<=>+=⎪⎩ 24. 求下面定解问题的一个解0, 0<(,0)0, (,)1, 0.xx yy u u y x u x u x x x +=<⎧⎨==>⎩ 25. 求下面定解问题的一个解0, , 0<(,)0, (,0)0, 0, (,0)1, 0.xx yy u u x R y u x x Ru x x u x x ππ+=∈<⎧⎪=∈⎨⎪=<=>⎩26. 设(0,)B R Ω=,1(0,)2R B Ω=,(,)u x y 在Ω内调和且在Ω上连续,在边界上非负,证明以下结果(1)(,),x y ∀∈Ω有(0,0)(,)(0,0),R r R r u u x y u R r R r-+≤≤+-其中r = (2)存在常数0M > 使得 11max (,)min (,).u x y M u x y ΩΩ≤。
数值计算:Legendre多项式

数值计算:Legendre多项式Legendre多项式的概念以及正交特性在此不多作描述,可以参考数学物理⽅程相关教材,本⽂主要讨论在数值计算中对于Legendre多项式以及其导数的计算⽅法。
Legendre多项式的计算递推公式(n+1)P n+1(x)=(2n+1)x⋅P n(x)−nP n−1(x)(n≥2)通式可以⽤幂级数表⽰为以下形式:P n(x)=[n2]∑k=0(−1)k(2n−2k)!2n k!(n−k)!(n−2k)!x n−2k n=0,1,2,⋯Legendre多项式前⼏项n01234P n(x)1x12(3x2−1)12(5x3−3x)18(35x4−30x2+3)在实际数值计算中,按照通项公式计算P n(x)会涉及到⼤量的乘除法运算,同时由于数据字节长度的限制,对于基数较⼤阶乘的运算会导致数据的溢出,因此,在实际的计算中,使⽤递推公式计算P n(x)更为合适。
%% Legendre多项式Pi(x)function L=myLegendre(N,x)% 返回Pi(x),i=0~N ,x为标量if N<2error("N is less then 2 !")endL=zeros(1,N);%⽣成Legendre多项式矩阵L(1,1)=1;L(1,2)=x;for i=2:N-1L(1,i+1)=((2*i-1)*x*L(1,i)-(i-1)*L(1,i-1))/i;endendLegendre多项式的导数对于x的⼀阶导数P′n(x)存在类似的递推公式(n+1)P′n+1(x)=(2n+1){x⋅P′n−1(x)+P n−1(x)}−nP′n−1(x)(n≥2)化简得到⼀阶导数项的递推公式与之前形式类似,但包含额外⼀项(2n+1)P n−1(x)(n+1)P′n+1(x)=(2n+1)x⋅P′n−1(x)−nP′n−1(x)+(2n+1)P n−1(x)(n≥2)Legendre多项式前⼏项n01234P′n(x)013x 12(15x2−3)12(35x3−15x)%% Legendre多项式Pi(x)对于x的⼀阶导数function L_x=myLegendre_x(N,x)L=myLegendre(N,x);L_x=zeros(1,N);%⽣成Legendre_x多项式矩阵L_x(1,1)=0;L_x(1,2)=1;for i=2:N-1L_x(1,i+1)=((2*i-1)*x*L_x(1,i)-(i-1)*L_x(1,i-1)+(2*i-1)*L(1,i))/i; endend计算结果%%demox=0.5;N=100;c=myLegendre(N,x);c_x=myLegendre_x(N,x);figure(2)subplot(2,1,1)plot(c,'-*');xlabel('N');title("P(x)") subplot(2,1,2)plot(c_x,'-o');xlabel('N');title("P'(x)") Processing math: 100%。
Legendre多项式及推广

定理 1只∽在实数域上有/ : / , 个零点 , 且全部是在 区间[1 1 - , 1
内部 相 异 的 。
15 年 , 9 6 莫叶曾用L gnr. 数 eede  ̄
∑ ( )
) 都有实根 1 和
() 9
证明 : 引入 函数g ∞= ( -)  ̄ m-,,, n 1 , x 1 , i - 1 …,一 时 各 Z " O 2
() 4
定 2 < , 触 E 则 L e r 数 乎 理: 若{ * 且 ) , 它的 e ne 几 处 g d级
处收敛。
证 明 见H 。 5 整 函 数
递 推 公 式嘲 (+)+ ) ( +) ( 一 一 ) =,… 5 : 1 . =2 1 1 ,)( ) 只( ) ( 2
( )已知函数如何求它的L gnr ̄数展开式 ; 1 eede
() 2 已知L gnr级数如何求它所表示的函数 的性质。 eede 设 函数 ) 实 数轴 的闭 区间 [11上 可积 , 在 - ,】 - 贝 L gnr级数展开式如下 : eede
() 2
l n , 正 数。 1 为 奇 1 二
e de n r ̄项 式 的运 算和 插 值提 供 方便 。 对 gn r级 数进 行 一 个概 括 的讲 述 , 结 关 于此 的 结 果 , 对 用L gn rt ̄ 表 示 ede 总 并 ee de :
整 函数 进 行 系统 的论 述 。 最后 介 . L g n r  ̄ e e de多项 式 的推 广 和 其 应 用 。
它 的 Lpae 分 表 示 是 alc积
确定 ( ) 即要 确定系数 c , 7, n 上式 两端 同乘尸 ∽并 在 区间 m
Legendre 多项式

o 當n=1時
當 n = 1 時,Legendre 方程式之通解的基底為:
8 4 2 y1 = 1 − x − x + L 24 y2和,以此為目標,討論 1 階之 Legendre 多項式 P1 ( x ) ,也就是 說,考慮:
P1 ( x ) = ky 2
上式可再整理為:
2a2 + n(n + 1)a0 = 0, 6a + (n + 2 )(n − 1)a = 0, 3 1 12a 4 + (n + 3)(n − 2 )a 2 = 0, 20a + (n + 4 )(n − 3)a = 0, 5 3 M
( 2a ) (2b) (2c) ( 2d )
6
1
=−
(n + 2)(n − 1) a
3!
1
a5 = −
(n + 4)(n − 3) a
20
3
=
(n − 3)(n − 1)(n + 2)(n + 4) a
5!
1
最後再將 a0 、 a1 、 a2 、…等代回原式,則:
y ( x ) = a0 + a1 x + a 2 x 2 + L = a0 + a1 x −
n n n +1
Pn′( x )
14. (n + 1)Pn+1 ( x ) = (2n + 1)xPn ( x ) − nPn−1 ( x ) , n = 1, 2, 3, L 15. Pnm ( x ) = (1 − x 2 )
m2
d m [Pn ( x )] dx m
- 38 -
legendre多项式

∂z2 = r5 =
r3
∂3V
15z3 − 9zr2
15 cos3 θ − 9 cos θ
∂z3 = −
r7
=−
r4
We notice that ∂nV /∂zn is an n-th degree polynomial of cos θ divided by rn+1. One formula which combines
Some equations relating (x, y, z) to (r, θ, φ) are,
x = r sin θ cos φ
(4)
y = r sin θ sin φ
(5)
z = r cos θ
(6)
r2 = x2 + y2 + z2
Байду номын сангаас
(7)
From (7) it follows that
∂r z
Legendre Polynomials
• Introduced in 1784 by the French mathematician A. M. Legendre(1752-1833).
• We only study Legendre polynomials which are special cases of Legendre functions. See sections 4.3, 4.7, 4.8, and 4.9 of Kreyszig.
1
Definition:
1 x2 + y2 + (z
− h)2
=
1 √
r2 − 2rh cos θ + h2
=
∞ n=0
数学物理方法_第5章 Legendre多项式(上)

在第三章中,我们已经求出了Legendre方 程(5.1.12)的通解,并且指出,Legendre方 程(5.1.12)加上自然条件 (5.1.13) 构成本征值问题,其本征值和本征函数依 次是
l n
1 ,
n 0,1, 2, ;
x Pn x .
2 m 2 (1 x )( x) 2 x( x) l (l 1) ( x) 0 (5.1.7c) 2 1 x
至此球坐标系下的Laplace方程(5.1.4) 分离变量的结果是得到三个常微分方程 (5.1.5)、(5.1.6)和(5.1.7). 方程 (5.1.5)加上周期性条件 (5.1.8) (0) (2 ) 构成本征值问题,解之得到 ( ) Am cos m Bm sin m, m 0,1, 2, .
n
的
(2n)! (2n)! n P2n ( x) P2n ( x); P2n (0) (1) n 2 (1) . 2 (2 n!) [(2n)!!]
P P2n1 (0) 0. 2n1 ( x) P 2n1 ( x);
2 Legendre多项式的积分表示 i. 施列夫利(Schlufli)积分 根据复变函数的Cauchy积分公式
(5.1.6a)
以及
1 d d m2 sin l (l 1) 2 0 sin d d sin
(5.1.7a)
对方程(5.1.7a)作变换
x cos 可得
即
2 (5.1.7b) d d m 2 0 (1 x ) l (l 1) 2 dx dx 1 x
5.1.1 5.1.2 5.1.3
legendre 谱方法

legendre 谱方法legendre谱方法是数学和物理学领域中一种用于解决某些特殊类型微分方程的方法。
起源于法国数学家Adrien-MarieLegendre,它也被称为Legendre多项式展开法。
Legendre谱方法最初用于解决常见的物理学问题,随着数学发展,它也用于解决更复杂的问题。
Legendre谱方法源于Legendre的基础理论,他在1830年发表的论文《Essai sur la theorie des fonctions analytiques》中提出了Legendre谱展开方法。
Legendre提出了一个方法,该方法将某种特定的函数f(x)从它的原始表达式展开为Legendre多项式的级数形式,这使他们能够以更容易解决的方式解决微分方程。
Legendre 多项式展开可用于解决积分方程,而Legendre谱方法则可用于解决带有因子的积分方程。
Legendre谱方法分为两个步骤:第一步是用Legendre多项式描述给定函数,第二步是将函数用Legendre谱表示法表示,即使用Legendre谱系数和Legendre多项式表示。
Legendre多项式可以表示为可以分解为两个Legendre多项式的形式,其中一个Legendre多项式用来表示函数的模数,另一个Legendre多项式用来表示函数的相位。
当一个函数以Legendre谱形式表示时,可以利用Legendre谱方法来解决不同类型的微分方程。
Legendre谱方法可用于解决波动方程,正弦方程,哈密尔顿方程,发散和收敛方程,HertzJacobi方程等。
在解决这类方程时,Legendre多项式是一种非常有用的工具,因为它可以用来表示XX X X函数的局部特征。
因此,Legendre谱方法可以被用来求解许多经典的物理问题,如有限深水内部波和非线性波的解析解等。
Legendre谱方法的优点之一是,它的解可以以数值的形式表示。
数值表示可以减少解的复杂性,因此可以更容易地解决问题。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
5.1.1 5.1.2 5.1.3
现在求解这个定解问题。在球坐标系 (r , , ) 2 中,Laplace方程 u(r, , ) 0 的表达式是
1 2 u 1 u 1 2u 0 r 2 sin 2 2 2 2 r r r r sin r sin
方程(5.1.7)叫做关联Legendre方程。在 m 0 时,它就退化为Legendre方程 (5.1.12) (1 x2 )( x) 2x( x) l (l 1)( x) 0,
(注意这里未知函数用 表示)这种退化, 有着真实的物理含义,它是物理问题具有
轴对称的反应。所谓轴对称问题即场量 u 与角度 无关,只是 r和 的函数。那么, 这会在什么情况下发生呢?重新考虑定解 问题(5.1.1)~(5.1.3)。非齐次边界 条件(5.1.2)是引起场量 u 发生变化的 唯一根源,如果这个非齐次函数不是角变量 的函数,则问题就具有轴对称性,我们 讨论的问题符合这个条件。
要确定其中的未知常数,需要进一步了解 Legendre多项式的性质,将在下一节中讨论。
当 m 0时,关联Legendre方程(5.1.7c) 的解也可用Legendre多项式表示。事实上, 对方程(5.1.7c)作变换
( x) (1 x )
2 m/ 2
Y ( x)
(5.1.17)
则函数 Y 满足方程
方程(5.1.6)是Euler方程,做变换 或直接用Maple求解得
r e
t
>dsolve(r^2*diff(R(r),r$2)+2*r* diff(R(r),r)-l*(l+1)*R(r)=0);
R( r )_C1 r _C2 r
l
( l1 )
_ C1和 _ C 2 是两个积分常数。或用 A , B 其中 l l 表示任意常数,将通解写成 l 1 l (5.1.10) R r Al r Bl r
P l ( x) 2l l ! dxl ( x2 1)l
2 l ( x 1) 按二项式 证明:将式(5.2.10)中的
定理展开,可得
k l/2 1 dl 2 ( 1) (2l 2k ) (l 2k 1) l 2 k l ( x 1) x l l l 2 l ! dx 2 k !(l k )! k 0 l/2 (2l 2k )! k l 2 k (1) l x 2 k !(l k )!(l 2k )! k 0
例1 在本来匀强的静电场 E 0 中,放置一个 导体球,球的半径为a,试研究导体球怎样 改变了匀强电磁场。 对比“大地与带电云的电场中导线”的例题, 这里是三维静电场问题,球外电势满足 Laplace方程,在距球无穷远处,电场保持 为原来的 E 0以球心为原点取球坐标系, 则定解问题是
2u 0 r a , u r E0 z E0 r cos , u r a 0.
可得
2
sin 2 d 2 dR sin d d 2 r sin m R dr dr d d
其中 m 是分离常数。由此可得两个微分方程 2 m 0 (5.1.5) 和 2
1 d 2 dR 1 d d m r sin 2 0 R dr dr sin d d sin
un r , An r n Bn r n1 Pn cos
n 0,1, 2,
(5.1.15)
将这些解叠加起来,得到级数解为
u r , An r n Bn r n1 Pn cos .
n 0
(5.1.16)
在我们讨论的问题中,自然条件(5.1.13) 是必需的。因为这里(即球坐标系下Laplace x cos , x 1对应 0和 方程的分离变量法中) 我们当然要求物理量 u 在各个方向上都有 有限值。
综上,定解问题(5.1.1)----(5.1.3)在具有 轴对称性质(即 u与 无关)的假设下, 具有如下所示的一般解称本征解
f
( n)
n! ( x) 2 i
(5.1.6a)
以及
1 d d m2 sin l (l 1) 2 0 sin d d sin
(5.1.7a)
对方程(5.1.7a)作变换
x cos 可得
即
2 (5.1.7b) d d m 2 0 (1 x ) l (l 1) 2 dx dx 1 x
用分离变量法求解,设 代入到方程中,得
u (r , , ) R(r )( )( )
将关于 的变量和关于 r , 的变量分离 r 2 sin 2 为此,用 遍乘上式,并适当移项,
R
d 2 dR R d d R d 2 0 r 2 sin 2 2 2 2 r dr dr r sin d d r sin d
Rodrigues公式由此得证。 利用Rodrigues公(5.2.10),可方便地 给出低阶的几个Legendre多项式的显式
P0 1 P 1 x cos 3x 2 1 3cos 2 1 P2 ( x) 2 23 3 5 x 3x 5cos 3cos P3 ( x) 2 2 1 1 4 2 P4 ( x) (35x 30 x 3) (35cos 4 20cos 2 9) 8 64 1 1 5 3 P5 ( x) (63x 70 x 15x) (63cos5 35cos3 30cos ) 8 128
对上面第二个方程,再将变量
r , 分离开来
1 d d m2 1 d 2 dR sin 2 r l (l 1) sin d d sin R dr dr
其中 l (l 1) 是第二次分离变量引入的常数, 它可以用一个实数 表示,为了讨论方便
d Pn ( x) x Pn ( x) (1 x ) (m 0,1, 2, m dx
mm 2 2m Nhomakorabea, n)
( 5.1.19)
一般解 u (r , , ) 应该是按本征函数序列的 叠加,并且是二重求和,求和指标为 n, m 首先本征解是
unm (r , , ) [ Ar n Br ( n1) ] [Cm cos m Dm sin m ]Pn m (cos )
利用Maple画出它们的图形如下
>plot({LegendreP(0,x),LegendreP(1,x),Leg endreP(2,x),LegendreP(3,x),LegendreP(4,x) ,LegendreP(5,x)},x=-1..1);
l
由图可见, P l (1) 1; P l ( x) 的奇偶性由l 奇偶性来决定。
在第三章中,我们已经求出了Legendre方 程(5.1.12)的通解,并且指出,Legendre方 程(5.1.12)加上自然条件 (5.1.13) 构成本征值问题,其本征值和本征函数依 次是
l n
1 ,
n 0,1, 2, ;
x Pn x .
而一般解就是
u (r , , ) unm (r , , )
n 0 m 0 n
§5.2 Legendre 多项式的性质 1、Legendre多项式的微分表示 我们已知,当 l为偶数时 Legendre多项式 l/2 是: (2l 2k )! k l 2 k
Pl ( x) (1)
2 m 2 (1 x )( x) 2 x( x) l (l 1) ( x) 0 (5.1.7c) 2 1 x
至此球坐标系下的Laplace方程(5.1.4) 分离变量的结果是得到三个常微分方程 (5.1.5)、(5.1.6)和(5.1.7). 方程 (5.1.5)加上周期性条件 (5.1.8) (0) (2 ) 构成本征值问题,解之得到 ( ) Am cos m Bm sin m, m 0,1, 2, .
令 l (l 1) [可以证明任意一个实数 都可表为 l (l 1) ,其中 l 为另一任意实数 或复数]。由此再得两个常微分方程 即 r R(r ) 2rR(r ) l (l 1) R(r) 0 (5.1.6b)
2
d 2 dR r l (l 1) R 0 dr dr
5.1 Legendre方程和Legendre多项 式的引出
在这一章中,我们讨论另一个重要的特殊函 数----Legendre多项式。首先将通过在球坐标 系中对Laplace方程分离变量,引入在第三章 中曾讨论过的Legendre方程,及其解 Legendre函数。并深入讨论Legendre方程在 区间 1,1 上有界解构成的另一类正交函数 系----Legendre多项式。 首先我们从一个实际例子出发,引出 Legendre方程。
n
的
(2n)! (2n)! n P2n ( x) P2n ( x); P2n (0) (1) n 2 (1) . 2 (2 n!) [(2n)!!]
P P2n1 (0) 0. 2n1 ( x) P 2n1 ( x);
2 Legendre多项式的积分表示 i. 施列夫利(Schlufli)积分 根据复变函数的Cauchy积分公式
k 0
当
l为奇数时,有
Pl ( x)
( l 1) / 2
x 2 k !(l k )!(l 2k )!
l
k 0