| 查看: 1141 | 回复: 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人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有27人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有3人回复
-
2027广东省杰青
已经有9人回复
-
朋友圈看到的
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有110人回复
-
明天放榜？
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有11人回复
-
重要消息，中午系统在维护
已经有10人回复
-
投稿咨询
已经有6人回复
找到一些相关的精华帖子，希望有用哦~
求助！光纤的相关问题~~
已经有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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
