ld = 50*10^-4;%%%%%
a1 = 28*10^-2;
a2 = 20*10^-1;
k0 = pi./ld;
ea = 3.7;
em = 2;
k1 = @(x) sqrt(ea*k0^2 - x.^2);
k2 = @(x) sqrt(em*k0^2 - x.^2);
f = @(x) cos(k1(x)*a1).*cos(k2(x)*a2)-(em^2*k1(x).^2 + ea^2*k2(x).^2)./(2*em*ea*k1(x).*k2(x)).*sin(k1(x)*a1).*sin(k2(x)*a2);
z = fzero(f,-1)
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ：24596024