应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】关于长周期光栅有效折射率求解，大家帮忙看一下
71/1返回列表
查看: 1819  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

shuanglei

铜虫 (小有名气)

[交流] 【求助】关于长周期光栅有效折射率求解，大家帮忙看一下 已有5人参与

以下是文献中求解光纤包层模有效折射率的方程，即neff为未知数，其他参量都已知。（文献中说neff有很多解）大家帮忙看一下，怎样才能求得neff的多个解，请各位大侠指点，万分感谢！
a1=2.625;                                                             %纤芯半径
a2=62.5;                                                               %包层半径
n1=1.458;                                                             %纤芯折射率
n2=1.450;                                                             %包层折射率
n3=1.0;                                                                  %空气折射率
lamda=1550;
z0=377;                                                             %真空电磁阻抗

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

u1=((2*pi/lamda)^2*(n1-neff^2))^(1/2);
u2=((2*pi/lamda)^2*(n2-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=diff(besselj(1,u1*a1))=besselj(0,u1*a1)-1/(u1*a1)*besselj(1,u1*a1)

m2=-besselk(0,e3*a2)-1/(e3*a2)*besselk(1,e3*a2);                        %m2=diff(besselk(1,e3*a2))=-besselk(0,e3*a2)-1/(e3*a2)*besselk(1,e3*a2)

m3=bessely(0,u2*a1)-1/(u2*a1)*bessely(1,u2*a1);                          %m3=diff(bessely(1,u2*a1))=bessely(0,u2*a1)-1/(u2*a1)*bessely(1,u2*a1)

m4=besselj(0,u2*a1)-1/(u2*a1)*besselj(1,u2*a1);                            %m4=diff(besselj(1,u2*a1))=besselj(0,u2*a1)-1/(u2*a1)*besselj(1,u2*a1)

m5=besselj(0,u2*a2)-1/(u2*a2)*besselj(1,u2*a2);                            %m5=diff(besselj(1,u2*a2))=besselj(0,u2*a2)-1/(u2*a2)*besselj(1,u2*a2)

m6=bessely(0,u2*a2)-1/(u2*a2)*bessely(1,u2*a2);                           %m6=diff(bessely(1,u2*a2))=bessely(0,u2*a2)-1/(u2*a2)*bessely(1,u2*a2)

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

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

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

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

s(a2)=m5*m3-m4*m6;                                                               
                                             
                                                                                                     

g1=1/t2*(u2*(j*k-t1*t2*u21*u32/(n2^2*a1*a2))*p(a2)-k*q(a2)+j*r(a2)-s(a2)/u2)/
      (-u2*((u32/n2^2*a2)*j-(u21/n1^2*a1)*k)*p(a2)+(u32/n2^2*a2)*q(a2)+(u21/n1^2*a1)*r(a2));
g0=t1*(u2*((u32/a2)*j-(n3^2*u21/n1^2*a1)*k)*p(a2)-(u32/a2)*q(a2)-(u21/a1)*r(a2))/
      (u2*((n3^2/n2^2)*j*k+t1*t2*u21*u32/(n1^2*a1*a2))*p(a2)-(n3^2/n1^2)*k*q(a2)+j*r(a2)-(n2^2/n1^2*u2)*s(a2));
g1=g0
其中n3 求neff
回复此楼

» 本帖已获得的红花（最新10朵）

sharinfang

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
硬着头皮钻下去？！
1楼 2010-07-13 08:32:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zjliu1265

至尊木虫 (正式写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
adu886886(金币+3):鼓励应助！欢迎常来仿真模拟板块！ 2010-07-14 07:54:48
我觉得你可以考虑用2分法解，我画了neff~g1-g0的曲线，发现只有一个解，
并非多解：
clc;close all;clear
a1=2.625;       %纤芯半径
a2=62.5;        %包层半径
n1=1.458;    %纤芯折射率
n2=1.450;    %包层折射率
n3=1.0;   %空气折射率
lamda=1550;
z0=377;     %真空电磁阻抗
neff=linspace(a1,a2,201);
t1=i*neff/z0;     %t1、t2都是复数   
t2=i*neff*z0;

u1=(2*pi/lamda)^2*sqrt(n1-neff.^2);
u2=(2*pi/lamda)^2*sqrt(n2-neff.^2);
e3=(2*pi/lamda)^2*sqrt(neff.^2-n3^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=diff(besselj(1,u1*a1))=besselj(0,u1*a1)-1/(u1*a1)*besselj(1,u1*a1)
m2=-besselk(0,e3*a2)-1./(e3*a2).*besselk(1,e3*a2);
%m2=diff(besselk(1,e3*a2))=-besselk(0,e3*a2)-1/(e3*a2)*besselk(1,e3*a2)
m3=bessely(0,u2*a1)-1./(u2*a1).*bessely(1,u2*a1);
%m3=diff(bessely(1,u2*a1))=bessely(0,u2*a1)-1/(u2*a1)*bessely(1,u2*a1)
m4=besselj(0,u2*a1)-1./(u2*a1).*besselj(1,u2*a1);
%m4=diff(besselj(1,u2*a1))=besselj(0,u2*a1)-1/(u2*a1)*besselj(1,u2*a1)
m5=besselj(0,u2*a2)-1./(u2*a2).*besselj(1,u2*a2);
%m5=diff(besselj(1,u2*a2))=besselj(0,u2*a2)-1/(u2*a2)*besselj(1,u2*a2)
m6=bessely(0,u2*a2)-1./(u2*a2).*bessely(1,u2*a2);
%m6=diff(bessely(1,u2*a2))=bessely(0,u2*a2)-1/(u2*a2)*bessely(1,u2*a2)
j=m1./(u1.*besselj(1,u1*a1));                                                  
k=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.*(j.*k-t1.*t2.*u21.*u32/(n2^2*a1*a2)).*pa2-k.*qa2+j.*ra2-sa2./u2)./...
      (-u2.*((u32/n2^2*a2).*j-(u21./n1^2*a1).*k).*pa2+(u32./n2^2*a2).*qa2+(u21/n1^2*a1).*ra2);
g0=t1.*(u2.*((u32/a2).*j-(n3^2*u21/n1^2*a1).*k).*pa2-(u32/a2).*qa2-(u21/a1).*ra2)./...
      (u2.*((n3^2/n2^2).*j.*k+t1.*t2.*u21.*u32/(n1^2*a1*a2)).*pa2-...
      (n3^2/n1^2)*k.*qa2+j.*ra2-(n2^2/n1^2*u2).*sa2);
FF=g1-g0;
%% 其中n3 subplot(121);plot(neff,real(FF));hold on;plot(xlim,[0,0],'r');
subplot(122);plot(neff,imag(FF));hold on;plot(xlim,[0,0],'r');

图形如下，可以看出只有一个解：
赞一下(2人) 回复此楼
2楼2010-07-14 07:51:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shuanglei

铜虫 (小有名气)
请问neff=linspace(a1,a2,201);中的201是什么意思。
赞一下 回复此楼
硬着头皮钻下去？！
3楼2010-08-29 15:43:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shuanglei

铜虫 (小有名气)
引用回帖:
Originally posted by zjliu1265 at 2010-07-14 07:51:23:
我觉得你可以考虑用2分法解，我画了neff~g1-g0的曲线，发现只有一个解，
并非多解：
clc;close all;clear
a1=2.625;       %纤芯半径
a2=62.5;        %包层半径
n1=1.458;    %纤芯折射率
n2=1.450;    %包 ...

请问neff=linspace(a1,a2,201);为什么只让他生成201维的行向量 ？
赞一下 回复此楼
硬着头皮钻下去？！
4楼2010-08-29 16:19:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sharinfang

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
送鲜花一朵
能不能告诉我的你的qq，我有一个超越方程不知道怎样编程
赞一下 回复此楼
5楼2012-10-16 18:00:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

littlewwld

铁虫 (小有名气)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
xiegangmai: 金币+1, 谢谢参与！ 2013-08-01 23:22:40
不同的模式的neff不同，要看你的光纤参数是单模还是多模。另外不同波长（不同模式）也有影响！
自己看一下光纤参数。
赞一下 回复此楼
拒接理所当然
6楼2013-07-28 21:34:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wawacgh

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
楼主，您好！这个问题您解决了吗？您知道如何算纤芯与同向包层模式的耦合系数吗？
赞一下 回复此楼
7楼2015-02-24 22:00:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 shuanglei 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿2110，化学学硕310分，本科重点双非求调剂 +17 努力奋斗112 2026-04-08 17/850 2026-04-09 23:24 by 大力水手力大无�
[考研] 一志愿华工085600 331分 +6 天下ww 2026-04-09 6/300 2026-04-09 18:59 by l_paradox
[考研] 电子信息270求调剂 +11 terminal469 2026-04-07 11/550 2026-04-09 18:28 by hy861222
[高分子] 环氧灌封胶 抗沉剂 +3 huagong159 2026-04-03 3/150 2026-04-09 17:11 by wei12shua
[考研] 368化学求调剂 +13 wwwwabcde 2026-04-07 14/700 2026-04-09 14:47 by heaven_jay
[考研] 302分求调剂 +4 凡语祈愿 2026-04-08 5/250 2026-04-08 22:03 by 土木硕士招生
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 材料调剂 +4 hzhahg 2026-04-06 4/200 2026-04-08 16:20 by luoyongfeng
[考研] 一志愿南昌大学，085600，344分求调剂 +11 调剂上岸玘 2026-04-05 12/600 2026-04-08 16:17 by luoyongfeng
[考研] 求调剂 +9 月@163.com 2026-04-07 11/550 2026-04-08 14:48 by qlm5820
[考研] 304求调剂 +16 c297914 2026-04-05 17/850 2026-04-08 13:00 by grayjzr
[考研] 336材料与化工085600求调剂 +19 水星记infp 2026-04-05 22/1100 2026-04-07 21:11 by yongzhesheng
[考研] 287求调剂 +3 通信学硕081000 2026-04-03 4/200 2026-04-06 21:03 by going home
[考研] 08专硕275调剂 +5 AaAa7420 2026-04-05 5/250 2026-04-05 18:01 by jkddd
[考研] 306分材料与化工求调剂 +7 黎吧啦啦你很有� 2026-04-03 7/350 2026-04-05 17:18 by Hdyxbekcb
[考研] 320求调剂 +3 一样圆 2026-04-04 3/150 2026-04-04 22:29 by 啵啵啵0119
[考研] 材料383求调剂 +5 郭阳阳阳成 2026-04-04 5/250 2026-04-04 19:06 by dongzh2009
[考研] 309求调剂 +4 快乐的小白鸽 2026-04-04 5/250 2026-04-04 15:55 by cql1109
[考研] 本科985，专业0812分336求调剂 +4 莫莫很行 2026-04-03 4/200 2026-04-03 21:31 by zhq0425
[考研] 一志愿重庆大学085404，总分314分，求调剂 +4 zf83hn 2026-04-03 4/200 2026-04-03 21:25 by 啵啵啵0119
信息提示
关闭
请填处理意见
关闭