matlab2.rar_resonator_ring_ring resonator_sorop

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)