Код: Выделить всё
import lcapy
#from lcapy import s, expr
#from IPython.display import display, Math
from sympy import Symbol
#from sympy import I, re, im
cct = Circuit("""
R 1 2
C 2 0""")
H = cct.transfer(1, 0, 2, 0)
print("H =",H)
H=H(jw).simplify()
print("H(jw) =",H)
H_num, H_den = H.as_numer_denom() # Extract numerator and denominator
H_num_real = re(H_num)
H_num_imag = im(H_num)
print("num_real:",H_num_real)
print("num_imag:", H_num_imag)
H_den_real = re(H_den)
H_den_imag = im(H_den)
print("den_real:",H_den_real)
print("den_imag:", H_den_imag)
num_phase = atan2(H_num_imag, H_num_real)
print("num_phase=",num_phase)
den_phase = atan2(H_den_imag,H_den_real)
Я попробовал функцию atan2() и получил ошибку. ожидая четко определенного выражения.
Подробнее здесь: https://stackoverflow.com/questions/790 ... hon-object