如何在Matplotlib中绘制同时包含X、Y误差的组合误差带？

要解决单独绘制X、Y误差带叠加效果杂乱的问题，我们可以通过生成所有误差边界顶点的凸包，绘制出同时包含X、Y误差及二者组合的单一闭合误差区域。

实现步骤 #

• 将列表格式的x、y及误差带数据转换为NumPy数组，方便批量计算
• 生成每个数据点对应的四个误差边界顶点（左下、左上、右上、右下）
• 使用凸包算法计算所有顶点的外包络线，得到闭合的误差区域轮廓
• 填充该轮廓区域，同时保留原始数据曲线

完整代码 #

import numpy as np
import matplotlib.pyplot as plt
from scipy.spatial import ConvexHull

# 生成数据
x = np.arange(-30, 31)
np.random.seed(0)
y = np.random.uniform(1, 10, len(x))

# 计算误差带上下限
x_lower = x - 1
x_upper = x + 1
y_lower = y * 0.98
y_upper = y * 1.02

# 生成所有误差边界的顶点
points = np.vstack([
    np.column_stack((x_lower, y_lower)),  # 左下顶点
    np.column_stack((x_lower, y_upper)),  # 左上顶点
    np.column_stack((x_upper, y_upper)),  # 右上顶点
    np.column_stack((x_upper, y_lower))   # 右下顶点
])

# 计算凸包，获取误差区域的闭合轮廓
hull = ConvexHull(points)
hull_points = points[hull.vertices, :]

# 绘图
plt.figure(figsize=(10, 6))
# 绘制原始数据曲线
plt.plot(x, y, label='随机数据', color='blue')
# 填充凸包形成的误差区域
plt.fill(hull_points[:, 0], hull_points[:, 1], color='lightblue', alpha=0.5, label='联合误差带')

# 添加图表元素
plt.xlabel('x')
plt.ylabel('y')
plt.title('包含X、Y联合误差的随机数据')
plt.grid(True)
plt.legend()
plt.show()

说明 #

• 凸包算法会自动计算所有误差顶点的最小闭合区域，完美覆盖X、Y单独误差及二者组合的所有可能范围
• 相比于分开绘制两个误差带，单一的凸包填充区域视觉上更整洁，也更准确地表达了数据的整体不确定性
• 如果不需要严格的最小凸包，也可以直接按顺序连接左右上下边界形成闭合区域，但凸包方案的视觉效果更优

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