Mathematica中如何遍历p、K数据集计算表达式得到sigma序列

步骤1：数据预处理 #

首先需定义公式中的常量R，同时将导入的数据集转为一维数值列表方便遍历，同时校验两个数据集长度是否匹配：

(* 替换为你实际使用的R的数值 *)
R = 1; 
(* 导入CSV数据并转为一维数值列表，如有表头可添加HeaderLines->1参数跳过 *)
p = Flatten@Import["C:\\Users\\Maryestein\\onlyp.csv", "Data"];
K = Flatten@Import["C:\\Users\\Maryestein\\onlyk.csv", "Data"];
(* 校验数据集长度是否一致，不匹配则终止运行 *)
If[Length[p] != Length[K], Print["错误：p和K数据集长度不匹配"]; Abort[]]

步骤2：遍历计算σ结果列表 #

推荐使用逻辑更直观的Table实现循环遍历，Mathematica的列表下标从1开始，K[[i]]、p[[i]]分别对应公式中的K_i、p_i。另外公式中i=1时的K_{i-1}默认设为0，你可以根据需求自行调整：

(* 获取数据集总长度 *)
n = Length[K];
(* 遍历计算所有σ平方的结果列表 *)
sigmaSqList = Table[
   (* 计算相邻K的差值，i=1时默认K0=0 *)
   deltaK = If[i == 1, K[[i]] - 0, K[[i]] - K[[i - 1]]];
   9/(2*Pi^2) * (Sin[K[[i]]*R] - K[[i]]*R*Cos[K[[i]]*R])^2 / (K[[i]]*R)^6 * p[[i]] * K[[i]]^2 * deltaK,
   {i, 1, n}
];
(* 如需得到σ的结果，直接对每个元素开平方即可 *)
sigmaList = Sqrt /@ sigmaSqList;

结果导出（可选） #

如果需要将计算结果保存为CSV文件，可使用以下代码：

Export["C:\\Users\\Maryestein\\sigma_result.csv", sigmaSqList, "CSV"]

注意事项 #

• 请提前将代码中的R替换为你实际使用的常量数值
• 若导入的CSV文件带有表头，可在Import中添加HeaderLines->1参数跳过表头行，示例：Import["xxx.csv", "Data", HeaderLines->1]
• 若计算出现无穷大结果，请先检查K数据集中是否存在等于0的项，公式分母包含K_i*R，需提前过滤K_i为0的行

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