微分方程

时间:2024-04-12 17:07:19

微分方程

微分方程

一、一阶微分方程的求解

1. 可分离变量型

  1. 能写成:
    y=f(x)g(y)y'=f(x)\cdot g(y)
    则:
    dyg(y)=f(x)dyg(y)=f(x)dx\frac{dy}{g(y)}=f(x)\\\int\frac{dy}{g(y)}=\int f(x)dx
  2. 能写成
    y=f(ax+by+c)y'=f(ax+by+c)
    令:
    u=ax+by+cu=ax+by+c
    则:
    u=a+bf(u)dua+bf(u)=dxdua+bf(u)=dxu'=a+bf(u)\\\frac{du}{a+bf(u)}=dx\\\int\frac{du}{a+bf(u)}=\int dx

2. 齐次型

  1. 能写成:
    y=f(yx)y'=f(\frac{y}{x})
    令:
    u=yxy=uxdydx=u+xdudxu=\frac{y}{x}\\y=ux\\\frac{dy}{dx}=u+x\frac{du}{dx}
    则:
    xdudx+u=f(u)duf(u)u=dxxduf(u)u=dxxx\frac{du}{dx}+u=f(u)\\\frac{du}{f(u)-u}=\frac{dx}{x}\\\int\frac{du}{f(u)-u}=\int\frac{dx}{x}
  2. 能写成
    1y=f(xy)\frac{1}{y'}=f(\frac{x}{y})
    令:
    u=xyx=uydxdy=u+ydudyu=\frac{x}{y}\\x=uy\\\frac{dx}{dy}=u+y\frac{du}{dy}
    则:
    ydudy+u=f(u)duf(u)u=dyyduf(u)u=dyyy\frac{du}{dy}+u=f(u)\\\frac{du}{f(u)-u}=\frac{dy}{y}\\\int\frac{du}{f(u)-u}=\int\frac{dy}{y}

3. 一阶线性型

能写成:
y+p(x)y=q(x)y'+p(x)y=q(x)
则:
ep(x)dxy+ep(x)dxp(x)y=ep(x)dxq(x)[ep(x)dxy]=ep(x)dxq(x)ep(x)dxy=ep(x)dxq(x)dx+Cy=ep(x)dx[ep(x)dxq(x)dx+C]e^{\int p(x)dx}\cdot y'+e^{\int p(x)dx}p(x)\cdot y=e^{\int p(x)dx}\cdot q(x)\\\left[e^{\int p(x)dx}\cdot y\right]'=e^{\int p(x)dx}\cdot q(x)\\e^{\int p(x)dx}\cdot y=\int e^{\int p(x)dx}\cdot q(x)dx+C\\y=e^{-\int p(x)dx}\left[\int e^{\int p(x)dx}\cdot q(x)dx+C\right]

二、二阶可降阶微分方程的求解

1. 赶尽杀绝yy

能写成:
y=f(x,y)y''=f(x,y')
yy,令:
y=p,y=pdpdx=f(x,p)y'=p,y''=p'\\\frac{dp}{dx}=f(x,p)
则解为:
p=φ(x,C1)p=\varphi(x,C1)
y=φ(x,C1)y'=\varphi(x,C1),通解为:
y=φ(x,C1)dx+C2y=\int\varphi(x,C1)dx+C2

2. 斩草除根xx

能写成:
y=f(y,y)y''=f(y,y')
xx,令:
y=p,y=dpdx=dpdydydx=dpdyppdpdy=f(y,p)y'=p,y''=\frac{dp}{dx}=\frac{dp}{dy}\cdot\frac{dy}{dx}=\frac{dp}{dy}\cdot p\\p\frac{dp}{dy}=f(y,p)
若求得其解为p=φ(y,C1)p=\varphi(y,C1),则:
p=dydx=φ(y,C1)p=\frac{dy}{dx}=\varphi(y,C1)
分离变量得:
dyφ(y,C1)=dx\frac{dy}{\varphi(y,C1)=dx}
两边积分:
dyφ(y,C1)=x+C2\int\frac{dy}{\varphi(y,C1)}=x+C2

三、高阶常系数线性微分方程的求解

1. 齐次线性微分方程的通解

  1. p24q>0p^2-4q>0,则λ1λ2\lambda_1\ne\lambda_2是特征方程的两个不等实根,则通解为:
    y=C1eλ1x+C2eλ2xy=C_1e^{\lambda_1x}+C_2e^{\lambda_2x}
  2. p24q=0p^2-4q=0,则λ1=λ2\lambda_1=\lambda_2是特征方程的两个相等实根,则通解为:
    y=(C1+C2x)eλxy=\left(C_1+C_2x\right)e^{\lambda x}
  3. p24q<0p^2-4q<0,设α±βi\alpha\pm\beta i是特征方程的一堆共轭复根,则通解为:
    y=eαx(C1cosβx+C2sinβx)y=e^{\alpha x}\left(C_1\cos{\beta x}+C_2\sin{\beta x}\right)

2. 非齐次线性微分方程的特解

  1. 当*项f(x)=Pn(x)eαxf(x)=P_n(x)e^{\alpha x}时,特解要设为:
    y=eαxQn(x)xky^*=e^{\alpha x}Q_n(x)x^k其中:
    {eαxQn(x)xnk={0,αλ1,αλ21,α=λ1α=λ22,α=λ1=λ2\left\{\begin{matrix} e^{\alpha x}照抄 \\ Q_n(x)为x的n次一般多项式 \\ k=\left\{\begin{matrix} 0,&\alpha\ne\lambda_1,\alpha\ne\lambda_2 \\ 1,&\alpha=\lambda_1或\alpha=\lambda_2 \\ 2,&\alpha=\lambda_1=\lambda_2 \end{matrix}\right. \end{matrix}\right.
  2. 当*项f(x)=eαx[Pm(x)cosβx+Pn(x)sinβx]f(x)=e^{\alpha x}\left[P_m(x)\cos{\beta x}+P_n(x)\sin{\beta x}\right]时,特解要设为:
    y=eαx[Ql(1)(x)cosβx+Ql(2)(x)sinβx]xky^*=e^{\alpha x}\left[Q_l^{(1)}(x)\cos{\beta x}+Q_l^{(2)}(x)\sin{\beta x}\right]x^k
    其中:
    {eαxl=max(m,n),Ql(1)(x),Ql(2)(x)xlk={0,α±βi1,α±βi \left\{\begin{matrix} e^{\alpha x}照抄 \\ l=\max(m,n),Q_l^{(1)}(x),Q_l^{(2)}(x)分别为x的两个不同的l次一般多项式 \\ k=\left\{\begin{matrix} 0, & \alpha\pm\beta i不是特征根 \\ 1, & \alpha\pm\beta i是特征根 \end{matrix}\right. \\ \end{matrix}\right.

3. 能写成y+py+qy=f(x)y''+py'+qy=f(x)

  1. λ2+pλ+q=0\lambda^2+p\lambda+q=0,写齐次方程的通解
  2. 设特解yy',代回方程求待定系数,写出特解
  3. 写出通解

4. 能写成y+py+qy=f1(x)+f2(x)y''+py'+qy=f_1(x)+f_2(x)

  1. λ2+pλ+q=0\lambda^2+p\lambda+q=0,写齐次方程的通解
  2. y+py+qy=f1(x)y''+py'+qy=f_1(x),写出特解y1y_1^*
  3. y+py+qy=f2(x)y''+py'+qy=f_2(x),写出特解y2y_2^*
  4. y1+y2y_1^*+y_2^*为特解
  5. 写出通解

5. 能写成x2y+pxy+qy=f(x)x^2y''+pxy'+qy=f(x)

6. y(n)(n3){y^{(n)}}'' (n\ge3)的情形