
本文目标:让没有任何微积分基础的读者,从零开始系统掌握微积分的核心概念——极限、导数、积分、微分方程的定义、性质、几何意义和历史背景,并理解它们如何应用于无人机飞控、光线追踪、深度学习等AI技术中。全文采用”概念→定义→性质→几何意义→计算方法→应用”的六步教学法,配合大量手算例子和工程实例。
第一部分:基础篇
一、微积分是什么——一个关于”变化”和”累积”的故事
1.1 微积分回答的两个核心问题
微积分(Calculus)只做两件事:
| 问题 |
分支 |
通俗理解 |
例子 |
| 变化的速率是多少? |
微分学 |
求”瞬间变化率” |
无人机此时的速度是多大? |
| 累积的总量是多少? |
积分学 |
求”累积总和” |
无人机飞了多远? |
这两个问题是互逆的——微分和积分互为逆运算,这就是微积分最核心的定理:微积分基本定理。
历史:微积分的创立是数学史上最重大的事件之一。17世纪,两位科学家几乎同时独立地创立了微积分:
- 艾萨克·牛顿(Isaac Newton,1643-1727):英国物理学家和数学家。他于1665-1666年(”奇迹年”)创立了”流数术”(fluxions),用于研究运动和引力。但他直到1687年的《自然哲学的数学原理》才公开发表。
- 戈特弗里德·莱布尼茨(Gottfried Leibniz,1646-1716):德国数学家。他于1684年发表了世界上第一篇微积分论文《一种求极大值和极小值的新方法》,比牛顿早三年。我们今天使用的记号 dy/dx 和 ∫ 都是莱布尼茨发明的。
两位科学家的优先权之争持续了数十年,但今天公认两人是独立发明的。现代微积分的基础是19世纪由柯西(Cauchy)和魏尔斯特拉斯(Weierstrass)用极限理论严格化的。
1.2 一个直观的引入——无人机的运动
假设一架无人机在飞行,时间是 t 秒,位置是 s(t) 米。
问题1:在 t=5 秒的那一刻,无人机的速度是多少?
- 不能简单地用”路程÷时间”,因为那算的是平均速度
- 需要知道瞬间速度——这就是导数
问题2:在 0 到 10 秒间,无人机飞了多远?
- 如果速度在变化,不能用”速度×时间”
- 需要将时间分割为无穷多小段,每段内速度近似不变,然后加起来——这就是积分
二、极限——微积分的基石
2.1 极限的直观理解
定义:极限描述的是:当一个变量无限接近某个值时,另一个变量会无限接近什么值。
历史:极限的思想可以追溯到古希腊的穷竭法(阿基米德,公元前3世纪)。但严格的极限定义直到19世纪才建立:
- 1821年,法国数学家柯西(Augustin-Louis Cauchy)在《分析教程》中给出了极限的初步定义
- 1850年代,德国数学家魏尔斯特拉斯(Karl Weierstrass)给出了现代使用的 ϵ-δ 极限定义,彻底消除了微积分中的”无穷小”含糊性
通俗理解:极限就像”永远逼近但可能永远达不到”——你离目标越来越近,误差可以任意小。
2.2 极限的数学定义
通俗定义:
x→alimf(x)=L
读作:”当 x 趋近于 a 时,f(x) 的极限是 L“。
含义:当 x 无限接近 a(但不等于 a)时,f(x) 无限接近 L。
严格的 ϵ-δ 定义(魏尔斯特拉斯):
对于任意给定的 ϵ>0,总存在一个 δ>0,使得当 0<∣x−a∣<δ 时,有 ∣f(x)−L∣<ϵ。
2.3 几个直观例子
例1:limx→2(3x+1)=?
- 当 x=1.9 时,3x+1=6.7
- 当 x=1.99 时,3x+1=6.97
- 当 x=1.999 时,3x+1=6.997
- 显然趋近于 7,所以 limx→2(3x+1)=7
例2:limx→0xsinx=1
这个极限是学习微积分遇到的第一个”意外结果”,它在信号处理、光波衍射中频繁出现。
例3:limx→0x1 — 不存在(因为趋近于无穷大)
2.4 极限的性质
| 性质 |
公式 |
含义 |
| 加法 |
lim(f+g)=limf+limg |
和的极限等于极限的和 |
| 减法 |
lim(f−g)=limf−limg |
差的极限等于极限的差 |
| 乘法 |
lim(f⋅g)=limf⋅limg |
积的极限等于极限的积 |
| 除法 |
limgf=limglimf |
商的极限等于极限的商(分母不为0) |
| 常数倍 |
lim(c⋅f)=c⋅limf |
常数可以提到极限外面 |
2.5 极限在深度学习中的一次出现
神经网络训练中的学习率衰减策略:
η(t)=t→∞limη0⋅1+αt1=0
学习率随训练步数增加而衰减到0,使模型在训练后期精细调整参数。
三、导数——变化的瞬间速率
3.1 导数的定义
定义:函数 f(x) 在 x=a 处的导数,记作 f′(a) 或 dxdf(a),定义为:
f′(a)=h→0limhf(a+h)−f(a)
几何意义:导数就是函数图像在 x=a 处的切线斜率。
物理意义:位置对时间的导数就是速度,速度对时间的导数就是加速度。
历史:导数的概念由牛顿(流数)和莱布尼茨(微商)分别独立创立。但”导数”这个名称和记号 dxdy 来自莱布尼茨,他认为导数就是两个无穷小量的商。”导数”(derivative)一词由拉格朗日在1797年引入。
3.2 一个完整的手算例子
求 f(x)=x2 在 x=2 处的导数:
f′(2)=h→0limh(2+h)2−22=h→0limh4+4h+h2−4=h→0limh4h+h2=h→0lim(4+h)=4
解释:y=x2 在 x=2 处的切线斜率为4,意味着在这一点,x 增加1个单位,y 大约增加4个单位。
3.3 常用函数的导数公式表
| 函数 |
导数 |
例子 |
| 常数 c |
0 |
f(x)=5⇒f′(x)=0 |
| xn |
nxn−1 |
f(x)=x3⇒f′(x)=3x2 |
| sinx |
cosx |
— |
| cosx |
−sinx |
— |
| ex |
ex |
(自身的导数等于自身) |
| lnx |
1/x |
— |
| ax |
axlna |
f(x)=2x⇒f′(x)=2xln2 |
3.4 导数的运算法则
| 法则 |
公式 |
助记 |
| 加法 |
(f+g)′=f′+g′ |
各求各的导,然后相加 |
| 乘法 |
(fg)′=f′g+fg′ |
“前导后不导 + 前不导后导” |
| 除法 |
(gf)′=g2f′g−fg′ |
“上导下不导减上不导下导,除以分母平方” |
| 链式法则 |
(f(g(x)))′=f′(g(x))⋅g′(x) |
“外层导×内层导” |
链式法则(最重要!) 是深度学习的数学基础。神经网络就是一层层的函数复合,反向传播就是反复应用链式法则。
链式法则的手算例子:
求 h(x)=(3x+1)2 的导数。
设 g(x)=3x+1,f(u)=u2。
则 h′(x)=f′(g(x))⋅g′(x)=2(3x+1)⋅3=6(3x+1)
3.5 偏导数——多变量函数的导数
定义:当一个函数有多个变量时(如 f(x,y)),对其中一个变量求导,将其他变量视为常数,称为偏导数。
∂x∂f(x,y)=h→0limhf(x+h,y)−f(x,y)
记号:偏导数用 ∂(读作”偏”或”partial”)表示,而不是 d。
例子:f(x,y)=x2y+3y2
- 对 x 求偏导:∂x∂f=2xy(将 y 视为常数)
- 对 y 求偏导:∂y∂f=x2+6y(将 x 视为常数)
梯度:将一个函数对所有变量的偏导数组成向量,称为梯度:
∇f=[∂x1∂f,∂x2∂f,…,∂xn∂f]T
梯度指向函数上升最快的方向,它的相反方向(负梯度)指向下降最快的方向——这就是梯度下降法的数学基础。
3.6 导数在无人机中的第一次出现
加速度计输出:无人机IMU中的加速度计测量的是加速度,但我们要的是速度和位置。
加速度 a(t)=dtdv⟹速度 v(t)=∫a(t)dt
角速度陀螺仪:陀螺仪测量的是角速度,但姿态需要的是角度:
角速度 ω(t)=dtdθ⟹姿态角 θ(t)=∫ω(t)dt
历史注记:IMU中的”积分”概念直接来自牛顿和莱布尼茨的微积分——每时每刻的角速度(导数)累积得到当前姿态(原函数)。这正是微积分基本定理的物理实现。
第二部分:积分篇
四、积分——累积的总和
4.1 积分的定义
定积分的定义(黎曼积分):
∫abf(x)dx=n→∞limi=1∑nf(xi∗)Δx
通俗理解:将区间 [a,b] 分成 n 个小区间,每个小区间用矩形近似,然后将所有矩形面积相加。当区间数无限增加时,这个和的极限就是定积分。
几何意义:定积分 ∫abf(x)dx 等于函数曲线与 x 轴之间的”有向面积”(x 轴上方为正,下方为负)。
历史:积分的概念可以追溯到古希腊的阿基米德(Archimedes,公元前287-212年),他用”穷竭法”计算抛物线弓形的面积,这是积分思想的雏形。17世纪,开普勒、卡瓦列里、费马等数学家拓展了积分方法。但现代积分的最终形式由:
- 牛顿:将积分视为”流数”的逆运算
- 莱布尼茨:发明了 ∫ 符号(拉长的S,表示Sum求和)
- 柯西(1823):给出了定积分的严格定义
- 黎曼(1854):给出了今天教科书上使用的黎曼积分定义
4.2 微积分基本定理——连接微分和积分
定理(牛顿-莱布尼茨公式):
如果 F′(x)=f(x)(即 F 是 f 的原函数),则:
∫abf(x)dx=F(b)−F(a)
这个定理的伟大之处:它将看似完全不同的两个概念——变化率(导数)和累积和(积分)——联系在了一起。计算积分不必再通过”无限求和”的繁琐过程,只需要找到原函数然后做减法。
历史:牛顿和莱布尼茨分别独立发现了这个定理。牛顿在1666年的手稿中已有表述,而莱布尼茨在1693年发表了这个定理。英国数学家巴罗(Isaac Barrow,牛顿的老师)在1669年已经认识到微分和积分的互逆关系,但没有明确提出这个定理。
4.3 常用积分公式表
| 导数公式 |
对应的积分公式 |
| dxd(xn+1)=(n+1)xn |
∫xndx=n+1xn+1+C |
| dxd(sinx)=cosx |
∫cosxdx=sinx+C |
| dxd(−cosx)=sinx |
∫sinxdx=−cosx+C |
| dxd(ex)=ex |
∫exdx=ex+C |
| dxd(lnx)=x1 |
∫x1dx=ln∣x∣+C |
其中 C 是积分常数——因为常数求导为0,所以原函数可以相差任意常数。
4.4 完整的定积分计算例子
计算 ∫02x2dx:
- 找到原函数:F(x)=3x3(因为 dxd3x3=x2)
- 代入上下限:
∫02x2dx=F(2)−F(0)=323−303=38
几何验证:y=x2 从 x=0 到 x=2 下方的面积是 38≈2.667。
4.5 积分在工程中的应用
| 应用 |
微分(测量) |
积分(计算) |
| 位置-速度-加速度 |
加速度计 a(t) |
v(t)=∫a(t)dt, s(t)=∫v(t)dt |
| 电流-电量 |
电流 I(t) |
电量 Q=∫I(t)dt |
| 热量-温度变化 |
加热速率 P(t) |
热量 H=∫P(t)dt |
| 概率密度-概率 |
PDF f(x) |
概率 P(a≤X≤b)=∫abf(x)dx |
第三部分:核心应用篇
五、微分方程——描述变化规律的数学语言
5.1 微分方程的定义
定义:包含未知函数及其导数的方程称为微分方程。
最简单的例子:牛顿冷却定律
dtdT=−k(T−Tenv)
含义:物体的温度变化率与温度差成正比。
历史:微分方程的起源与微积分本身一样早。
- 1671年,牛顿在他的著作中将许多物理问题表述为微分方程
- 1690年代,雅各布·伯努利(Jacob Bernoulli)和约翰·伯努利(Johann Bernoulli)兄弟系统研究微分方程
- 18世纪,欧拉、丹尼尔·伯努利、拉普拉斯等数学家将微分方程应用于物理、力学和天文学
5.2 无人机飞行动力学中的微分方程
六自由度运动方程就是一组微分方程:
力的方程:
mdtdv=Faero+Fgravity+Fthrust
含义:质量 × 加速度(速度的导数)= 所有力的和
力矩方程:
Idtdω=τaero+τcontrol
含义:惯性张量 × 角加速度(角速度的导数)= 所有力矩的和
姿态的方程:
dtdϕθψ=100sinϕtanθcosϕsinϕ/cosθcosϕtanθ−sinϕcosϕ/cosθpqr
含义:欧拉角的变化率(左)由当前姿态和角速度(右)共同决定。
这些微分方程的意义:只要知道当前的状态(位置、速度、姿态、角速度)和控制输入(油门、舵面偏角),就可以通过求解微分方程预测无人机未来的状态——这就是仿真的数学基础。
5.3 如何求解微分方程——数值方法
绝大多数微分方程没有解析解(写不出公式),所以工程中都用数值方法近似求解。
欧拉方法——最简单的数值解法:
xn+1=xn+Δt⋅f(xn,tn)
其中 f(x,t)=dtdx,Δt 是时间步长。
从零开始的例子:
假设 v(t) 是无人机速度,加速度 a=9.8m/s2,Δt=0.1s,初速度 v0=0。
计算前3步:
- 第0步:v0=0
- 第1步:v1=0+0.1×9.8=0.98m/s
- 第2步:v2=0.98+0.1×9.8=1.96m/s
- 第3步:v3=1.96+0.1×9.8=2.94m/s
更精确的方法:四阶龙格-库塔法(Runge-Kutta 4, RK4)——这是PX4固件和Gazebo仿真中使用的标准方法,精度远高于欧拉法。
六、梯度与优化——深度学习的学习算法
6.1 梯度下降法
定义:梯度下降是机器学习中最核心的优化算法。它利用一阶导数(梯度)信息,沿着函数下降最快的方向更新参数:
θn+1=θn−η∇L(θn)
其中:
- θ 是模型的参数(权重)
- η 是学习率(步长)
- ∇L 是损失函数对参数的梯度
几何理解:想象你站在一座山上(损失函数曲面),目标是下到山谷(找到最小值)。你环顾四周,找到最陡的下降方向,然后迈一步。重复这个过程——这就是梯度下降。
历史:梯度下降法由柯西在1847年首次提出,但当时没有计算机,只能用于小规模问题。2010年代,随着深度学习的兴起,梯度下降成为训练神经网络的核心算法。2014年,Kingma和Ba提出的Adam优化器结合了动量法和自适应学习率,是目前最广泛使用的变体。
6.2 反向传播——链式法则的工程奇迹
反向传播(Backpropagation)是训练神经网络的方法,其数学本质就是反复应用链式法则。
以一个3层神经网络为例:
L=Loss(f3(f2(f1(x))))
要求 ∂W1∂L(第1层的权重梯度):
∂W1∂L=∂f3∂L⋅∂f2∂f3⋅∂f1∂f2⋅∂W1∂f1
从输出层开始,逐层向前求导——这就是”反向”传播名称的由来。
一个极小但完整的例子:
假设 f(x)=wx+b(线性层),损失 L=21(ypred−ytrue)2。
给定 x=2,ytrue=5,w=1,b=0:
- 前向传播:ypred=1×2+0=2,L=21(2−5)2=4.5
- 反向传播:
- ∂ypred∂L=ypred−ytrue=2−5=−3
- ∂w∂ypred=x=2
- ∂w∂L=∂ypred∂L⋅∂w∂ypred=−3×2=−6
- 更新参数:wnew=w−η⋅(−6)=w+6η
这就是YOLO、端到端自动驾驶、所有神经网络训练的数学本质——链式法则 × 梯度下降。
6.3 导数在光线追踪中的应用
光线追踪中的重要性采样需要从概率分布中采样方向。而概率分布的累积分布函数(CDF)的反函数的导数决定采样密度。
微表面BRDF模型中的几何项 G 涉及对入射角和出射角的导数——这些导数决定了材质在光线追踪中的表现。
七、微积分知识体系总览
7.1 完整概念地图
1 2 3 4 5 6 7 8 9 10 11
| 极限 (极限的ε-δ定义) ↓ 导数 (瞬间变化率) ←────→ 积分 (累积和) ↓ ↓ 导数公式 积分公式 ↓ ↓ 链式法则 微积分基本定理 (连接两者) ↓ 偏导数 → 梯度 → 梯度下降法 → 深度学习训练 ↓ 微分方程 → 数值解法 → 物理仿真/无人机飞控
|
7.2 各概念在AI/无人机中的出现频率
| 概念 |
🛩️飞行动力学 |
🎨光线追踪 |
🤖YOLO |
🔗端到端 |
| 极限 |
⭐ |
⭐⭐⭐ |
— |
⭐ |
| 导数 |
⭐⭐⭐ |
⭐⭐⭐ |
⭐⭐ |
⭐⭐⭐ |
| 偏导数 |
⭐⭐ |
⭐⭐ |
⭐⭐⭐ |
⭐⭐⭐ |
| 梯度 |
⭐ |
⭐⭐ |
⭐⭐⭐ |
⭐⭐⭐ |
| 链式法则 |
— |
— |
⭐⭐⭐ |
⭐⭐⭐ |
| 积分 |
⭐⭐⭐ |
⭐⭐⭐ |
— |
⭐ |
| 微分方程 |
⭐⭐⭐ |
— |
— |
⭐⭐ |
| 数值积分 |
⭐⭐⭐ |
⭐⭐ |
— |
⭐⭐ |
| 梯度下降 |
— |
— |
⭐⭐⭐ |
⭐⭐⭐ |
7.3 给自学者的推荐阅读顺序
- 先理解极限(第2章)——微积分的基础,但不必陷在 ϵ-δ 中
- 再理解导数(第3章)——理解”瞬间变化率”的几何意义
- 理解积分(第4章)——理解”累积和”的几何意义
- 理解微积分基本定理(第4.2节)——微积分的精髓
- 学习微分方程(第5章)——看它如何描述真实世界
- 最后理解反向传播(第6.2节)——理解深度学习为什么能工作
推荐资源:
- 3Blue1Brown《微积分的本质》(YouTube系列)——最好的可视化教程
- Gilbert Strang《微积分》(MIT 18.01)——经典教材
- James Stewart《微积分》(Calculus)——最全面的入门教材
- Khan Academy 微积分——免费、循序渐进
7.4 核心公式速查卡
| 公式 |
含义 |
| f′(a)=limh→0hf(a+h)−f(a) |
导数的定义 |
| dxdxn=nxn−1 |
幂函数求导 |
| dxd[f(g(x))]=f′(g(x))g′(x) |
链式法则(深度学习的数学基础) |
| ∫abf(x)dx=F(b)−F(a) |
微积分基本定理 |
| θn+1=θn−η∇L(θn) |
梯度下降法 |
| dtdv=m1∑F |
牛顿第二定律(微分方程形式) |
| ∇f=[∂f/∂x1,…,∂f/∂xn]T |
梯度(指向最陡上升方向) |
参考文献与推荐资源
- Stewart, J. (2015). Calculus: Early Transcendentals. 8th ed., Cengage Learning. — 全球使用最广的微积分教材
- Strang, G. (2017). Calculus. 2nd ed., Wellesley-Cambridge Press. — MIT教授的经典之作
- 3Blue1Brown. Essence of Calculus (YouTube系列). — 最好的微积分可视化教程
- Goodfellow, I., Bengio, Y., & Courville, A. (2016). Deep Learning. MIT Press. Chapter 4-6. — 深度学习的数值计算基础
- Kingma, D. P., & Ba, J. (2015). Adam: A Method for Stochastic Optimization. ICLR. — Adam优化器的原始论文
- Rumelhart, D. E., Hinton, G. E., & Williams, R. J. (1986). Learning representations by back-propagating errors. Nature, 323, 533-536. — 反向传播算法的原始论文
- Burden, R. L., & Faires, J. D. (2010). Numerical Analysis. 9th ed., Brooks/Cole. — 数值微分与数值积分
- Stevens, B. L., Lewis, F. L., & Johnson, E. N. (2015). Aircraft Control and Simulation. 3rd ed., Wiley. — 飞行动力学中的微分方程应用