Что-то типа такого должно сработать. Нужно чтобы диапазон

был достаточно большим (

, где

- количество полученных точек). Фаза у комплексных чисел будет сбита на константу.
Код:
dx = % шаг дискретизации по x
n = % кол-во точек дискретизации
E = % значения ф-и E(x)
dp = (2*pi)/(n*dx)
P = 0:dp:1
k = length(P)
Z = fft(E, k)
T = asin(P)
plot(T, abs(Z))
Получилось. Спасибо Вам.