D7_8常系数非齐次线性微分方程

目录上页下页返回结束),(0为常数qpyqypy,02qrpr特征方程:xrxrCCy21ee21实根xrxCCy1e)(21)sincos(e21xCxCyx特征根通解目录上页下页返回结束常系数非齐次线性微分方程第八节型)(e)(xPxfmxxxPxflxcos)([e)(型]sin)(~xxPn一、二、第七章目录上页下页返回结束)(xfyqypy),(为常数qp二阶常系数线性非齐次微分方程:根据解的结构定理,其通解为Yy*y非齐次方程特解齐次方程通解求特解的方法根据f(x)的特殊形式,的待定形式,代入原方程比较两端表达式以确定待定系数.①—待定系数法目录上页下页返回结束)([exQx)()2(xQp])()(2xQqp)(exPmx一、型)(e)(xPxfmx为实数,)(xPm设特解为,)(e*xQyx其中为待定多项式,)(xQ])()([e*xQxQyx])()(2)([e*2xQxQxQyx代入原方程,得为m次多项式.)(xfyqypy(1)若不是特征方程的根,则取从而得到特解形式为.)(e*xQymxQ(x)为m次待定系数多项式目录上页下页返回结束(2)若是特征方程的单根,为m次多项式,故特解形式为(3)若是特征方程的重根,,02p)(xQ则是m次多项式,故特解形式为xmxQxye)(*2小结对方程①,)2,1,0(e)(*kxQxyxmk)(xQ)(xPm)()(2xQqp即即当是特征方程的k重根时,可设特解目录上页下页返回结束例1.的一个特解.解:本题而特征方程为不是特征方程的根.设所求特解为代入方程:比较系数,得31,110bb于是所求特解为0,0目录上页下页返回结束例2.的通解.解:本题特征方程为,0652rr其根为对应齐次方程的通解为设非齐次方程特解为xbxbxy210e)(*比较系数,得1,2110bb因此特解为.e)1(*221xxxy代入方程得xbbxb01022所求通解为.e)(2221xxx,2目录上页下页返回结束二、型xxPxxPxfnlxsin)(~cos)(e)(xmxPxf)i(e)()(xmxP)i(e)(第二步求出如下两个方程的特解xmxPyqypy)i(e)(yqypy分析思路:第一步将f(x)转化为第三步利用叠加原理求出原方程的特解第四步分析原方程特解的特点xmxP)i(e)(目录上页下页返回结束第一步利用欧拉公式将f(x)变形xxfe)(i2)(~2)(xPxPnlx)i(ei2)(~2)(xPxPnlx)i(exmxPxf)i(e)()(xmxP)i(e)(xmxP)i(e)(xmxP)i(e)(则令,,maxlnm)(xPl2eeiixx)(~xPni2eeiixx目录上页下页返回结束第二步求如下两方程的特解i是特征方程的k重根(k=0,1),xmkxQxy)i(1e)())((次多项式为mxQm故xmxPyqypy)i(111e)()()(等式两边取共轭:xmxPyqypy)i(111e)(1y这说明为方程③的特解.xmxPyqypy)i(e)(②xmxPyqypy)i(e)(③设则②有特解:目录上页下页返回结束第三步求原方程的特解利用第二步的结果,根据叠加原理,原方程有特解:11*yyyxkxexmxmQQiiee原方程yqypyxxPxxPnlxsin)(~cos)(exkxe)sini(cosxxQm)sini(cosxxQmxkxexRmcosxRmsin~mmRR~,其中均为m次多项式.xmxP)i(e)(xmxP)i(e)(目录上页下页返回结束第四步分析的特点yxRxRxyyymmxksin~cose11因11yy*yy所以mmRR~,因此均为m次实多项式.11yyy本质上为实函数,11yy目录上页下页返回结束小结:xxPxxPnlxsin)(~cos)(e对非齐次方程yqypy),(为常数qpxRxRxymmxksin~cose*则可设特解:其中为特征方程的k重根(k=0,1),i上述结论也可推广到高阶方程的情形.目录上页下页返回结束例3.的一个特解.解:本题特征方程,2,0故设特解为不是特征方程的根,代入方程得xxxadxcxcbxa2cos2sin)433(2cos)433(012r,)(xxPl,0)(~xPn比较系数,得9431,da于是求得一个特解13a043cb03c043ad0cb目录上页下页返回结束例4.的通解.解:特征方程为,092r其根为对应齐次方程的通解为比较系数,得因此特解为)3sin33cos5(*xxxy代入方程:xaxb3sin63cos6所求通解为为特征方程的单根,)3sin33cos5(xxx因此设非齐次方程特解为目录上页下页返回结束例5.解:(1)特征方程有二重根所以设非齐次方程特解为(2)特征方程有根xxyyxsin3e)2()4(利用叠加原理,可设非齐次方程特解为)sincos(xkxdx设下列高阶常系数线性非齐次方程的特解形式:目录上页下页返回结束例6.求物体的运动规律.解:问题归结为求解无阻尼强迫振动方程tphxktxsindd222•当p≠k时,齐次通解:tkCtkCXcossin21)(sintkAtpbtpaxcossin非齐次特解形式:0,22bpkha因此原方程④之解为若设物体只受弹性恢复力f,sin的作用pthF和铅直干扰力Oxx代入④可得:④目录上页下页返回结束当干扰力的角频率p≈固有频率k时,)(sintkAxtppkhsin22自由振动强迫振动!22将很大振幅pkh•当p=k时,)cossin(tkbtkatx非齐次特解形式:代入④可得:khba2,0方程④的解为Oxxtphxktxsindd222④目录上页下页返回结束若要利用共振现象,应使p与k尽量靠近,或使)(sintkAxtktkhcos2随着t的增大,强迫振动的振幅这时产生共振现象.可无限增大,若要避免共振现象,应使p远离固有频率k;p=k.自由振动强迫振动对机械来说,共振可能引起破坏作用,如桥梁被破坏,电机机座被破坏等,但对电磁振荡来说,共振可能起有利作用,如收音机的调频放大即是利用共振原理.Oxx目录上页下页返回结束内容小结xmxPyqypye)(.1为特征方程的k(＝0,1,2)重根,xmkxQxye)(*则设特解为]sin)(~cos)([e.2xxPxxPyqypynlx为特征方程的k(＝0,1)重根,ixkxye*则设特解为]sin)(~cos)([xxRxxRmm3.上述结论也可推广到高阶方程的情形.目录上页下页返回结束思考与练习时可设特解为xy*xbxacos)(*yxdxcxbxa2sin)(2cos)(xk2e时可设特解为提示:xdcxsin)(1.(填空)设]sin)(~cos)([xxRxxRmm目录上页下页返回结束2.求微分方程xyyye44的通解(其中为实数).解:特征方程,0442rr特征根:221rr对应齐次方程通解:2时,,exAy令代入原方程得,2)2(1A故原方程通解为2时,,e2xxBy令代入原方程得,21B故原方程通解为目录上页下页返回结束3.已知二阶常微分方程xcybyaye有特解2(1e),exxyx求微分方程的通解.解:将特解代入方程得恒等式xxxxcxbaabaee)1(e)2(e)1(比较系数得01baca201ba0a1b2c故原方程为对应齐次方程通解:xxCCYee21xxxyee原方程通解为xxCCyee21xxe目录上页下页返回结束作业P3651(3),(4),(7),(8);2(2),(4);3习题课2第九节