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

[复制链接]
旅行666666 发表于 2022-5-14 11:10:21
本帖最后由 旅行666666 于 2022-5-14 20:02 编辑

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]的图像,也就是图中蓝色的那条线

波长lambda与有效折射率实部neff的图像

波长lambda与有效折射率实部neff的图像

代码中的em,要带入到函数f中的

代码中的em,要带入到函数f中的

隐函数f(lambda,neff)

隐函数f(lambda,neff)

2 条回复


谢中华 发表于 2022-5-15 08:48:31
不知道 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])
复制代码

旅行666666 发表于 2022-5-17 12:18:26
谢中华 发表于 2022-5-15 08:48
不知道 j 是什么?虚数单位吗?

j是虚部单位,我运行后好像没有图像,我想画的是lambda与neff的实部的图像,那这个取实部的real加在哪里了?
您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
相关文章
热门教程
站长推荐
快速回复 返回顶部 返回列表