Sympy代码为何无法识别RC电路传递函数中的符号变量s？

编写RC电路传递函数计算代码时，LaTeX输出显示传递函数包含s，但用h.has(s)检测时却抛出错误，提示未包含s。

问题原因 #

你手动创建的sympy.symbols('s')和LCapy内部使用的拉普拉斯变量s是不同的符号对象。SymPy的has()方法判断的是对象身份（即是否为同一个实例），而非符号名称，因此即使名字都是s，也会检测失败。

解决方案 #

方案1：使用LCapy内置的符号s

直接调用LCapy提供的s变量，而非自行创建，这样和LCapy内部使用的是同一个对象，has()方法可以正确识别：

import lcapy
from IPython.display import display, Math

# 使用LCapy内置的拉普拉斯变量s
s = lcapy.s

# 定义RC电路网表
netlist = """
R1 1 2
C1 2 0
"""

circuit = lcapy.Circuit(netlist)
h = circuit.transfer(0, 1, 0, 2).simplify()
display(Math(f"H(s) = {h.latex()}"))

# 现在可以正确检测到s
if not h.has(s):
    raise ValueError("Transfer function h does not contain 's'.")

方案2：转换为SymPy表达式后检查符号名称

如果需要使用自己定义的s，可以将LCapy的传递函数转换为SymPy表达式，通过遍历自由符号的名称来判断：

import lcapy
import sympy as sp
from IPython.display import display, Math

netlist = """
R1 1 2
C1 2 0
"""

circuit = lcapy.Circuit(netlist)
h = circuit.transfer(0, 1, 0, 2).simplify()
display(Math(f"H(s) = {h.latex()}"))

# 将LCapy对象转换为SymPy表达式
h_sym = h.expr
# 检查是否存在名为's'的自由符号
has_s = any(sym.name == 's' for sym in h_sym.free_symbols)
if not has_s:
    raise ValueError("Transfer function h does not contain 's'.")

补充说明 #

LCapy的TransferFunction对象内部维护了自己的符号系统，直接使用其内置的s是最稳妥的方式。如果需要和SymPy的符号系统交互，通过.expr属性转换为SymPy表达式后再操作即可。

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