关于f(lambda,neff)=0的隐函数，想得到横轴为lambda，纵轴为Re(neff)之间的关系。（即波长与有效折射率实部之间的关系）

clc, clear all;
c=3e17;    %常量光速c
e0=1.6e-19;    %常量
e00=3.7;    %常量
ed=1;      %常量
w=50e-9;    %常量
r=0.018*e0;    %常量
wp=9.1*e0;   %常量
lambda=600:0.3:1200; %横轴的波长
omiga=2*pi*c./lambda;  %频率
em=e00-(wp.^2)./(omiga.^2+j.*omiga*r);  %em与频率之间的关系
f=em.*sqrt(neff.^2-ed).*tanh(w*pi.*sqrt(neff.^2-ed)./lambda)+ed.*sqrt(neff.^2-em);  %关于横轴波长，纵轴有效折射率的函数f(lambda,neff)=0



下面的程序怎么写，我想画出横轴为波长lambda[600,1200]，纵轴为有效折射率的实部neff[0,5]的图像，也就是图中蓝色的那条线

不知道 j 是什么？虚数单位吗？

1. clc, clear all;
   
2. c = 3e17;    %常量光速c
   
3. e0 = 1.6e-19;    %常量
   
4. e00 = 3.7;    %常量
   
5. ed = 1;      %常量
   
6. w = 50e-9;    %常量
   
7. r = 0.018*e0;    %常量
   
8. wp = 9.1*e0;   %常量
   
9. omiga = @(lambda)2*pi*c./lambda;  %频率
   
10. em = @(lambda)e00-(wp^2)./(omiga(lambda).^2+j*omiga(lambda)*r);  %em与频率之间的关系
    
11. fun = @(lambda,neff)em(lambda).*sqrt(neff.^2-ed).*tanh(w*pi*sqrt(neff.^2-ed)./lambda)+ed*sqrt(neff.^2-em(lambda));  %关于横轴波长，纵轴有效折射率的函数f(lambda,neff)=0
    
12. fimplicit(fun,[600,1200])

复制代码

谢中华 发表于 2022-5-15 08:48
不知道 j 是什么？虚数单位吗？

j是虚部单位，我运行后好像没有图像，我想画的是lambda与neff的实部的图像，那这个取实部的real加在哪里了？