常微分方程
常微分方程(ODEs)是涉及一个或多个自变量以及其导数的方程。术语“常微分”用于与“偏微分”相对,后者用于指偏微分方程(PDEs)。ODEs 在多种科学领域中都很重要,因为它们通常描述了系统随时间的变化。这些方程的解可以揭示出所研究系统的行为和特性。
基本概念和定义
常微分方程是包含函数一个或多个导数的方程。主要目的是确定未知函数,通常表示为y(x)
。对于 ODEs,仅包括对一个变量的导数,通常是时间或空间。ODEs 通常可以表示为:
F(x, y, y', y'', ..., y (n) ) = 0
其中y'
、y''
、...、y (n)
表示y
对自变量x
的导数。
微分方程的阶数
微分方程的阶数是方程中最高阶导数的阶数。例如:
y' + y = 0
是一阶微分方程。y'' + y' - y = x
是二阶微分方程。
线性与非线性 ODEs
如果一个 ODE 可以表示为未知函数及其导数的线性多项式,则称为线性 ODE。线性 ODE 看起来像这样:
a n (x)y (n) + a n-1 (x)y (n-1) + ... + a 1 (x)y' + a 0 (x)y = g(x)
其中a i (x)
和g(x)
仅是x
的函数。如果 ODE 不能以这种形式表示,则是非线性的。
求解常微分方程
有多种方法可以求解常微分方程,所选方法取决于 ODE 的类型和阶数。让我们讨论一些基本方法:
可分离方程
如果一个方程可以重新排列成一种形式,使变量可以在方程的两边分离,则称为可分离方程。例如:
dy/dx = g(y)h(x)
可重新排列为:
(1/g(y)) dy = h(x) dx
现在可以分别对两边进行积分。
示例:
dy/dx = yx
这可以重写为:
(1/y) dy = x dx
对两边进行积分:
ln|y| = (1/2)x 2 + C
求解y
,得到:
y = ± e (1/2)x 2 + C
一阶线性方程
一阶线性微分方程的形式为:
dy/dx + P(x)y = Q(x)
这里通常使用积分因子法:
IF = e ∫P(x) dx
乘以积分因子将方程转换为可以直接积分的形式。
示例:
dy/dx + y = e x
首先,确定积分因子,IF = e ∫1 dx = e x
。
两边乘以e x
:
e x dy/dx + e x y = e 2x
变为:
d/dx(e x y) = e 2x
对两边积分:
e x y = (1/2)e 2x + C
因此,y = (1/2)e x + Ce -x
。
高阶线性微分方程
高阶线性微分方程在数学建模和分析中非常重要。它们通常可以表示为:
a n (x)y (n) + a n-1 (x)y (n-1) + ... + a 1 (x)y' + a 0 (x)y = g(x)
示例:
通常讨论的二阶线性 ODE 是:
y'' + p(x)y' + q(x)y = g(x)
如果g(x) = 0
,则方程是齐次的。
解法
特征方程
如果g(x) = 0
,假设解的形式为y = e rx
。代入 ODE 得到特征方程,通常是r
的多项式:
ar 2 + br + c = 0
求解r
以找到基于不同、重复或复根的通解。解的形式通常是:
- 不同根
r 1 , r 2
:y = c 1 e r 1 x + c 2 e r 2 x
- 重复根
r
:y = (c 1 + c 2 x)e rx
- 复根
α ± βi
:y = e αx (c 1 cos(βx) + c 2 sin(βx))
特殊解
当考虑方程g(x) ≠ 0
时,解由齐次解y h
和通过如常数变异或参数变差等方法获得的特解y p
组成。
常微分方程的应用
ODEs 广泛应用于现实中的建模物理系统和现象。一些常见的例子有:
人口动力学
人口增长常常可以用 ODE 来建模。最简单的模型假设增长率与当前人口规模成正比:
dy/dt = ry
其中r
是增长率。
电路
基尔霍夫定律应用于包含电阻、电容和电感的电路会导致 ODE 生成。例如,一个 RLC 电路可以被描述为:
L(di/dt) + Ri + (1/C)∫i dt = E(t)
牛顿冷却定律
一个物体的热量损失速率与物体和周围环境的温差成正比:
dT/dt = k(T - T env )
质量-弹簧-阻尼系统
这些系统可以通过二阶 ODE 模型化。阻尼谐振子的基本方程为:
m(d 2 x/dt 2 ) + c(dx/dt) + kx = 0
其中m
是质量,c
是阻尼系数,k
是弹簧的刚度。
通过实例理解
变量分离的例子
考虑:
dy/dx = xy
分离变量,我们得到:
dy/y = x dx
对两边积分:
ln|y| = (1/2)x 2 + C
求解y
,得到:
y = Ce (1/2)x 2
一阶线性方程的例子
考虑:
dy/dx + 2y = e x
积分因子为:
IF = e ∫2 dx = e 2x
乘以该因子:
e 2x dy/dx + 2e 2x y = e 3x
这使之更简单:
d/dx(e 2x y) = e 3x
积分右侧,我们得到:
e 2x y = (1/3)e 3x + C
所以y = (1/3)e x + Ce -2x
。
理解常微分方程对分析物理、工程及其它各类系统的行为至关重要。方法和方程类型的多样性为建模大量现象提供了一个稳健的框架。