sesan_lamda.m

本软件使用matlab预言编写

clc
clear
format long
load HE_Lamda.mat N_si N_he N_cl KJ lamda f f_eff
c=3.*10^8;
n=length(KJ);
m=length(lamda);
L=[];
K=[];
for i=0:n-1
    [L1,D1]=ChaDao(lamda,N_he(i.*m+1:i.*m+m));
    [L2,D2]=ChaDao(L1,D1);
    D=-10.^12.*L2.*D2./c;
    L=[L,L2(FindD(D))];
    K=[K,KJ(i+1)]
    plot(L2,D),hold on
end
hold off
axis([0.5,3.5,-100,100])
figure
%plot(BoChang,Dw)
plot(K,L)
%plot(BoChang,0*BoChang),hold off
axis([0,4,0,4])
xlabel('kongjing')
ylabel('bochang')
title('普通单模光纤色散曲线')
grid on