使用numpy.savetxt和numpy.loadtxt编写和读取复数

时间:2014-04-22 22:48:46

标签: python arrays numpy complex-numbers

我需要编写和读取复数。我想使用numpy.savetxtnumpy.loadtxt来执行此操作。由于我编写的代码相当大,我创建了一个测试文件来尝试编写和读取复数。

到目前为止,我已经能够使用numpy.savetxt编写复数。代码如下:

import numpy

d1 = -0.240921619563 - 0.0303165074169j
d2 = -0.340921619563 - 0.0403165074169j
d3 = -0.440921619563 - 0.0503165074169j
d4 = -0.540921619563 - 0.0603165074169j

array = numpy.array([d1, d2, d3, d4])

save = open("test.dat", "w")
numpy.savetxt(save, array.reshape(1, array.shape[0]), newline = "\r\n", fmt = "%.10f")

save.close()

这给出了以下输出:

 (-0.2409216196+-0.0303165074j)  (-0.3409216196+-0.0403165074j)  (-0.4409216196+-0.0503165074j)  (-0.5409216196+-0.0603165074j)

我现在想要做的就是读取/加载数据。我的脚本是:

import numpy

d = numpy.loadtxt("test.dat")

这段代码还不够,我目前无法加载数据。我的问题类似于this one。但是,通过+-手动替换-,我仍然无法加载数据。我认为解决方案在于dtype的{​​{1}}选项。我虽然无法弄明白。

非常感谢您的帮助!

1 个答案:

答案 0 :(得分:4)

感谢Warren Weckesser!你建议的链接对我帮助很大。我现在有两个工作脚本:一个用于使用numpy.savetxt编写复数,另一个用于使用numpy.loadtxt从文件中读取/加载复数。

以下列出了代码。

写作:

import numpy

d1 = -0.240921619563-0.0303165074169j
d2 = -0.340921619563-0.0403165074169j
d3 = -0.440921619563-0.0503165074169j
d4 = -0.540921619563-0.0603165074169j

array = numpy.array([d1, d2, d3, d4])

save = open("test.dat","w")
numpy.savetxt(save, array.reshape(1, array.shape[0]), newline = "\r\n", fmt = '%.4f%+.4fj '*4)

save.close()

读/加载:

import numpy

coeffs = numpy.loadtxt("test.dat", dtype = numpy.complex128)