使用gnuplot绘制黑体函数（普朗克曲线）

我试图用gnuplot绘制一条标准黑体曲线（普朗克函数）和维恩近似来匹配一组给定的恒星数据

数据在0到1.4个单位的（任意）y范围内精细绘制。然而，无论我做什么，普朗克曲线和维恩曲线几乎处处为0–所有点都在1e-27范围内

我是否定义了错误的函数，或者gnuplot是否在我不知道的情况下将非常小的数字设置为零

下面是我正在使用的gnuplot代码：

设置标题“Planck和Wien函数（黑体曲线）”字体“Helvetica，14”\
设置“ps3g0ivplot.eps”
设置xrange[“0”：“1.1e4”]\
#设置Y范围[“0”：“1.4”]\
设置标签“任意流量”字体“Helvetica，14”\
设置xlabel“波长（埃）”字体“Helvetica，14”\
设置mxtics\
设置粘液\
#常数
h=6.62606876e-34#普朗克常数
c=2.99792458e8#光速
k=1.3806503e-23#玻尔兹曼常数
普朗克（x，t）=\
（2*h*c*c/x**5）*（exp（-1*h*c/（x*k*t））-1）
维也纳（x，t）=\
（2*h*c**2/x**5）*exp（-1*h*c）/（x*k*t））
#设定温度，t
t=6400开氏度，与ps3g0iv.dat中给定的星相匹配
#普朗克印刷机（4500吨）
#印刷维恩（4500吨）
密谋\
Wien（x，t）lt 2 lw 3标题“Wien”\
普朗克（x，t）lt 1 lw 3标题“普朗克”\
#“ps3g0iv.dat”lt 9 pt 1 pointsize 1标题“ps3g0iv.dat”带行
退出

如果x的单位为埃，则必须将x替换为

x*1e-10

才能得到正确的公式。只需提示一下，您不需要；\谢谢你，克里斯托夫，我会试试的。感谢Jonatan提供的风格提示，这些都是前一段时间需要的东西的残羹剩饭，但我把它们拿出来了。普朗克函数是错误的，应该是

（2*h*c**2/x**5）/（exp（h*c/（x*k*t））-1）

。x的缩放只是x范围和grange选择的问题。谢谢Christoph，我想当你试图在凌晨2点完成作业时就会发生这种情况。我想也需要一些缩放。我将在一点时间内尝试一下。