當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

shuanglei

超級版主

應(yīng)助: 0 (幼兒園)
金幣: 56.5
散金: 8
紅花: 1
帖子: 130
在線: 37.7小時(shí)
蟲號: 1051296
注冊: 2010-07-03
性別: GG
專業(yè): 光學(xué)

[交流] 【求助】關(guān)于長周期光柵有效折射率求解，大家?guī)兔匆幌? 已有5人參與

以下是文獻(xiàn)中求解光纖包層模有效折射率的方程，即neff為未知數(shù)，其他參量都已知。（文獻(xiàn)中說neff有很多解）大家?guī)兔匆幌拢鯓硬拍芮蟮胣eff的多個(gè)解，請各位大俠指點(diǎn)，萬分感謝！
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都是復(fù)數(shù)
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

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

sharinfang

» 猜你喜歡

292分，材料與化工，申請調(diào)劑已經(jīng)有13人回復(fù)
316求調(diào)劑已經(jīng)有7人回復(fù)
材料專碩322分已經(jīng)有4人回復(fù)
286分調(diào)劑已經(jīng)有16人回復(fù)
337求調(diào)劑已經(jīng)有8人回復(fù)
一志愿北交大材料工程總分358 已經(jīng)有7人回復(fù)
362求調(diào)劑已經(jīng)有13人回復(fù)
332求調(diào)劑已經(jīng)有7人回復(fù)
342求調(diào)劑已經(jīng)有11人回復(fù)
319求調(diào)劑已經(jīng)有4人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

大家?guī)兔砜匆幌逻@個(gè)反應(yīng)，按照文獻(xiàn)的方法怎么都做不出... 已經(jīng)有50人回復(fù)
請大家?guī)兔匆幌潞舜抛V吧！已經(jīng)有5人回復(fù)
ABAQUS安裝成功后，運(yùn)行出現(xiàn)問題，請大家?guī)兔匆幌?/a> 已經(jīng)有14人回復(fù)

大家?guī)兔匆幌�，這樣的回復(fù)是什么意思呢？已經(jīng)有16人回復(fù)
【急】請大家?guī)兔匆幌逻@個(gè)offer是否合乎要求，新手散光所有金幣求助~~ 已經(jīng)有12人回復(fù)
大家?guī)兔匆幌? 這個(gè)反應(yīng)為什么不能進(jìn)行已經(jīng)有6人回復(fù)
求助啊，SCI修稿后一月有余，詢問編輯，給了一封回信，大家?guī)兔匆幌率鞘裁匆馑及?/a> 已經(jīng)有24人回復(fù)

大家?guī)兔匆幌翫SC測試圖已經(jīng)有8人回復(fù)
初寫論文菜鳥，大家?guī)兔匆幌逻@幾張圖是用什么軟件畫的或生成的？已經(jīng)有10人回復(fù)
大家?guī)兔匆幌�，好奇怪，我的電泳圖，三個(gè)加樣孔在分離膠合在一起了已經(jīng)有7人回復(fù)
請大家?guī)兔匆幌�，套磁的教授回信，該怎么回�?fù) 已經(jīng)有18人回復(fù)
審稿意見回來了，大家?guī)兔匆幌��！?/a> 已經(jīng)有9人回復(fù)

【求助】請求大家?guī)兔匆幌聦徃迦说幕貜?fù) 已經(jīng)有6人回復(fù)
【求助/交流】大家?guī)兔匆幌逻@是什么菌已經(jīng)有55人回復(fù)
【求助】大家?guī)兔匆幌挛易龅募ぐl(fā)態(tài) 已經(jīng)有8人回復(fù)
【求助】大家?guī)兔匆幌逻@個(gè)LC/MS/MS的問題已經(jīng)有9人回復(fù)
【分享】大家?guī)兔匆幌逻@個(gè)Ni的峰歸屬于NiO還是Ni2O3 已經(jīng)有7人回復(fù)

硬著頭皮鉆下去？！

1樓 2010-07-13 08:32:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wawacgh

主管區(qū)長

應(yīng)助: 0 (幼兒園)
金幣: 33
帖子: 3
在線: 5.1小時(shí)
蟲號: 3693859
注冊: 2015-02-22
專業(yè): 光學(xué)

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

樓主，您好！這個(gè)問題您解決了嗎？您知道如何算纖芯與同向包層模式的耦合系數(shù)嗎？

贊一下

回復(fù)此樓

7樓2015-02-24 22:00:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 7 個(gè)回答

zjliu1265

實(shí)習(xí)版主

應(yīng)助: 18 (小學(xué)生)
金幣: 12051.2
散金: 200
帖子: 642
在線: 679.6小時(shí)
蟲號: 585320
注冊: 2008-08-08
性別: GG
專業(yè): 光學(xué)

★ ★ ★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
adu886886(金幣+3):鼓勵(lì)應(yīng)助！歡迎常來仿真模擬板塊！ 2010-07-14 07:54:48

我覺得你可以考慮用2分法解，我畫了neff~g1-g0的曲線，發(fā)現(xiàn)只有一個(gè)解，
并非多解：
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都是復(fù)數(shù)
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');

圖形如下，可以看出只有一個(gè)解：

贊一下(2人)

回復(fù)此樓

2樓2010-07-14 07:51:23

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

shuanglei

兌換貴賓

應(yīng)助: 0 (幼兒園)
金幣: 56.5
散金: 8
紅花: 1
帖子: 130
在線: 37.7小時(shí)
蟲號: 1051296
注冊: 2010-07-03
性別: GG
專業(yè): 光學(xué)

請問neff=linspace(a1,a2,201);中的201是什么意思。

贊一下

回復(fù)此樓

硬著頭皮鉆下去？！

3樓2010-08-29 15:43:56

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

shuanglei

版主

應(yīng)助: 0 (幼兒園)
金幣: 56.5
散金: 8
紅花: 1
帖子: 130
在線: 37.7小時(shí)
蟲號: 1051296
注冊: 2010-07-03
性別: GG
專業(yè): 光學(xué)

引用回帖:

Originally posted by zjliu1265 at 2010-07-14 07:51:23:
我覺得你可以考慮用2分法解，我畫了neff~g1-g0的曲線，發(fā)現(xiàn)只有一個(gè)解，
并非多解：
clc;close all;clear
a1=2.625; %纖芯半徑
a2=62.5; %包層半徑
n1=1.458; %纖芯折射率
n2=1.450; %包 ...

請問neff=linspace(a1,a2,201);為什么只讓他生成201維的行向量？

贊一下

回復(fù)此樓

硬著頭皮鉆下去？！

4樓2010-08-29 16:19:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 273求調(diào)劑 +16	李芷新1 2026-03-31	16/800	2026-04-01 17:01 by cq2548
[考研] 0703一志愿南師大334求調(diào)劑 +4	seven7yu 2026-03-30	4/200	2026-04-01 16:10 by oooqiao
[考研] 材料科學(xué)與工程339求調(diào)劑 +10	hyz0119 2026-03-31	11/550	2026-04-01 16:00 by hyz0119
[考研] 0856初試324分求調(diào)劑 +5	想上學(xué)求調(diào) 2026-04-01	5/250	2026-04-01 15:28 by 考研上岸呀一定?/a>
[考研] 086000生物與醫(yī)藥初試274求調(diào)劑 +5	小叮當(dāng)來了 2026-03-30	6/300	2026-04-01 14:25 by mxj1203
[考研] 279求調(diào)劑 +6	莫xiao 2026-04-01	6/300	2026-04-01 11:11 by 逆水乘風(fēng)
[考研] 0710生物學(xué)求調(diào)劑！ +6	敘述文 2026-03-31	6/300	2026-04-01 09:39 by JourneyLucky
[考研] 080500-315分復(fù)試調(diào)劑 +9	上岸3821 2026-03-31	9/450	2026-03-31 17:29 by 唐沐兒
[考研] 考研調(diào)劑求助 +7	13287130938 2026-03-31	7/350	2026-03-31 16:39 by 690616278
[考研] 085601 材料工程 313分求調(diào)劑 +6	Ong3 2026-03-27	6/300	2026-03-31 15:18 by yanflower7133
[考研] 求調(diào)劑 +8	11ggg 2026-03-30	8/400	2026-03-31 13:56 by nanaliuyun
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +6	星沉uu 2026-03-26	7/350	2026-03-31 10:19 by GdShizy
[考研] 本科211總分289，08工學(xué)真心求調(diào)劑 +3	utopiaE 2026-03-30	3/150	2026-03-30 23:42 by ms629
[考研] 哈爾濱工業(yè)大學(xué)材料與化工專碩378求調(diào)劑 +3	塔比烏斯 2026-03-30	3/150	2026-03-30 22:55 by 無際的草原
[有機(jī)交流] 考研調(diào)劑 +8	watb 2026-03-26	8/400	2026-03-30 18:40 by 544594351
[考研] 356求調(diào)劑 +3	gysy?s?a 2026-03-28	3/150	2026-03-29 00:33 by 544594351
[考研] 本科新能源科學(xué)與工程，一志愿華理能動(dòng)285求調(diào)劑 +7	AZMK 2026-03-28	11/550	2026-03-28 21:01 by xxxsssccc
[考研] 312，生物學(xué)求調(diào)劑 +3	小譯同學(xué)abc 2026-03-28	3/150	2026-03-28 15:32 by 落�？伤�
[考研] 081200-314 +3	LILIQQ 2026-03-27	4/200	2026-03-28 09:41 by 保護(hù)地球你我做?/a>
[考研] 341求調(diào)劑 +7	青檸檬1 2026-03-26	7/350	2026-03-27 00:19 by wxiongid