如何将scipy.integrate.odeint生成的(t,y)数据按列导出至Excel？

嘿，这个问题我太懂了！你遇到的情况是因为直接把t和y数组传给pd.DataFrame时，它默认把每个数组当成了一行数据，而不是我们需要的一列。下面给你两种简单的解决方法，保证导出后t在第一列，y在第二列：

方法一：用numpy.column_stack合并数组 #

先把t和y合并成一个二维数组（每一行对应一组t_i, y_i），再创建DataFrame：

import numpy as np
import pandas as pd
from scipy.integrate import odeint

# 先模拟你的ODE求解过程（替换成你自己的代码）
def ode_model(y, t):
    return -y  # 示例ODE：dy/dt = -y

y_initial = 1.0
t = np.linspace(0, 10, 50)  # 时间点数组
y = odeint(ode_model, y_initial, t)  # 求解得到的y数组（二维）

# 关键步骤：合并t和y为二维数组，创建DataFrame
combined_data = np.column_stack((t, y))
df = pd.DataFrame(combined_data, columns=["时间t", "解y"])

# 导出到Excel，记得关闭索引（index=False）
df.to_excel("ode_result.xlsx", index=False)

方法二：用字典构造DataFrame #

直接用列名作为键，对应的数组作为值来创建DataFrame，注意如果y是odeint返回的二维数组，需要先展平成一维：

# 承接上面的t和y数据
y_flattened = y.flatten()  # 把二维的y转成一维数组
df = pd.DataFrame({
    "时间t": t,
    "解y": y_flattened
})

df.to_excel("ode_result.xlsx", index=False)

为什么原来的方法会出错？ #

如果你之前写的是pd.DataFrame([t, y])，Pandas会把整个t数组当成第一行的唯一元素，y数组当成第二行的唯一元素，所以导出后就是两行长长的字符串。而我们需要的是每一行对应一对t和y的数值，所以必须按列的方式组织数据~

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