Matlab阶跃折射率分布均匀光纤的光场分布模拟-abbr_e825a256c50d9e59443513624838742e.rar

Matlab阶跃折射率分布均匀光纤的光场分布模拟-abbr_e825a256c50d9e59443513624838742e.rar
本设计基于matlab语言，对阶跃折射率分布均匀光纤（圆柱型光纤）模型中的光场分布进行了模拟。

clear all
n1=1.5075;
n2=1.5;
r=9.8*10^（-6）;
dn=（n1-n2）/n1;
for ii=1:101
    r1（ii）=（r/100）*ii;
end
for ii=1:60
    n（ii）=n1*sqrt（（1-2*dn*（r1（ii）/r1（60））））;
end
for ii=61:101
    n（ii）=n2;
end
rh=9.8*10^（-8）;
k0=2*pi/（1.55*10^（-6））;
m=1;
zh=20*10^（-8）;
n0=1.503;
for ii=1:101
    b（ii1）=（k0^2）*（rh^2）*（n（ii）^2-（n0）^2）-2*（2-m^2/2）;
end%系数b
for ii=1:101
        f（ii1）=（rh^2）*（k0^2）*（n（ii）^2-（n0）^2）-（rh^2）*（m^2）/r^2-2;
end%系数f
a=4*（k0^2）*（n0^2）*（rh^2）/（1+i*k0*n0*zh）;
c=（1+exp（i*m*pi））*（2-m^2/2）;
d=（1-k0*n0*i*zh）/（1+k0*n0*i*zh）;
e=zeros（1011）;
g=zeros（1011）;
for ii=2:101
    e（ii1）=1-1/（2*（ii-1））;
end
for ii=2:101
      g（ii1）=1+1/（2*（ii-1））;  
end%系数整理结束
 %x=0;%光开始输入
 %h1=0.03; 
 %for jj=1:101
  %  u（1jj）=exp（-（（x+jj）*h1）^2）;
 %end%光结束
 r0=2.9815*10^（-6）;
for ii=1:101
    r（ii）=rh*（ii-1）;
end
for jj=1:101
    u（1jj）=0.25*10^6*r（jj）*exp（-3.33*10^5*（r（jj）-r0））^2;
end
for m1=1:101%大程序开始
A=zeros（101101）;
A（11）=a+b（11）;
A（12）=c;
for ii=2:101
    A（iiii-1）=e（ii1）;
end
for ii=2:100
    A（iiii）=f（ii1）+a;
end
for ii=2:99
    A（iiii+1）=g（ii1）;
end%左面的系数矩阵A结束
B=zeros（101101）;
B（11）=a+b（11）*d;
B（12）=c*d;
for ii=2:101
    B（iiii-1）=e（ii1）*d;
end
for ii=2:100
    B（iiii）=f（ii1）*d+a;
end
for ii=2:99
    B（iiii+1）=g（ii1）*d;
end%右面的系数矩阵B结束
d1=ones（101101）;
y=ones（101101）;
%z=ones（101101）;
d1（:m1）=B*u.‘;%转置
  d2=ones（101101）;
for ii=1:101
    d2（m1ii）=d1（iim1）;%变成一行一行的
end
    %下面用追赶法
  h2（m11）=c./（a+b（11））;
    for jj=2:100
    h2（m1jj）=g（jj1）./（f（11）+a-e（jj1）*h2（m1jj-1））;
    end
    y（m11）=d2（m11）/（a+b（11））;
    for jj=2:101
    y（m1jj）=（d1（jjm1）-e（jj1）.*y（m1jj-1））./（f（11）+a-e（jj1）.*h2（m1jj-1））;%problem
    end
    z（m1101）=y（m1101）;
    for jj=100:-1:1
    z（m1jj）=y（m1jj）-h2（m1jj）.*z（m1jj+1）;
    end
    for jj=1:101
        u（1jj）=z（m1jj）;
    end
    for jj=1:101
    q（m1jj）=z（m1jj）;
    end
    end%最大循环的end
    [xy]=meshgrid（0:1:1000:zh:99*zh）;
    mesh（abs（q））;
    surf（abs（q））;
    shading interp 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2144  2006-12-09 17:16  yuanzhus.m

     文件       2144  2006-12-09 17:16  是matlab语言，对阶跃折射率分布均匀光纤（圆柱型光纤）模型中的光场分布进行了模拟\yuanzhus.m

     目录          0  2010-06-03 22:38  是matlab语言，对阶跃折射率分布均匀光纤（圆柱型光纤）模型中的光场分布进行了模拟

----------- ---------  ---------- -----  ----

                 4724                    5