matlab求光纤直径,MATLAB计算微纳光纤有效折射率和论文不符

仿照欧老师109页的程序对微纳光纤的有效折射率进行了计算，TE01模特征方程简化后比较简单，算出来和论文相符，但计算HE1m模式和EH模时特征方程复杂，init参量不知道怎么设置，和论文中的曲线不符，还有一个问题想请教，HE模式和EH模式满足的特征方程相同，那是如何区别两种模式的有效折射率的？特征方程如下：下面是我的程序，

atemp=a(j);

Vtemp=k_0*atemp*(n_1^2-n_2^2)^0.5; %归一化频率V

Utemp=NaN;

i=0;

while (isnan(Utemp) && i

init=3.5*(N-i)/N;

try

Utemp=fzero(@(Utemp) ... %调用fzero求解本征方程

((besselj(0,Utemp)-besselj(1,Utemp)/Utemp)/(Utemp*besselj(1,Utemp))+ ...

(besselk(0,sqrt(Vtemp^2-Utemp^2))-besselk(1,sqrt(Vtemp^2-Utemp^2)))/(sqrt(Vtemp^2-Utemp^2)*besselk(1,sqrt(Vtemp^2-Utemp^2))))* ...

((besselj(0,Utemp)-besselj(1,Utemp)/Utemp)/(Utemp*besselj(1,Utemp))+ ...

(besselk(0,sqrt(Vtemp^2-Utemp^2))-besselk(1,sqrt(Vtemp^2-Utemp^2)))/(n_1^2*sqrt(Vtemp^2-Utemp^2)*besselk(1,sqrt(Vtemp^2-Utemp^2))))- ...

(1-Utemp^2/(atemp^2*k_0^2*n_1^2))*(1/Utemp^2+1/(Vtemp^2-Utemp^2))^2,init);

catch

end

i=i+1;

end

U(j)=Utemp; %归一化横向相位参数U

end

beta=sqrt(k_0.^2.*n_1.^2-(U./a).^2); %传播常数beta

n_eff=beta./k_0; %模式有效折射率

figure

hold on

axis([0.2 2 1 1.45])

plot(a,n_eff,'r:');

xlabel('半径(um)')

ylabel('有效折射率')

title('不同半径微纳光纤的模式有效折射率')

grid on

toc

论文中的曲线如下：

2016-5-5 21:26 上传

特征方程

2016-5-5 21:29 上传

论文中曲线

2016-5-5 21:30 上传

我画出来的

matlab求光纤直径,MATLAB计算微纳光纤有效折射率和论文不符相关推荐

大气复折射率matlab,MATLAB计算微纳光纤有效折射率和论文不符
仿照欧老师109页的程序对微纳光纤的有效折射率进行了计算,TE01模特征方程简化后比较简单,算出来和论文相符,但计算HE1m模式和EH模时特征方程复杂,init参量不知道怎么设置,和论文中的曲线不符, ...
浙大计算机学院博士张磊,浙江大学张磊、童利民教授团队《ACS AMI》：用于感知组织硬度的微纳光纤触觉传感器...
随着微创精密手术和手术机器人的不断发展,病人康复质量显著提高.然而,微创手术中的触觉缺失或触觉失真导致微创手术存在偶发性器官破坏风险.准确获取组织器官的触觉信息不仅能够降低微创手术风险,还能够显著提升 ...
matlab求刚度,matlab直接刚度法计算结构频率
matlab直接刚度法计算结构频率 syms E I K r l x cja1=-E*I*K*r^3*(cos(x)*sinh(x)+sin(x)*cosh(x))jc1=-E*I*K*r^2*sin ...
matlab求近似解,matlab求近似解
. 实验一方程根的近似计算一.问题求非线性方程的根二.实验目的 1.学会使用 matlab 中内部函数 roots.solve.fsolve.fzero 求解方程,并用之解决实际...... ...
matlab求两点间距离,matlab如何求一个N*2的矩阵的任意两点间的距离？
matlab如何求一个N*2的矩阵的任意两点间的距离? mip版关注:64 答案:3 悬赏:0 解决时间 2021-02-23 20:55 已解决 2021-02-23 16:37 假设有个矩 ...
matlab求点,MATLAB求两点中点
已知两点坐标,连接这两点,求线段的中点坐标.有没有什么公式啊? A(x1,y1).B(x2,y2),则AB中点是M((x1+x2)/2,(y1+y2)/2) 如何用MATLAB做出一个求两点间最短距离 ...
matlab求心率,matlab心率分析
[实例简介] 国外大学的一款matlab开发软件,用于心率分析.有源代码,是学习的好参考. [实例截图] [核心代码] matlab心率分析 └── matlab心率分析 └── HRVAS_v1.0 ...
matlab 求概率密度,MATLAB如何使用pdf函数计算指定分布的概率密度函数
MATLAB如何使用pdf函数计算指定分布的概率密度函数 [语法说明] Y=pdf('name',X,A) Y=pdf('name',X,A,B) Y=pdf('name',X,A,B,C) 字符串n ...
怎么用matlab求特征向量,MATLAB用eig()函数求【特征值】【特征向量】【归一化
在MATLAB中,计算矩阵A的特征值和特征向量的函数是eig(A),常用的调用格式有 5种: (1) E=eig(A):求矩阵A的全部特征值,构成向量E. 想求最大特征值用:max(eig(A))就好 ...

matlab求光纤直径,MATLAB计算微纳光纤有效折射率和论文不符

matlab求光纤直径,MATLAB计算微纳光纤有效折射率和论文不符相关推荐

最新文章

热门文章