# 常微分方程

在数学分析中，**常微分方程**（英语：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
$$

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


---

# Agent Instructions: Querying This Documentation

If you need additional information that is not directly available in this page, you can query the documentation dynamically by asking a question.

Perform an HTTP GET request on the current page URL with the `ask` query parameter:

```
GET https://zzy.gitbook.io/mybook/calculus/ordinary_differential_equation.md?ask=<question>
```

The question should be specific, self-contained, and written in natural language.
The response will contain a direct answer to the question and relevant excerpts and sources from the documentation.

Use this mechanism when the answer is not explicitly present in the current page, you need clarification or additional context, or you want to retrieve related documentation sections.