In [1]:
import matplotlib.pyplot as plt
In [2]:
from numpy import *
In [3]:
def optical_force(mpb_data):
    wavg = 0.5*(mpb_data[:-1,1]+mpb_data[1:,1])
    dw = mpb_data[1:,1]-mpb_data[:-1,1]
    vg = 0.5*(mpb_data[:-1,2]+mpb_data[1:,2])
    ds = mpb_data[1,0]-mpb_data[0,0]
    return multiply(multiply(-1/wavg, dw/ds), 1/vg)
In [4]:
odd_data = genfromtxt('parallel-wvgs-yodd.dat', delimiter=",");
In [5]:
force_odd = optical_force(odd_data)
In [6]:
s_odd = odd_data[:-1,0]
In [7]:
even_data = genfromtxt('parallel-wvgs-yeven.dat', delimiter=",");
In [8]:
force_even = optical_force(even_data)
In [9]:
s_even = even_data[:-1,0]
In [10]:
plt.plot(s_odd,force_odd,'b-',label='antisymmetric');
In [11]:
plt.plot(s_even,force_even,'r-',label='symmetric');
In [12]:
plt.xlabel("waveguide separation s/a");
In [13]:
plt.ylabel("optical force (F/L)(ac/P)");
In [14]:
plt.legend(loc='upper right');
In [15]:
plt.xticks(arange(0,1.2,0.2));
In [16]:
plt.yticks(arange(-1.5,1.0,0.5));
In [17]:
plt.show()