Python实现多次Lotka-Volterra仿真批量生成规范命名PDF图表

核心实现逻辑 #

• 命名规则完全依托Python原生字符串格式化实现，无需手动做前导零判断：对循环的仿真序号直接采用{sim_round:02d}的零填充格式，自动生成01、02...10的标准化文件名，文件排序也能和仿真顺序严格对应，不会出现chart_1、chart_10排在最前的乱序问题。
• os模块仅做最必要的轻量操作，不引入冗余逻辑：
  • 用os.makedirs(输出目录, exist_ok=True)提前创建存放PDF的文件夹，省略手动建目录的操作，目标目录已存在时也不会抛出报错
  • 用os.path.join()拼接文件全路径，自动适配Windows、Linux/macOS的路径分隔符，避免手动拼接斜杠导致的跨平台路径兼容问题
• 每轮仿真的PDF生成完全独立，从结构上规避所有图堆进同一个文件的常见问题：
  • PdfPages实例必须放在单次循环内部，通过上下文管理器初始化，绝不放在循环外做全局复用
  • 每轮PDF写完后调用plt.close('all')清空所有matplotlib画布，避免多轮绘图时的内存泄漏、跨轮串图问题
• 整体无冗余IO开销：不需要跑完所有仿真再遍历文件列表批量重命名，循环过程中直接按既定规则生成目标文件，整体运行效率和单次仿真的开销几乎没有差别。

可直接复用的核心代码片段 #

import os
import matplotlib.pyplot as plt
from matplotlib.backends.backend_pdf import PdfPages

# 提前创建独立的输出目录，集中存放所有仿真PDF
os.makedirs("lotka_volterra_outputs", exist_ok=True)

# 循环执行10次独立仿真，序号从1开始匹配文件命名规则
for sim_round in range(1, 11):
    # 👇 以下为已实现的单次仿真逻辑，不需要做任何ODE层面的修改
    # 随机生成alpha、beta、delta、gamma参数
    # 调用scipy.integrate.odeint完成方程求解
    # 绘制种群时间变化曲线、捕食者-猎物种群相图
    # 👆 原有单次运行的代码保持不变即可

    # 逐轮生成独立PDF
    file_name = f"chart_{sim_round:02d}.pdf"
    full_save_path = os.path.join("lotka_volterra_outputs", file_name)
    # 每轮单独初始化PdfPages，自动完成文件打开、写入、关闭流程
    with PdfPages(full_save_path) as pdf:
        # 将当前轮生成的所有画布依次写入PDF
        for fig_id in plt.get_fignums():
            pdf.savefig(plt.figure(fig_id))
    
    # 清空当前轮所有画布，避免状态残留
    plt.close('all')

常见注意事项 #

• 切勿在循环外部初始化PdfPages实例并全局复用，这类写法会导致10次仿真的所有图表全部写入同一个PDF文件，完全无法实现独立单文件的需求
• 不要直接用原始数字序号拼接文件名（如chart_1.pdf、chart_2.pdf），这类命名在文件管理器中会出现chart_1、chart_10排在chart_2前面的排序乱序问题，带前导零的两位数字命名可以完全规避这类问题
• 每轮仿真写完PDF后记得清理matplotlib画布：如果持续缓存绘图对象不释放，连续运行多轮时会出现内存占用过高、跨轮绘图内容串线的问题

内容的提问来源于stack exchange，提问作者Gromulus-Romulus