极坐标下图形的面积（数学系常用表白曲线）

极坐标下图形的面积（数学系常用表白曲线）

小象实战讲义 · 人工智能数学基础

在学习了直角坐标系下用定积分求面积后，本节我们将目光投向另一种强大的坐标系——极坐标系。许多在直角坐标系下表达复杂的曲线（如三叶玫瑰线、心形线），在极坐标下却有着简洁优美的方程。掌握极坐标下面积的求法，不仅能拓宽我们处理几何问题的工具箱，也为后续学习向量分析、场论等知识打下基础。学完本节，你将能计算多种极坐标曲线的面积，并深刻理解微元法的普适性。

💡 核心导读

本节我们将沿着以下路径，掌握极坐标下求面积的核心思想与方法：

1. 极坐标回顾：理解极坐标 $(r, \theta)$ 如何唯一确定平面上的点，建立与直角坐标系的联系。
2. 面积公式推导：运用微元法，将曲边扇形面积近似为无数个小三角形面积之和，从而导出核心公式。
3. 经典案例解析：通过计算“三叶玫瑰线”和“心形线”的面积，实战演练公式应用，并学习利用对称性简化计算。
4. 计算验证：使用 Python 的 SymPy 库进行符号积分，验证心形线面积公式的正确性。

极坐标系回顾

我们熟悉的平面直角坐标系 $(x, y)$ 用一对垂直的坐标轴来定位点。极坐标系则采用另一种视角：用距离和角度来定位。

在平面上选定一个极点 $O$（相当于原点）和一条极轴（通常取为 $x$ 轴正方向）。对于平面上任意一点 $P$：

• $r$：表示点 $P$ 到极点 $O$ 的距离，$r \ge 0$。
• $\theta$：表示以极轴为始边，逆时针旋转到射线 $OP$ 所成的角。

这样，有序对 $(r, \theta)$ 就唯一确定了点 $P$ 的位置。从线性代数的观点看，直角坐标系和极坐标系只是为同一个二维空间选用了两组不同的“基”，而坐标就是向量在这组基下的系数。

在极坐标系中，一条曲线通常由方程 $r = r(\theta)$ 给出，其中 $\theta$ 在某个区间 $[\alpha, \beta]$ 内变化。

极坐标下面积的微元法推导

我们的目标是求由曲线 $r = r(\theta)$、射线 $\theta = \alpha$ 和 $\theta = \beta$ 所围成的曲边扇形的面积 $S$。

微积分的思想是“化整为零，再积零为整”。我们无法直接计算曲边扇形的面积，但可以将其无限细分。

1. 细分自变量：将角度区间 $[\alpha, \beta]$ 进行无限细分。考虑其中一个微小的角度区间 $[\theta, \theta + d\theta]$。
2. 寻找面积微元：对应于这个微小角度区间的是一块窄窄的曲边扇形。当 $d\theta$ 非常小时，这段曲线近似为直线，该曲边扇形可以近似看作一个圆心角为 $d\theta$，半径为 $r(\theta)$ 的扇形。
3. 计算微元面积：扇形的面积公式是 $\frac{1}{2} \times (\text{半径})^2 \times (\text{圆心角弧度})$。因此，面积微元 $dS$ 为： $$dS = \frac{1}{2} [r(\theta)]^2 d\theta$$
4. 积分求和：将所有这些从 $\theta = \alpha$ 到 $\theta = \beta$ 的微小面积 $dS$ “累加”起来，就得到了整个曲边扇形的面积，即对 $dS$ 求定积分：

$$S = \int_{\alpha}^{\beta} dS = \int_{\alpha}^{\beta} \frac{1}{2} [r(\theta)]^2 d\theta$$

这就是极坐标下平面图形面积的计算公式。其核心在于，用无数个顶角为 $d\theta$、腰长为 $r(\theta)$ 的微小等腰三角形面积之和，来逼近曲边扇形的面积。

经典案例解析

案例一：三叶玫瑰线

三叶玫瑰线的极坐标方程为 $r(\theta) = a \sin(3\theta)$，其中 $a > 0$。其图像关于原点对称，有三片形状相同的“叶子”。

问题：求一整条三叶玫瑰线的总面积。

分析与求解：

1. 利用对称性：由于三片叶子完全相同，我们只需计算其中一片的面积，再乘以 3 即可。进一步观察，一片叶子本身关于其中心线 $\theta = \frac{\pi}{6}$ 对称。因此，我们可以只计算从 $\theta = 0$ 到 $\theta = \frac{\pi}{6}$（此时 $r$ 从 0 增加到最大值 $a$，对应半片叶子的上升沿）这“半片”叶子的面积，再乘以 6。
2. 套用公式：设一片叶子的面积为 $S_{\text{leaf}}$，则总面积 $S = 3 S_{\text{leaf}} = 6 \times S_{\text{半叶}}$。 $$S_{\text{半叶}} = \frac{1}{2} \int_{0}^{\pi/6} [a \sin(3\theta)]^2 d\theta = \frac{a^2}{2} \int_{0}^{\pi/6} \sin^2(3\theta) d\theta$$
3. 计算积分：利用三角恒等式 $\sin^2 x = \frac{1 - \cos 2x}{2}$ 进行降幂。 $$\begin{aligned} \int_{0}^{\pi/6} \sin^2(3\theta) d\theta &= \int_{0}^{\pi/6} \frac{1 - \cos(6\theta)}{2} d\theta \ &= \frac{1}{2} \left[ \theta - \frac{1}{6}\sin(6\theta) \right]_{0}^{\pi/6} \ &= \frac{1}{2} \times \frac{\pi}{6} = \frac{\pi}{12} \end{aligned}$$
4. 得出结果： $$S_{\text{半叶}} = \frac{a^2}{2} \times \frac{\pi}{12} = \frac{\pi a^2}{24}$$ $$S = 6 \times S_{\text{半叶}} = 6 \times \frac{\pi a^2}{24} = \frac{\pi a^2}{4}$$

结论：三叶玫瑰线的总面积为 $\frac{\pi a^2}{4}$。

案例二：心形线（表白曲线）

心形线的极坐标方程为 $r(\theta) = a(1 - \sin\theta)$，其中 $a > 0$。这是一个关于 $y$ 轴对称的图形。

问题：求心形线的面积。

分析与求解：

1. 利用对称性：心形线关于 $y$ 轴对称。因此，我们可以计算右半部分（例如 $\theta$ 从 $-\frac{\pi}{2}$ 到 $\frac{\pi}{2}$）的面积，再乘以 2。
2. 套用公式： $$S = 2 \times \frac{1}{2} \int_{-\pi/2}^{\pi/2} [a(1 - \sin\theta)]^2 d\theta = a^2 \int_{-\pi/2}^{\pi/2} (1 - \sin\theta)^2 d\theta$$
3. 展开并化简：展开被积函数，并利用对称区间上奇函数的积分为零的性质。 $$\begin{aligned} (1 - \sin\theta)^2 &= 1 - 2\sin\theta + \sin^2\theta \ \int_{-\pi/2}^{\pi/2} (1 - \sin\theta)^2 d\theta &= \int_{-\pi/2}^{\pi/2} (1 + \sin^2\theta) d\theta - \underbrace{\int_{-\pi/2}^{\pi/2} 2\sin\theta d\theta}$${0 \text{ (奇函数)}} \ &= \int{-\pi/2}^{\pi/2} (1 + \sin^2\theta) d\theta \end{aligned}
4. 计算积分：再次使用 $\sin^2\theta = \frac{1 - \cos 2\theta}{2}$。 $$\begin{aligned} \int_{-\pi/2}^{\pi/2} (1 + \sin^2\theta) d\theta &= \int_{-\pi/2}^{\pi/2} \left(1 + \frac{1 - \cos 2\theta}{2}\right) d\theta \ &= \int_{-\pi/2}^{\pi/2} \left(\frac{3}{2} - \frac{1}{2}\cos 2\theta\right) d\theta \ &= \left[ \frac{3}{2}\theta - \frac{1}{4}\sin 2\theta \right]_{-\pi/2}^{\pi/2} \ &= \left( \frac{3\pi}{4} - 0 \right) - \left( -\frac{3\pi}{4} - 0 \right) = \frac{3\pi}{2} \end{aligned}$$
5. 得出结果： $$S = a^2 \times \frac{3\pi}{2} = \frac{3\pi a^2}{2}$$

结论：心形线的面积为 $\frac{3\pi a^2}{2}$。

📝 动手练一练

1. 双纽线面积：双纽线的极坐标方程为 $r^2 = a^2 \cos(2\theta)$。求其所围成图形的总面积。 提示：注意 $\cos(2\theta) \ge 0$ 时 $r$ 才有定义，并利用对称性。

2. 阿基米德螺线：阿基米德螺线的方程为 $r = a\theta$ ($a>0$, $\theta \ge 0$)。求该螺线第一圈 ($0 \le \theta \le 2\pi$) 与极轴所围图形的面积。

参考答案：

1. 双纽线关于极点和坐标轴都对称。由 $r^2 \ge 0$ 得 $\cos(2\theta) \ge 0$，解得 $\theta \in [-\frac{\pi}{4}, \frac{\pi}{4}] \cup [\frac{3\pi}{4}, \frac{5\pi}{4}]$。图形由两片相似的“叶子”组成。计算一片的面积： $$S_{\text{leaf}} = \frac{1}{2} \int_{-\pi/4}^{\pi/4} a^2 \cos(2\theta) d\theta = \frac{a^2}{2} \left[ \frac{1}{2}\sin(2\theta) \right]$${-\pi/4}^{\pi/4} = \frac{a^2}{2}Sleaf​=21​∫−π/4π/4​a2cos(2θ)dθ=2a2​[21​sin(2θ)]−π/4π/4​=2a2​ 因此总面积 $S = 2S${\text{leaf}} = a^2。

2. 直接套用公式： $$S = \frac{1}{2} \int_{0}^{2\pi} (a\theta)^2 d\theta = \frac{a^2}{2} \int_{0}^{2\pi} \theta^2 d\theta = \frac{a^2}{2} \times \frac{(2\pi)^3}{3} = \frac{4\pi^3 a^2}{3}$$

计算验证：使用 Python 验证心形线面积

我们可以使用 SymPy 符号计算库来验证心形线面积公式 $\frac{3\pi a^2}{2}$ 的正确性。

import sympy as sp
import numpy as np

# 定义符号
theta, a = sp.symbols('theta a', real=True, positive=True)

# 定义心形线极坐标方程 r = a*(1 - sin(theta))
r = a * (1 - sp.sin(theta))

# 定义面积微元表达式 (1/2)*r^2 d(theta)
area_expr = (1/sp.Integer(2)) * r**2

# 计算心形线的总面积（积分区间为 0 到 2π）
total_area_sym = sp.integrate(area_expr, (theta, 0, 2*sp.pi))

# 化简结果
total_area_simplified = sp.simplify(total_area_sym)
print(f"心形线 (r = a*(1-sinθ)) 的符号积分总面积: {total_area_simplified}")

# 与理论公式对比
theoretical_area = (3*sp.pi/2) * a**2
print(f"理论公式给出的总面积: {theoretical_area}")
print(f"两者是否相等: {total_area_simplified.equals(theoretical_area)}")

# 数值验证：取 a=2，分别用符号结果和数值积分计算
a_val = 2
# 1. 代入符号结果求值
area_from_sym = total_area_simplified.subs(a, a_val).evalf()
print(f"\n当 a = {a_val} 时，由符号积分得到的面积: {area_from_sym}")

# 2. 数值积分验证（使用梯形法则）
def r_func(theta_val):
    return a_val * (1 - np.sin(theta_val))

theta_vals = np.linspace(0, 2*np.pi, 10000) # 细分区间
r_vals = r_func(theta_vals)
area_element_vals = 0.5 * r_vals**2 # 面积微元值
area_numerical = np.trapezoid(area_element_vals, theta_vals) # 梯形积分
print(f"当 a = {a_val} 时，数值积分（梯形法则）得到的面积: {area_numerical}")
print(f"数值误差: {abs(area_from_sym - area_numerical)}")

运行上述代码，你将看到符号积分的结果直接输出为 $\frac{3\pi a^2}{2}$，与我们的推导完全一致。数值积分的结果在足够细分的区间下也与理论值高度吻合，验证了公式的正确性。

本章小结

本节我们深入探讨了在极坐标系下如何计算平面图形的面积，这是定积分几何应用的又一重要篇章。

要点回顾：

• 核心公式：由曲线 $r = r(\theta)$ 与射线 $\theta=\alpha$, $\theta=\beta$ 所围曲边扇形的面积为 $S = \frac{1}{2}\int_{\alpha}^{\beta} [r(\theta)]^2 d\theta$。
• 推导精髓：微元法。将图形无限细分为圆心角为 $d\theta$ 的微小扇形，其面积近似为 $\frac{1}{2}r^2 d\theta$，积分即得总面积。
• 解题关键：善于观察图形的对称性（如轴对称、中心对称、周期性），可以极大简化积分区间和计算过程。
• 微积分思想：无论是直角坐标还是极坐标，求面积、体积、弧长等问题的本质都是“寻找正确的微元，然后积分求和”。这是将复杂、不规则的整体问题转化为简单、规则的局部问题的通用方法论。

行动清单：

1. 推导一遍：合上讲义，尝试独立推导极坐标面积公式，理解用三角形近似扇形的几何意义。
2. 动手编程：运行或改写提供的 Python 代码，尝试验证三叶玫瑰线或其他极坐标曲线的面积，感受符号计算与数值计算的能力。
3. 拓展思考：极坐标方程 $r = \sin(2\theta)$ 和 $r = \sin(4\theta)$ 对应的图形分别是什么样子？你能根据对称性，猜想并尝试计算它们的面积吗？

— 小象教研组