数学物理方法8863cab8a843a136efd636c74773f6a5第9章_第1页
数学物理方法8863cab8a843a136efd636c74773f6a5第9章_第2页
数学物理方法8863cab8a843a136efd636c74773f6a5第9章_第3页
数学物理方法8863cab8a843a136efd636c74773f6a5第9章_第4页
数学物理方法8863cab8a843a136efd636c74773f6a5第9章_第5页
已阅读5页,还剩101页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

第九章二阶常微分方程级数解法变换法本征值问题,9.2特殊函数常微分方程,9.1正交曲线坐标系,9.3常点邻域的级数解法,9.4正则奇点邻域上的级数解法,9.5施图姆刘维本征值问题,对于圆的Dirichlet问题,其边界条件,若分离变量,则,但若选极坐标,9.1正交曲线坐标系,边界条件分离不出来,则,边界条件能分离出来,若以q1,q2,q3正交坐标,则它与直角坐标的相互关系为,如,(1)、柱坐标,1、正交曲线坐标系,(1)、柱坐标,其中,(2)、球坐标,其中,(1)、柱坐标,2、正交曲线坐标系中的u,直角坐标系中,类似有,上面第一式两边除以2再加上第二式得,在柱坐标系中,直角坐标系中,在柱坐标系中,直角坐标系中,在极坐标系中,(2)、球坐标系中,9.2特殊函数常微分方程,(1)、球坐标系,(一)、Laplace方程u=0,首先试图将此变量变r与和分离,代入,令,化简为两个方程,两边除以R,Y,乘以r2,上边第一式化为,是欧拉型常数方程,令,称为球函数方程,称为球函数方程,令,接着试图将变量和分离,代入,用,除以两边,用,除以两边,令,令,令,令,令,该方程称为连带勒让德方程,如m=0,称为勒让德方程,球坐标系,连带勒让德方程,(2)、柱坐标系,试图将变量变与和z分离,代入,用,除以两边,代入,令,令,用,除以两边,代入,令,令,令,即,即,称为贝塞尔方程,即,称为虚宗量贝塞尔方程,柱坐标系,考虑三维波动方程,(二)、波动方程,首先试图将时间变量t与空间变量r,代入,令,简化为,令,分解为,第一个方程的解为,称为亥姆霍兹方程,考虑三维输运方程,(三)、输运方程,首先试图将时间变量t与空间变量r,代入,令,简化为,令,分解为,第一个方程的解为,为亥姆霍兹方程,(1)、球坐标系,(三)、亥姆霍兹方程,首先试图将此变量变r与和分离,代入,两边除以R,Y,乘以r2,令,化简为两个方程,上边第一式化为,这称为l阶球贝塞尔方程,称为球函数方程,令,若k=0,l阶球贝塞尔方程退化为欧拉型方程,l阶球贝塞尔方程,l阶球贝塞尔方程,为l+1/2阶贝塞尔方程,l阶球贝塞尔方程,连带勒让德方程,(2)、柱坐标系,试图将变量变与和z分离,代入,用,除以两边,代入,令,令,代入,令,令,为m阶贝塞尔方程,令,代入,若,为m阶贝塞尔方程,m阶球贝塞尔方程退化为欧拉型方程,为m阶贝塞尔方程,m阶贝塞尔方程退化为欧拉型方程,9.3常点邻域的级数解法,考虑二阶常微分方程,初始条件为,可以用级数求,更一般,对于复变函数w(z),初始条件为,z为复变函数,z0为选定点,C0,C1为复常数,若p(z)和q(z)在点z0的邻域内解析,z0称为方程的常点,若z0是p(z)和q(z)的奇点,z0称为方程的奇点,对于常点邻域内解析的情形,可用级数解法求解,其中ak有待确定,(1)、勒让德方程的级数解,即,在x0=0的邻域内解析,可以用级数求,勒让德方程,或,令,代入有,即,比较各次幂系数有,从而可得,外推,外推,将以上系数代入得,其中,将以上系数代入得,其中,(2)、解的收敛性,利用达朗贝尔判别法有,说明,收敛,说明,发散,可以证明在,(3)、勒让德多项式,由于,勒让德方程的级数解若能退化为多项式,则发散问题解决,考虑到,若取,则y0(x)只到x2n项,因为,被称为l阶勒让德多项式,若取,则y0(x)只到x2n项,因为,但可取,则y1(x)仍发散,从而,将l限制在0或正整数,使“解在x=1保持有限“,称为勒让德方程的自然边界条件,考虑二阶常微分方程,9.4正则奇点邻域上的级数解法,若z0是p(z)和q(z)的奇点,z0称为方程的奇点,解也以方程为奇点,则在,存在两个线性独立解,或,(一)、正则奇点邻域上的级数解,存在两个线性独立解,或,在,若存在的两个线性独立解为有限个负幂项,称z0为方程的正则奇点,考虑p(z)以z0为不高于一阶极点,q(z)以z0为不高于二阶极点,则,可以证明z0为方程的正则奇点,有,或,或,其中s1和s2是如下判定方程的两个根,取,若s1s2整,则取第一个w2(z)等式,否则取第二个,的证明,代入,的证明,考虑z的最低幂次zs-2,考虑阶贝塞尔方程,(二)、阶贝塞尔方程,p(z)以x0=0为一阶极点,q(z)以x0=0为二阶极点,判定方程,代入,(1)、阶贝塞尔方程,考虑阶贝塞尔方程,代入,先取s1=,代入阶贝塞尔方程,先取s1=,代入阶贝塞尔方程,即,即,比较各次幂系数有,得,其中,而,故只要x有限,级数收敛,(2)、解的收敛性,利用达朗贝尔判别法有,故只要x有限,级数收敛,(3)、贝塞尔函数,有,取,称为贝塞尔函数,表示为阶贝塞尔函数,同理有,用,阶贝塞尔方程的通解为,J和J-称为第一类柱函数,当m时,J和J-线性无关,当=m时,J和J-线性相关,证明如下:,(4)、整数m阶贝塞尔方程,令k-m=l,因为k-m+10,故当=m时,J和J-线性相关,需要寻找另一与Jm无关的解,取,阶贝塞尔方程的一个特解为,称为阶诺伊曼函数,阶诺伊曼函数为第二类柱函数,阶贝塞尔方程的通解可取为,但当=m时,C为欧拉常数,考虑(l+1/2)阶贝塞尔方程,取l=0,代入=1/2,(5)、(l+1/2)阶贝塞尔方程,1/2阶贝塞尔函数,1/2阶贝塞尔函数,s1s2=1第二个特解应为,但可试着用=-1/2代入贝塞尔函数,与,线性无关,故J1/2和J-1/2可作为=1/2贝塞尔方程的线性独立解,=1/2贝塞尔方程的通解为,l+1/2贝塞尔函数为,s1s2=2l+1第二个特解应为,但可试着用=-(l+1/2)代入贝塞尔函数,与,线性无关,=l+1/2贝塞尔方程的通解为,当x0时,,(6)、x=0处的自然边界条件,剩下,若研究区域含x=0,要去掉,称x=0处的具有自然边界条件,考虑阶虚宗量贝塞尔方程,(三)、虚宗量贝塞尔方程,(1)、阶虚宗量贝塞尔方程,令,为阶贝塞尔方程,和,阶诺伊曼函数,令阶虚宗量贝塞尔函数为实数,阶虚宗量贝塞尔方程的一般解为,考虑m阶虚宗量贝塞尔方程,(2)、整数m阶虚宗量贝塞尔方程,令,为m阶贝塞尔方程,m阶虚宗量贝塞尔函数为实数,故要寻找另一个独立解,而,m阶虚宗量贝塞尔函数为实数,9.5施图姆刘维本征值问题,(一)、施图姆刘维本征值问题,考虑形式为,的带参量的二阶常微分方程,只有某些非零方程才有解,称为本征值,对应的非零解称为本征函数,前面介绍的方程都属于施图姆刘维型方程,(i)考虑贝塞尔方程,(i)贝塞尔方程,注意:贝塞尔方程是在柱坐标系中代换得到的,自然边界条件,(ii)勒让德方程,自然边界条件,(iii)连带勒让德方程,自然边界条件,(iv)一般二阶常微分方程,边界条件,连带勒让德方程,对于有自然边界条件的情况,贝塞尔方程,勒让德方程,(二)、施图姆刘维本征值问题的性质,(1)、若k(x),k(x),q(x)连续,或最多以x=a,x=b为一阶极点,则存在无限多本征值,(2)、所有本征值均大于零,对应有非零本征函数,(3)、有带全重(x)的正交关系,本征函数yn和ym满足,证:

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论