数学物理方法_第5章 Legendre多项式(上)
legendre多项式推导
![legendre多项式推导](https://img.taocdn.com/s3/m/ddfdfbf564ce0508763231126edb6f1aff007197.png)
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 多项式作为基函数。
数学物理方程-第五章格林函数法
![数学物理方程-第五章格林函数法](https://img.taocdn.com/s3/m/5f427104cc175527072208bb.png)
第五章 格林函数法在第二章中利用分离变量法求出了矩形区域和圆域上位势方程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多项式](https://img.taocdn.com/s3/m/a9cf0fe6bb0d4a7302768e9951e79b8968026860.png)
数值计算: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多项式及推广
![Legendre多项式及推广](https://img.taocdn.com/s3/m/f43334cf58f5f61fb73666e7.png)
定理 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 多项式
![Legendre 多项式](https://img.taocdn.com/s3/m/2a7d2849cf84b9d528ea7a6d.png)
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多项式
![legendre多项式](https://img.taocdn.com/s3/m/b7733ade3186bceb19e8bbba.png)
∂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多项式(上)
![数学物理方法_第5章 Legendre多项式(上)](https://img.taocdn.com/s3/m/621d28cb8bd63186bcebbce3.png)
在第三章中,我们已经求出了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 谱方法](https://img.taocdn.com/s3/m/8ee7ec6e68eae009581b6bd97f1922791788be44.png)
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谱方法的优点之一是,它的解可以以数值的形式表示。
数值表示可以减少解的复杂性,因此可以更容易地解决问题。
数学物理第五章-格林函数法
![数学物理第五章-格林函数法](https://img.taocdn.com/s3/m/c3de6e7030126edb6f1aff00bed5b9f3f80f7249.png)
则有
(P[, P(0P),P0)(P, P(1P),P10),]
P
(P,
P0
),
P
即
G(P, P0) (P, P0) (P, P1)
为上半空间的格林函数,且有
G(P, P0) (P, P0) (P, P1)
1(
4
1 (1 1)
4 r0 r1
1
(x )2 (y )2 (z )2
( x,
y)
1
0
x2 x2
0 y2 y2
1 1
的解并求出u(0, 0, a)的值(常数a 0)?
解 由上半空间的泊松公式
1 )
(x )2 (y )2 (z )2
G 1 (
1
1
)
4 (x )2 ( y )2 (z )2 (x )2 ( y )2 (z )2
直接计算可得
G
G
1
n
z
z0
2
3
[(x )2 ( y )2 2 ]2
又
u(
,,
)
G n
ds
GfdV
1
2
(x, y)
3
R2 [(x )2 ( y )2 2 ]2
因此
u
B
n
ds
B
u
4
2
ds
u(x,
y,
z
)
(10)
u
u
B
n
ds
B
4
2
ds
u(x,
y,
z
)
(10)
其中 P(x, y, z ) B.------积分中值定理
同理可得
B
u n
legendre多项式递推公式推导
![legendre多项式递推公式推导](https://img.taocdn.com/s3/m/ec93c7c40342a8956bec0975f46527d3240ca6ad.png)
legendre多项式递推公式推导以Legendre多项式递推公式推导为标题的文章一、引言在数学中,Legendre多项式是一类具有重要应用的正交多项式。
它是法国数学家Adrien-Marie Legendre在18世纪末通过研究椭圆函数而引入的。
Legendre多项式在物理学、工程学和数学分析等领域中都有广泛的应用。
本文将详细介绍Legendre多项式的递推公式推导过程。
二、Legendre多项式的定义Legendre多项式是通过一种特殊的递推关系得到的。
它是通过以下方式定义的:1. Legendre多项式的零次项为P0(x) = 1;2. Legendre多项式的一次项为P1(x) = x;3. 对于n≥2,Legendre多项式可以通过如下递推公式得到:Pn(x) = (2n-1)xPn-1(x) - (n-1)Pn-2(x)三、推导过程为了推导出Legendre多项式的递推公式,我们将利用以下两个性质:1. 递推关系的推导过程中,我们使用的是首项系数为1的Legendre多项式;2. 我们将递推关系中的Pn(x)展开成幂级数形式。
我们可以用幂级数形式表示Legendre多项式Pn(x):Pn(x) = a0 + a1x + a2x^2 + ...其中,a0、a1、a2等为常数。
然后,我们将递推关系Pn(x) = (2n-1)xPn-1(x) - (n-1)Pn-2(x)代入上述幂级数展开式中:a0 + a1x + a2x^2 + ... = (2n-1)x(a0 + a1x + a2x^2 + ...) - (n-1)(a0 + a1x + a2x^2 + ...)我们可以对比幂级数展开式中各个幂次的系数,得到递推关系:a0 = - (n-1)a0a1 = (2n-1)a0 - (n-1)a1a2 = (2n-1)a1 - (n-1)a2...可以发现,递推关系中的每一项都与前两项有关。
Legendre多项式
![Legendre多项式](https://img.taocdn.com/s3/m/715bbcf5aef8941ea76e056c.png)
1 / 21
Legendre 方程 L-多项式的表示和性质 L-多项式的应用 连带 Legendre 方程 连带 L-多项式的表示和性质 球函数
.
.....
......
.
...
...
目录
1. Legendre 方程
2. L-多项式的表示和性质
3. L-多项式的应用
4. 连带 Legendre 方程
Legendre 方程 L-多项式的表示和性质 L-多项式的应用 连带 Legendre 方程 连带 L-多项式的表示和性质 球函数
.
.....
......
.
...
...
《数学物理方法:特殊函数》 § 2.6 球函数: Legendre 多项式
主讲人:卢东强
上海大学 上海市应用数学和力学研究所
/portals/419/dqlu /dqlu
Pnl (x) 的微分表示 (Rodrigues 公式):
Pnl (x)
=
(1
−
x2)n/2
1 2ll!
dl+n dxl+n
(x2
−
1)l.
Pnl (x) 的复积分表示 (Schl¨afli 公式):
Pnl (x)
=
(1
−
x2)n/2 (l
+ n)! l!
1 2πi
∮
C
2l
(z2 − 1)l (z − x)l+n+1
.Fourier 级数.
解: f (x) 为三次多项式, 应表示成 P0(x), P1(x), P2(x), P3(x) 的线性组合:
∑3 2x3 + 3x + 4 = flPl(x)
第5章格林函数法
![第5章格林函数法](https://img.taocdn.com/s3/m/dd7963d9d15abe23482f4dff.png)
第5章格林函数法格林(Green)函数,又称为点源影响函数,是数学物理中的一个重要概念.格林函数代表一个点源在一定的边界条件下和初始条件下所产生的场.知道了点源的场,就可以用叠加的方法计算出任意源所产生的场.格林函数法是解数学物理方程的常用方法之一.5.1 格林公式TΣ上具有连续一阶导数,在区域及其边界中具有连续二阶导数,应用矢量分析的高斯定理d d T T div =∇∫∫∫∫∫∫i A V =A V (5.1.1)单位时间内流体流过边界闭曲面S 的流量单位时间内V 内各源头产生的流体的总量将对曲面Σ的积分化为体积分d ()d d d T T Tu u V u V u V Σ∇=∇∇=Δ+∇∇∫∫∫∫∫∫∫∫∫∫∫i i i S v v v v (5.1.2)()uv u v u v∇=∇⋅+∇以上用到公式称上式为第一格林公式.同理有d ()d d d T T T u u V u V u V Σ∇=∇∇=Δ+∇∇∫∫∫∫∫∫∫∫∫∫∫i i i S v v v v (5.1.3)上述两式相减得到()d ()d Tu u u u V Σ∇−∇=Δ−Δ∫∫∫∫∫i S v v v v的外法向偏导数.5.1.4)为第二格林公式.进一步改写为()d ()d Tu S u u V n Σ∂∂−=Δ−Δ∂∂∫∫∫∫∫ v u v v v n (5.1.4)5.2 泊松方程的格林函数法讨论具有一定边界条件的泊松方程的定解问题.泊松方程()() u f Δ=−r r (5.2.1)(5.2.2)是区域边界Σ上给定的函数.是第一、第二、第三类边界条件的统一描述典型的泊松方程(三维稳定分布)边值问题()()[]()u f u u n αβϕΣΣΔ=−⎧⎪∂⎨+=⎪∂⎩r r r (5.2.3)上沿界面外法线方向的偏导数格林函数的引入及其物理意义引入:为了求解定解问题(5.2.3),我们必须定义一个与此定解问题相应的格林函数0(,)G r r 它满足如下定解问题,边值条件可以是第一、二、三类(,)()[]0G G G n δαβΣΔ=−−⎧⎪∂⎨+=⎪∂⎩00r r r r (5.2.4)()δ−0r r 代表三维空间变量的δ函数,在直角坐标系中其形式为0()()()()x x y y z z δδδδ−=−−−r r 函数前取负号是为了以后构建格林函数方便格林函数的物理意义【2】:在物体内部(T 内)0r 处放置一个单位点电荷,而该物体的界面保持电位为零, 那么该点电荷在物体内产生的电势分布,就是定解问题(5.2.4)的解――格林函数.由此可以进一步理解通常人们为什么称格林函格林函数互易定理:因为格林函数0(,)G r r 代表0r 处的脉冲(或点源)在r 处所产生的影响(或所产生的场),所以它只能是距离0||−r r 的函数,故它应该遵守如下的互易定理:(,)()G G ,=r r r r (5.2.5))得到())d (()())d T u S u G G u V n ∂⋅=Δ−Δ∂∫∫∫r r r (5.2.6)0()]d (()())d ())()()]d T G u S G u u G Vf u V δ∂−⋅=Δ−Δ∂−+−∫∫∫r r r r r r r n (5.2.7)根据δ函数性质有:00()()]d ()T u V u δ−=∫∫∫r r r r (5.2.8)故有0(,)()]d G u S ∂−∂r r r)r n (5.2.9)泊松方程的基本积分公式.00000000((,))d [(,)()]d u G V G u S n Σ∂∂+−∂∂∫∫ r )r r r r r n 格林函数满足互易定理并利用格林函数的对称性则得到(5.2.10)解的基本思想:通过上面解的形式(5.2.9)我们容易观察出引用格林函数的目的:主要就是为了使一个非齐次方程(5.2.1)与任意边值问题(5.2.2)所构成的定解问题转化为求解一个特定的边值问题(5.2.4). 一般后者的解容易求得,通(5.2.9)即可求出(5.2.1)和(5.2.2)定解问题的解.考虑格林函数所满足的边界条件讨论如下:1.第一类边值问题:()()|()u f u ϕΣΣΔ=−⎧⎨=⎩r r r (5.2.11)相应的格林函数0(,)G r r 是下列问题的解:000(,)(-)(,)|0 G G δΣΔ=−⎧⎨=⎩r r r r r r (5.2.12)考虑到格林函数的齐次边界条件,由公式(5.2.9)可得第一类边值问题的解000(,)()(,)()d ()d T G u G f V S ϕΣ∂=−∂∫∫∫∫∫ nr r r r r r r (5.2.13)另一形式的第一类边值问题的解000(,)()d G S ∂∂0n r r r (5.2.5)2.第二类边值问题()()|()p u f unϕΣΔ=−⎧⎪∂⎨=⎪∂⎩r r r 是下列问题的解:(5.2.15)00,)|0n Σ=r (5.2.16)5.2.9)可得第二类边值问题解00(,)()d ()(,)d G f V G SϕΣ+∫∫ r r r r r r (5.2.17)3.第三类边值问题()() []()p u f u u n αβϕΣΔ=−⎧⎪∂⎨+=⎪∂⎩r r r 是下列问题的解:(5.2.18)0(,)]0G G n βΣ∂+=∂r r (5.2.19)边值条件,两边同乘以格林函数G(5.2.19)的边值条件的两边同乘以函数u得[]0Gu G nαβΣ∂+=∂G ϕ[]()p uG u G nαβϕΣ∂+=∂r )得到第三类边值问题的解001,)()d ((,)d f V G S ϕβΣ+∫∫ r r r r)r r (5.2.20)格林函数的互易性则得到000001)()d ()(,)d 0f V G S ϕβΣ+∫∫r r r r r (5.2.21)这就是第三边值问题解的积分表示式.右边第一个积分表示区域T 中分布的源0()f r 在r点产生的场的总和.第二个积分则代表边界上的状况对r点场的影响的总和.两项积分中的格林函数相同.这说明泊松方程的格林函数是点源在一定的边界条件下所产生的对于拉普拉斯方程0()0f ≡r 第一边值问题的解为0000(,)()()]d G u S ϕΣ∂=−∂∫∫ r r r r n (5.2.22)第三边值问题的解为1()()(,)d u G S ϕβΣ=∫∫ r r r r (5.2.23)5.3 无界空间的格林函数基本解无界区域这种情形公式(5.2.10)中的面积分应为零,故有000()(,)()d T u G f V =∫∫∫r r r r (5.3.1)选取()u r 和0(,)G r r 分别满足下列方程()()u f Δ=−r r (5.3.2)00(,)(-)G δΔ=−r r r r (5.3.3)5.3.1 三维球对称对于三维球对称情形,我们选取00=r 对(5.3.3)式两边在球内积分)d V(5.3.4)T∫∫∫(5.3.5)5.1.1)得到2(,0)d (,0)d sin d d S S G G V G r r θθϕ∂⋅∇=∇⋅=∂∫∫∫∫ r r S (5.3.6)故有2sin d d (,0)d 1S T G r G V r θθϕ∂=Δ=−∂∫∫∫∫∫ r 使上式恒成立,有2(,0)4π1G r r∂=−∂r 14πcr=+0G →因此0c =,,故得到1(,0)4πG r=r对于三维无界球对称情形的格林函数可以选取为001(,)4π||G =−r r r r (5.3.7)代入(5.3.1)得到三维无界区域问题的解为0(5.3.8)上式正是我们所熟知的静电场的电位表达式5.3.2 二维轴对称情形用单位长的圆柱体来代替球.积分在单位长的圆柱体内进行,即(,0)d ()d TTG V VδΔ=−∫∫∫∫∫∫r r ()d 1V δ=∫∫∫r ,0)d (,0)d SV G =∇∫∫i r SG只是垂直于轴,且向外的分量,所以上式在圆柱体上、下底的面积分为零,只剩下沿侧面的积分,即d d ()d 1T Gr z V r ϕδ=−=−∫∫∫r选取的圆柱的高度为单位长,则很容易得到下面的结果12πG r r∂=−∂11(,0)ln 2πG c r =+r 令积分常数为0,得到11(,0)ln 2πG r=r 0011(,)ln 2π||G =−r r r r (5.3.9))代入式(5.3.1)得到二维无界区域的解为000011()()ln d 2π|S u f S |=−∫∫r r r r。
legendre 多项式
![legendre 多项式](https://img.taocdn.com/s3/m/e7d621c289eb172dec63b707.png)
m = 0,1, 2,L l l = 0,1, 2,L
复数形式的球函数:
Yl
m
(θ , ϕ ) = Pl ( cos θ ) e
m
imϕ
m = 0, ±1, ±2,L ± l , l = 0,1, 2,L
8
独立的 l 阶球函数有2l+1个。
球函数的正交关系与模
Clm = ∫
2π
0
∫
π
0
* Ylm (θ , ϕ ) f (θ , ϕ ) sin θ dθ dϕ .
10
第十一章 柱函数
11
在柱坐标系下对 Laplace 方程和 Helmholtz 方 程分离变量得到 Bessel 方程:
d R dR 2 2 x +x ± ( x m m ) R = 0. x = µ ρ 2 dx dx
由 Legendre 多项式的微分表示立即可得:
Pl
m
( x)
(1 − x ) =
2 l!
l
m 2 2
l d 2 x − 1) . l +m ( dx
l +m
这一关系又称作 Rodrigues 公式 公式。
根据 Cauchy 公式,连带 Legendre 函数还可以表 示为如下积分形式 积分形式 积分形式:
作变换: Θ ( x ) = (1 − x 作变换:
m 2 2
) y ( x)
代入方程并整理可得: 代入方程并整理可得:
1 − x 2 ) y ''− 2 ( m + 1) xy '+ l (l + 1) − m ( m + 1) y = 0 (
勒让德多项式
![勒让德多项式](https://img.taocdn.com/s3/m/572691efaa00b52acec7ca14.png)
六、结论与心得
在调试的过程中,一开始我发现原来的程序进行到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]);//输出
n阶legendre 多项式的根
![n阶legendre 多项式的根](https://img.taocdn.com/s3/m/64cf6b4300f69e3143323968011ca300a6c3f6aa.png)
n阶legendre 多项式的根
Legendre多项式是数学中的一种特殊函数,它在物理学、工程学、计算机科学等领域都有广泛的应用。
其中,n阶Legendre多项式的根是一个重要的概念,它在数学和物理学中都有着重要的作用。
我们来了解一下Legendre多项式。
Legendre多项式是一类正交多项式,它的定义如下:
Pn(x) = (1/2^n * n!) * d^n/dx^n [(x^2 - 1)^n]
其中,n为多项式的阶数,x为自变量。
Legendre多项式的根是指多项式Pn(x)的n个实数解,它们分别称为第1个、第2个、第3个……第n个Legendre多项式的根。
Legendre多项式的根在数学和物理学中都有着广泛的应用。
在数学中,它们被用于求解微分方程、计算积分等问题。
在物理学中,它们被用于描述电磁场、量子力学等现象。
Legendre多项式的根具有一些特殊的性质。
首先,它们是实数,且在区间[-1,1]内均匀分布。
其次,它们具有对称性,即第k个根和第n-k+1个根互为相反数。
此外,它们还满足一些重要的关系式,如下所示:
Pn(xk) = 0
Pn(xk) / Pn'(xk) = (-1)^k / (n+1)
其中,xk为第k个根,Pn(x)为n阶Legendre多项式,Pn'(x)为Pn(x)的导数。
n阶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