我正在学习使用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做了什么?
答案 0 :(得分:9)
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])
另请注意尽可能,转置将返回一个新视图(数据不会被复制),因此这是一个非常有效的操作。
答案 1 :(得分:0)
代表转置。因此,X被对齐,以便将兔子分配给第一个元素,将狐狸分配给第二个元素。这些是数组(相当于线性代数中的矩阵)而不是列表,因此对齐确实很重要。