24小时热门版块排行榜

返回列表

shuanglei

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 56.5
散金: 8
红花: 1
帖子: 130
在线: 37.7小时
虫号: 1051296
注册: 2010-07-03
性别: GG
专业: 光学

[交流] 【求助】MATLAB求解光纤包层模有效折射率,请教matlab语法问题

clc;
clear;

a1=2.625;                                                          %纤芯半径
a2=62.5;                                                             %包层半径
n1=1.458;                                                          %纤芯折射率
n2=1.450;                                                          %包层折射率
n3=1.0;                                                             %空气折射率
lamda=1.55;
z0=377;                                                             %真空电磁阻抗

t1=1i*neff/z0;                                                       %t1、t2都是复数 ?????
t2=1i*neff*z0;

u1=((2*pi/lamda)^2*(n1^2-neff^2))^(1/2);
u2=((2*pi/lamda)^2*(n2^2-neff^2))^(1/2);                            %??????
e3=((2*pi/lamda)^2*(neff^2-n3^2))^(1/2);

u21=1/u2^2-1/u1^2;
u32=1/e3^2+1/u2^2;

m1=besselj(0,u1*a1)-1/(u1*a1)*besselj(1,u1*a1);                         %m1为besselj(1,u1*a1)的一阶微商

m2=-besselk(0,e3*a2)-1/(e3*a2)*besselk(1,e3*a2);                         %m2为besselk(1,e3*a2)的一阶微商

m3=bessely(0,u2*a1)-1/(u2*a1)*bessely(1,u2*a1);                         %m3为bessely(1,u2*a1)的一阶微商

m4=besselj(0,u2*a1)-1/(u2*a1)*besselj(1,u2*a1);                         %m4为besselj(1,u2*a1)的一阶微商

m5=besselj(0,u2*a2)-1/(u2*a2)*besselj(1,u2*a2);                         %m5为besselj(1,u2*a2)的一阶微商

m6=bessely(0,u2*a2)-1/(u2*a2)*bessely(1,u2*a2);                         %m6为bessely(1,u2*a2)的一阶微商

bej=m1/(u1*besselj(1,u1*a1));
bek=m2/(e3*besselk(1,e3*a2));

pa2=besselj(1,u2*a2)*bessely(1,u2*a1)-besselj(1,u2*a1)*bessely(1,u2*a2);

qa2=besselj(1,u2*a2)*m3-m4*bessely(1,u2*a2);

ra2=m5*bessely(1,u2*a1)-besselj(1,u2*a1)*m6;

sa2=m5*m3-m4*m6;


g1=(1/t2)*(u2*(bej*bek+t1*t2*u21*u32/(n2^2*a1*a2))*pa2-bek*qa2+bej*ra2-sa2/u2)/(-u2*(u32/(n2^2*a2)*bej-u21/(n1^2*a1)*bek)*pa2+u32/(n1^2*a2)*qa2+u21/(n1^2*a1)*ra2);
g0=t1*(u2*(u32/a2*bej-n3^2*u21/(n2^2*a1)*bek)*pa2-u32/a2*qa2-u21/a1*ra2)/(u2*(n3^2/n2^2*bej*bek+t1*t2*u21*u32/(n1^2*a1*a2))*pa2-n3^2/n1^2*bek*qa2+bej*ra2-n2^2/(n1^2*u2)*sa2);

g1=g0;

求解neff（光纤各阶包层模的有效折射率，其中neff有约一百多个解）。我对MATLAB语法不懂，只把方程（即光纤包层模的本征方程）用MATLAB写了出来。请大虾帮忙解一下。（精确到小数点后8位）

[ Last edited by GrasaVampiro on 2010-8-27 at 08:39 ]

回复此楼