最新消息:

【乞助】MATLAB求解光纤包层模有效折射率,叨教

365bet admin 浏览

  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 ]

    发表我的评论
    取消评论

    表情

    您的回复是我们的动力!

    • 昵称 (必填)

    网友最新评论