| 查看: 1249 | 回复: 4 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】求解光纤包层模的有效折射率 已有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 ] |
回复此楼» 猜你喜欢
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有61人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有8人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有12人回复
-
桂林理工大学物理学专业招收调剂
已经有11人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助！光纤的相关问题~~
已经有9人回复
- 如何去除AAO模板
已经有5人回复
- 用DPM模型时，solution里面会不会有 discrece model方程
已经有4人回复
- 请问谁知道康宁单模光纤SMF-28e的折射率是多少？
已经有3人回复
- 咨询光纤外部（周围）折射率与光纤本身的有效折射率之间的联系与区别，
已经有3人回复
- 【求助】求解N-烷基化反应中遇见的问题
已经有8人回复
- 昨晚第一次梦到了基金，今天竟得到了点点模棱两可的消息！
已经有171人回复
- 【求助】关于长周期光栅有效折射率求解，大家帮忙看一下
已经有6人回复
- 【求助】求材料的有效折射率？
已经有8人回复
struggle2011
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1094.5
- 散金: 200
- 帖子: 116
- 在线: 50.4小时
- 虫号: 1026276
- 注册: 2010-05-22
- 性别: GG
- 专业: 光学
4楼2018-03-03 14:22:57
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
小木虫(金币+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 ] |
2楼2010-08-26 11:35:51
struggle2011
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1094.5
- 散金: 200
- 帖子: 116
- 在线: 50.4小时
- 虫号: 1026276
- 注册: 2010-05-22
- 性别: GG
- 专业: 光学
5楼2018-03-03 14:36:59
