sesan.m

本软件使用matlab预言编写

clc
clear
format long
load HE_0.7.mat N0 N1 N2 KJ lamda f f_eff
plot(lamda,N0,'-b',lamda,N1,':',lamda,N2,'-b')
c=3.*10^8;
b_a=polyfit(lamda,N1,10)
tmp=polyval(b_a,lamda)
figure
plot(lamda,tmp,lamda,N1,'o')
b1=polyder(b_a);
b2=polyder(b1);
D=-10.^12.*lamda.*polyval(b2,lamda)./c;


figure
%plot(BoChang,Dw)
plot(lamda,D)
%plot(BoChang,0*BoChang),hold off
axis([1.2,1.6,-1000,100])
xlabel('波长（um）')
ylabel('色散D（ps/km*nm）')
title('普通单模光纤色散曲线')
grid on