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

shuanglei

銅蟲 (小有名氣)

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

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

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

» 猜你喜歡

362求調(diào)劑已經(jīng)有14人回復(fù)
求調(diào)劑 302分初試 0854 已經(jīng)有4人回復(fù)
299求調(diào)劑已經(jīng)有4人回復(fù)
266分，一志愿電氣工程，本科材料，求材料專業(yè)調(diào)劑已經(jīng)有3人回復(fù)
312求調(diào)劑已經(jīng)有4人回復(fù)
315求調(diào)劑已經(jīng)有7人回復(fù)
土木304求調(diào)劑已經(jīng)有3人回復(fù)
316求調(diào)劑已經(jīng)有16人回復(fù)
一志愿華東理工大學(xué)，080500學(xué)碩，317分，求調(diào)劑已經(jīng)有13人回復(fù)
材料與化工306分找調(diào)劑已經(jīng)有9人回復(fù)

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

大家?guī)兔砜匆幌逻@個反應(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ī)兔匆幌逻@個offer是否合乎要求，新手散光所有金幣求助~~ 已經(jīng)有12人回復(fù)
大家?guī)兔匆幌? 這個反應(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ī)兔匆幌�，好奇怪，我的電泳圖，三個加樣孔在分離膠合在一起了已經(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ī)兔匆幌逻@個LC/MS/MS的問題已經(jīng)有9人回復(fù)
【分享】大家?guī)兔匆幌逻@個Ni的峰歸屬于NiO還是Ni2O3 已經(jīng)有7人回復(fù)

硬著頭皮鉆下去？！

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

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

sharinfang

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
帖子: 11
在線: 29.3小時
蟲號: 1727628
注冊: 2012-03-30
性別: GG
專業(yè): 光子與光電子器件

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

送鮮花一朵

能不能告訴我的你的qq，我有一個超越方程不知道怎樣編程

贊一下

回復(fù)此樓

5樓2012-10-16 18:00:16

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

查看全部 7 個回答

zjliu1265

至尊木蟲 (正式寫手)

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

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

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

圖形如下，可以看出只有一個解：

贊一下(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小時
蟲號: 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小時
蟲號: 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)只有一個解，
并非多解：
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ā)表

[考研] 環(huán)境工程297分求調(diào)劑一志愿杭高院 +13	GENJIOW 2026-03-31	14/700	2026-04-02 16:51 by chran16
[考研] 346求調(diào)劑 +5	鄭誠樂 2026-04-02	5/250	2026-04-02 16:38 by SZW_UJN
[考研] 286分調(diào)劑 +20	Faune 2026-03-30	22/1100	2026-04-02 13:24 by clyblh
[考研] 08工科求調(diào)劑290分 +5	1314捧花 2026-04-02	8/400	2026-04-02 13:16 by 喬噠噠噠
[考研] 085900土木水利336分求調(diào)劑 +4	Zhangjiangj 2026-03-31	6/300	2026-04-02 11:40 by 1753564080
[考研] 0856材料與化工調(diào)劑，339 +13	10213207 2026-03-31	13/650	2026-04-02 11:35 by michael2011
[考研] 085600，材料與化工321分，求調(diào)劑 +13	大饞小子 2026-03-27	13/650	2026-04-02 10:48 by sanrepian
[考研] 266求調(diào)劑 +10	陽陽哇塞 2026-04-01	10/500	2026-04-02 09:07 by vgtyfty
[考研] 材料求調(diào)劑 +8	呢呢妮妮 2026-04-01	8/400	2026-04-02 07:13 by yjolah
[考研] 285求調(diào)劑 +11	AZMK 2026-04-01	11/550	2026-04-01 22:40 by peike
[考研] 309求調(diào)劑 +19	誰不是少年 2026-03-29	19/950	2026-04-01 15:47 by jp9609
[考研] 一志愿鄭大085600，310分求調(diào)劑 +6	李瀟可 2026-03-26	6/300	2026-04-01 14:44 by chenqifeng666
[考研] 求調(diào)劑 +5	零八# 2026-03-27	5/250	2026-04-01 14:40 by yulian1987
[碩博家園] 博一被送出聯(lián)培感覺不適應(yīng)怎么辦 +3	全村的狗 2026-03-31	3/150	2026-04-01 10:44 by 328838485
[考研] 土木304求調(diào)劑 +3	兔突突突， 2026-03-31	3/150	2026-04-01 09:42 by JourneyLucky
[考研] 求調(diào)劑生物學(xué) 377分 +6	zzll03 2026-03-31	6/300	2026-03-31 17:33 by 唐沐兒
[考研] 080500-315分復(fù)試調(diào)劑 +9	上岸3821 2026-03-31	9/450	2026-03-31 17:29 by 唐沐兒
[考博] 材料專業(yè)申博 +5	杜雨婷dyt 2026-03-29	5/250	2026-03-31 11:19 by oooqiao
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +5	小張zxy張 2026-03-26	10/500	2026-03-31 00:29 by jp9609
[考研] 一志愿上海理工能源動力（085800）310分求調(diào)劑 +3	zhangmingc 2026-03-27	4/200	2026-03-27 19:01 by 給你你注意休息