07 伪彩图 pcolor

07　伪彩图 pcolor

用途

等距网格的标量场着色，与 pcolormesh 同源，网格等距时更省事。

示例（5 张）

07-1　开关柜温度场的伪彩图

07-2　接地网地表电位伪彩图

07-3　供电分区负荷密度伪彩图

07-4　风电场风功率密度伪彩图

07-5　电缆终端电场强度伪彩图

代码

以下 5 个函数取自 code/第一批_分布与矩阵261107.py，直接调用即可出图。样式（宋体 + Times、10 磅、刻度朝内四边、留白自检）由 code/公共配置261107.py 的 setup_cn_style() 和 save() 统一提供，本节不再重复。

# 摘自 code/第一批_分布与矩阵261107.py

def pc_01():
    """开关柜温度场的伪彩图。"""
    fig, ax = new_fig(left=0.19, bottom=0.21, right=0.86)
    x = np.arange(0, 0.61, 0.03)
    y = np.arange(0, 0.91, 0.03)
    X, Y = np.meshgrid(x, y)
    Z = 30 + 26 * np.exp(-((X - 0.3) / 0.16) ** 2 - ((Y - 0.45) / 0.15) ** 2) + 5 * Y
    m = ax.pcolor(X, Y, Z, cmap=CMAP_SEQ, shading="auto")
    ax.set_xlabel("柜宽 " + IT["x"] + " / m")
    ax.set_ylabel("柜高 " + IT["y"] + " / m")
    cbar(fig, ax, m, "温度 " + IT["T"] + " / °C")
    save(fig, "07_伪彩图", "07_伪彩图_开关柜温度场")

def pc_02():
    """接地网地表电位伪彩图。"""
    fig, ax = new_fig(left=0.19, bottom=0.21, right=0.86)
    x = np.arange(-30, 31, 2.0)
    y = np.arange(-30, 31, 2.0)
    X, Y = np.meshgrid(x, y)
    V = np.zeros_like(X)
    for gx, gy, k in [(-8, -6, 1.0), (8, -6, 0.85), (0, 8, 0.92)]:
        V += k / np.sqrt((X - gx) ** 2 + (Y - gy) ** 2 + 1.2) * 62
    m = ax.pcolor(X, Y, V, cmap=CMAP_SEQ, shading="auto")
    ax.set_xlabel("东向距离 " + IT["x"] + " / m")
    ax.set_ylabel("北向距离 " + IT["y"] + " / m")
    ax.set_aspect("equal")
    cbar(fig, ax, m, "地表电位 " + IT["U"] + " / V")
    save(fig, "07_伪彩图", "07_伪彩图_接地网地表电位")

def pc_03():
    """供电分区负荷密度伪彩图。"""
    fig, ax = new_fig(left=0.19, bottom=0.21, right=0.86)
    rng = np.random.default_rng(703)
    x = np.arange(0, 12.5, 0.5)
    y = np.arange(0, 9.5, 0.5)
    X, Y = np.meshgrid(x, y)
    Z = (1200 * np.exp(-((X - 4.2) / 2.2) ** 2 - ((Y - 4.6) / 1.8) ** 2)
         + 900 * np.exp(-((X - 9.0) / 1.6) ** 2 - ((Y - 3.0) / 1.4) ** 2))
    m = ax.pcolor(X, Y, Z, cmap=CMAP_SEQ, shading="auto")
    ax.set_xlabel("东向距离 " + IT["x"] + " / km")
    ax.set_ylabel("北向距离 " + IT["y"] + " / km")
    ax.set_aspect("equal")
    cbar(fig, ax, m, "负荷密度 / (kW/km²)")
    save(fig, "07_伪彩图", "07_伪彩图_供电分区负荷密度")

def pc_04():
    """风电场风功率密度伪彩图。"""
    fig, ax = new_fig(left=0.19, bottom=0.21, right=0.86)
    x = np.arange(0, 6.1, 0.25)
    y = np.arange(0, 4.1, 0.25)
    X, Y = np.meshgrid(x, y)
    Z = (320 + 180 * np.exp(-((X - 2.0) / 1.4) ** 2 - ((Y - 2.4) / 1.0) ** 2)
         + 140 * np.exp(-((X - 4.6) / 1.0) ** 2 - ((Y - 1.3) / 0.9) ** 2))
    m = ax.pcolor(X, Y, Z, cmap=CMAP_SEQ, shading="auto")
    ax.set_xlabel("东向距离 " + IT["x"] + " / km")
    ax.set_ylabel("北向距离 " + IT["y"] + " / km")
    ax.set_aspect("equal")
    cbar(fig, ax, m, "风功率密度 / (W/m²)")
    save(fig, "07_伪彩图", "07_伪彩图_风功率密度")

def pc_05():
    """电缆终端电场强度伪彩图。"""
    fig, ax = new_fig(left=0.19, bottom=0.21, right=0.86)
    x = np.arange(-40, 41, 2.0)
    y = np.arange(0, 121, 3.0)
    X, Y = np.meshgrid(x, y)
    E = (42 * np.exp(-(np.abs(X) / 14) ** 1.3) * (1 - 0.004 * Y)
         + 12 * np.exp(-((np.abs(X) - 22) / 7) ** 2) * np.exp(-Y / 60))
    m = ax.pcolor(X, Y, E, cmap=CMAP_SEQ, shading="auto")
    ax.set_xlabel("轴向位置 " + IT["x"] + " / mm")
    ax.set_ylabel("径向位置 " + IT["y"] + " / mm")
    cbar(fig, ax, m, "场强 " + IT["E"] + " / (kV/mm)")
    save(fig, "07_伪彩图", "07_伪彩图_电缆终端电场")

说明

• 图片尺寸 1181×944 px（100 mm × 80 mm @ 300 dpi），三维图为 1181×1122 px。
• 图内不放标题，图题由本文件的二级/三级标题承担。
• 数据为贴近广东实际量级的仿真算例，随机种子固定，重跑结果一致。