使用Matplotlib绘制极坐标方程r=2sinθ+2cosθ时未生成预期圆形的技术问题

嘿，这个问题我太熟悉了！我来帮你拆解清楚为什么会出现这种“心形线”的错觉，以及怎么修复它。

首先，先确认你的数学推导没错：方程 r = 2sinθ + 2cosθ 确实对应一个圆。我们把它转换成直角坐标就能看明白：
两边同乘r得 r² = 2rsinθ + 2rcosθ，代入 r² = x²+y²、rsinθ=y、rcosθ=x，整理后得到：
(x-1)² + (y-1)² = 2
这是一个圆心在(1,1)、半径√2的圆，和Wolfram Alpha的结果一致。

那为什么你的Matplotlib代码画出了心形线？问题出在Matplotlib处理极坐标负r的方式，以及绘制顺序：
Matplotlib的极坐标绘图有个默认行为：当r < 0时，它会自动把角度θ加上π，然后用|r|作为距离来绘制这个点。从极坐标的定义来说，(r, θ)和(|r|, θ+π)是同一个点，但问题在于，你的代码是按θ从0到2π递增的顺序绘制的——当r由正变负时，Matplotlib会突然把点的位置跳到相反的角度，导致连线直接穿过圆心，形成了类似心形的交叉轨迹。

举个具体的例子：当θ接近3π/4时，r趋近于0；当θ超过3π/4后，r变成负数，Matplotlib会把这部分的θ加上π，也就是映射到7π/4到11π/4（即7π/4到2π再到3π/4）的区间，这时候绘制的点会从圆的右上角突然跳到左下角，然后继续连线，看起来就像心形线了。

你尝试用abs(r)的思路是对的！当你把r取绝对值后，所有r都是非负的，Matplotlib不会再调整角度，每个θ对应的点都在圆的正确位置上，而且绘制顺序是连续的，这样就能画出完整的圆。你说“上半部分符合预期”可能是个小误解，其实完整运行abs(r)的代码应该能得到正确的圆形。

下面给你几种可行的修复方案：

方案一：确保r非负（最简单） #

直接对r取绝对值，让Matplotlib按正常顺序绘制：

import numpy as np
from matplotlib import pyplot as plt

theta = np.arange(0, 2*np.pi, 0.01)
# 取绝对值确保r始终为正
r = abs(2 * np.cos(theta) + 2 * np.sin(theta))

fig = plt.figure()
ax = fig.add_subplot(polar=True)
ax.plot(theta, r)
plt.show()

方案二：拆分theta区间，避免跳跃 #

把theta分成r≥0和r<0两部分，分别绘制，这样每部分的连线都是连续的：

import numpy as np
from matplotlib import pyplot as plt

theta = np.arange(0, 2*np.pi, 0.01)
r = 2 * np.cos(theta) + 2 * np.sin(theta)

# 拆分r正和r负的区间
mask_pos = r >= 0
mask_neg = r < 0

fig = plt.figure()
ax = fig.add_subplot(polar=True)
# 分别绘制两部分
ax.plot(theta[mask_pos], r[mask_pos])
ax.plot(theta[mask_neg], r[mask_neg])
plt.show()

方案三：直接用直角坐标绘图（最直观） #

既然我们已经知道这是个圆，也可以直接用直角坐标绘制，完全避开极坐标的坑：

import numpy as np
from matplotlib import pyplot as plt

# 生成圆上的点
theta_circle = np.arange(0, 2*np.pi, 0.01)
x = 1 + np.sqrt(2)*np.cos(theta_circle)
y = 1 + np.sqrt(2)*np.sin(theta_circle)

fig = plt.figure()
ax = fig.add_subplot(aspect='equal')
ax.plot(x, y)
plt.show()

这样三种方法都能得到你预期的圆形啦！

内容的提问来源于stack exchange，提问作者matrucious