python数据分析，难懂cookbook代码答案

【问题标题】：python data analysis, difficulty understanding cookbook codepython数据分析，难懂cookbook代码
【发布时间】：2012-10-27 01:26:49
【问题描述】：

我正在学习使用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 是做什么的？

【问题讨论】：

我猜odeint 函数将返回一个具有T 属性的对象，该属性是一个touple 或2 个成员的列表，即(rabbits, foxes)，但我不知道更多细节
这是一个矩阵转置操作，见here，更多细节见here
交互式地，您可以随时输入help(X.T)，以显示文档。

标签： python numpy matplotlib scipy

【解决方案1】：

X.T 是X 的转置。所以，在那一行中，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])

还要注意只要有可能，转置会返回一个新的视图（数据不会被复制），所以这是一个非常高效的操作。

【讨论】：

【解决方案2】：

T 代表转置。所以 X 是对齐的，兔子被分配到第一个元素，狐狸被分配到第二个元素。这些是数组（相当于线性代数中的矩阵）而不是列表，因此对齐确实很重要。

【讨论】：

我怀疑它返回一个矩阵——它可能返回一个ndarray，在某些情况下它的行为与矩阵不同（例如与另一个ndarray相乘）
我想说的是，这些都是特殊类型，这就是为什么 alginment 很重要。数组有对齐，列表没有。