如何为对数图创建每十年的点间距

Question

into

我有一些频率从freq_start_hz = X到freq_stop_hz = Y。

我试图基于每十年的许多样本(num_samp_per_decade(，在范围[freq_start_hz, freq_stop_hz]之间进行对数(基数10(，包括端点。

我注意到numpy具有logspace方法(链接(，它使您能够根据示例的总数num创建某个范围base ** start到base ** stop的对数分隔。

您可以帮助我创建python代码，该代码每十年都会创建甚至对数间距？

示例

freq_start_hz = 10，freq_stop_hz = 100，num_samp_per_decade = 5

这很容易，因为它只有十年。因此，可以使用以下内容来创建它：

import numpy as np
from math import log10

freq_start_hz = 10
freq_stop_hz = 100
num_samp_per_decade = 5
freq_list = np.logspace(
    start=log10(freq_start_hz), 
    stop=log10(freq_stop_hz), 
    num=num_samp_per_decade, 
    endpoint=False,
    base=10,
)
freq_list = np.append(freq_list, freq_stop_hz)  # Appending end
print(freq_list.tolist())

输出是[10.0, 17.78279410038923, 31.622776601683793, 56.23413251903491, 100.0]

注意：这很好，因为我是这样设计的。如果freq_start_hz = 8，此方法现在不再起作用，因为它现在跨越了数十年。

结论

我希望在某个地方，math，numpy，另一个scipy库中有一个预制方法或我的互联网搜索未出现的其他库。

Answer 1

根据数十年的数量计算点数。

from math import log10
import numpy as np
start = 10
end = 1500
samples_per_decade = 5
ndecades = log10(end) - log10(start)
npoints = int(ndecades) * samples_per_decade
#a = np.linspace(log10(start), log10(end), num = npoints)
#points = np.power(10, a)
points = np.logspace(log10(start), log10(end), num=npoints, endpoint=True, base=10)
print(points)