Skip to content

第一章 绪论 ​

第一章 绪论 ​

来源:课堂讲义《第一章 绪论1.pdf》· 主讲:周盾(njust.edu.cn) 笔记由课堂 PDF 整理,公式按数学含义重建;模糊处标 [?] 待校对。

课程信息 ​

  • 教材:王高雄等《常微分方程》(第四版),高等教育出版社
  • 参考书:
    1. 朱思铭《常微分方程学习辅导与习题解答》
    2. 丁同仁《常微分方程教程》
    3. 王柔怀、伍卓群《常微分方程讲义》
    4. 庞特里亚金《常微分方程》(倪明康译)
    5. 张芷芬等《微分方程定性理论》
  • 学时:第 5–14 周,共 48 学时
  • 考核:总成绩 = 平时 20% + 期中 10% + 期末 70%;每周二交作业
  • 课程内容:第一章 绪论 → 第二章 一阶微分方程的初等积分法 → 第三章 解的存在定理 → 第四章 高阶微分方程 → 第五章 线性微分方程组 → 第六章 非线性微分方程 → 第七章 一阶线性偏微分方程* → 第八章 边值问题
  • 课程目的:掌握各类可求解方程(组)的类型与解法;熟悉解的基本性质(存在性、唯一性);了解稳定性分析、定性分析等基本方法

微分方程的背景 ​

  • 微分方程及相关学科几乎贯穿各类科学研究,是现代数学的重要分支,也是解决实际问题的工具
  • 与代数方程不同:代数方程寻求未知量与已知量的关系;微分方程还涉及未知量的变化率(导数),描述"变化问题"
    • 例:自由落体关心"下落距离如何随时间变化";火箭关心"飞行轨道如何随时间演化"
  • 基本思想:找出已知函数与未知函数的关系,从含未知函数及其导数的方程中求得未知函数表达式
  • 历史:与微积分同时产生——牛顿用级数求解简单微分方程;雅各布·贝努利、欧拉、克雷洛、达朗贝尔、拉格朗日丰富其理论;Poincaré 与 Lyapunov 奠定现代微分方程理论(动力系统)基础

1.1 常微分方程模型 ​

RL / RLC 电路 ​

  • 基尔霍夫(Kirchhoff)第二定律:闭合回路中所有支路电压的代数和为零

数学摆 ​

人口模型(Malthus) ​

  • 假设:人口自然增长中,净相对增加率(单位时间净增长数 / 人口总数)为常数 rr

人口模型的改进:logistic 模型(Verhulst) ​

  • 引入环境最大容纳量 NmN_m,假设净相对增长率为

r(1−N(t)Nm)r\left(1 - \frac{N(t)}{N_m}\right)

SI 传染病模型 ​

  • 易感染者(Susceptible)与已感染者(Infective)两类

Volterra 被捕食–捕食模型 ​

  • 时刻 tt 被食鱼总数 x(t)x(t),捕食鱼总数 y(t)y(t);单位时间相遇次数 bxybxy;捕食鱼自然减少率与 yy 成正比
  • 另有两种群竞争模型

Lorenz 方程(蝴蝶效应) ​

  • Edward Lorenz(1963,MIT)给出第一个混沌现象——蝴蝶效应,为大气流体动力学模型的简化常微分方程组

建模总结 ​

  • 微分方程反映量与量之间的关系,与时间有关,是动态系统
  • 从自然规律出发,考虑主要因素,构造含自变量、未知函数及其导数的关系式 → 数学模型
  • 研究解和解结构的性质,与实际对照,不断改进模型;可由微分方程发现或预测新规律

1.2 微分方程的基本概念 ​

微分方程的定义 ​

联系**自变量、未知函数及未知函数的导数(或微分)**的关系式称为微分方程。

例 1:下列关系式都是微分方程:

dydx=x2(1)\frac{dy}{dx} = x^2 \tag{1}

x dy−y dx=0(2)x\,dy - y\,dx = 0 \tag{2}

d2xdt2+tx(dxdt)3+x=0(3)\frac{d^2x}{dt^2} + tx\left(\frac{dx}{dt}\right)^3 + x = 0 \tag{3}

d4xdt4+5d2xdt2+3x=sin⁡t[?](4)\frac{d^4x}{dt^4} + 5\frac{d^2x}{dt^2} + 3x = \sin t \tag{4}[?]

∂z∂x+∂z∂y=z(5)\frac{\partial z}{\partial x} + \frac{\partial z}{\partial y} = z \tag{5}

∂2u∂x2+∂2u∂y2+xy−u=0[?](6)\frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} + xy - u = 0 \tag{6}[?]

常微分方程与偏微分方程 ​

  • 常微分方程(ODE):自变量个数只有一个
    • 例:(1)–(4)
  • 偏微分方程(PDE):自变量个数为两个或两个以上
    • 例:(5)(6)
  • 注:本课程主要研究常微分方程,简称微分方程或方程

微分方程的阶 ​

微分方程中出现的未知函数的最高阶导数(或微分)的阶数称为微分方程的阶数。

  • (1) 一阶;(2) 一阶;(3) 二阶;(4) 四阶

nn 阶微分方程的一般形式:

F(x,y,dydx,…,dnydxn)=0F\left(x, y, \frac{dy}{dx}, \ldots, \frac{d^n y}{dx^n}\right) = 0

其中 xx 是自变量,yy 是未知函数,且一定含有 dnydxn\dfrac{d^n y}{dx^n}。

线性与非线性 ​

如果方程 F(x,y,dydx,…,dnydxn)=0F\left(x, y, \dfrac{dy}{dx}, \ldots, \dfrac{d^n y}{dx^n}\right) = 0 的左端是 y,y′,…,y(n)y, y', \ldots, y^{(n)} 的一次有理式,则称为 nn 阶线性方程;否则为非线性方程。

nn 阶线性微分方程的一般形式:

an(x)dnydxn+an−1(x)dn−1ydxn−1+⋯+a1(x)dydx+a0(x)y=f(x)a_n(x)\frac{d^n y}{dx^n} + a_{n-1}(x)\frac{d^{n-1} y}{dx^{n-1}} + \cdots + a_1(x)\frac{dy}{dx} + a_0(x)y = f(x)

其中 an(x),…,a1(x),a0(x),f(x)a_n(x), \ldots, a_1(x), a_0(x), f(x) 为 xx 的已知函数。

  • 线性:(1)(2)(4);非线性:(3)(含 (dxdt)3\left(\frac{dx}{dt}\right)^3 项)

微分方程的解 ​

定义:若函数 y=φ(x)y = \varphi(x) 在区间 II 上有直到 nn 阶的连续导数,且对 ∀x∈I\forall x \in I 有

F(x,φ(x),φ′(x),…,φ(n)(x))≡0F\left(x, \varphi(x), \varphi'(x), \ldots, \varphi^{(n)}(x)\right) \equiv 0

则称 y=φ(x)y = \varphi(x) 为方程在 II 上的一个解,称为显式解。

例:验证 y=sin⁡xy = \sin x 与 y=cos⁡xy = \cos x 都是 y′′+y=0y'' + y = 0 在 (−∞,+∞)(-\infty, +\infty) 上的解。 (y=sin⁡x⇒y′=cos⁡x,y′′=−sin⁡x⇒y′′+y=−sin⁡x+sin⁡x=0y = \sin x \Rightarrow y' = \cos x, y'' = -\sin x \Rightarrow y'' + y = -\sin x + \sin x = 0 ✓)

显式解与隐式解 ​

若关系式 Φ(x,y)=0\Phi(x, y) = 0 确定的隐函数 y=φ(x)y = \varphi(x) 是方程的解,则称 Φ(x,y)=0\Phi(x, y) = 0 为方程的隐式解。

例:一阶方程 dydx=−xy\dfrac{dy}{dx} = - \dfrac{x}{y} 有显式解 y=−1−x2y = -\sqrt{1-x^2} 与 y=1−x2y = \sqrt{1-x^2},隐式解 x2+y2=1x^2 + y^2 = 1。

注:显式解与隐式解统称为微分方程的解。

通解与特解 ​

定义:若微分方程的解中含有任意常数,且所含相互独立的任意常数的个数与方程的阶数相同,则称为方程的通解。

例:y=c1sin⁡x+c2cos⁡xy = c_1 \sin x + c_2 \cos x(c1,c2c_1, c_2 为任意常数)是 y′′+y=0y'' + y = 0 的通解。

nn 阶微分方程通解的一般形式:

y=φ(x,c1,c2,…,cn)y = \varphi(x, c_1, c_2, \ldots, c_n)

其中 c1,…,cnc_1, \ldots, c_n 为相互独立的任意常数(独立性的行列式判别条件:Wronski 型行列式 ≠0\neq 0)。

例:验证 y=c1ex+c2e−x+c3e2x+12y = c_1 e^x + c_2 e^{-x} + c_3 e^{2x} + \dfrac{1}{2} 是 y′′′−2y′′−y′+2y=6y''' - 2y'' - y' + 2y = 6 的通解(代入计算可证,且 Wronski 行列式非零)。

  • 隐式通解(通积分):隐式解中含有任意常数且个数与阶数相同
  • 以后不区分显式通解与隐式通解,统称通解

特解:在通解中给任意常数以确定的值而得到的解。

例:通解 y=c1sin⁡x+c2cos⁡xy = c_1 \sin x + c_2 \cos x 中取 c1=0,c2=1c_1 = 0, c_2 = 1 得 y=cos⁡xy = \cos x;取 c1=1,c2=0c_1 = 1, c_2 = 0 得 y=sin⁡xy = \sin x,均为 y′′+y=0y'' + y = 0 的特解。

定解条件与初值问题(Cauchy 问题) ​

  • 定解条件:根据实际问题附加的条件;求满足定解条件的解的问题称定解问题
  • 初始条件:nn 阶微分方程的初始条件指 nn 个条件(当 x=x0x = x_0 时):

y(x0)=y0,y′(x0)=y0(1),…,y(n−1)(x0)=y0(n−1)y(x_0) = y_0, \quad y'(x_0) = y_0^{(1)}, \quad \ldots, \quad y^{(n-1)}(x_0) = y_0^{(n-1)}

其中 x0,y0,y0(1),…,y0(n−1)x_0, y_0, y_0^{(1)}, \ldots, y_0^{(n-1)} 为给定的 n+1n+1 个常数。

  • 定解条件为初始条件时,相应定解问题称为初值问题(Cauchy 问题)

例(P9):验证 y=c1ex+c2e−4xy = c_1 e^x + c_2 e^{-4x} 是 y′′+5y′+4y=0y'' + 5y' + 4y = 0 的通解,并求满足 y(0)=2,y′(0)=1y(0) = 2, y'(0) = 1 的特解。

解:由初始条件得方程组

{c1+c2=2c1−4c2=1⇒c1=95,  c2=15[?]\begin{cases} c_1 + c_2 = 2 \\ c_1 - 4c_2 = 1 \end{cases} \Rightarrow c_1 = \frac{9}{5}, \; c_2 = \frac{1}{5}[?]

(PDF 原文计算为 c1=3,c2=−1c_1 = 3, c_2 = -1,与 y(0)=2, y′(0)=1y(0)=2,\ y'(0)=1 核对不符,疑为讲义笔误,待校对)

校对点

讲义原文解出 c1=3, c2=−1c_1 = 3,\ c_2 = -1,代回 y′(0)=3−4(−1)=7≠1y'(0) = 3 - 4(-1) = 7 \neq 1,与初始条件矛盾;按 y(0)=2,y′(0)=1y(0)=2, y'(0)=1 应得 c1=9/5,c2=1/5c_1 = 9/5, c_2 = 1/5。

积分曲线与方向场 ​

  • 积分曲线:一阶方程 dydx=f(x,y)\dfrac{dy}{dx} = f(x, y) 的解 y=φ(x)y = \varphi(x) 在 xyxy 平面上表示的曲线;通解对应一族积分曲线
  • 方向场:在 f(x,y)f(x, y) 的定义域 DD 内每一点画一个以 f(x,y)f(x, y) 为斜率、中心在该点的直线段,带有这种直线段的区域称为方程的方向场
  • 等斜线:方向场中方向相同的点的几何轨迹。方程 dydx=f(x,y)\dfrac{dy}{dx} = f(x, y) 的等斜线为 f(x,y)=kf(x, y) = k(kk 为参数)

微分方程组 ​

驻定(自治)与非驻定(非自治) ​

  • 驻定(自治)系统:右端不显含自变量
  • 非驻定(非自治)系统:右端显含自变量

相空间、轨线与奇点 ​

  • 相空间:不含自变量、仅由未知函数组成的空间
  • 轨线:积分曲线在相空间中的投影
  • 奇点(平衡点):平衡解(驻定解、常数解)对应的点

例:{u˙=vv˙=−u\begin{cases} \dot{u} = v \\ \dot{v} = -u \end{cases},u(0)=1,v(0)=0u(0) = 1, v(0) = 0 的解为 (u(t),v(t))=(cos⁡t,−sin⁡t)(u(t), v(t)) = (\cos t, -\sin t),轨线为相平面中的圆(u2+v2=1u^2 + v^2 = 1)。

垂直等倾线、水平等倾线 ​

作业 ​

  • p12:1(1、3、6)
  • p13:3(1、5、8),4,6,7

相关笔记 ​

  • ODE常微分方程

用 VitePress 搭建 · 托管在 GitHub Pages