北京大学数学物理方法上课件_6二阶线性常微分方程的幂级数解法_.DOCVIP

6二阶线性常微分方程的幂级数解法 d2wdw+p(z)+q(z)w=0dz2dz(1)二阶常微分方程的标准形式是 6.1二阶线性常微分方程的常点和奇点 设p(z)和q(z)为二阶线性常微分方程 dwd2w+p(z)+q(z)w=0dzdz方程的常点和奇点 的系数.如果p(z),q(z)在z0点解析,则z0点称为方程的常点.如果p(z),q(z)中至少有一个在z0点不解析,则z0称为方程的奇点. Example6.1超几何方程 z(1?z) Solution写成dwd2w+[γ?(1+α+β)z]?αβw=0dzdzd2wγ?(1+α+β)zdwαβ+?w=0dzz(1?z)dzz(1?z)其系数是 p(z)=γ?(1+α+β)z z(1?z) αβq(z)=?z(1?z) p(z),q(z)有两个奇点:z=0,1.所以z=0,1是方程的奇点,其它点都是方程的常点 HYPERLINK /retype/zoom/197ab78aa0116c175f0e48a2?pn=1x=0y=77raww=12rawh=11o=png_6_0_0_0_0_0_0_918_1188type=picaimh=11md5sum=e2edd3796ae0e7c4b7d9e38f103c5c3csign=34e5491830zoom=png=0-333jpg=0-306 \t _blank .Example6.2Legendre方程 d2ydy(1?x)2?2x+l(l+1)y=0dxdx2 Solution写成d2y2xdyl(l+1)?+y=0dx1?xdx1?x由系数知:x=±1为方程的奇点 HYPERLINK /retype/zoom/197ab78aa0116c175f0e48a2?pn=1x=0y=128raww=12rawh=12o=png_6_0_0_0_0_0_0_918_1188type=picaimh=12md5sum=e2edd3796ae0e7c4b7d9e38f103c5c3csign=34e5491830zoom=png=0-333jpg=0-306 \t _blank . 6.2方程常点邻域内的解 方程常点邻域内的解 Theorem6.1.如果p(z)和q(z)在圆|z?z0|R内解析,则在此圆内常微分方程初值问题 d2wdw+p(z)+q(z)w=0dz2dz w(z0)=C0,w??(z0)=C1 C0,C1为任意常数,有唯一的解w(z),并且w(z)在这个圆内单值解析. 根据这个定理,可以用幂级数解法求方程常点邻域内的解. 首先将p(z)和q(z)在|z?z0|R内展开成Taylor级数 p(z)= q(z)=∞ k=0∞ k=0ak(z?z0)kbk(z?z0)k(2)(3) 根据定理,又可以把w(z)在相同的邻域|z?z0|R内展开为Taylor级数 w(z)= 由初值条件可得 c0=w(z0)=C0, Example6.3求Legendre方程 2xdyl(l+1)d2y?+y=0dx1?xdx1?x在x=0点邻域内的解,其中l是一个参数. Solutionx=0为常点.因为x=±1为方程奇点,所以p(x),q(x)在|x|1内解析,于是可以在|x|1内找到方程的解析解.可令解为∞ ckxky(x)= k=0∞ k=0ck(z?z0)k(4)c1=w (z0)=C1将幂级数展开式一起代入微分方程,即可用待定系数法求出解的各个系数ck. 下一步若将方程系数p(x),q(x)展开成Taylor级数,则解的过程比较繁.我们将方程仍写成 d2ydy(1?x)?2x+l(l+1)y=0dxdx2 将幂级数解代入,有 (1?x) +l(l+1) 左边第一项可写成 ∞ k=2 ∞ k=02∞ ckk(k?1)xk?2?2x∞ k=0ckkxk?1k=0∞ k=0ckxk=0ckk(k?1)xk?2?∞ k=0ckk(k?1)xk∞ k=0= 代入,整理、和并后,得 ∞ ck+2(k+2)(k+1)x?kckk(k?1)xk{(k+2)(k+1)ck+2?[k(k+1)?l(l+1)]ck}xk=0 k=0 比较x各幂次的系数 (k+2)(k+1)ck+2?[k(k+1)?l(l+1)]ck=0 得到系数之间的递推关系 c(k?l)(k+l+1) k+2=ck 反复利用递推关系 cc2n=0 (2n)!(2n?l?2)(2n?l?4)···(?l) (2n+l?1)(2n+l?3)···(l+1) c2n+1=c1 (2n+1)!(2n?l?1)(2n?l?3)···(?l+1) (2n+l)(2n+l?2)···(l+2) 利用Γ函数性质 Γ(z+1)=zΓ(z) Γ(z+n)=(z+n?1)(z+n?2)···(z