###########################################################################
# Scriptfile: mod_response.lsf
#
# Description:
# This file calculates the response of a modulator
# 
# 
#
# Copyright 2012, Lumerical Solutions, Inc.
###########################################################################

V = linspace(0.5,1,21);

neff=[	2.50652+1.177e-006i,  
	2.50652+1.17854e-006i,  
	2.50652+1.18042e-006i,  
	2.50652+1.1828e-006i,  
	2.50652+1.18591e-006i,  
	2.50652+1.19e-006i,  
	2.50652+1.19538e-006i,  
	2.50652+1.20238e-006i,  
	2.50652+1.21148e-006i,  
	2.50652+1.22368e-006i,  
	2.50652+1.24235e-006i,  
	2.50651+1.28087e-006i,  
	2.50651+1.38715e-006i,  
	2.5065+1.69323e-006i,  
	2.50647+2.52136e-006i,  
	2.50641+4.65729e-006i, 
	2.5063+9.65742e-006i,  
	2.50611+1.85249e-005i,  
	2.5059+2.92719e-005i,  
	2.50571+3.98431e-005i,  
	2.50554+4.9736e-005i];

# Calculate the modulation response
L = 1000e-6;
la0 = 1.55e-6;
arg = 2*pi*real(neff)/la0*L;
ang1 = arg - floor(real(neff)/la0*L)*2*pi;
ang2 = ang1(1);
T = 0.25*abs(exp(-1i*ang1 - 2*pi*imag(neff)/la0*L) + exp(-1i*ang2 - 2*pi*imag(neff(1))/la0*L))^2;
plot(V,T,"Voltage(V)","Transmission (a.u.)");

?"Extinction ratio: " + num2str(10*log10(max(T)/min(T))) + "dB";
?"Insertion loss: " + num2str(-10*log10(max(T))) + "dB";