请求解析MATLAB代码运算顺序并协助转换为Python

原MATLAB代码分步拆解 #

把核心代码拆成5个独立运算步骤：

1. 计算起始直径的常用对数：log10_first = log10(first_diameter)
2. 计算终止直径的常用对数：log10_last = log10(last_diameter)
3. 计算对数区间的步长：step = (log10_last - log10_first) / lkdiam
4. 生成对数尺度下的等间隔序列：log_seq = log10_first:step:log10_last（共74个元素，73个间隔+1个起点）
5. 将对数序列还原为原始尺度，再转置为列向量：logdp = (10.^log_seq)'

Python等价实现 #

使用numpy复刻该逻辑，代码如下：

import numpy as np

first_diameter = 5e-9
last_diameter = 15e-6
lkdiam = 73

# 分步实现
log10_first = np.log10(first_diameter)
log10_last = np.log10(last_diameter)
step = (log10_last - log10_first) / lkdiam
# 加step确保包含终点（规避浮点精度问题）
log_seq = np.arange(log10_first, log10_last + step, step)
# 转换为原始尺度并转成列向量
logdp = np.power(10, log_seq).reshape(-1, 1)

核心逻辑说明 #

这段代码是在对数尺度上生成等间隔点，再映射回原始尺度，最终得到的数值序列呈几何分布（相邻元素比值近似恒定），而非均匀的算术分布。

内容的提问来源于stack exchange，提问作者Ajgar Dev