人民币奖励,答出来2天内款到拜托各位了3Q方程如下f所示,这个方程要求的是在wl知道的情况下求解neff,而对于任意一个波长wl,neff有多个复数解.具体参数如下：n2=1.4628;n3=1.48;a2=62.5e-6;m=1;wl=1550e-

Question

人民币奖励,答出来2天内款到拜托各位了 3Q
方程如下f所示,这个方程要求的是在wl知道的情况下求解neff,而对于任意一个波长wl,neff有多个复数解.具体参数如下：n2=1.4628;n3=1.48;a2=62.5e-6;m=1;wl=1550e-9;这个方程的关键在于n3大于n2,如果n3小于n2,可以用二分法求出neff的多个实数根,但当n3>n2后有多个复数根,想了几天没有想出来,书也翻了N本,还是没有头绪,希望高手们都来讨论一下.function f=HE_n3dayun2gs(wl,neff)format longk0=2*pi/wl;beta_cl=neff*k0;u2=a2*k0*sqrt(n2^2-neff^2);w2=a2*k0*sqrt(neff^2-n3^2);v=sqrt(u2^2+w2^2);left=besselj(m-1,u2)/(u2*besselj(m,u2));besselk_differential=-besselk(m-1,w2)-(m/w2)*besselk(m,w2); %besselk的微分right1=-(n2^2+n3^2)*besselk_differential/(2*n2^2*w2*besselk(m,w2));right2=m/u2^2;right3=(n2^2-n3^2)*besselk_differential/(2*n2^2*w2*besselk(m,w2));right4=m*v^2*beta_cl/(u2^2*w2^2*n2*k0);right=right1+right2-sqrt(right3^2+right4);f=left-right;如有高手答出来的,人民币奖励,2天内到账,不过由于目前还是学生,故数额不多,但至少够一个人去吃一个普通的大餐!同时也欢迎大家一起进来讨论!

Answer 1

也期待matlab的结果.如果用Forcal解,要用到IMSL库的贝塞尔函数,不知与matlab的贝塞尔函数的函数名是如何对应的?说说贝塞尔函数的意义,如：besslj：第一类0阶Bessel函数.或者：修正的第一类1阶Bessel函数.或者：阶数为整数,参数为复数(x,y).或者其他besselk：... ... 查看原帖>>