数理方程第6讲.pptVIP

第六章 贝塞尔函数 §6.1 贝塞尔方程的导出 考虑固定边界的圆膜振动问题：设有半径为l、边界固定的薄圆盘, 且初始横向位移和速度已知, 求圆膜各点的振动规律.这归结为求解下述定解问题: 或 为了求出这个方程满足条件 由于我们讨论的是圆形薄膜问题，则 (q)应该是以2p为周期的周期函数, 这就决定了 只能等于如下的数: 0,1,2,…,n,…对应于 n=n, 有 0(q)=a0/2 (为常数), n(q)=ancos nq + bnsin nq (n=1,2,…).以 n=n代入(6.10)得 则得 §6.2 贝塞尔方程求解 用x,y来表示自变量和函数值, 则 阶贝塞尔方程为 化简后写成 1o a0(c2- 2)=0; 2o a1[(c+1)2- 2]=0; 3o [(c+k)2- 2]ak+ak-2=0 (k=2,3,…)由1o 得c=? , 代入2o 得a1=0. 现暂取c= , 代入3o 得 4o 由此知(6.14)的一般项为 这样(6.14)中一般项的系数变成 β阶第一类贝塞尔函数: 取c=-β时, 用同样方法可得(6.13)的另一特解 当β不为整数时, 这两个特解Jβ(x)与J-β(x)是线性无关的, 由齐次线性常微分方程的通解的结构定理知道, (6.13)的通解为 y=AJβ(x)+BJ-β(x), (6.19)其中A,B为两个任意常数.当然, 在β不为整数的情况, 方程(6.13)的通解除了可以写成(6.19)式以外还可写成其他的形式, 只要能够找到该方程的另一个与Jβ(x)线性无关的特解, 它与Jβ(x)就可构成(6.13)的通解, 这样的特解是容易找到的. 例如, 在(6.19)中取A=cotβp, B=-cscβp, 则得到(6.13)的一个特解 显然, Yβ(x)与Jβ(x)是线性无关的, 因此, (6.13)的通解可写成 y=AJβ(x)+BYβ(x). (6.21)由(6.20)式所确定的函数Yβ(x)称为第二类贝塞尔函数, 或称牛曼函数. §6.3 当β为整数时贝塞尔方程的通解 这时JN(x)与J-N(x)已不能构成贝塞尔方程的通解了. 为了求出贝塞尔函数的通解, 还要求出一个与JN(x)线性无关的特解.定义第二类贝塞尔函数为 求上式的极限可得 根据这个函数的定义, 它确是贝塞尔方程的一个特解, 而且与Jn(x)是线性无关的(因为当x=0时, Jn(x)为有限值, 而Yn(x)为无穷大).综上所述, 不论n是否为整数, 贝塞尔方程(6.13)的通解都可表示为 y=AJn(x)+BYn(x).其中A,B为任意常数, n为任意实数. §6.4 贝塞尔函数的母函数及递推公式 考虑零阶与一阶贝塞尔函数之间的关系.在(6.17)中令n=0及n=1得: 取出第一个级数的第k+2项求导数, 得 即 以上结果可以推广, 现将Jn(x)乘以xn求导数, 得 将(6.26)和(6.27)两式左端的导数求出来, 并经过化简, 则分别得 第二类贝塞尔函数也具有与第一类贝塞尔函数相同的递推公式 6.5 函数展成贝塞尔函数的级数 6.5.1 贝塞尔函数的零点1o Jn(x)有无穷多个单重实零点, 且这无穷多个零点在x轴上关于原点是对称分布着的, 因而Jn(x)必有无穷多个零点;2o Jn(x)的零点与Jn+1(x)的零点是彼此相间分布的, 即Jn(x)的任意两个相邻零点之间必存在一个且仅有一个Jn+1(x)的零点; 函数f(r)在(0,l)上具有一阶连续导数及分段连续的二阶导数,且f(r)在r=0处有界, 在r=l处等于零, 则f(r)在(0,l)上可展开为绝对且一致收敛的级数 即 §6.6 贝塞尔函数应用举例 例1 设有半径为1的薄均匀圆盘, 边界上温度为零摄氏度, 初始时刻圆盘内温度分布为1-r2, 其中r是圆盘内任一点的极半径, 求圆内温度分布规律.解 采用极坐标系, 定解条件与q无关, 温度只能是r,t的函数, 可归结为求解下列定解问题 由物理意义, 还有条件|u|?, 且当t?+?时, u?0.令 u(r,t)=F(r)T(t),代入方程(6.44)得 方程(6.48)的解为 b=mn(0) (n=1,2,3,?),综合以上结果可得 Fn(r)=J0(mn(0)r), 从而 另外 从而 作业 习题六 第101页第3、5、8、10、12题 即 同理可得 将这两式相减及相加, 分别得到 用(6.28)式可以用较低阶的贝塞尔函数把较高阶的贝塞尔函数表示出来, 如果我们已有零阶和一阶贝塞尔函数, 就可以做到. (6.26) (6.27) 再考虑一个特殊阶的贝塞尔函数的递推公式，即