如何在Python中基于一般式椭圆代码实现椭圆上任意点绘制？

问题描述 #

用户提供了一段绘制一般式旋转椭圆的Matplotlib代码，希望实现类似Geogebra中Point on Object的功能——能在椭圆上任意位置绘制或移动点。原代码修正缺失依赖后的版本如下：

import numpy as np
import matplotlib as mpl
import matplotlib.pyplot as plt

x_range = np.linspace(-6, +6)
y_range = np.linspace(-10, +1)
x, y = np.meshgrid(x_range, y_range)

A, B, C, D, E, F = 0.25, 0.40, 0.25, 1.50, 2.00, 1.50

f = lambda x, y, A, B, C, D, E, F: A*x**2 + B*x*y +C*y**2 + D*x + E*y + F
equation = f(x, y, A, B, C, D, E, F)

fig, ax = plt.subplots(1, 1, figsize = (4, 4), tight_layout = True)

ax.contour(x, y, equation, levels = [0], colors = 'blue')

plt.show()

实现方法 #

方法1：参数化椭圆方程（精准控点） #

将一般式椭圆转换为标准参数形式，通过角度参数直接生成椭圆上的任意点，还能轻松调整点的位置：

import numpy as np
import matplotlib.pyplot as plt
from matplotlib.patches import Ellipse

# 原椭圆参数
A, B, C, D, E, F = 0.25, 0.40, 0.25, 1.50, 2.00, 1.50

# 计算椭圆中心(h,k)
delta = B**2 - 4*A*C
h = (2*C*D - B*E)/delta
k = (2*A*E - B*D)/delta

# 计算椭圆旋转角度（转换为角度制）
theta_rad = 0.5 * np.arctan2(B, A - C)
theta_deg = np.degrees(theta_rad)

# 计算长半轴a和短半轴b
n = 2*(A*h**2 + C*k**2 + B*h*k - F)
a = np.sqrt(n / (2*(A*np.cos(theta_rad)**2 + C*np.sin(theta_rad)**2 + B*np.cos(theta_rad)*np.sin(theta_rad))))
b = np.sqrt(n / (2*(A*np.sin(theta_rad)**2 + C*np.cos(theta_rad)**2 - B*np.cos(theta_rad)*np.sin(theta_rad))))

# 生成椭圆上的目标点（修改theta_point即可移动点，范围0~2π）
theta_point = np.pi/3  # 示例：60度位置的点
x_point = h + a*np.cos(theta_point)*np.cos(theta_rad) - b*np.sin(theta_point)*np.sin(theta_rad)
y_point = k + a*np.cos(theta_point)*np.sin(theta_rad) + b*np.sin(theta_point)*np.cos(theta_rad)

# 绘制椭圆与点
fig, ax = plt.subplots(1, 1, figsize=(4,4), tight_layout=True)
ellipse = Ellipse((h,k), 2*a, 2*b, angle=theta_deg, edgecolor='blue', facecolor='none')
ax.add_patch(ellipse)
ax.scatter(x_point, y_point, color='red', s=50, zorder=5)

ax.set_xlim(-6,6)
ax.set_ylim(-10,1)
ax.set_aspect('equal')  # 保持坐标轴比例，避免椭圆变形
plt.show()

只需修改theta_point的数值，就能让点在椭圆上任意移动。

方法2：交互式点击添加点 #

通过Matplotlib的事件处理，实现类似Geogebra的手动点击添加点功能：

import numpy as np
import matplotlib.pyplot as plt

A, B, C, D, E, F = 0.25, 0.40, 0.25, 1.50, 2.00, 1.50
# 定义椭圆方程判断函数
is_on_ellipse = lambda x,y: abs(A*x**2 + B*x*y + C*y**2 + D*x + E*y + F) < 1e-3

fig, ax = plt.subplots(1,1, figsize=(4,4), tight_layout=True)
x_range = np.linspace(-6,6)
y_range = np.linspace(-10,1)
x,y = np.meshgrid(x_range, y_range)
ax.contour(x,y, A*x**2 + B*x*y + C*y**2 + D*x + E*y + F, levels=[0], colors='blue')
ax.set_aspect('equal')

# 点击事件处理
def onclick(event):
    if event.inaxes != ax:
        return
    x_click, y_click = event.xdata, event.ydata
    if is_on_ellipse(x_click, y_click):
        ax.scatter(x_click, y_click, color='red', s=50, zorder=5)
        fig.canvas.draw()

fig.canvas.mpl_connect('button_press_event', onclick)
plt.show()

运行后点击椭圆附近位置，只要误差在允许范围内，就会在椭圆上添加红点。

补充说明 #

• 参数化方法适合需要精准控制点位置、批量生成或动态移动点的场景；
• 交互式方法更贴近手动操作的体验，适合临时添加点；
• 必须设置ax.set_aspect('equal')，否则椭圆会因坐标轴比例失衡变形，影响点的位置判断。

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