简单微分方程的求解

一、一阶微分方程1.线性齐次方程'y()0pxy①分离变量法求解②两边同时乘以()pxdxe，积分因子法通解：()pxdxyCe2.线性非齐次方程'y()()pxygx①常数变易法②两边同时乘以()pxdxe，积分因子法通解：()()(())pxdxpxdxyeCgxedx线性微分方程的解有一些很好的性质，例如（1）齐次方程的解或者恒等于零，或者恒不等于零（2）齐次方程任何解的线性组合仍是它的解（3）齐次方程的任一解与非齐次方程任一解之和仍是非齐次方程的解（4）非齐次方程任意两解之差必是对应齐次方程的解（5）非齐次方程的任一解与对应齐次方程的通解之和是非齐次方程的通解。3.Bernoulli方程'()()ypxygxy（1）0时，该方程为线性非齐次方程（2）1时，该方程为线性齐次方程（3）0,1时，作变量替换1zy，该方程转化为(1)()(1)()dzpxzgxdx，这是关于未知函数z的一阶线性方程4.Riccati方程2()()()dypxyqxyfxdxRiccati方程在一般情况下无法用初等积分求出解，只是对一些特殊情况或者事先知道了它的一个特解，才能求出其通解。（1）当()px、()qx、()fx都是常数时，是可分离变量方程，用分离变量法求解。（2）当()0px时，是线性方程。（3）当()0fx时，是Bernoulli方程。当()fxr，设已有一特解1()yx命1()()()zxyxyx，代得211(2)dzdydypzpyqzdxdxdx这是一个关于z的Bernoulli方程。（4）当Riccati方程的形式为22dylbayydxxx，可利用变量替换zxy，将方程化为可分离变量方程2(1)dzxazlzbdx当Riccati方程的一个特解()yx已知时，我们利用变换()yzx，代入方程后可得：22()()(2()())()(())()dzdxpxzzxxqxzxfxdxdx由于()yx是方程的解，从上式消去相关的项后得：2(2()()())()dzpxxqxzpxzdx，这是一个Bernoulli方程。（5）当Riccati方程的形式为2mdyaybxdx，其中a、b、m都是常数，且设0a，又设0x和0y，则当440,2,,,(1,2,)2121kkmkkk时，方程可通过适当的变换化为变量可分离方程。5.可分离变量方程'()()yfxgy()()dyfxdxgy，通解为()()dyfxdxCgy6.齐次方程()dyygdxx作变量替换yzx，则dydzzxdxdx，即()dzgzzdxx通解为ln()dzxCgzz。7.全微分方程与积分因子设(,)uFxy是一个连续可微的二元函数，则它的全微分为：(,)(,)(,)FxyFxydudFxydxdyxy若有函数使得：(,)(,)(,)dFxyMxydxNxydy则称(,)(,)0MxydxNxydy为全微分方程，此时，微分方程的解就是(,)FxyC微分方程的成立条件：设函数(,)Mxy和(,)Nxy在一个矩形区域R中连续且有连续的一阶偏导数，则(,)(,)0MxydxNxydy是全微分方程的充要条件是(,)(,)MxyNxyyx微分方程的解为000(,)(,)(,)xyxyFxyMsydsNxsds（线积分法）此时还可应用偏积分法与凑微分法如：22(cossin)(1)0xxxydxyxdy重新分组整理为22cossin()0xxdxxydxyxdyydy如果有函数(,)xy，使得方程(,)(,)(,)(,)0xyMxydxxyNxydy是全微分方程（恰当方程），则(,)xy称为方程(,)(,)0MxydxNxydy的一个积分因子积分因子一般很难求解，但有如下情况可求：（1）微分方程(,)(,)0MxydxNxydy有一个依赖于x的积分因子的充要条件是(,)(,)()/(,)MxyNxyNxyyx仅于x有关，则积分因子可求：(,)(,)()/(,)()MxyNxyNxydxyxxe（2）微分方程(,)(,)0MxydxNxydy有一个依赖于y的积分因子的充要条件是(,)(,)()/(,)NxyMxyMxyxy仅于y有关，则积分因子可求：(,)(,)()/(,)()MxyNxyMxydyyxxe积分因子是求解微分方程的一个极为重要的办法，绝大多数方程的求解都可以通过寻找到一个合适的积分因子来解决。但求一个微分方程的积分因子十分困难，需要灵活运用各种微分法的技巧与经验。例如，当一个微分方程中出现xdyydx的项时，函数1xy、21x、21y和221xy都有可能成为其积分因子，可以根据方程中其他项进行适当的选择。下面的几个方程和对应的积分因子分别为：10,xyxdyydxxydx2210,xxdyydxxdx2210,yxdyydxydx22221)(0,xyxdyydxxydx另外，若有微分方程：1122((,)(,))((,)(,))0MxydxNxydyMxydxNxydy其中第一组和第二组各有积分因子1(,)xy和2(,)xy，使得1111(,)((,)(,))(,)xyMxydxNxydydFxy2222(,)((,)(,))(,)xyMxydxNxydydFxy由于对任意可微函数1()uG和2()uG，111(,)((,))xyFxyG是第一组的积分因子，222(,)((,))xyFxyG是第二组的积分因子。如果能选取的1()uG和2()uG，使得：111222(,)((,))(,)((,))xyFxyxyFxyGG则111(,)((,))(,)xyFxyxyG就是该微分方程的一个积分因子。8.变量替换法（1）形如()faxbycdydx的方程对于这种类型的方程，引入新变量zaxbyc则abdydzdxdx，于是原方程就化为()abfzdzdx这是一个变可分离方程，它的通解为()dzxCabfz此时注意：形如111222axbycaxbycdydx的微分方程，若上下二元一次方程组有解，则利用齐次解法依靠解的坐标点化简此式，若无解则利用变量替换法求解。（2）形如()()0yfxydxxgxydy的方程对于这类方程，引入新变量zxy，则zxy，2xdzzdxxdy原方程可以化为(()())()0zfzgzdxgzdzx，这是一个可分离变量方程。（3）用变量替换法求解微分方程是十分灵活的，依赖于方程的形式和求导的经验，在学习过程中要多积累。9.一阶隐式微分方程解法10.近似解法（1）逐次迭代法逐次迭代法是利用证明初始值问题解的存在唯一性时所构造的Picard迭代序列的前若干项来近似初始值问题的解，其近似序列为：00()yxy0100()(,())xxyxyftytdt001()(,())xnnxyxyftytdt当初始值问题满足解的存在唯一性定理的条件时，上面的迭代序列在一个区间一致收敛到它的解。故当n较大时，()nyx就是初始值问题解的一个较好的近似。（2）Taylor级数法设初始值问题的解可以在0x的邻域内展开为收敛幂级数：00()()nnnyxaxx由Taylor级数理论知，na是由yx的n阶导数确定的，即：()000()()()!nnnyxyxxxn于是，级数形式的解实际上就是要求出()yx在0x点的各阶导数值。如果我们能计算出()yx前面一些导数值()0()nyx时，就可以利用函数()000()()()!nNnNnyxyxxxn来近似初始值问题的解()yx。由复合链导法则和方程初始值得：(0)000()()yxyxy0(1)000()(,)xxdyyxfxydx002(2)''00000002()(,())(,)(,)(,)xyxxxxdydyxfxyxfxyfxyfxydxdx0(3)''''''0000000()((,)(,)(,)(,)2(,)(,)xyxxxyxxdyxfxyfxyfxyfxyfxyfxydx''2'''2000000000000(,)(,)(,)(,)((,))(,)yyyxyfxyfxyfxyfxyfxyfxy根据需要，当函数(,)fxy已知时，我们可以计算出解()yx在0x点直到N阶导数值从而得出()yx的近似表达式。从另一观点看，近似解()000()()()!nNnNnyxyxxxn实质上就是要确定()yx的级数表达式00()()nnnyxaxx中的前面若干个系数na，我们可以将00()()nnnyxaxx代入初始值问题00(,),()dyfxyyxydx，比较所得等式两边0()xx的同次幂的系数即可，这种方法即为待定系数法。（3）数值方法如Runge-Kutta微分方程数值求解方法。