Annexe
Annexe 1 : représentation de l'indice
réfraction de la silice
>> B1=0.6961663; >> B2=0.4079426; >>
B3=0.8974794; >> C1=0.0684043; >> =0.1162414; >> C3=9.896161;
>> lambda=0.8:0.1:1.8;
>> ne=sqrt(1 + (B1./(1-C1./lambda.^2)) +
(B2./(1-./lambda.^2)) + (B3./(1- >> plot(lambda,ne)
Annexe 2 : représentation de l'indice de
groupe
>> B1=0.6961663; >> B2=0.4079426; >>
B3=0.8974794; >> C1=0.0684043; >> =0.1162414; >> C3=9.896161;
>> lambda=0.8:0.1:1.8;
>> ne=sqrt(1 + (B1./(1-C1./lambda.^2)) +
(B2./(1-./lambda.^2)) + (B3./(1- C3./lambda.^2)));
>> n=ne.^2;
>> n1=1./sqrt(n);
>> x1=(2.*B1.*C1.^2.*lambda)./(lambda.^2-C1.^2).^2;
>> x2=(2.*B2.*.^2.*lambda)./(lambda.^2-.^2).^2;
>> x3=(2.*B3.*C3.^2.*lambda)./(lambda.^2-C3.^2).^2;
>> m=x1+x2+x3;
>> z=(lambda./2).*m.*n1;
>> Ng=ne+z;
>> plot(lambda,Ng)
Annexe 3 : Programme de représentation de la
dispersion du guide, la dispersion du matériau et la dispersion
chromatique (dispersion totale).
%Dmat
x=0.8:0.1:1.8;
c=30000000;
a1=0.0696404;
a2=0.1162412;
a3=9.896161;
b1=0.696163;
b2=0.4019426;
b3=0.8974794;
f11=1+(b1*x)./((x."2)-(a1."2));
f12=1+(b2*x)./((x."2)-(a2."2)); f13=1+(b3*x)./((x."2)-(a3."2));
f=0.5./sqrt(f11+f12+f13);
g11=(2.*b1.*a1.*a1.*(3.*x.*x+a1.*a1))./((x.*x-a1.*a1)."3);
g12=(2.*b2.*a2.*a2.*(3.*x.*x+a2.*a2))./((x.*x-a2.*a2)."3);
g13=(2.*b3.*a3.*a3.*(3.*x.*x+a3.*a3))./((x.*x-a3.*a3)."3); g=g11+g12+g13;
k=-(1./4).*(f11+f12+f13)."-1.5;
d11=(2.*b1.*a1.*a1.*x)./((x.*x-a1.*a1)."2);
d12=
|
(2
|
.*b2.*a2.*a2.*x)./((x.*x-a2.
|
*a2)
|
."2);
|
d13=
|
(2
|
.*b3.*a3.*a3.*x)./((x.*x-a3.
|
*a3)
|
."2);
|
d=(d11+d12+d13)."2; Dmat=-10."12.*(x./c).*(f.*g+k.*d);
plot(x,Dmat)
%Dguid
Dguid=-54.36./x; %Dchr
Dchr=Dmat+Dguid; hold on
plot(x,Dmat,x,Dchr,x,Dguid)
grid
|