python数据分析,难以理解食谱代码

时间:2012-10-26 17:20:39

标签: python numpy matplotlib scipy

我正在学习使用python进行数据分析等等。我对scipy cookbook中此代码中发生的事情感到有些困惑。

当食谱描述整合然后绘图过程时,通过matplotlib,它首先有:

t = linspace(0, 15, 1000)
X0 = array([10,5])
X = scipy.integrate.odeint(dX_dt, X0, t)
rabbits, foxes = X.T

这段代码有什么作用?

rabbits, foxes = X.T

具体来说,X.T做了什么?

2 个答案:

答案 0 :(得分:9)

X.TX的转置。因此,在该行中,X必须是形状为(N,2)的数组。当你转置它时,你得到一个形状(2,N)的数组,可以解压缩。

考虑:

>>> import numpy as np
>>> a = np.arange(10).reshape((5,2))
>>> a
array([[0, 1],
       [2, 3],
       [4, 5],
       [6, 7],
       [8, 9]])
>>> B,C = a
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
ValueError: too many values to unpack
>>> B,C = a.T
>>> B
array([0, 2, 4, 6, 8])
>>> C
array([1, 3, 5, 7, 9])

另请注意尽可能,转置将返回一个新视图(数据不会被复制),因此这是一个非常有效的操作。

答案 1 :(得分:0)

代表转置。因此,X被对齐,以便将兔子分配给第一个元素,将狐狸分配给第二个元素。这些是数组(相当于线性代数中的矩阵)而不是列表,因此对齐确实很重要。