D7_8常系数非齐次

常系数非齐次线性微分方程机动目录上页下页返回结束第八节型)()(xPexfmxxxPexflxcos)([)(型]sin)(~xxPn一、二、第七章)(xfyqypy),(为常数qp二阶常系数线性非齐次微分方程:根据解的结构定理,其通解为Yy*y非齐次方程特解齐次方程通解求特解的方法根据f(x)的特殊形式,的待定形式,代入原方程比较两端表达式以确定待定系数.①—待定系数法机动目录上页下页返回结束)([xQex)()2(xQp])()(2xQqp)(xPemx一、型)()(xPexfmx为实数,)(xPm设特解为,)(*xQeyx其中为待定多项式,)(xQ])()([*xQxQeyx])()(2)([*2xQxQxQeyx代入原方程,得(1)若不是特征方程的根,则取从而得到特解形式为.)(*xQeymx为m次多项式.Q(x)为m次待定系数多项式机动目录上页下页返回结束(2)若是特征方程的单根,为m次多项式,故特解形式为(3)若是特征方程的重根,,02p)(xQ则是m次多项式,故特解形式为xmexQxy)(*2小结对方程①,)2,1,0()(*kexQxyxmk此结论可推广到高阶常系数线性微分方程.)(xQ)(xPm)()(2xQqp即即当是特征方程的k重根时,可设特解机动目录上页下页返回结束例1.的一个特解.解:本题而特征方程为不是特征方程的根.设所求特解为代入方程:比较系数,得31,110bb于是所求特解为0,0机动目录上页下页返回结束例2.的通解.解:本题特征方程为,0652rr其根为对应齐次方程的通解为设非齐次方程特解为xebxbxy210)(*比较系数,得1,2110bb因此特解为.)1(*221xexxy代入方程得xbbxb01022所求通解为.)(2221xexx,2机动目录上页下页返回结束例3.求解定解问题0)0()0()0(123yyyyyy解:本题特征方程为其根为设非齐次方程特解为代入方程得故21322CC故对应齐次方程通解为1CYxeC2xeC23原方程通解为1CyxeC2xeC23由初始条件得,0机动目录上页下页返回结束于是所求解为xeeyxx2141432解得41143321CCC机动目录上页下页返回结束二、型xxPxxPexfnlxsin)(~cos)()(ximexPxf)()()(ximexP)()(第二步求出如下两个方程的特解ximexPyqypy)()(yqypy分析思路:第一步将f(x)转化为第三步利用叠加原理求出原方程的特解第四步分析原方程特解的特点ximexP)()(机动目录上页下页返回结束第一步利用欧拉公式将f(x)变形xexf)(ixPxPnl2)(~2)(xie)(ixPxPnl2)(~2)(xie)(ximexPxf)()()(ximexP)()(ximexP)()(ximexP)()(则令,,maxlnm)(xPl2xixiee)(~xPnieexixi2机动目录上页下页返回结束第二步求如下两方程的特解i是特征方程的k重根(k=0,1),ximkexQxy)(1)())((次多项式为mxQm故ximexPyqypy)(111)()()(等式两边取共轭:ximexPyqypy)(111)(1y这说明为方程③的特解.ximexPyqypy)()(②ximexPyqypy)()(③设则②有特解:机动目录上页下页返回结束第三步求原方程的特解利用第二步的结果,根据叠加原理,原方程有特解:11*yyyxkexximximeQeQ原方程yqypyxxPxxPenlxsin)(~cos)(xkex)sin(cosxixQm)sin(cosxixQmxkexxRmcosxRmsin~mmRR~,其中均为m次多项式.机动目录上页下页返回结束第四步分析的特点yxRxRexyyymmxksin~cos11因11yy*yy所以mmRR~,因此均为m次实多项式.11yyy本质上为实函数,11yy机动目录上页下页返回结束小结:xxPxxPenlxsin)(~cos)(对非齐次方程yqypy),(为常数qpxRxRexymmxksin~cos*则可设特解:其中为特征方程的k重根(k=0,1),i上述结论也可推广到高阶方程的情形.机动目录上页下页返回结束例4.的一个特解.解:本题特征方程,2,0故设特解为不是特征方程的根,代入方程得xxxadxcxcbxa2cos2sin)433(2cos)433(012r,)(xxPl,0)(~xPn比较系数,得9431,da于是求得一个特解13a043cb03c043ad0cb机动目录上页下页返回结束例5.的通解.解:特征方程为,092r其根为对应齐次方程的通解为比较系数,得因此特解为)3sin33cos5(*xxxy代入方程:xaxb3sin63cos6所求通解为为特征方程的单根,)3sin33cos5(xxx因此设非齐次方程特解为机动目录上页下页返回结束例6.解:(1)特征方程有二重根所以设非齐次方程特解为(2)特征方程有根xexyyxsin3)2()4(利用叠加原理,可设非齐次方程特解为)sincos(xkxdx设下列高阶常系数线性非齐次方程的特解形式:机动目录上页下页返回结束例7.求物体的运动规律.解:问题归结为求解无阻尼强迫振动方程tphxktxsindd222•当p≠k时,齐次通解:tkCtkCXcossin21)(sintkAtpbtpaxcossin非齐次特解形式:0,22bpkha因此原方程④之解为第七节例1(P294)中若设物体只受弹性恢复力f,sin的作用ptHF和铅直干扰力xox代入④可得:④机动目录上页下页返回结束当干扰力的角频率p≈固有频率k时,)(sintkAxtppkhsin22自由振动强迫振动!22将很大振幅pkh•当p=k时,)cossin(tkbtkatx非齐次特解形式:代入④可得:khba2,0方程④的解为机动目录上页下页返回结束若要利用共振现象,应使p与k尽量靠近,或使)(sintkAxtktkhcos2随着t的增大,强迫振动的振幅这时产生共振现象.可无限增大,若要避免共振现象,应使p远离固有频率k;p=k.自由振动强迫振动xox对机械来说,共振可能引起破坏作用,如桥梁被破坏,电机机座被破坏等,但对电磁振荡来说,共振可能起有利作用,如收音机的调频放大即是利用共振原理.机动目录上页下页返回结束内容小结xmexPyqypy)(.1为特征方程的k(＝0,1,2)重根,xmkexQxy)(*则设特解为]sin)(~cos)([.2xxPxxPeyqypynlx为特征方程的k(＝0,1)重根,ixkexy*则设特解为]sin)(~cos)([xxRxxRmm3.上述结论也可推广到高阶方程的情形.机动目录上页下页返回结束思考与练习时可设特解为xy*xbxacos)(*yxdxcxbxa2sin)(2cos)(xek2时可设特解为提示:xdcxsin)(1.(填空)设]sin)(~cos)([xxRxxRmm机动目录上页下页返回结束2.求微分方程xeyyy44的通解(其中为实数).解:特征方程,0442rr特征根:221rr对应齐次方程通解:2时,,xeAy令代入原方程得,2)2(1A故原方程通解为2时,,2xexBy令代入原方程得,21B故原方程通解为机动目录上页下页返回结束3.已知二阶常微分方程xecybyay有特解,)1(2xxexey求微分方程的通解.解:将特解代入方程得恒等式xxxxecexbaeaeba)1()2()1(比较系数得01baca201ba0a1b2c故原方程为对应齐次方程通解:xxeCeCY21xxexey原方程通解为xxeCeCy21xex机动目录上页下页返回结束