如果有问题的话，根据图像中的公式解，我就根据它们写了程序，但是很难解。 感谢大神看问题出在哪里。 步骤如下。
clear clc
syms Bta
a=3.04e-6，b=62.5e-6
n1=1.451，n2=1.444，n3=1.333
k0=2*pi/(400e-9 )、k1=2*pi/) 1550e-9 )、k2=2*pi/) 1551e-9 ) )。
v1=a*sqrt((BTA^2)-)-(n3^2) ) *(k0^2) )
u=a*sqrt((N1^2) ) *(k0^2)-(Bta^2) )
w=b*sqrt () BTA^2)-)-(n2^2) ) *(k0^2) )
v=b*sqrt () BTA^2)-)-(n3^2) ) *(k0^2) )
s=a/b
t=(besselj(0，u )u*besselj ) 1，u ) ) besselk ) 0，) w*s ) *besselk(1) )，) w*s ) )
T1=(besselk(0，v )/) v*besselk ) 1，v ) ) besseli(0) 0，w )w*besseli ) 1，w ) )
T2=(besselj(0，u )u*besselj ) 1，u ) ) besseli(0(0) 0，(w*s ) besseli )，) w*s ) )
T3=(besselk(0，v )/) v*besselk ) 1，v ) )-besselk(0) 0，w )/(w * ) besselk ) 1，w ) )
T4=(besselk(1，w ) *besseli ) 1，) w*s ) )/) besselk ) 1，) w*s ) *besseli ) 1，w ) )
DJ=diff(Besselj(1，u )，Bta ) ) ) ) ) ) ) ) ) )。
dk=diff(Besselk(1，v1 )，Bta ) ) ) ) ) ) ) ) ) ) ) )。
T5=(DJ/) u*besselj(1，u ) ) dk/) V1* ) besselk ) 1，v1 ) )
T6=() k2^2) ) DJ/) u*besselk(1，v1 ) ) ) dk/) V1*besselk ) ) )
T7=(BTA )2) ) ) ()1/(v1 ^2) )2) )
BTA=solve () (t*T1 )/(t2*t3 )=t4 )、) t5*t6=t7 ) )
11.PNG
(66.87 KB，下载次数： 0 )
2018-4-30 14:58上传
12.PNG
(118.82 KB，下载次数： 0 )
2018-4-30 14:58上传