绘制4类潮汐图表及计算潮汐比例时遇格式错误求助

问题原因 #

你触发的TypeError: unsupported format string passed to numpy.ndarray.__format__，是因为计算得到的(K1+O1)/(M2+S2)结果是numpy数组类型，而Python的格式化字符串无法直接处理数组对象，必须先将其转换为标量值（如float）或采用数组专用的格式化方式。

解决方案 #

根据你的数组情况，选择以下对应处理方式：

情况1：结果是单元素数组 #

如果计算出的比例是只有一个值的数组，直接提取标量即可：

# 假设ratio是单元素numpy数组
ratio = (K1 + O1) / (M2 + S2)
# 方法1：用item()提取标量
print(f"潮汐比例: {ratio.item():.2f}")
# 方法2：强制转换为float
print(f"潮汐比例: {float(ratio):.2f}")

情况2：结果是多元素数组 #

如果比例是包含多个值的数组（比如对应不同时间点或潮汐类型），可以选择遍历输出或用numpy的格式化工具：

# 遍历每个元素逐个输出
for r in ratio:
    print(f"潮汐比例: {r:.2f}")

# 或者用numpy.array2string统一格式化输出
print(f"潮汐比例: {np.array2string(ratio, formatter={'float_kind':lambda x: f'{x:.2f}'})}")

完整修正示例代码 #

import numpy as np
import matplotlib.pyplot as plt

# 生成48小时时间序列
t = np.linspace(0, 48, 1000)

# 4种潮汐类型的参数（振幅+相移）
tide_types = [
    {"name": "半日潮主导", "K1_amp": 1, "O1_amp": 1, "M2_amp": 3, "S2_amp": 2, "phase": 0},
    {"name": "全日潮主导", "K1_amp": 3, "O1_amp": 2, "M2_amp": 1, "S2_amp": 1, "phase": np.pi/4},
    {"name": "混合潮（半日倾向）", "K1_amp": 2, "O1_amp": 1.5, "M2_amp": 2.5, "S2_amp": 2, "phase": np.pi/2},
    {"name": "混合潮（全日倾向）", "K1_amp": 2.5, "O1_amp": 2, "M2_amp": 1.5, "S2_amp": 1, "phase": 3*np.pi/4},
]

# 绘制潮汐曲线并计算比例
plt.figure(figsize=(12, 8))
for idx, tide in enumerate(tide_types, 1):
    # 计算潮汐高度
    height = (
        tide["K1_amp"] * np.cos(2*np.pi*(1/24)*t + tide["phase"]) +
        tide["O1_amp"] * np.cos(2*np.pi*(1/25.82)*t + tide["phase"]) +
        tide["M2_amp"] * np.cos(2*np.pi*(2/24)*t + tide["phase"]) +
        tide["S2_amp"] * np.cos(2*np.pi*(2/25)*t + tide["phase"])
    )
    
    # 子图绘制
    plt.subplot(2, 2, idx)
    plt.plot(t, height)
    plt.title(tide["name"])
    plt.xlabel("时间(小时)")
    plt.ylabel("潮汐高度")
    
    # 计算并打印潮汐比例（已修正）
    ratio = (tide["K1_amp"] + tide["O1_amp"]) / (tide["M2_amp"] + tide["S2_amp"])
    print(f"{tide['name']}: (K1+O1)/(M2+S2) = {ratio:.2f}")

plt.tight_layout()
plt.show()

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