计算饱和水汽密度函数并绘图时遇AttributeError的解决方法

先梳理可运行的基础代码框架 #

先把你的代码补全为标准结构，方便排查问题：

# 定义饱和水汽压计算系数
c0 = -2.915127
c2 = 0.001124901
c4 = -0.000000158311
c6 = 0.0000000000123218

# 饱和水汽压计算函数
def psat(temp):
    return c0 + c2*(temp + 100)**2 + c4*(temp + 100)**4 + c6*(temp + 100)**6

# 生成-10到40℃的温度序列
import numpy as np
temps = np.arange(-10, 41, 1)
psat_values = psat(temps)

常见AttributeError原因及对应解决办法 #

• 绘图库调用错误
  比如用matplotlib时，没正确导入pyplot模块就直接调用plt.plot，或者导入方式错误。
  解决：确保导入和调用步骤正确：

  import matplotlib.pyplot as plt
  plt.plot(temps, psat_values)
  plt.xlabel('温度(℃)')
  plt.ylabel('饱和水汽压')
  plt.title('温度-饱和水汽压关系')
  plt.show()
  

  注意：不能只导入matplotlib就直接用matplotlib.plot，必须调用pyplot子模块。

• 变量名冲突
  如果代码中把psat、plt或np这类标识符重新赋值成了普通变量（比如psat = 1013），会导致调用时找不到函数/模块的属性。
  解决：检查代码，改掉所有和库名、函数名冲突的变量名。

• 对非对象错误调用属性
  比如直接对numpy数组调用plot方法（psat_values.plot()），但numpy数组本身没有plot属性，必须通过matplotlib的函数来绘图。
  解决：始终用plt.plot(自变量, 因变量)的格式调用绘图方法，不要直接给数组加属性调用。

• 温度序列类型错误
  如果用普通Python列表（比如temps = list(range(-10,41))）传入函数，会导致矢量化运算失败，进而引发连锁报错。
  解决：用numpy.arange生成numpy数组类型的温度序列，确保函数能正常做批量计算。

额外提示：你的函数计算的是水汽压，不是密度 #

如果目标是饱和水汽密度，需要用理想气体状态方程转换：

# 水汽摩尔质量(kg/mol)、气体常数(J/(mol·K))
M = 0.01801528
R = 8.314
# 转换为热力学温度，同时将水汽压单位转为Pa（假设原psat输出为hPa）
rho = (psat_values * 100) * M / (R * (temps + 273.15))
# 绘图
plt.plot(temps, rho)
plt.ylabel('饱和水汽密度(kg/m³)')
plt.show()

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