如何使用Sympy存储方程，以便后续代入参数进行数值求解？

你的思路完全可行，只需要修正两处细节就能实现需求：

一、修正方程保存的代码 #

Pickle需要二进制模式读写文件，原代码用文本模式打开文件会报错，修改后保存方程的代码如下（Jupyter Notebook中）：

import sympy as sp
import pickle

# 定义符号变量
a, b, c = sp.symbols("a b c")
# 创建并简化方程（示例方程）
eqn = sp.Eq(b + c, a)

# 用二进制写入模式保存方程
with open("eqn.pkl", "wb") as f:
    pickle.dump(eqn, f)

注意：建议用.pkl作为pickle文件的后缀，更易识别文件类型。

二、加载方程并代入参数求解 #

加载方程后，不能直接给变量赋值（比如b=1）——这会覆盖SymPy的符号变量，需要用SymPy的subs()方法替换符号，或者在solve()中传入参数字典：

方法1：先替换参数再求解 #

import sympy as sp
import pickle

# 加载保存的方程
with open("eqn.pkl", "rb") as f:
    eqn = pickle.load(f)

# 定义参数值（键为SymPy符号对象）
params = {sp.symbols('b'): 1, sp.symbols('c'): 2}
# 替换方程中的参数
eqn_substituted = eqn.subs(params)
# 求解目标变量a
result = sp.solve(eqn_substituted, sp.symbols('a'))
print(result)  # 输出: [3]

方法2：直接在solve中传入参数 #

import sympy as sp
import pickle

with open("eqn.pkl", "rb") as f:
    eqn = pickle.load(f)

# 直接将参数字典传入subs后求解
result = sp.solve(eqn.subs({'b':1, 'c':2}), 'a')
print(result)  # 输出: [3]

关键说明 #

• 始终用二进制模式(wb/rb)操作pickle文件，避免编码错误；
• 替换参数时，确保引用的是SymPy的符号对象（而非普通变量），如果怕混淆，可在加载后重新获取符号：b_sym, c_sym, a_sym = sp.symbols("b c a")；
• 对于复杂物理方程，可先用SymPy的simplify()方法优化方程后再保存，提升后续求解效率。

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