第十三讲：常微分方程

概念

微分方程的解是一个函数，微分方程可以含有微分，一阶或者高阶导数。未知数是一元函数的称为常微分方程，微分方程的未知数函数最高阶数n称为方程的阶，方程的阶数和方程的通解的未知数个数相等(因为要求出原函数必须求积分n次，每次求积分都会出现一个未知常数C)，，方程的初始条件通常会给定$y(x_0)=A,y^{\prime }(x_0)=A_1......y^n(x_0)=A_n$,每次给出一个就可以解出一个未知常数，对于解的结构，常微分非齐次方程的解=齐次方程的通解+一个特解(这个概念其实在线性代数里也出现过，很好理解，通解部分的运算结果得0，不会影响特解的结果，特解本身就是符合方程的，所以合理),由于括号里的性质易得，齐次部分相同的两个二阶非齐次方程的特解相加，也是两个非齐次部分拼成的新方程的特解。

方程类型

可分离变量型

$\frac{dy}{dx}=f(x)g(y)==>\int \frac{dy}{g(y)}=\int f(x)dx$
分离变量然后两边同时积分即可，最好写出显示解y=f(x)+C(每次积分常数不要漏掉).

可化为分离变量型 (换元法)

• $\frac{dy}{dx}=f(ax+by+c)=f(u)$
  $\frac{du}{dx}=a+b\frac{dy}{dx}$

• $\frac{dy}{dx}=\phi (\frac{y}{x})=\phi (u)$
  $\frac{dy}{dx}=\frac{d(ux)}{dx}=\frac{xdu+udx}{dx}=x\frac{du}{dx}+u$

一阶线性微分方程

• $y^{\prime }+p(x)y=q(x)$
  $y=e^{-\int p(x)dx}[\int e^{\int p(x)dx}q(x)dx+C]$
  推导：
  $y^{\prime }+p(x)y$可以写成$e^{-\int p(x)dx}\ast [y{e}^{\int p(x)dx}{]}^{\prime }$,将未知函数y分离出来就可以得到上面的式子,两边同时积分时，记得加上常数C。
  补充：p(x)=ln|φ(x)|时，求解x时可以不加绝对值，否则都加。

• $y^{\prime }+p(x)y=q(x)y^n$(伯努利方程)
  换元法，两边同除以$y^n$，令$z=y^{1-n}$,$\frac{dz}{dx}=(1-n)\ast y^{-n}\frac{dy}{dx}$，式子变为$\frac{1}{1-n}\frac{dz}{dx}+p(x)z=q(x)$,之后按照一阶常微分方程解出得数即可。

二阶可降阶微分方程

• $y^{\prime \prime }=f(x,y^{\prime })$
  由于没有未知函数y，直接把y’看成未知函数用一阶方程的解法来解即可，解出y‘后再用一阶方程的解法解一次，解两次一阶方程就可以解出来y函数，注意由于积分了两次这个方程会有两个未知数C1,C2。

• $y^{\prime \prime }=f(y,y^{\prime })$

  $y^{\prime }=p,y^{\prime \prime }=\frac{dp}{dy}\ast \frac{dy}{dx}=\frac{dp}{dy}\ast p$，分离变量再积分求出p的方程(记得别漏写常数C1)，然后$p=\frac{dy}{dx}$分离变量再积分求出y的方程，记得此时要添加常数C2.

高阶线性微分方程

二阶常系数齐次

• $y^{\prime \prime }+p{y}^{\prime }+qy=0$(p,q为常数)
  该方程判别式：$△=p^2-4q$,特征值为λ1，λ2。通解为y。
  1. $△>0,y=C_1{e}^{{\lambda }_{1}x}+C_2{e}^{{\lambda }_{2}x}$
  2. $△=0,y=(C_1+C_{2}x)e^{\lambda x}$
  3. $△<0,y=e^{\alpha x}(C_{1}cos\beta x+C_{2}sin\beta x),(\alpha \pm \beta i是复数特征根)$
     注意：C1，C2在给定初始条件(函数在某点的值)的情况下可以解出。

二阶常系数非齐次

齐次方程的通解+非齐次的特解=非齐次通解

$P_n(x),P_m(x)$为x的n次和m次多项式。

• $y^{\prime \prime }+p{y}^{\prime }+qy=e^{\alpha x}{P}_n(x)$
  设特解$y^{\ast }=e^{\alpha x}{Q}_n(x)x^k$

  1. $e^{\alpha x}$照抄
  2. $Q_n(x)$为n次x的多项式
  3. $k=(\alpha =={\lambda }_1)+(\alpha =={\lambda }_2)$
     注意：Qn(x)中的多项式系数要把特解带入原方程去解出来。

• $y^{\prime \prime }+p{y}^{\prime }+qy=e^{\alpha x}[P_n(x)cos\beta x+P_m(x)sin\beta x]$
  设特解$y^{\ast }=e^{\alpha x}[Q_l^1(x)cos\beta x+Q_l^2(x)sin\beta x]x^k$

  1. $e^{\alpha x}$照抄
  2. $l=max(n,m),Q^1{Q}^2$为两个系数不同的L次x的多项式
  3. $k=(\alpha \pm \beta i=={\lambda }_1)||(\alpha \pm \beta i=={\lambda }_2)$
     注意：Q1(x),Q2(x)中的多项式系数要把特解带入原方程去解出来。

• 欧拉方程