python - 在python中生成对数刻度

Question

我需要生成一个列表，其中包含 10^a 和 10^b 之间的 logscale 刻度，并且 a < b，但是到目前为止我找不到任何方便的函数。所以我继续手动完成:

 # List containing 0.1, 0.2, 0.3, ..., 800, 900, 1000 
 x = np.arange(0.1, 1.0, 0.1).tolist() +
     np.arange(1, 10, 1).tolist() + 
     np.arange(10, 100, 10).tolist() + 
     np.arange(100, 1000, 100).tolist()

有没有一种方法可以做到这一点？

Answer 1

这可以通过 broadcasting 实现——

N = 4 # number of "levels"
out = ((10**np.arange(N)[:,None])*np.arange(0.1, 1.0, 0.1)).ravel()

另一个与 np.linspace ——

np.linspace([0.1,1,10,100],[1,10,100,1000],9, endpoint=False, axis=1).ravel()

用 N 概括——

s = [10**i for i in range(-1,N)]
np.linspace(s[:-1],s[1:],9, endpoint=False, axis=1).ravel()