第七章1.在00x的领域上求解常微分方程)(02''是常数wywy.解:kkxkxxaxy2420)()!2(1)1()(!41)(!211)(12531)()!12(1)1()(!51)(!31)(kkxkxxxaxaxasincos10.2.在00x的领域上求解雅可比方程0)1()2()1('''2yyx.解:00x是方程的常点.设0nnnxay(1)则0)1()1(kkkxa(2)11')()(nnnxnay01)1()(kkkxak(3)1')2()2(nnnnxayx1)2(kkkxka(4)0222'')2)(1()1(kkknnnxakkxanny(5)22''2)1()1(nkknnnxkakxnanyx(6)把(2)(6)代入雅可脾方程可得系数递推公式:0)1()2()1()1)(()2)(1(12kkkakkkakakk122kkakakakkkkk)2)(1()1()2()1(,即122kkakakakkkk)2)(1()1)((可以写出前面几个系数,但难以写出一般的系数公式。3.在00x的领域上求解0'''yxyxy。解:这里00x是方程的正则奇点。令1)(kkkxaxy,代入方程,把各幂次项集合如下表。1sxsxksx''xysass)1(1)1(ssas1))(1(ksaksks'xyssaksaks)(ysaksa由最低幂项系数为零,得0)1(sass,0,1,021ssas(1)取11ss,这时系数递推公式是:12)2)(1(kkakkka(2)当0k时,,0,0,02112aaaa推知),1(,0kak从而得到方程的一个特解:xaxy11)((3)由于10121ss为整数,所以对于判定方程较小根2s对应的另一个特解的形式是:02ln)(kkkxbxAxxy且11'2ln)(kkkxkbxAAxy22''2)1()(kkkxbkkxAxy代入原方程,集合如下表.0xxxlnxkx''2xyA22b1)1(kkbk'2xyA1bAkkb2y0bA1bkb令上表各幂次项系数为零,有0bA0AA,知A为任意02112bbAb,知22Ab,1b为任意)2(,)1(11kbkkkbkk推之有:AAbb!31!212132132123Abb!4!3!243234Akkkbk!)!1()!2(2212!)!1()!2(1ln!)!1()!2(ln)(kkkkxkkkxxxAAxkkkxbAxAxxy(4)故21021!)!1()!2(1ln)()()(kkxkkkxxaxaxyxyxy。欧勒型常微分方程在00x的邻域上求解)(0222为常数mymyxyx解:把方程写成标准形式,系数xxp/1)(,22/)(xmxq,则00x为)(xp的一阶极点,为)(xq的二阶极点。因此0x为方程的正则奇点.把kskssksksskskssxaksksxsasxassxyxxaksxasxsaxyxxamxamxamxym)1)(()1()1()()()1()()(11021102112022代入常微分方程,然后相加,令各个幂次合并后的系数分别为零,得到一系列方程为222002221122222[(1)]0[]0[(1)(1)]0[(1)]0[(2)(1)(2)]0[(2)]0sssmasmasssmasmasssmasma考虑到第一个系数00a,则第一个方程为判定方程022ms由此得到第一项的幂次为ms1和ms2先取ms1,那么第二个方程成为00])1[(1122aamm,第三个方程成为00])2[(2222aamm,同理得到)1(0kak。因此,常微分方程的一个特解为mxxy)(1再取ms2,同理得)1(0kak。那么,常微分方程的另一个特解为mxxy)(2因此常微分方程的通解为mmDxCxxy)(。9.在00x的邻域上求解阶贝塞尔方程0)(222yxyxyx(为常数)解:把该方程写成标准形式,则系数xxp/1)(和22/1)(xxq在00x点分别为一阶奇点和二阶奇点,因此00x点为贝塞尔方程的正则奇点.把ksksskskssksksskskssxaksksxsasxassxyxxaksxasxsaxyxxaxaxaxyxaxaxaxyx)1)(()1()1()()()1()()()(11021102112022221202代入常微分方程,并合并各项,令各个幂次合并后的系数分别为零,得到一系列方程为0])[(0])2[(0])1[(0][2220222122022kkaaksaasasas考虑到第一个系数00a,则第一个方程为判定方程022s由此得到第一项的幂次为1s和2s第二个方程00])1[(1122aa。利用后面的各式进行系数的递推,得到递推公式为0])[(222kkaaks也即2222))((1)(1kkkaksksaksa先取1s,递推公式成为2)2(1kkakka,因此20023142042022111(22)21!(1)210(23)311(24)42!(1)(2)211(1)!(1)(2)()20kkkkaaaaaaaaaakka因此贝塞尔方程的一个特解为2410211()11!(1)22!(1)(2)21(1)!(1)(2)()2kkxxyxaxxkk该级数的收敛半径为)2(lim|/|lim2kkaaRkkkk因此只要x有限,级数就是收敛的。通常取)1(210a式中)(x为函数,若为整数,那么!)1(。此时把这个解称为阶贝塞尔函数,记为)(xJ022)1(!1)1()(kkkxkkxJ再取2s,递推公式成为2)2(1kkbkkb因此20023142042022111(22)21!(1)210(32)311(42)42!(1)(2)211(1)!(1)(2)()20kkkkbbbbbbbbbbkkb因此,贝塞尔方程的另一个特解为kkxkkxxxbxy2422022)()2)(1(!1)1(2)2)(1(!2122)1(!111)(该级数的收敛半径为)2(lim|/|lim2kkbbRkkkk因此只要x有限,级数就是收敛的。通常取)1(210b把这个解称为阶贝塞尔函数,记为)(xJ022)1(!1)1()(kknkxkkxJ因此,阶贝塞尔方程的通解就是两个特解的线性叠加.)()()(21xJCxJCxy式中1C和2C为常数。若取ctg1C、csc2C,并代入通解中可以得到一个特解,该特解可以作为阶贝塞尔方程的第二个线性独立的特解,称为阶诺依曼函数,也即sin)(cos)()(xJxJxN因此阶贝塞尔方程的通解还可以写成)()()(21xNCxJCxy10.在00x的邻域上求解21阶贝塞尔方程0])([22122yxyxyx解:由判定方程解出两个根为211s和212s。对应于较大的根211s,有一个特解为贝塞尔函数为xxxkxxkxxkkkkxxkkkxkkkxkkxJxykkkkkkkkkkkkkkkkkkksin2)()!12(1)1(2)()!12(1)1(2)(135)12)(12(24)22(2)1(2)(21135)12)(12(!)1(2)()1)((!)1(2)(!1)1()()(0120202020221212121022312121212121由于判定方程两根之差121ss为整数,第二个特解的形式为2/12/12lnsin2ln)()(21kkkkkkxaxxxAxaxxAJxy该特解代入贝塞尔方程中,可以证明0A令各个幂次的项分别等于零,可以得到xxaxycos1)(212若令/221a,则xxxJcos2)(21因此方程的通解为)()()(212121xJCxJCxy由此可以推广到半奇数)(21l阶贝塞尔方程的求解。判定方程为211ls和)(211ls,两根之差为1221lss为整数。对应于211ls的一个特解为21l阶的贝塞尔函数,也即022121212)1(!)1()(kklklxklkxJ第二个线性独立的特解可以用)2/1(2ln)()(21lkkklxbxxAJxy代入方程中,可以证明0A,因此第二个特解仍可以用0221)(21212)1(!)1()(kklklxklkxJ因此)(21l阶贝塞尔方程的通解为)()()()(212121xJCxJCxyll11.在00x的邻域上求解整数m阶贝塞尔方程0)(222ymxyxyx(m为整数)解:贝塞尔方程对应ms1的一个特解为02022)!(!1)1(2)1(!1)1()(kkmkkkmkmxkmkxkmkxJ而对应ms2的另一个特解为022)1(!1)1()(kkmkmxkmkxJ如果m是整数,只要mk,则1km是负整数,而负整数的函数为无限大,因此这个级数实际上只从mk开始,也即21()(1)!(1)2mkkmkmxJxkmk令mkl,则20201()(1)()!(1)21(1)(1)()!!2(1)()mllmmlmlmllmmxJxlmlxlmlJx因此,第二个特解实际上就是第一个特解.对于