SciPy 中 10j 的值是多少？

Question

我正在学习 Python 和 SciPy。我遇到了以下两种表达方式：

a = np.concatenate(([3], [0]*5, np.arange(-1, 1.002, 2/9.0)))

和

b = np.r_[3,[0]*5,-1:1:10j]

两个表达式输出相同的数组。我不明白第二个表达式中的 10j。它的价值是什么？非常感谢您的帮助。

Answer 1

这是创建np.linspace 的简写。

根据docs for np.r_：

如果使用切片表示法，语法start:stop:step 等价于括号内的np.arange(start, stop, step)。但是，如果 step 是一个虚数（即100j），则它的整数部分被解释为所需的点数，并且包含 start 和 stop。换句话说，start:stop:stepj 被解释为括号内的np.linspace(start, stop, step, endpoint=1)。

因此，对于您的具体情况，-1:1:10j 将导致 (1 - (-1)) / 9 = 0.222222... 的步长大小为以下数组：

>>> np.r_[-1:1:10j]
array([-1.        , -0.77777778, -0.55555556, -0.33333333, -0.11111111,
        0.11111111,  0.33333333,  0.55555556,  0.77777778,  1.        ])

虽然这发生给您的答案与np.arange(-1, 1.002, 2/9.0) 相同，但请注意arange 通常不是创建这样一个数组的好方法，因为在aranges 中使用非整数步长是bad idea：

使用非整数步长（例如 0.1）时，结果通常会不一致。对于这些情况，最好使用 linspace。

Answer 2

给你：引用https://docs.scipy.org/doc/numpy/reference/generated/numpy.r_.html

但是，如果 step 是一个虚数（即 100j），则它的整数部分被解释为所需的点数，并且包含 start 和 stop。