www.pudn.com > waveguide_dispersion.zip > waveguide.m, change:2011-12-13,size:728b


clc 
j=1; 
for wl=1.2:0.4/39:1.6%波长 
d=10*10^(-6);%厚度 
dn=0.003; 
 
c=3*10^5; %km/s 
l=1;%每1km 
a=d/2; 
wl1=wl.*10^(-6); 
wl2=wl.*wl; 
n1=sqrt(1+0.6961663.*wl2./(wl2-0.068403^2)+0.4079426.*wl2/(wl2-0.1162414^2)+0.8974794.*wl2/(wl2-9.896161^2)); 
n2=n1-n1.*dn; 
k=2*pi./wl1; 
 
n_=neff(n1,n2,n2,d,wl1);%有效折射率 
V=k.*a.*n1.*sqrt(2.*dn); 
b=(n_.^2-n2.^2)/(n1.^2-n2.^2); 
 
t=l.*n_./c; 
 
%N=10^2; 
dw=10^(-6);%wl1./N;%小量 
wl11=wl1+dw; 
n_1=neff(n1,n2,n2,d,wl11); 
t1=l.*n_1./c; 
dt=t1-t;%时延差 
Dw1(j)=dt./dw; 
Dw(j)=Dw1(j)*10^3;%单位转换s/(m.km)--ps/(nm.km) 
j=j+1; 
end 
wl=1.6:-0.4/39:1.2; 
j=1; 
plot(wl,Dw); 
title('波导色散'); 
xlabel('wavelength(um)'); 
ylabel('dispersion(ps/(nm-km))'); 
 
figure(1); 
grid on