In [1]:

    

%pylab inline
from scipy import signal
import math


    

    




Populating the interactive namespace from numpy and matplotlib

In [10]:

    

s1 = np.array([1.0, 1.0j, 1+1j])
s2 = np.array([0.5, 0.5j, 0])
p1 = np.angle(s2/s1)
p2 = np.angle(s2) - np.angle(s1)
print p1
print p2


    

    




[ 0.  0.  0.]
[ 0.          0.         -0.78539816]

In [9]:

    

o1 = zeros(p2.size, dtype=complex)
o1.real, o1.imag = cos(p2), sin(p2)
o2 = np.exp(1j * p2)
print o1
print o2


    

    




[ 1.00000000+0.j          1.00000000+0.j          0.70710678-0.70710678j]
[ 1.00000000+0.j          1.00000000+0.j          0.70710678-0.70710678j]

In [ ]: