| 查看: 1165 | 回复: 0 | |||
[交流]
【求助】MATLAB求解光纤包层模有效折射率,请教matlab语法问题
|
|
clc; clear; a1=2.625; %纤芯半径 a2=62.5; %包层半径 n1=1.458; %纤芯折射率 n2=1.450; %包层折射率 n3=1.0; %空气折射率 lamda=1.55; 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); g1=g0; 求解neff（光纤各阶包层模的有效折射率，其中neff有约一百多个解）。我对MATLAB语法不懂，只把方程（即光纤包层模的本征方程）用MATLAB写了出来。请大虾帮忙解一下。（精确到小数点后8位） [ Last edited by GrasaVampiro on 2010-8-27 at 08:39 ] |
回复此楼» 猜你喜欢
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
申请科研项目纯靠实力能成功吗？
已经有7人回复
-
极限拉扯，如何抉择
已经有14人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
找到一些相关的精华帖子，希望有用哦~
求助！光纤的相关问题~~
已经有9人回复
- 请教2-氨基噻唑与2,4-二氯嘧啶的4位氯缩合
已经有7人回复
- 关于matlab线性规划 背包问题的 解法 求救呀~~方程知道 但不知道如何解答····
已经有4人回复
- 请教matlab反应动力学参数估计遇到的问题，谢谢
已经有15人回复
- 请教关于matlab矩阵转换问题
已经有6人回复
- 怎么样设置matlab在excel中的输出位置
已经有10人回复
- 请教Matlab图形保存问题
已经有5人回复
- matlab的函数调用问题
已经有3人回复
- 咨询光纤外部（周围）折射率与光纤本身的有效折射率之间的联系与区别，
已经有3人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 【求助】利用MATLAB求解的外围轮廓包络线为什么不与图形外围轮廓重合？怎样才可以精确
已经有8人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
10
