2

我在 MATLAB 中实现了对 216x216 矩阵进行操作的代码，这些矩阵包含数字数据和有时是字符串。我对这些矩阵所做的操作大多类似于过滤某个阈值以上的矩阵，找到所有高于某个值的矩阵索引，找到一个高于 X 的值列表，然后找到它们之间的连续差异，一些字符串替换操作。做矩阵点积等。我需要访问数千个文件来生成这些矩阵（我在 MATLAB 中使用的 dlmread）。

现在我需要用通常与操作系统捆绑在一起的任何其他语言来实现上述项目，比如 Perl、c 或 python 或开源语言。

我做了一个简短的搜索，发现 python 是一个很好的研究工具。python 是否有一些用于矩阵运算的 MATLAB 等价物（例如将文件直接读入数组、查找、dlmwrite 等）

因为我的代码已经有很多没有这些 MATLAB 函数的循环，所以代码会变得更加混乱且难以维护。

或者你能指出任何其他的选择。我熟悉小 Perl，但不熟悉 python 或 R。

4

2 回答 2

6

从这个页面开始比较 NumPy 和 Matlab。

以下是有关您的帖子的一些示例：

In [5]: import scipy

In [6]: X = scipy.randn(3,3)

In [7]: X
Out[7]: 
array([[-1.16525755,  0.04875437, -0.91006082],
       [ 0.00703527,  0.21585977,  0.75102583],
       [ 1.12739755,  1.12907917, -2.02611163]])

In [8]: X>0
Out[8]: 
array([[False,  True, False],
       [ True,  True,  True],
       [ True,  True, False]], dtype=bool)

In [9]: scipy.where(X>0)
Out[9]: (array([0, 1, 1, 1, 2, 2]), array([1, 0, 1, 2, 0, 1]))

In [10]: X[X>0] = 99

In [11]: X
Out[11]: 
array([[ -1.16525755,  99.        ,  -0.91006082],
       [ 99.        ,  99.        ,  99.        ],
       [ 99.        ,  99.        ,  -2.02611163]])

In [12]: Y = scipy.randn(3,2)

In [13]: scipy.dot(X, Y)
Out[13]: 
array([[-124.41803568,  118.42995937],
       [-368.08354405,  199.67131528],
       [-190.13730231,  161.54715769]])

（无耻插件：我曾经在 Python 和 Matlab 之间做的比较。）

于 2012-05-25T03:12:50.407 回答
3

你可能想看看NumPy / SciPy，它们可能会帮助你做你想做的事。此外，还有大量用户，这将使在需要时更容易获得帮助。两个库的一般介绍在这里

于 2012-05-25T02:59:35.857 回答