同感想象麻木

Question

我在忙着写一些同情代码来处理虚数符号表达式。

首先，我想让它把x和y作为实数，找到x=y的解。所以我可以这样做。

x, y = sympy.symbols("x y", real=True)  
print(sympy.solve([x-sympy.I*y]))

（SymPy solve接受值列表，所有值都必须为0。所以x-iy=0=>；x=iy）。SymPy会正确地告诉我

[{x: 0, y: 0}]

但是，如果我这样做（理论上相同）：

x, y = sympy.symbols("x y")
print(sympy.solve([x-sympy.I*y, sympy.im(y), sympy.im(x)]))

现在辛皮告诉我

[{re(y): y, re(x): I*y, im(x): 0, x: I*y, im(y): 0}]

这在技术上是正确的，但并没有为我做任何事。这只是SymPy中的一个限制吗，或者我可以通过这样约束复数x和y来得到x=y=0吗？

Answer 1

下面是一个例子，说明如何更一般地解决问题。我使用了一个提示，即使用实际的虚构字符将是一个麻烦，并使用collect_const（）函数来执行缩减。

'''
Converts T to Pi Circuit Topology Symbolically
'''

from sympy import simplify, Symbol, pprint, collect_const

 # Use this for your imaginary symbol
j = Symbol('j')

# Circuit symbols
R = Symbol('R')
w = Symbol('w')
L = Symbol('L')
C = Symbol('C')

# Arbitrary Circuit Element Equations
R1 = 1 / (j*w*C)
R2 = R + j*w*L
R3 = R + j*w*L

# Circuit conversion equations
RN = R1 * R2 + R2 * R3 + R1 * R3 
RA = RN / R1
RB = RN / R2
RC = RN / R3

#Print the original circuit element equations
pprint(R1)
print("\n")
pprint(R2)
print("\n")
pprint(R3)

#Print the original solved equation, followed by the appropriately reduced equation
print("\nOriginal\t")
pprint(RA)
print("\nReduced\t")
pprint(collect_const(simplify(RA),j))

print("\nOriginal\t")
pprint(RB)
print("\nReduced\t")
pprint(collect_const(simplify(RB),j))

print("\nOriginal\t")
pprint(RC)
print("\nReduced\t")
pprint(collect_const(simplify(RC),j))

Answer 2

因为SymPy比复数更擅长于简化实数对，所以下面的策略有帮助：为实部/虚部设置实数变量，然后从中形成复数。

from sympy import *
x1, x2, y1, y2 = symbols("x1 x2 y1 y2", real=True)  
x = x1 + I*x2
y = y1 + I*y2

现在x和y可以作为复杂的变量在一个方程中使用，比如你的

sol = solve([x-I*y, im(y), im(x)])
print(x.subs(sol[0]), y.subs(sol[0])) 

输出：0 0。