sesan_f.m

本软件使用matlab预言编写

clc
clear
format long
load HE_f.mat N_si N_he N_cl KJ lamda f f_eff
c=3.*10^8;
n=length(f);
m=length(lamda);
L=[];
F=[];
for i=1: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(FindZ(D))+L2(FindZ(D)-1))./2];
    L=[L,L2(FindD(D))];
    F=[F,f(i+1)]
    plot(L2,D),hold on
end
hold off
axis([0.5,3.5,-150,100])
figure
%plot(BoChang,Dw)
plot(F,L)
%plot(BoChang,0*BoChang),hold off
axis([0,0.7,0,4])
xlabel('kongjing')
ylabel('bochang')
title('普通单模光纤色散曲线')
grid on