clc
clear all
a=1;
r1=0.9;
r2=0.9;
r0=5.025e-6; %ring1
l=2*pi*r0;
n1=3.191;
%n1=3.191044517;
n2=3.193525;%10k
n3=3.19653481;%20k
n21=3.18874; %-10k
n31=3.18546519; %-5k
lambda=1.54e-6:0.00000000001:1.56e-6;
%lambda=1.5493e-6:0.0000000000000001:1.5496e-6; %FSR
k0=2.*pi./lambda;
phi=k0.*n1.*l;
phi1=k0.*n2.*l;
phi2=k0.*n3.*l;
phi21=k0.*n21.*l;
phi31=k0.*n31.*l;
Tp=((r2.*a).^2-2.*a.*r1.*r2.*cos(phi)+r1.^2)./(1-2.*r1.*a.*r2.*cos(phi)+(r1.*r2.*a).^2);
%Tp1=((r2.*a).^2-2.*a.*r1.*r2.*cos(phi1)+r1.^2)./(1-2.*r1.*a.*r2.*cos(phi1)+(r1.*r2.*a).^2);
Td=0.54*(1-r1.^2).*(1-r2.^2).*a./(1-2.*r1.*a.*r2.*cos(phi)+(r1.*r2.*a).^2);
Td1=0.54*(1-r1.^2).*(1-r2.^2).*a./(1-2.*r1.*a.*r2.*cos(phi1)+(r1.*r2.*a).^2);
Td2=0.54*(1-r1.^2).*(1-r2.^2).*a./(1-2.*r1.*a.*r2.*cos(phi2)+(r1.*r2.*a).^2);
Td21=0.54*(1-r1.^2).*(1-r2.^2).*a./(1-2.*r1.*a.*r2.*cos(phi21)+(r1.*r2.*a).^2);
Td31=0.54*(1-r1.^2).*(1-r2.^2).*a./(1-2.*r1.*a.*r2.*cos(phi31)+(r1.*r2.*a).^2);
%Td=(1-r1.^2).*(1-r2.^2).*a./(1-2.*r1.*a.*r2.*cos(phi+phi1)+(r1.*r2.*a).^2);
plot(lambda,Td31,'-.g','Linewidth',2)
hold on
plot(lambda,Td21,'-.r','Linewidth',2)
hold on
plot(lambda,Td,'Linewidth',2)
hold on
plot(lambda,Td1,'r','Linewidth',2)
hold on
plot(lambda,Td2,'g','Linewidth',2)
hold off
legend('P=-200kPa','P=-100kPa','P=0Pa','P=100kPa','P=200kPa')
xlabel('Wavelength (�m)','fontsize',20)
ylabel('Normalized Power','fontsize',20)