| 查看: 1338 | 回复: 2 | |||
[交流]
利用Matlab处理Lumerical FDTD的三角纳米片电场分布仿真结果 已有2人参与
|
|
FDTD计算得到的电场分布，但是FDTD通过另存为jpg或者截屏所得到的图片分辨率很低，得到的图片往往不能直接使用。因此，可以通过脚本输入到Maltab，然后再利用Matlab处理图片并输出。 但是将数据从FDTD输出到Matlab中，并不是想象中那么简单，经历了好几次坑，反复摸索之后，得到了一种比较可行的方案，介绍如下。 1. FDTD原始结果 这里我们选用三角纳米片的电场分布仿真结果进行举例。图1是FDTD直接输出的结果（截图），可以明显看出，x方向和y方向的比例并不相同，而且不容易调节成比例尺相同，我目前有两个可行的方案可以解决这个问题，一种是先建立一个方形的图片，利用签字笔在电脑屏幕上画出方形的边界，然后再反复调节FDTD的图片，使其边界和画出的边界重合；第二种方案相对更精准也更方便，借用Snipaste截图软件创建一个方形的贴图，这个好处是这个贴图可以一直置于顶层，然后再调节FDTD中图片的边界即可。这两种方案都是调节好之后进行截图，因为直接另存，FDTD输出的图片更加模糊，而且也没有在FDTD Solutions软件中找到可以设置分辨率的选项（FDTD Solutions版本为2018版），因此可以使用高分辨率截图软件或者较高分辨率的软件，然后将图片放到PS或者AI中进行分辨率的进一步调节。FDTD中能够调节的着实比较少，很有必要继续调整。  2. FDTD结果导出到Matlab FDTD数据导入到Matlab主要参考Lumerical官网的介绍文档matlabsave。 可以通过如下脚本实现电场的读取： E=getresult("z=0","E"  ;E2=getelectric("z=0" ;E1=sqrt(E2); matlabsave("tri120.mat",E,E1); 其中E得到的是一个struct结果: 包含以下字段的 struct: E: [48841×3 double] Lumerical_dataset: [1×1 struct] f: 4.2827e+14 lambda: 7.0000e-07 x: [221×1 double] y: [221×1 double] z: 0 struct中E是一个包含三个x, y, z三个方向的电场分布，是矢量信息，而通过getelectric并开根号得到的是电场的标量数值。值得一说的是，对结构体struct中内容的调用可以采用E.x等命令进行直接读取。 值得注意的是，FDTD和Matlab采用的矩阵规则不同，输出后的结果需要进行非共轭转置，才能够和FDTD的结果进行对应。 x，y， E 等的信息除了利用getdata或getresult获取以外，也可以通过右键单击，将结果输出到脚本，可以得到相同的效果。  3. Matlab作图 将电场结果，x，y，z数据导出到matlab中之后，作图可以有多种方案，contour，contourf，pcolor等均可以实现，pcolor更为方便。  contourf和pcolor的语句如下： [C,h] = contourf(E,M) %counterf, M越大，越精细 h = pcolor(x,y,E) %pcolor 如图3为pcolor输出的结果。作图后，仍需进行如下调整： set(h,'LineStyle','none'); %线型设置为无 colorbar %色标 colormap jet %色标模式选用jet axis equal %x, y 等比例尺 最后，如果有仿真相关需求，欢迎通过微信公众号联系我们。 微信公众号：320科技工作室。 |
回复此楼» 猜你喜欢
焦虑
已经有4人回复
-
材料与化工（0856）304求B区调剂
已经有6人回复
-
深圳大学硕士招生（2026秋，传感器方向，仅录取第一志愿）
已经有7人回复
-
（081700）化学工程与技术-298分求调剂
已经有33人回复
-
0703化学求调剂
已经有7人回复
-
一志愿西南交大，材料专硕317求调剂
已经有3人回复
-
315求调剂
已经有9人回复
-
求调剂
已经有4人回复
-
293求调剂
已经有3人回复
-
求b区学校调剂
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
hustoe2011
至尊木虫 (文坛精英)
- 应助: 178 (高中生)
- 金币: 22855.1
- 散金: 9962
- 红花: 47
- 帖子: 22053
- 在线: 4425.9小时
- 虫号: 1333598
- 注册: 2011-06-29
- 性别: GG
- 专业: 光学和光电子材料
2楼2020-10-13 08:55:26
3楼2021-03-06 01:12:06
