%做關(guān)于Sapphire(n_o)、Sapphire(n_e)、Fuse silica、BK7 (K9)、 SF10晶體二階色散曲線圖
% Sapphire(n_o)
B1=1.5039759;                     
B2=0.55069141;
B3=6.5927379;
C1=0.0054804113;
C2=0.014799428;
C3=402.89514;
c=3*10^8;
x=0.4:0.001:2.0;
n=(1+(B1.*x.^2)./(x.^2-C1)+(B2.*x.^2)./(x.^2-C2)+(B3.*x.^2)./(x.^2-C3)).^0.5;
X1=gradient(n);                                                                    %dn/dx
X2=gradient(X1);                                                                   %d2n/dx2
D1=-x.^2/(2*pi.*c).*X1;                                                            %dn/dw
D2=2*x.^3/(4*pi.*pi.*c.*c).*X1+x.^4/(4*pi.*pi.*c.*c).*X2;                          %d2n/dw2
DD2=(2*D1+2*pi.*c.*D2./x)./c;
plot(x,DD2.*10^30,'b')
hold on;
% Sapphire(n_e)
B1=1.4313493;
B2=0.65054713;
B3=5.3414021;
C1=0.0052799261;
C2=0.014238265;
C3=325.01783;
x=0.4:0.001:2.0;
n=(1+(B1.*x.^2)./(x.^2-C1)+(B2.*x.^2)./(x.^2-C2)+(B3.*x.^2)./(x.^2-C3)).^0.5;
X1=gradient(n);                                                                    %dn/dx
X2=gradient(X1);                                                                   %d2n/dx2
D1=-x.^2/(2*pi.*c).*X1;                                                            %dn/dw
D2=2*x.^3/(4*pi.*pi.*c.*c).*X1+x.^4/(4*pi.*pi.*c.*c).*X2;                          %d2n/dw2
DD2=(2*D1+2*pi.*c.*D2./x)./c;
plot(x,DD2.*10^30,'r')
hold on;
% Fuse silica
B1=0.6961663;
B2=0.4079426;
B3=0.8974794;
C1=0.0046791483;
C2=0.013512063;
C3=97.934003;
x=0.4:0.001:2.0;
n=(1+(B1.*x.^2)./(x.^2-C1)+(B2.*x.^2)./(x.^2-C2)+(B3.*x.^2)./(x.^2-C3)).^0.5;
X1=gradient(n);                                                                    %dn/dx
X2=gradient(X1);                                                                   %d2n/dx2
D1=-x.^2/(2*pi.*c).*X1;                                                            %dn/dw
D2=2*x.^3/(4*pi.*pi.*c.*c).*X1+x.^4/(4*pi.*pi.*c.*c).*X2;                          %d2n/dw2
DD2=(2*D1+2*pi.*c.*D2./x)./c;
plot(x,DD2.*10^30,'g')
hold on;
% BK7 (K9)
B1=1.0396121;
B2=0.23179234;
B3=1.0104695;
C1=0.0060006987;
C2=0.020017914;
C3=103.56065;
x=0.4:0.001:2.0;
n=(1+(B1.*x.^2)./(x.^2-C1)+(B2.*x.^2)./(x.^2-C2)+(B3.*x.^2)./(x.^2-C3)).^0.5;
X1=gradient(n);                                                                    %dn/dx
X2=gradient(X1);                                                                   %d2n/dx2
D1=-x.^2/(2*pi.*c).*X1;                                                            %dn/dw
D2=2*x.^3/(4*pi.*pi.*c.*c).*X1+x.^4/(4*pi.*pi.*c.*c).*X2;                          %d2n/dw2
DD2=(2*D1+2*pi.*c.*D2./x)./c;
plot(x,DD2.*10^30,'m')
hold on;
% SF10
B1=1.6162598;
B2=0.25922933;
B3=1.1749087;
C1=0.013606860;
C2=0.061596046;
C3=121.92271;
x=0.4:0.001:2.0;
n=(1+(B1.*x.^2)./(x.^2-C1)+(B2.*x.^2)./(x.^2-C2)+(B3.*x.^2)./(x.^2-C3)).^0.5;
X1=gradient(n);                                                                    %dn/dx
X2=gradient(X1);                                                                   %d2n/dx2
D1=-x.^2/(2*pi.*c).*X1;                                                            %dn/dw
D2=2*x.^3/(4*pi.*pi.*c.*c).*X1+x.^4/(4*pi.*pi.*c.*c).*X2;                          %d2n/dw2
DD2=(2*D1+2*pi.*c.*D2./x)./c;
plot(x,DD2.*10^30,'c')
title('二階色散-波長曲線');
xlabel('wavelength');
ylabel('二階色散');
legend('Sapphire(n_o)','Sapphire(n_e)','Fuse silica','BK7','SF10')
hold on;