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