24小时热门版块排行榜

>论坛更新日志 (3935)
>虫友互识 (949)
>论文投稿 (284)
>休闲灌水 (284)
>导师招生 (178)
>文献求助 (164)
>基金申请 (149)
>硕博家园 (106)
>考博 (87)
>公派出国 (76)
>博后之家 (75)
>教师之家 (70)
>考研 (62)
>论文道贺祈福 (41)
>招聘信息布告栏 (40)
>找工作 (36)

返回列表

shuanglei

铜虫 (小有名气)

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

[交流] 【求助】求解光纤包层模的有效折射率已有3人参与

% function fm=zheshe(neff)

neff=1.23;

a1=2.625;                                                          %纤芯半径
a2=62.5;                                                             %包层半径
n1=1.458;                                                          %纤芯折射率
n2=1.450;                                                          %包层折射率
n3=1.0;                                                             %空气折射率
lamda=1550;
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);

fm=g1-g0;

return
求解neff 。我调试好久没有弄出来请各位大侠帮帮指教一下。

[ Last edited by luoqiquan on 2010-10-21 at 05:20 ]

回复此楼

» 猜你喜欢

【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有144人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求助！光纤的相关问题~~ 已经有9人回复
如何去除AAO模板已经有5人回复
用DPM模型时，solution里面会不会有 discrece model方程已经有4人回复
请问谁知道康宁单模光纤SMF-28e的折射率是多少？已经有3人回复
咨询光纤外部（周围）折射率与光纤本身的有效折射率之间的联系与区别，已经有3人回复
【求助】求解N-烷基化反应中遇见的问题已经有8人回复
昨晚第一次梦到了基金，今天竟得到了点点模棱两可的消息！已经有171人回复
【求助】关于长周期光栅有效折射率求解，大家帮忙看一下已经有6人回复
【求助】求材料的有效折射率？已经有8人回复

硬着头皮钻下去？！

1楼 2010-08-24 15:56:27

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lhfx_313

至尊木虫 (文坛精英)

应助: 1 (幼儿园)
贵宾: 0.155
金币: 29747.1
散金: 2852
红花: 20
沙发: 102
帖子: 11130
在线: 1568.7小时
虫号: 550265
注册: 2008-04-23
专业: 通信理论与系统

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
luoqiquan(金币+2):thank you 2010-08-30 17:34:09

求解neff?? 是求解fm吧？
试试下面的

clc;
clear;
neff=1.23;

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

t1=1*i*neff/z0;                                                       %t1、t2都是数 ?????
t2=1*i*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);
%m1besselj(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);

fm=g1-g0

结果：fm =

      0 - 0.0191i

[ Last edited by lhfx_313 on 2010-8-26 at 11:38 ]