комплексный экспоненциальный множитель от квадрата root - PullRequest
0 голосов
/ 03 мая 2020

Я играю с гауссовым волновым пакетом в sympy, а после и интегралом я получаю следующее выражение

enter image description here

Я хочу извлечь E ** (I * pi / 2) из ​​квадрата root, так что это упрощается с помощью комплексной экспоненты выше, но мне не удалось сделать это с моим (ограниченным) знанием симпати. Я хочу продолжать работать с выражением, и все, что я делаю после этого, приводит только к некомпактным выражениям.

Вот код, который я использую (я пытался избавиться от некоторых символов, но у sympy были проблемы с интеграцией без них ...):

from sympy import *
init_printing(use_latex=True, fontsize='18pt')

x = symbols('x', real=True)
k = Symbol('k', real=True)

a = Symbol('a', real=True, positive=True)
l = Symbol('l', real=True, positive=True)
A = Symbol('A', real=True, positive=True)

psi = E**(-a*x**2)*E**(I*l*x)

ft_func = psi*E**(-I*k*x)

phi = integrate(ft_func, (x,-oo,oo))

psit_func = phi*E**(I*(k*x-A*k**2))

psit = simplify(integrate(powsimp(expand(psit_func)), (k,-oo,oo), conds='none'))

Я использую Python 3.6.9 и Sympy 1.4.

Спасибо

...