Book3 Ch09 深入圆锥曲线 数学要素 从加减乘除到机器学习
Book3 Ch09 深入圆锥曲线 数学要素 从加减乘除到机器学习
Book3 Ch09 深入圆锥曲线 数学要素 从加减乘除到机器学习
9 深入圆锥曲线
探寻和数据科学、机器学习之间联系
地球是人类的摇篮,但我们不能永远生活在摇篮里。
Earth is the cradle of humanity, but one cannot live in a cradle forever.
—— 康斯坦丁·齐奥尔科夫斯基 (Konstantin Tsiolkovsky) | 俄罗斯火箭专家 | 1857 ~ 1935
◄ matplotlib.patches.Rectangle() 绘制通过定位点,以及设定宽度和高度的矩形
◄ matplotlib.pyplot.contour() 绘制等高线图
◄ matplotlib.pyplot.contourf() 绘制填充等高线图
◄ numpy.cosh() 双曲余弦函数
◄ numpy.isinf() 判断是否存在无穷
◄ numpy.maximum() 计算最大值
◄ numpy.sinh() 双曲正弦函数
◄ numpy.tanh() 双曲正切函数
◄ sympy.Eq() 定义符号等式
◄ sympy.evalf() 将符号解析式中未知量替换为具体数值
◄ sympy.plot_implicit()绘制隐函数方程
◄ sympy.symbols() 定义符号变量
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]
Page 2 | Chapter 9 深入圆锥曲线 | Book 3《数学要素》 | 鸢尾花书:从加减乘除到机器学习
正圆
分类
椭圆
双曲线
抛物线
离心率
有趣的圆锥曲线
双曲线和椭圆的叠加
和矩形相切
0<p<1
取值对应形状
p=1
圆锥曲线
1<p<2
超椭圆
p=2
p=
p和q不同值
超椭球
双曲正弦
双曲余弦
双曲函数
双曲正切
和自然指数函数关系
圆锥曲线的一般形式
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]
Page 3 | Chapter 9 深入圆锥曲线 | Book 3《数学要素》 | 鸢尾花书:从加减乘除到机器学习
9.1 圆锥曲线:探索星辰大海
虽然正圆、椭圆、抛物线、双曲线这样的数学概念现在见诸于中学课本,但是它们现如今依
旧展现着巨大能量。比如,在星辰大海的征途中,圆锥曲线扮演重要角色。
Hyperbolic trajectory
Unbound orbits
图 1. 航天器的几种轨道
探索火星约每 26 个月有一个发射窗口,这是因为地球在低轨道绕太阳运行,而火星在高轨道
绕行。地球和火星的公转周期不同,两个行星大约每 26 个月“相遇”一次,也就是说地球与火星之
间的距离最近。
航天器首先进入绕太阳圆周运动的低轨道。
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]
Page 4 | Chapter 9 深入圆锥曲线 | Book 3《数学要素》 | 鸢尾花书:从加减乘除到机器学习
太空船在低轨道 A 点处上瞬间加速后,进入一个椭圆形的转移轨道。注意,加速瞬间火星位
于 B。太空船由此椭圆轨道的近拱点开始,抵达远拱点后再瞬间加速,进入火星所在的目标轨
道。反过来,霍曼转移轨道亦可将太空船送往较低的轨道,不过是两次减速而非加速。
Earth at arrival, C
Mars at arrival, D
Earth at launch, A
图 2. 探索火星的霍曼轨道
9.2 离心率:联系不同类型圆锥曲线
不同类型圆锥曲线可以通过同离心率 (eccentricity) e 联系起来:
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]
Page 5 | Chapter 9 深入圆锥曲线 | Book 3《数学要素》 | 鸢尾花书:从加减乘除到机器学习
图 3. 离心率连续变化条件下一组圆锥曲线
9.3 一组有趣的圆锥曲线
本节介绍一组有趣的圆锥曲线,解析式如下:
x12 x22 xx
+ − 2 1 2 = 1 (2)
m2 n2 mn
Ellipse Hyperbola
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]
Page 6 | Chapter 9 深入圆锥曲线 | Book 3《数学要素》 | 鸢尾花书:从加减乘除到机器学习
(a), 0 ρ 2 (b), 0 ρ 1
图 4. m = n = 1,圆锥曲线随 ρ 变化,ρ 非负
(a), 2 ρ 0 (b), 1 ρ 0
图 5. m = n = 1,圆锥曲线随 ρ 变化,ρ 非正
当 m = n = 1 时,且 ρ = 1 时,(2) 为:
( x − y) =1
2
(3)
以上解析式对应两条直线:
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]
Page 7 | Chapter 9 深入圆锥曲线 | Book 3《数学要素》 | 鸢尾花书:从加减乘除到机器学习
x − y = 1, x − y = −1 (4)
(a), 0 ρ 2 (b), 2 ρ 0
9.4 特殊椭圆:和给定矩形相切
这一节,我们要在特殊条件约束下绘制椭圆。
给定如图 7 所示的三类矩形,假定它们的中心都位于原点。本节绘制和矩形四个边相切的椭
圆。椭圆可以是正椭圆,也可以是旋转椭圆。
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]
Page 8 | Chapter 9 深入圆锥曲线 | Book 3《数学要素》 | 鸢尾花书:从加减乘除到机器学习
2n 2n
2n
2m 2m
2m
图 7. m、n 大小关系不同的矩形
x12 x22 2 x1 x2
+ − = 1− 2 (5)
m2 n2 mn
(5) 可以进一步写成:
1 x12 x22 2 x1 x2
+ − =1 (6)
1 − 2 m2 n2 mn
如图 8 所示,以矩形的中心为原点构造平面直角坐标系,容易计算得到矩形和椭圆相切的切
点 A、B、C、D 的坐标为:
A ( m, n ) , B ( m, n ) , C ( −m, − n ) , D ( − m, −n ) (7)
请大家格外注意 AC 连线,我们将在本系列丛书的条件概率和线性回归话题中谈到这条直
线。
x2
(n, 0) B
( m, 0) (m, 0)
x1
D ( n, 0)
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]
Page 9 | Chapter 9 深入圆锥曲线 | Book 3《数学要素》 | 鸢尾花书:从加减乘除到机器学习
图 8. 四个切点的位置
正椭圆
当 ρ = 0 时,椭圆为正椭圆,即,
x12 x22
+ =1 (8)
m2 n2
x2
C A
x1
图 9. 当 ρ = 0 时,四个切点的位置
线段
当 ρ = 1 时,椭圆退化为一条线段,对应解析式为:
x1 x2
− =0 (10)
m n
当 ρ = −1 时,椭圆也是一条线段:
x1 x2
+ =0 (11)
m n
两种情况对应的图像为图 10。
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]
Page 10 | Chapter 9 深入圆锥曲线 | Book 3《数学要素》 | 鸢尾花书:从加减乘除到机器学习
x2
ρ=1
x1
ρ=
图 10. 当 ρ = ±1 时,椭圆退化成线段
旋转椭圆
图 11. m = n 时,和给定正方形相切椭圆
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]
Page 11 | Chapter 9 深入圆锥曲线 | Book 3《数学要素》 | 鸢尾花书:从加减乘除到机器学习
二元高斯分布
我们之所以讨论这种特殊形态的椭圆,是因为它和二元高斯分布的概率密度函数直接相关。
二元高斯分布 (bivariate Gaussian distribution) 的概率密度函数 fX,Y(x,y) 解析式如下:
x − 2 2
− − − −
exp
1 1 1 x y y
f X ,Y ( x, y ) = − 2 X ,Y + (12)
2 (1 − 2 ) X
X X Y Y
2π X Y 1 − X2 ,Y X ,Y X Y Y
Ellipse
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]
Page 12 | Chapter 9 深入圆锥曲线 | Book 3《数学要素》 | 鸢尾花书:从加减乘除到机器学习
9.5 超椭圆:和范数有关
超椭圆 (superellipse) 是对椭圆的拓展,最常见的超椭圆的解析式为:
p p
x1 x
+ 2 =1 (13)
a b
还有两个特殊的情况,当 p = 1 时,超椭圆图形为菱形:
x1 x2
+ =1 (14)
a b
当 p = +∞时,超椭圆图形为长方形,对应的解析式为:
x x2
max 1 , =1 (15)
a b
第一个例子
当 a = 2,b = 1 时,超椭圆的解析式为:
p p
x1 x
+ 2 =1 (16)
2 1
图 14 所示为 p 取不同值时,超椭圆的形状。
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]
Page 13 | Chapter 9 深入圆锥曲线 | Book 3《数学要素》 | 鸢尾花书:从加减乘除到机器学习
2 2 2
1 1 1
x2 0 x2 0 x2 0
1 1 1
2 2 2
2 1 0 1 2 2 1 0 1 2 2 1 0 1 2
x1 x1 x1
2 2 2
1 1 1
x2 0 x2 0 x2 0
1 1 1
2 2 2
2 1 0 1 2 2 1 0 1 2 2 1 0 1 2
x1 x1 x1
2 2 2
1 1 1
x2 0 x2 0 x2 0
1 1 1
2 2 2
2 1 0 1 2 2 1 0 1 2 2 1 0 1 2
x1 x1 x1
第二个例子
当 a = 1,b = 1 时,超椭圆的解析式为:
x1 + x2 = 1
p p
(17)
图 15 所示为 p 取不同值时,超椭圆的形状。
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]
Page 14 | Chapter 9 深入圆锥曲线 | Book 3《数学要素》 | 鸢尾花书:从加减乘除到机器学习
2 2 2
1 1 1
x2 0 x2 0 x2 0
1 1 1
2 2 2
2 1 0 1 2 2 1 0 1 2 2 1 0 1 2
x1 x1 x1
2 2 2
1 1 1
x2 0 x2 0 x2 0
1 1 1
2 2 2
2 1 0 1 2 2 1 0 1 2 2 1 0 1 2
x1 x1 x1
2 2 2
1 1 1
x2 0 x2 0 x2 0
1 1 1
2 2 2
2 1 0 1 2 2 1 0 1 2 2 1 0 1 2
x1 x1 x1
p 和 q 两个参数
将 (13) 解析式进一步推广,得到如下二维平面的超椭圆解析式:
p q
x1 x
+ 2 =1 (18)
a b
其中,p 和 q 为正数。
x1 + x2 = 1
p q
(19)
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]
Page 15 | Chapter 9 深入圆锥曲线 | Book 3《数学要素》 | 鸢尾花书:从加减乘除到机器学习
超椭球
x1 r x2 r r x3 t
+ + =1 (20)
a b c
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]
Page 16 | Chapter 9 深入圆锥曲线 | Book 3《数学要素》 | 鸢尾花书:从加减乘除到机器学习
t = 0.8
t=1
t=2
t=3
( )
1 p
D p
= xi
p 1/ p
= x1 + x2 + + xD
p p
x (21)
i =1
p
其中,
x = x1 xD
T
x2 (22)
本系列丛书将在《矩阵力量》一册系统讲解范数。
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]
Page 17 | Chapter 9 深入圆锥曲线 | Book 3《数学要素》 | 鸢尾花书:从加减乘除到机器学习
x2
x1
9.6 双曲函数:基于单位双曲线
当 a = 1 和 b = 1 时,双曲线为单位双曲线 (unit hyperbola):
类似前文提到过的三角函数和单位圆之间关系,单位双曲线可以用来定义双曲函数
(hyperbolic function)。
双曲正切 tanh(),可以通过如下比例计算得到:
sinh
tanh = (24)
cosh
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]
Page 18 | Chapter 9 深入圆锥曲线 | Book 3《数学要素》 | 鸢尾花书:从加减乘除到机器学习
x2 − y 2 = 1
coshθ
sinhθ
Area = θ
1
图 19. 单位双曲线和双曲函数的关系
4 coshθ
tanhθ
2
4 sinhθ
8
8 6 4 2 0 2 4 6 8
θ
表 1. 用英文表达双曲函数
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]
Page 19 | Chapter 9 深入圆锥曲线 | Book 3《数学要素》 | 鸢尾花书:从加减乘除到机器学习
和指数函数关系
exp ( ) − exp ( − )
sinh =
2
exp ( ) + exp ( − )
cosh = (25)
2
sinh exp ( ) − exp ( − )
tanh = =
cosh exp ( ) + exp ( − )
exp(θ)
4 8
exp( θ) exp(θ)
2 6 exp( θ)
0 4
coshθ
2 2
sinhθ
4 0
2 0 2 2 0 2
θ θ
9.7 圆锥曲线一般式
圆锥曲线的一般形如下:
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]
Page 20 | Chapter 9 深入圆锥曲线 | Book 3《数学要素》 | 鸢尾花书:从加减乘除到机器学习
满足下列条件,圆锥曲线为正圆:
满足下列条件,圆锥曲线为正椭圆,即没有旋转:
满足下列条件,圆锥曲线为正双曲线:
满足下列任一等式,圆锥曲线为正抛物线:
矩阵运算
把 (26) 写成如下矩阵运算式:
T T
1 x1 2 A B x1 D x1
+ +F =0
2 x2 B 2C x2 E x2
(31)
进一步写成:
1 T
x Qx + w T x + F = 0 (32)
2
其中,
2 A B D
Q= , w = E (33)
B 2C
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]
Page 21 | Chapter 9 深入圆锥曲线 | Book 3《数学要素》 | 鸢尾花书:从加减乘除到机器学习
正如牛顿所言,“我不知道世人看我的眼光。依我看来,我不过是一个在海边玩耍的孩子,不
时找到几个光滑卵石、漂亮贝壳,而惊喜万分。而展现在我面前的是,真理的浩瀚海洋,静候探
索。”
人类何尝不是在宇宙某个角落玩耍的一群孩子,手握的知识不过沧海一粟,却雄心万丈一心
要去探索星辰大海。
但也正是这群孩子将无数的不可能变成了可能,现在他们已经在地月系、甚至太阳系的边缘
跃跃欲试。
今人不见古时月,今月曾经照古人。宇宙的星辰大海一直都在人类眼前,它从未走远。路漫
漫其修远兮,吾将上下而求索。
地球不过是人类的摇篮,我们的征途是星辰大海。这句话含蓄而浪漫。刘慈欣《三体》中则
说的更为露骨而冷酷——“我们都是阴沟里的虫子,但总还是得有人仰望星空。”
本 PDF 文件为作者草稿,发布目的为方便读者在移动终端学习,终稿内容以清华大学出版社纸质出版物为准。
版权归清华大学出版社所有,请勿商用,引用请注明出处。
代码及 PDF 文件下载:https://github.com/Visualize-ML
本书配套微课视频均发布在 B 站——生姜 DrGinger:https://space.bilibili.com/513194466
欢迎大家批评指教,本书专属邮箱:[email protected]