§4.3高阶微分方程的降阶和幂级数解法一、可降阶的一些方程类型n阶微分方程的一般形式:0),,,,()('nxxxtF1不显含未知函数x,或更一般不显含未知函数及其直到k-1(k1)阶导数的方程是)57.4(0),,,,()()1()(nkkxxxtF阶方程的则可把方程化为若令knyyxk,)()58.4(0),,,,()('knyyytF解得),,,(1knccty积分即),,,(1)(knkcctx为任常数这里nncccctx,,),,,,(11解题步骤:则方程化为令,)(yxk第一步:0),,,,()('knyyytF第二步:求以上方程的通解),,,(1knccty即),,,(1)(knkcctx第三步:对上式求k次积分,即得原方程的通解为任常数这里nncccctx,,),,,,(11)57.4(0),,,,()()1()(nkkxxxtF解令,44ydtxd则方程化为01ytdtdy这是一阶方程,其通解为,cty即有,44ctdtxd对上式积分4次,得原方程的通解为,54233251ctctctctcx例1.014455的通解求方程dtxdtdtxd2不显含自变量t的方程,一般形式:)59.4(,0),,,()('nxxxF,,'作为新的自变量而把作为新的未知函数此时以xxy,ydtdx因为dtdy22dtxddxdydtdx,dxdyy3232dxddxdtdtdtdtd)(dxdyydxdxdyyd)(dtdx,222dxydy2)(dxdyy用数学归纳法易得:来表达可用)(,,,)1()1()(nkdxyddxdyyxkkk将这些表达式代入(4.59)可得:2222(,,,(),)0dydydyFxyyyydxdxdx即有新方程0),,,,()1()1(nndxyddxdyyxG它比原方程降低一阶解题步骤:第一步:原方程化为自变量为新的为新的未知函数并令,,,'xyxy0),,,,()1()1(nndxyddxdyyxG第二步:求以上方程的通解),,,(11nccxy第三步:解方程),,,(11nccxdtdx即得原方程的通解解令,,作为新的自变量并以xydtdx则方程化为02ydxdyxy从而可得,0y及,xydxdy这两方程的全部解是,1xcy例2.0)(222的通解求方程dtdxdtxdx再代回原来变量得到,1xcdtdx所以得原方程的通解为12,ctxce3已知齐线性方程的非零特解,进行降阶1(1)0xx设是二阶齐线性方程22()()0,(4.69)dxdxptqtxdtdt的非零解令1xxy则'''11xxyxy''''''''1112xxyxyxy代入(4.69)得'''''''111111[2()][()()]0xyxptxyxptxqtxy即''''111[2()]0xyxptxy''''111[2()]0xyxptxy引入新的未知函数',zy方程变为'111[2()]0dzxxptxzdt是一阶线性方程,解之得()21,ptdtczex因而()112211[],(4.70)ptdtxxccedtx12,cc这里是任常数.则()21211,ptdtycedtcx因此(4.69)的通解为1x因它与之比不等于常数,12,xx故线性无关120,cc令=1得(4.69)的一个解:()21211,ptdtxxedtx()112211[],(4.70)ptdtxxccedtx12,cc这里是任常数.22()()0,(4.69)dxdxptqtxdtdt解题步骤:第一步:1xxy令方程变为''''1112[()]0xyxptxy第二步:'zy令方程变为'1112[()]0dzxxptxzdt解之得()21,ptdtczex即()112211[],(4.70)ptdtxxccedtx第三步:1210,ccx令=1得与线性无关一个解:()21211,ptdtxxedtx第四步:(4.69)的通解为()112211[],(4.70)ptdtxxccedtx12,cc这里是任常数.注一般求(4.69)的解直接用公式(4.70)解这里12sin(),tptxtt由(4.70)得22122[]sindtttccedtt例322sin20.tdxdxxxtdttdt已知是方程的解,试求方程的通解2122sin[sinttcctt12sin[tcct121[sincos]ctctt12,cc这里是任常数.sintxt21]dttcot]t(2)一般已知齐线性方程111()()0(4.2)nnnnndxdxatatxdtdt2,,,,kxxx1的k个线性无关的解0,1,2,,,ixik显然,kxxy令则'''kkxxyxy''''''''2kkkxxyxyxy()()'(1)''(2)()(1)2nnnnnkkkknnxxynxyxyxy代入(4.2)得()'(1)1[()]nnkkkxynxatxy()(1)1[()]0nnkknkxatxaxy(4.2)kx因为的解,y故的系数恒为零,y即化为不含的方程,',zyxk令则在0的区间上方程变为(1)(1)11()()0,(4.67)nnnzbtzbtz'(),1,2,,1(4.67)1iikxzikkx且是的个线性无关的解事实上21,,,(4.2),kxxx1由为的解及以上变换知'()kkxzxxzdtx或21,,,(4.67),kzzz1因此是的解若12211kkzzz10则121kkkkkxxxxxx12k-1即12211kkkkxxxx102,,,,kxxx1由线性无关知121,,kk全为021,,,,kzz1故z线性无关因此,对(4.67)仿以上做法,1,kzudt令z-2un则可把方程化为关于的阶线性方程(2)(3)12()()0,(4.68)nnnuctuctu,k且可(4.68)的-2个线性无关的解'1(),1,2,,2iikzuikz-nk以上做法一直下去,可降低阶.二、二阶线性方程的幂级数解法对二阶变系数齐线性方程22()()0(4.72)dydypxqxydxdx其求解问题,归结为寻求它的一个非零解.下面考虑该方程及初始条件'(1)0000(),()yxyyxy的情况用级数表示解?00)x(不失一般性,可设定理10,(4.72)pxqxxxR若方程(4.72)中系数()和()都可展成的幂级数,且收敛区间为则方程有形如0,(4.73)nnnyax=.xR的特解,也以为级数的收敛区间定理112()()()(),(4.72)pxqxxpxxqxxR若方程(4.72)中系数和都具有这样的性质,即和均可展成x的幂级数,且收敛区间为则方程有形如00,(4.75)nnnnnnyxaxax00,.axR的特解,这里是一个待定常数,级数(4.75)也以为收敛区间例4''240(0)0,.xyyyy求方程y满足初始条件(0)=1的解解设级数1nnaaxax0y=为方程的解,(1,2,)aii这里是一个待定常数,由初始条件得:10,1;aa0因而22nnyxaxax=122nnyaxnax=1223232(1)nnyaaxnnax=将它代入方程,合并同类项,并令各项系数等于零,得220a332240a42243440aaa22(1)2(2)40nnnnnanaa即20,a31,a40,a,22,1nnaan因而51,2!a60,a71,3!a80,a91,4!a也即211,!kak20,ka;k对一切正整数成立故方程的解为52132!!kxxxxky=422(1)2!!kxxxxk2xxe例522222()0(4.74).dydyxxxnydxdxn求解n阶Bessel方程这里为非负常数解将方程改写为2222210dydyxnydxxdxx易见,它满足定理11条件,且222()1,()xpxxqxxn,11xx按展成的幂级数收敛区间为由定理方程有形如0,(4.75)knkyax00,a的解,这里是一个待定常数,将(4.75)代入(4.74)中,得220kkxaxk(+k)(+k-1)10kkxaxk(+k)220()0kkxnaxkx比较的同次幂系数得220()0an221[(1)]0an222[()]0,2,3,kkaknak(4.76)a0因为0,220,n则有,n从而,n为确定起见暂令由(4.76)得10,a2,2,3,(2)kkaakknk即2121,(21)(221)kkaaknk222,2(22)kkaaknk1,2,k从而可得210,1,2,kak022(1),2!(1)(2)()kkkaaknnnk1,2,k0,nBessel因此在时得到方程的一个解201021(1),(4.77)2!(1)()nkknkkayaxxknnk0a若将任常数取为012(1)nan-10()(1)()xppexdxppp这里,注意到时.因此(4.77)变为2101(1)()(),(4.77)!(1)2kknnkxyJxknk()(4.74),.nJxBesselnBessel是由方程定义的特殊函数称为阶函数,n当时完全类可得210,1,2,kak022(1),2!(1)(2)()kkkaaknnnk1,2,k若取012(1)nan则可得(4.74)的另一个特解2201(1)()(),(4.78)!(1)2kknnkxyJxknk()(4.74),.nJxBesselnBessel是由方程定义的特殊函数称为-阶函数由达朗贝尔判别法,对任x值(4.77),(4.78)收敛.-,()()nnJxJx因此,当不等于非负整数时和都是(4.74)的解,且线性无关.因而(4.74)的通解为12()(),nnycJxcJx12,.cc这里为任常数()nnakn2k当等于正整数,而,不能从(4.76)确定因此,不能象上面一样求得通解;()nJx但可用一3介绍的降阶法,求出与线性无关的解,因此,(4.74)的通解为11221()[],()dxxnnyJxccedxJx1221()[],()nnJxccdxxJx12,.cc为任常数例62'29(4)0.25xyxyxy求方程的通解解2tx引入新变量我们有2dydydtdydxdtdxdt2222(2)4dyddydtdydxdtdtdxdt代入方程得22229()025dydytttydtdt3,5nBessel这是的方程故方程的通解为132355()(),ycJtcJt代回原来的变量得原方程的通解为132355(2)(2),ycJxcJx12,.cc为任常数作业P1652,5,P1658,10