我正在尝试理解并修改以下代码:
import numpy as np
def f(x,y):
return (1 - x / 2 + x**5 + y**3) * np.exp(-x**2 -y**2)
n = 256
x = np.linspace(-3, 3, n)
y = np.linspace(-3, 3, n)
X,Y = np.meshgrid(x, y)
pl.axes([0.025, 0.025, 0.95, 0.95])
pl.contourf(X, Y, f(X, Y), 8, alpha=.75, cmap=pl.cm.hot)
C = pl.contour(X, Y, f(X, Y), 8, colors='black', linewidth=.5)
pl.clabel(C, inline=1, fontsize=10)
pl.xticks(())
pl.yticks(())
pl.show()
这里我们有一组点(x,y)
和每个点的值,用f(x,y)
计算
现在,我在txt输出文件中以x;y;output
的形式提供了一组计算结果,例如我用csv
模块读取。关键是我不理解这里的数据类型,可能是meshgrid。假设每个点都是字典key
中的键FH_DICT
,因此FH_DICT[key]
将在上面的代码中扮演f(x,y)
的角色。但我不知道如何实现它,因为每个点的输出值不容易表达为数学函数。
感谢您的时间。