传输线输入阻抗公式ZL复数解求解代码TypeError报错排查与修正

问题描述 #

基于传输线理论的输入阻抗公式，给定频率、特征阻抗Z0、传输线长度、输入阻抗Zin，求解复数形式的负载阻抗ZL时，使用SymPy编写代码出现TypeError: Cannot convert expression to float错误，原代码如下：

import sympy as sy
import math

#calculate beta
pi = sy.pi.evalf()
freq = 20*10**6 #float(input('frequency[MHz]='))
lamb = 2.99792458*10**8/freq  #2.99792458*10**8[ms-1]
beta = 2*pi/lamb

#input values
z0 = 50 #float(input('Z0[ohm]='))
leng = 0.5 #float(input('length[m]='))
zin = 50 #complex(input('Zin[ohm]='))

#Zin equation
zl = sy.Symbol('zl')
tan = math.tan(beta*leng)
print('tan=',tan)
a = sy.solveset(zin,z0*(complex(zl,z0*tan))/(complex(z0,zl*tan)))
print(a)

错误原因 #

1. 数值与符号函数混用：math.tan仅支持浮点数运算，无法兼容SymPy的符号变量zl，触发类型转换错误，应使用SymPy内置的sy.tan支持符号运算。
2. 复数构造方式错误：complex()仅能接收数值参数，不能直接包装SymPy符号zl，传输线公式的复数形式需用SymPy的虚数单位sy.I来表达。
3. 求解函数参数顺序错误：solveset的正确调用格式是solveset(方程, 待求解变量)，原代码未构建有效等式，且参数顺序颠倒。

修正后的代码 #

import sympy as sy

# 输入参数（可根据需求修改或改为输入交互）
freq = 20 * 10**6
c = 2.99792458 * 10**8  # 光速
z0 = 50
leng = 0.5
zin = 50  # 支持复数输入，例如 zin = complex(50, 10)

# 计算相位常数beta
beta = 2 * sy.pi * freq / c

# 定义复数类型的负载阻抗符号
zl = sy.Symbol('zl', complex=True)

# 构建传输线输入阻抗等式
input_impedance_eq = sy.Eq(
    zin,
    z0 * (zl + z0 * sy.I * sy.tan(beta * leng)) / (z0 + zl * sy.I * sy.tan(beta * leng))
)

# 求解负载阻抗ZL
solution = sy.solveset(input_impedance_eq, zl)

print("负载阻抗ZL的解：")
print(solution)

关键说明 #

• 全程使用SymPy的符号运算工具（sy.pi、sy.tan、sy.I），避免数值与符号类型冲突。
• 通过sy.Eq构建标准等式，匹配solveset的求解要求。
• 为zl指定complex=True，明确其复数属性，提升求解精度。

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