docs/math/coordinate.md
author: Ir1d, HeRaNO, Chrogeek, abc1763613206, mxdyzmx
在小学或初中已经学习过角的 静态定义：具有公共端点的两条射线组成的图形叫做角．
但是该定义将角度限制在了 $[0, 360^\circ]$，这给深入研究带来了一定的困难，还有其他的问题无法解释清，比如：旋转 $720^\circ$ 是什么意思？
在高中数学，讲了角的 动态定义：平面内一条射线绕其端点从一个位置旋转到另一个位置形成的图形叫做角．
开始的位置称为 始边，结束的位置称为 终边．并规定：
这样就把角的概念推向了 任意角．
???+ note "注意" 零角始边和终边重合，但始边和终边重合的角并不都是零角，如以 $360^\circ$ 为倍数的角．
实际应用中经常有角度到各种参数的转换，而使用弧度制描述角可以减少系数的使用．所以接下来，介绍 弧度制：
把长度等于半径长的弧所对的圆心角称为 $1$ 弧度的角，用符号 $\text{rad}$ 表示，读作：弧度．
根据前面的规定，正角的弧度为正，负角的弧度为负，零角的弧度为 $0$，如果半径为 $r$ 的圆的圆心角 $\alpha$ 所对弧长为 $l$，则：
$$ |\alpha|=\dfrac{l}{r} $$
利用这个公式还可以写出弧长和扇形面积公式，在此略过．
于是，$360^\circ$ 角的弧度为 $2\pi$，这样有了对应关系之后就可以进行角度值和弧度制的转化了：
$$ k \operatorname{rad} = \frac{\pi}{180^\circ} n^\circ $$
考虑一个角，将其终边再旋转一周，甚至多周，始边位置不动，那么终边位置永远是相同的，称这些角为终边位置相同的角．
与角 $\alpha$ 终边位置相同的角的集合很容易得出，为 ${\varphi \mid \varphi = \alpha + 2k\pi, k \in \mathbf{Z}}$．
可以理解为：给这个角的边不停加转一圈，终边位置不变．
???+ note "$\pi$ 和 $\tau$ 两个数学常数" 目前西方数学界有一些观点认为，「真正的圆周率」应为 $2\pi$，将这个值记为希腊字母 $\tau$．新圆周率的支持者们选择在 6 月 28 日庆祝「真正的」圆周率日．
比如，在弧度制下，一个周角是 $2\pi$，直接对 $2\pi$ 进行等分可以得到周角的等分．又例如，在复变函数中频繁出现 $2\pi$ 的组合，等等．
为了迎合中国各地区约定俗成的习惯，在 **OI Wiki**，采用参数 $\pi$ 表示圆周率．
???+ note "编程中圆周率的习惯写法"
在 C/C++ 语言中，一般取 $\pi$ 为 acos(-1)，只有这个值是最接近 $\pi$ 的浮点数．使用 acos(-1) 或者 4 * atan(1) 写出来的 $\pi$ 是 $3.14159265358979310000$．
采用其他值，例如 `acos(-1.0/2.0)`，`acos(1.0/2.0)`，`asin(1.0/2.0)` 等等，写出来的 $\pi$ 是 $3.14159265358979360000$，这就不是最接近 $\pi$ 的浮点数了．
如果你背得下来，也可以直接写 $3.1415926535897932$．
在同一个平面上互相垂直且有公共原点的两条数轴构成平面直角坐标系（Rectangular Coordinates）．
通常，两条数轴分别置于水平位置与垂直位置，取向右与向上的方向分别为两条数轴的正方向．水平的数轴叫做 $x$ 轴（x-axis）或横轴，垂直的数轴叫做 $y$ 轴（y-axis）或纵轴，$x$ 轴 $y$ 轴统称为坐标轴，它们的公共原点 $O$ 称为平面直角坐标系的原点（origin），以点 $O$ 为原点的平面直角坐标系记作平面直角坐标系 $xOy$．
$x$ 轴 $y$ 轴将坐标平面分成了四个象限（quadrant），右上方的部分叫做第一象限，其他三个部分按逆时针方向依次叫做第二象限、第三象限和第四象限．象限以数轴为界，横轴、纵轴上的点及原点不在任何一个象限内．一般情况下，$x$ 轴 $y$ 轴取相同的单位长度，但在特殊的情况下，也可以取不同的单位长度．
在平面直角坐标系中，对于平面上的任意一点，都有唯一的一个有序数对（即点的坐标（coordinates））与它对应；反过来，对于任意一个有序数对，都有平面上唯一的一点与它对应．
对于平面内任意一点 $C$，过点 $C$ 分别向 $x$ 轴、$y$ 轴作垂线，垂足在 $x$ 轴、$y$ 轴上的对应点 $a, b$ 分别叫做点 $C$ 的横坐标、纵坐标，有序数对（ordered pair）$(a, b)$ 叫做点 $C$ 的直角坐标．一个点在不同的象限或坐标轴上，其坐标都不一样．
考虑实际情况，比如航海，说「点 $B$ 在点 $A$ 的北偏东 $30^\circ$ 方向上，距离为 $100$ 米」，而不是「以 $A$ 为原点建立平面直角坐标系，$B(50,50\sqrt 3)$」．
这样：
就建立了 极坐标系．
设 $A$ 为平面上一点．
那么有序数对 $(\rho,\varphi)$ 即为 $A$ 的 极坐标．
由终边相同的角的定义可知，$(\rho,\varphi)$ 与 $(\rho,\varphi + 2k\pi)\ (k\in \mathbf{Z})$ 其实表示的是一样的点．特别地，极点的极坐标为 $(0,\varphi)\ (\varphi \in \mathbf{R})$，于是平面内的点的极坐标表示有无数多种．
如果规定 $\rho \ge 0,0 \le \varphi < 2\pi$，那么除极点外，其他平面内的点可以用唯一有序数对 $(\rho,\varphi)$ 表示，而极坐标 $(\rho,\varphi)$ 表示的点是唯一确定的．
当然，有时候研究极坐标系下的图形有些不方便．要想转到直角坐标系下研究，有互化公式．点 $A(\rho,\varphi)$ 的直角坐标 $(x,y)$ 可以如下表示：
$$ \begin{aligned} x &= \rho \cos \varphi \ y &= \rho \sin \varphi \end{aligned} $$
进而可知：
$$ \begin{aligned} \rho^2 &= x^2 + y^2\ \tan \varphi &= \frac{y}{x}\ \ \ \ (x\not =0) \end{aligned} $$
于是有 $\rho = \sqrt{x^2+y^2}$．
但具有相同 $\dfrac{y}{x}$ 的 $\tan\varphi$ 有两个可能的 $\varphi$ 的值，此时还需要根据 $x, y$ 的值来确定方向．具体地，定义函数：
$$ \operatorname{atan2}(y, x) = \begin{cases} \arctan(\frac{y}{x}) & \text{if } x > 0 \ \arctan(\frac{y}{x}) + \pi & \text{if } y \ge 0, x < 0 \ \arctan(\frac{y}{x}) - \pi & \text{if } y < 0, x < 0 \ \pi/2 & \text{if } y > 0, x = 0 \ -\pi/2 & \text{if } y < 0, x = 0 \ \text{any} & \text{if } y = 0, x = 0 \end{cases} $$
则 $\varphi = \operatorname{atan2}(y, x)$．注意上述函数的值域为 $(-\pi, \pi]$．
在 C/C++ 语言的 <math.h> 或 <cmath> 库里定义了 该函数，调用 atan2(y, x) 即可．
使用如下方法建立空间直角坐标系：
这样就构成了一个空间直角坐标系，称为空间直角坐标系 $O-xyz$．定点 $O$ 称为该坐标系的原点．
任意两条坐标轴确定一个平面，这样可确定三个互相垂直的平面，统称为坐标面．其中 $x$ 轴与 $y$ 轴所确定的坐标面称为 $xOy$ 面，类似地有 $yOz$ 面和 $zOx$ 面．三个坐标面把空间分成八个部分，每一部分称为一个卦限．
取定空间直角坐标系 $O-xyz$ 后，就可以建立空间的点与三元组之间的一一对应关系．
设点 $M$ 为空间的一点，过点 $M$ 分别作垂直于 $x$ 轴、$y$ 轴和 $z$ 轴的平面．设三个平面与 $x$ 轴、$y$ 轴和 $z$ 轴的交点依次为 $P, Q, R$，点 $P, Q, R$ 分别称为点 $M$ 在 $x$ 轴、$y$ 轴和 $z$ 轴上的投影．又设点 $P, Q, R$ 在 $x$ 轴、$y$ 轴和 $z$ 轴上的坐标依次为 $x, y, z$，于是点 $M$ 确定了一个三元组 $(x, y, z)$．
反之，如果给定一个三元组 $(x, y, z)$，可以在 $x$ 轴上取坐标为 $x$ 的点 $P$，在 $y$ 轴上取坐标为 $y$ 的点 $Q$，在 $z$ 轴上取坐标为 $z$ 的点 $R$，然后点 $P, Q, R$ 分别作垂直于 $x$ 轴、$y$ 轴和 $z$ 轴的三个平面，它们相交于空间的一点 $M$，点 $M$ 就是由三元组 $(x, y, z)$ 所确定的点．
这样一来，空间的点 $M$ 与三元组 $(x, y, z)$ 之间就建立了一一对应的关系．把三元组 $(x, y, z)$ 称为点 $M$ 的坐标，记作 $M(x, y, z)$，其中 $x$ 称为横坐标、$y$ 称为纵坐标、$z$ 称为竖坐标．
空间柱坐标系，将极坐标扩展为三维的方式：从应用于平面工作中的极坐标系开始，然后过极点 $O$ 添加垂直于该平面的 $z$ 轴，方向朝上．
为了找到由柱坐标 $(\rho, \varphi, z)$ 所描述的点，可以首先在极坐标系下处理 $\rho$ 和 $\varphi$，然后根据 $z$ 坐标沿着 $z$ 轴「向上」或「向下」移动．
两坐标系下 $z$ 的值是相同的．
$(x,y)$ 与 $(\rho, \varphi)$ 的相互转换参见上文 平面直角坐标系与极坐标系的相互转换．
球坐标可以通过以下方法确定：
这样即可到达球坐标 $(r,\vartheta,\varphi)$ 所描述的点．其中 $\vartheta$ 称为 天顶角，$\varphi$ 称为 方位角．
???+ warning "Warning" 由于诸多原因，有的地方使用 $\phi$ 表示天顶角，用 $\theta$ 表示方位角．阅读文章遇到球坐标系时请务必注意这一点．
同时，在写文章时，如果用到了球坐标系，建议提前声明清楚使用什么符号表示天顶角和方位角．
两坐标系下 $\varphi$ 的值是相同的．
从柱坐标系到球坐标系：
$$ \begin{aligned} r &= \sqrt{\rho^2 + z^2} \ \vartheta &= \begin{cases} \arctan\left(\frac{\rho}{z}\right) & \text{if }z > 0 \ \pi/2 & \text{if }z = 0, \rho \not= 0 \ \arctan\left(\frac{\rho}{z}\right) + \pi & \text{if }z < 0 \ \end{cases} \end{aligned} $$
注意对于柱坐标系下的点 $(0,0,0)$，其球坐标的 $\vartheta$ 不明确．
从球坐标系到柱坐标系：
$$ \begin{aligned} \rho &= r \sin \vartheta \ z &= r \cos \vartheta \end{aligned} $$
可以结合上文 平面直角坐标系与极坐标系的相互转换 和上文 柱坐标系与球坐标系的相互转换 一起使用，或直接使用下面的公式：
从空间直角坐标系到球坐标系：
$$ \begin{aligned} r &= \sqrt{x^2 + y^2 + z^2} \ \vartheta &= \arccos\left(\frac{z}{\sqrt{x^2 + y^2 + z^2}}\right) \ \varphi &= \operatorname{atan2}(y, x) \end{aligned} $$
其中 $\operatorname{atan2}$ 的定义见 平面直角坐标系与极坐标系的相互转换．
注意对于空间直角坐标系下的点 $(0,0,0)$，其球坐标的 $\vartheta$ 和 $\varphi$ 取值不明确．
从球坐标系到空间直角坐标系：
$$ \begin{aligned} x &= r \sin \vartheta \cos \varphi \ y &= r \sin \vartheta \sin \varphi \ z &= r \cos \vartheta \end{aligned} $$