figuren_clb.m

pcf project dds sdfsd sdcsc sdcsc sdxcs gh fgb dfv fdgbvfg b fg fb fgbv gbfbf s bgtb fgbfv b fbv

clear
clc
na=1;% air hole refractive index
pitch=6.41e-6; %distance between to air hole
d = pitch*.375;
wavelength = 1.55e-6;
nsi = get_nsi(wavelength);
n_s = nsi;
ns = nsi;
n_clx = get_nclx(wavelength,d,pitch,na,n_s);
n_effx=get_neffx(wavelength,pitch,n_clx,ns);
w_G= get_WpcfGx(pitch,ns,n_clx,wavelength);
w = w_G
R0 = 500*(10).^(-3);
% i =0;
% for y = -1*(1+.02)*w:.1e-6:(1+.02)*w;
%     i = i+1;
%     neff(i) =  sqrt((n_effx.^2) .* (1 + 2.*y ./ R0));
% end
% y = -1*(1+.02)*w:.1e-6:(1+.02)*w;
% nm = neff;
% plot(y,nm,'Color',[0 0 1],'LineStyle','--','LineWidth',2,'Marker','o');
% hold on
% line([-1*w w],[n_effx n_effx])
y = [w,-w]
n_clx =  sqrt((n_clx.^2) .* (1 + 2.*y ./ R0))
w_G= get_WpcfGx(pitch,ns,n_clx,wavelength)
% legend('d= 0.50 \mu m ; cladding index','d= 0.50 \mu m ; mode index',...
%        'd= 1,00 \mu m ; cladding index','d= 1,00 \mu m ; mode index',...
%        'Location','SouthWest' )
% xlabel('\lambda')
% ylabel('Cladding index for \Lambda = 2.3 \mu m')
% % creat tex box
% annotation1=annotation('textbox','Position',[0.3839 0.8087 0.3256 0.0746],...
%     'EdgeColor','none','String',{'refractive index for pure silica'},...
%     'FitHeightToText','on');
% hold off