以下是一个使用 `scipy` 库中的 `quad` 函数进行数值积分的 Python 程序。该程序会针对用户指定范围（如 0 到 1、1 到 2、...）对某个函数逐段求积，并输出结果。 ```python from scipy.integrate import quad # 定义要计算积分的目标函数 def target_function(x): return x**2 + 2 * x + 1 # 示例：目标函数为 f(x) = x^2 + 2x + 1 # 计算多个区间的积分值 def calculate_integrals(func, intervals): results = [] for interval in intervals: a, b = interval integral_result, _ = quad(func, a, b) results.append(integral_result) return results if __name__ == "__main__": # 指定区间列表 [(start_1, end_1), (start_2, end_2), ...] intervals = [(i, i+1) for i in range(5)] # 如 [0,1], [1,2], ..., [4,5] integrals = calculate_integrals(target_function, intervals) print("各个区间的积分结果:") for idx, result in enumerate(integrals): print(f"区间 [{idx}, {idx+1}] 上的积分值: {result:.6f}") ``` ### 给出解释： 1. **目标函数定义 (`target_function`)** - 我们定义了一个示例函数 \( f(x) = x^2 + 2x + 1 \)，你可以根据需求替换为目标函数。 2. **核心逻辑部分 (`calculate_integrals`)** - 使用了 SciPy 的 `quad()` 方法来进行数值积分操作。 - 遍历传入的每个区间 `(a,b)` 并调用 `quad` 对应地返回积分结果及误差估计值 `_` （我们只关心结果，所以忽略误差信息）。 3. **生成并处理多个连续区间** - 在主模块中利用列表推导式创建从 `[0,1]`, `[1,2]`, 直至所需上限的一系列子区间。 - 将这些区间传递给 `calculate_integrals` 来获取各片段上对应的积分值数组。 4. **打印最终结果** - 最后按顺序列出每一个分片内得到的具体积分答案供查阅分析之需。 #### 注意事项： - 如果需要改变积分上下限或调整其他参数，则可以直接修改 `intervals` 变量内容即可适配新的场景要求；同时确保导入必要的科学运算库安装完成后再执行脚本代码运行环境准备妥当。