数理方程第8讲（副本）.ppt

此文件可在网址下载 第五章 贝塞尔函数 §5.1 贝塞尔方程的引出 设有半径为R的薄圆盘, 侧面绝缘, 边界温度保持为零摄氏度, 初始温度已知, 求圆盘内的瞬时温度分布规律.这归结为求解下述定解问题: 或 为了求出这个方程满足条件 由于u(x,y,t)是单值函数, 所以V(x,y)也必是单值的, 因此Q(q)应该是以2p为周期的周期函数, 这就决定了m只能等于如下的数: 0,12,22,…,n2,…对应于mn=n2, 有 Q0(q)=a0/2 (为常数), Qn(q)=ancos nq + bnsin nq (n=1,2,…).以mn=n2代入(5.10)得 则得 §5.2 贝塞尔方程求解 用x,y来表示自变量和函数值, 则n阶贝塞尔方程为 化简后写成 1o a0(c2-n2)=0; 2o a1[(c+1)2-n2]=0; 3o [(c+k)2-n2]ak+ak-2=0 (k=2,3,…)由1o 得c=?n, 代入2o 得a1=0. 现暂取c=n, 代入3o 得 4o 4o 由此知(5.14)的一般项为 这样(5.14)中一般项的系数变成 n阶第一类贝塞尔函数: 取c=-n时, 用同样方法可得(5.13)的另一特解 当n不为整数时, 这两个特解Jn(x)与J-n(x)是线性无关的, 由齐次线性常微分方程的通解的结构定理知道, (5.13)的通解为 y=AJn(x)+BJ-n(x), (5.19)其中A,B为两个任意常数.当然, 在n不为整数的情况, 方程(5.13)的通解除了可以写成(5.19)式以外还可写成其他的形式, 只要能够找到该方程的另一个与Jn(x)线性有关的特解, 它与Jn(x)就可构成(5.13)的通解, 这样的特解是容易找到的. 例如, 在(5.19)中取A=cot np, B=-csc np, 则得到(5.13)的一个特解 显然, Yn(x)与Jn(x)是线性无关的, 因此, (5.13)的通解可写成 y=AJn(x)+BYn(x). (5.21)由(5.20)式所确定的函数Yn(x)称为第二类贝塞尔函数, 或称牛曼函数. §5.3 当n为整数时贝塞尔方程的通解 这时JN(x)与J-N(x)已不能构成贝塞尔方程的通解了. 为了求出贝塞尔函数的通解, 还要求出一个与JN(x)线性无关的特解.定义第二类贝塞尔函数为 求上式的极限可得 根据这个函数的定义, 它确是贝塞尔方程的一个特解, 而且与Jn(x)是线性无关的(因为当x=0时, Jn(x)为有限值, 而Yn(x)为无穷大).综上所述, 不论n是否为整数, 贝塞尔方程(5.13)的通解都可表示为 y=AJn(x)+BYn(x).其中A,B为任意常数, n为任意实数. §5.4 贝塞尔函数的递推公式 考虑零阶与一阶贝塞尔函数之间的关系.在(5.17)中令n=0及n=1得: 取出第一个级数的第k+2项求导数, 得 即 以上结果可以推广, 现将Jn(x)乘以xn求导数, 得 将(5.26)和(5.27)两式左端的导数求出来, 并经过化简, 则分别得 第二类贝塞尔函数也具有与第一类贝塞尔函数相同的递推公式 5.5 函数展成贝塞尔函数的级数 5.5.1 贝塞尔函数的零点1o Jn(x)有无穷多个单重实零点, 且这无穷多个零点在x轴上关于原点是对称分布着的, 因而Jn(x)必有无穷多个零点;2o Jn(x)的零点与Jn+1(x)的零点是彼此相间分布的, 即Jn(x)的任意两个相邻零点之间必存在一个且仅有一个Jn-1(x)的零点; 5.5.2 贝塞尔函数的正交性 通常把定积分 利用§2.6中关于特征函数系的完备性可知, 任意在[0,R]上具有一阶连续导数及分段连续的二阶导数的函数f(r), 只要它在r=0处有界, 在r=R处等于零, 则它必能展开成如下形式的绝对且一致收敛的级数 即 §5.6 贝塞尔函数应用举例 例1 设有半径为1的薄均匀圆盘, 边界上温度为零摄氏度, 初始时刻圆盘内温度分布为1-r2, 其中r是圆盘内任一点的极半径, 求圆内温度分布规律.解 采用极坐标系, 定解条件与q无关, 温度只能是r,t的函数, 可归结为求解下列定解问题 由物理意义, 还有条件|u|?, 且当t?+?时, u?0.令 u(r,t)=F(r)T(t),代入方程(5.44)得 方程(5.48)的解为 b=mn(0) (n=1,2,3,?),综合以上结果可得 Fn(r)=J0(mn(0)r), 从而 另外 从而 *§5.8 贝塞尔函数的渐近公式当x很大时, 近似有 作业 习题五 第127页开始第6题, 第7题 或 由此得 因为t?+?时, u?0, 所以l0, 令l=b2 则 此时方程(5.47)的通