> For the complete documentation index, see [llms.txt](https://zzy.gitbook.io/mybook/llms.txt). Markdown versions of documentation pages are available by appending `.md` to page URLs; this page is available as [Markdown](https://zzy.gitbook.io/mybook/calculus/ordinary_differential_equation.md).

# 常微分方程

在数学分析中，**常微分方程**（英语：ordinary differential equation，简称ODE）是未知函数只含有一个自变量的微分方程。

很多科学问题都可以表示为常微分方程，例如根据牛顿第二运动定律，物体在力的作用下的位移$s$和时间$t$的关系就可以表示为如下常微分方程：

$$
m\frac{d^2s}{dt^2}=f(s)
$$

其中 $m$是物体的质量，$f(s)$是物体所受的力，是位移的函数。所要求解的未知函数是位移 $s$，它只以时间 $t$ 为自变量。

## 可分离方程

### 一阶，变量 $x$ 和 $y$ 均可分离（一般情况, 下面有特殊情况）

$$
P\_{1}(x)Q\_{1}(y)+P\_{2}(x)Q\_{2}(y) {\frac {dy}{dx}}=0 \\
P\_{1}(x)Q\_{1}(y) dx+P\_{2}(x)Q\_{2}(y) dy=0
$$

#### 解法

分离变量（除以$ P*{2}Q*{1}$）。

#### 通解

$$
\int ^{x}{\frac {P\_{1}(\lambda )}{P\_{2}(\lambda )}} d\lambda +\int ^{y}{\frac {Q\_{2}(\lambda )}{Q\_{1}(\lambda )}} d\lambda =C
$$

### 一阶，可化为变量可分离型

$$
\frac {dy}{dx}=F(ax+by+c)
$$

其中$a$，$b$不全为零

#### 解法

令

$$
u=ax+by+c
$$

，则

$$
{\frac{du}{dx}} = a+b{\frac{dy}{dx}}
$$

#### 通解

$$
x=\int ^{ax+by+c}{\frac {d\lambda }{a+bf(\lambda )}}+C
$$

## 一般一阶微分方程

### 一阶，齐次

$$
{\frac {dy}{dx}}=F({\frac {y}{x}})
$$

#### 解法

令 $y=ux$，然后通过分离变量 $u$ 和 $x$ 求解。

#### 通解

$$
\ln(Cx)=\int ^{\frac {y}{x}}{\frac{d\lambda }{F(\lambda )-\lambda }}
$$

### 一阶，可分离变量

$$
yM(xy)+xN(xy) {\frac {dy}{dx}}=0 \\

yM(xy) dx+xN(xy) dy=0
$$

#### 解法

分离变量（除以 $xy$）。

#### 通解

$$
\ln(Cx)=\int ^{xy}{\frac {N(\lambda ) d\lambda }{\lambda \[N(\lambda )-M(\lambda )]}}
$$

如果$N=M$, 解为$xy=C$。

### 恰当微分, 一阶

$$
M(x,y){\frac {dy}{dx}}+N(x,y)=0 \\

M(x,y),dy+N(x,y) dx=0
$$

其中 ${\frac {\partial M}{\partial x}} = {\frac {\partial N}{\partial y}}$

#### 解法

整个积分。

#### 通解

$$
F(x,y)=\int ^{y}M(x,\lambda ) d\lambda +\int ^{x}N(\lambda ,y) d\lambda +Y(y)+X(x)=C
$$

其中 $Y(y)$ 和 $X(x)$ 是积分出来的函数而不是常数，将它们列在这里以使最终函数$F(x,y)$满足初始条件。

### 反常微分, 一阶

$$
M(x,y){\frac {dy}{dx}}+N(x,y)=0 \\

M(x,y) dy+N(x,y) dx=0
$$

其中${\frac {\partial M}{\partial x}}\neq {\frac {\partial N}{\partial y}}\\,!$

#### 解法

积分变量 $μ(x,y)$满足

$$
{\frac {\partial (\mu M)}{\partial x}}={\frac {\partial (\mu N)}{\partial y}}
$$

#### 通解

如果可以得到 $μ(x,y)$:

$$
F(x,y)=\int ^{y}\mu (x,\lambda )M(x,\lambda ) d\lambda +\int ^{x}\mu (\lambda ,y)N(\lambda ,y) d\lambda +Y(y)+X(x)=C
$$

## 一般二阶微分方程

### 二阶, 自治

$$
{\frac {\mathrm {d} ^{2}y}{\mathrm {d} x^{2}}}=F(y)
$$

#### 解法

原方程乘以 ${\frac {2\mathrm {d} y}{\mathrm {d} x}}$ ，代换

$$
2{\frac {\mathrm {d} y}{\mathrm {d} x}}{\frac {\mathrm {d} ^{2}y}{\mathrm {d} x^{2}}}={\frac {\mathrm {d}}{\mathrm {d} x}} ({\frac {\mathrm {d} y}{\mathrm {d} x}})^{2}
$$

，然后两次积分。

#### 通解

$$
x=\pm \int ^{y}{\frac {\mathrm {d} \lambda }{\sqrt {2\int ^{\lambda }F(\epsilon ),\mathrm {d} \epsilon +C\_{1}}}}+C\_{2}
$$

### 二阶可降阶，${\frac{d^2y}{dx^2}=f(x, {\frac{dy}{dx}})}$

#### 解法

令${\frac{dy}{dx} = p(x)}$，$\frac{d^2y}{dx^2}={\frac{dp}{dx}}$，则原方程变为

$$
{\frac{dp}{dx}} = f(x, p)
$$

### 二阶可降阶，${\frac{d^2y}{dx^2}=f(y, {\frac{dy}{dx}})}$

#### 解法

令${\frac{dy}{dx} = p(x)}$，$\frac{d^2y}{dx^2}={\frac{dp}{dx}}={\frac{dp}{dy}} \cdot {\frac{dy}{dx}} = {\frac{dp}{dy}} \cdot p$，则原方程变为

$$
{\frac{dp}{dy}} = f(y, p)
$$

## 线性方程 (最高到$n$阶)

### 一阶线性，非齐次的函数系数

$$
{\frac {dy}{dx}}+P(x)y=Q(x)
$$

#### 解法

积分因子: $e^{\int ^{x}P(\lambda )\\,d\lambda }$。

#### 通解

$$
y=e^{-\int ^{x}P(\lambda ) d\lambda }\[\int ^{x}e^{\int ^{\lambda }P(\epsilon ) d\epsilon }Q(\lambda ) {d\lambda }+C]
$$

### 二阶线性，非齐次的常系数

$$
{\frac {d^{2}y}{dx^{2}}}+b{\frac {dy}{dx}}+cy=r(x)
$$

#### 解法

余函数 $y*{c}$: 设 $y*{c}=\mathrm {e} ^{\alpha x}$，代换并解出 $\alpha$ 中的多项式，求出线性无关函数 $e^{\alpha *{j}x}$。 特解 $y*{p}$：一般运用常数变易法，虽然对于非常容易的 $r(x)$ 可以直观判断。

#### 通解

$$
y=y\_{c}+y\_{p}
$$

如果 $b^{2}>4c$, 则:

$$
y\_{c}=C\_{1}e^{(-b+{\sqrt {b^{2}-4c}}){\frac {x}{2}}}+C\_{2}e^{(b+{\sqrt {b^{2}-4c}}){\frac {x}{2}}}
$$

如果 $b^{2}=4c$, 则:

$$
y\_{c}=(C\_{1}x+C\_{2})e^{-{\frac {bx}{2}}}
$$

如果 $b^{2}<4c$, 则:

$$
y\_{c}=e^{-{\frac {bx}{2}}} \[C\_{1}\sin {({\sqrt {|b^{2}-4c|}}{\frac {x}{2}})}+C\_{2}\cos {({\sqrt {|b^{2}-4c|}}{\frac {x}{2}})}]
$$

### $n$ 阶线性，非齐次常系数

$$
\sum *{j=0}^{n}b*{j}{\frac {d^{j}y}{dx^{j}}}=r(x)
$$

#### 解法

余函数 $&#x79;*{c}$：设 $y*{c}=\mathrm {e} ^{\alpha x}$，代换并解出 $\alpha$ 中的多项式，求出线性无关函数 $e^{\alpha *{j}x}$。 特解 $y*{p}$：一般运用常数变易法，虽然对于非常容易的 $r(x)$ 可以直观判断。

#### 通解

$$
y=y\_{c}+y\_{p}
$$

由于 $\alpha *{j}$ 为 n 阶多项式的解：$\prod* {j=1}^{n} (\alpha -\alpha \_{j})=0$，于是：

对于各不相同的 $\alpha \_{j}$，

$$
y\_{c}=\sum *{j=1}^{n}C*{j}e^{\alpha \_{j}x}
$$

每个根 $\alpha *{j}$ 重复 $k*{j}$ 次，

$$
y\_{c}=\sum *{j=1}^{n} (\sum *{\ell =1}^{k*{j}}C*{\ell }x^{\ell -1})e^{\alpha \_{j}x}
$$

对于一些复数值的 $\alph&#x61;*{j}$，令 $\alpha*{j} = \ch&#x69;*{j} + i\gamma*{j}$，使用欧拉公式，前面结果中的一些项就可以写成

$$
e^{\alpha \_{j}x}=e^{\chi \_{j}x}\cos(\gamma \_{j}x+\phi \_{j})
$$

的形式，其中 $\phi \_{j}$ 为任意常量（相移）。

### 伯努利方程

$$
{\frac {dy}{dx}}+P(x)y=Q(x){\frac {d^{n}y}{dx^{n}}}, n \neq 1, 2
$$

#### 解法

变形为$y^{-n} \cdot {\frac{dy}{dx}} + P(x)y^{1-n}=Q(x)$

令$z = y^{1-n}$，得${\frac{dz}{dx}} = (1-n)y^{-n}{\frac{dy}{dx}}$， 则${\frac{1}{1-n}}{\frac{dz}{dx}}+P(x)z=Q(x)$

#### 通解

按照一阶线性微分方程求解即可。

### 欧拉方程

$$
x^{n}{\frac {d^{n}y}{dx^{n}}}+p\_{1}x^{n-1}{\frac {d^{n-1}y}{dx^{n-1}}}+...+p\_{n-1}x{\frac {dy}{dx}}+p\_{n}y=f(x)
$$

#### 解法

令$x= e^t$或$t=lnx$，则

$$
{\frac{dy}{dx}} = {\frac{dy}{dt}}{\frac{dt}{dx}}= {\frac{1}{x}}{\frac{dy}{dt}} \\{\frac{d^{2}y}{dx^{2}}} = {\frac{1}{x^2}}({\frac{d^2y}{dt^2}} - {\frac{dy}{dt}}) \\
{\frac{d^{3}y}{dx^{3}}} = {\frac{1}{x^3}}({\frac{d^3y}{dt^3}} - 3{\frac{d^2y}{dt^2}} +2 {\frac{dy}{dt}}) \\
...
$$

用$D$表示${\frac{d}{dt}}$，转化为

$$
x^{k}{\frac{d^{k}y}{dx^{k}}}=D(D-1)...(D-k+1)y
$$

转化为常系数的微分方程求解。