FDTD Script 远场图-farfield3dintegrate/farfieldspherical

最新推荐文章于 2025-03-26 17:04:00 发布

zhdy56789

最新推荐文章于 2025-03-26 17:04:00 发布

阅读量642

点赞数 2

分类专栏： FDTD 文章标签：笔记

原文链接：https://blog.csdn.net/jxxl_1314/article/details/130329297

版权

FDTD 专栏收录该内容

16 篇文章

订阅专栏

FDTD Script命令学习-farfield3dintegrate/farfieldspherical

本文介绍了FDTD仿真中的两个关键命令：farfield3dintegrate用于在锥形区域内对3D远场投影数据进行积分，而farfieldspherical则将远场数据转换为球坐标系。这两个命令涉及的角度参数以度为单位但在计算中转换为弧度。示例代码展示了如何使用这些命令来计算特定锥形区域内的能量比例以及生成球坐标系下的电场分布图像。
摘要由CSDN通过智能技术生成
目录

一、farfield3dintegrate命令介绍
（1）语法
（2）示例
二、farfieldspherical命令介绍
（1）语法

（2）示例

一、参考FDTD官网farfield3dintegrate命令介绍：

https://optics.ansys.com/hc/en-us/articles/360034410174-farfield3dintegrate-Script-command

1. farfield3dintegrate命令功能:在3D模拟中，在以立体角theta0(与z轴的夹角)和方位角phi0(与x轴的夹角) 为中心、在指定的半角宽度的锥形区域下的远场投影积分。远场电场是方向余弦（ux,uy）的函数，但是farfileld3dintegrate会自动对变量进行变更（根据下一句话，应该是进行单位换算）。类似的，角度是指定以度为单位，但是进行积分计算前会转化为弧度。不同监视器取向ux,uy,na,nb的介绍可以参考farfield3d命令介绍。

（1）语法：out=farfield3dintegrate(E2,ux,uy,halfangle,theta0,phi0)；见上述积分公式，该语法是对3D远场投影数据的积分。

（2）示例

计算光源达到中心角theta=phi=0,30°的锥形远场区域的能量比例。

代码如下：
m="monitor1";
res = 201;
E2 = farfield3d(m,1,res,res); #返回远场投影的电场强度数据，结果为201*201阶矩阵；
ux = farfieldux(m,1,res,res);#返回结果为201*201阶矩阵，方向矢量；
uy = farfielduy(m,1,res,res);
halfangle=30;#半角大小为30°
theta0=0;
phi0=0;
cone_30 = farfield3dintegrate(E2, ux, uy, halfangle, theta0, phi0); # integrate over 30 degree cone，采用上述公式进行积分计算；
total = farfield3dintegrate(E2, ux, uy); # integrate over entire hemisphere，对半角大小为90°的半球电场强度进行积分计算；
T = transmission(m); # fraction of source power transmitted into far field，返回监视器monitor1中的透过率数据；
?cone_30/total; # fraction of far field power within a 30 degree cone，显示投影到远场中半角大小为30°电场强度与整个半球远场区域的电场强度比例；
?cone_30/total*T; # fraction of source power transmitted into the far field within a 30 degree cone，显示投影到远场中半角大小为30°电场强度与整个光源的强度比例；
附：transmission命令可以参考：https://optics.ansys.com/hc/en-us/articles/360034405354-transmission-Script-command

二、参考FDTD官网farfieldspherical命令介绍：

https://optics.ansys.com/hc/en-us/articles/360034410194-farfieldspherical-Script-command

1. farfieldspherical命令功能：将3D模拟的远场电场E（ux,uy)数据转换到球坐标系下一位阵列的电场E（theta,phi）中。远场投影函数一般会返回关于ux,uy（方向余弦）函数的投影。farfieldspherical用来将这个远场投影数据转化为到更加普遍的theta,phi单位中（功能即坐标系的转化：三维笛卡尔坐标系转换为球坐标系）。不同监视器取向ux,uy,na,nb的介绍可以参考farfield3d命令介绍。

（1）语法：out=farfieldspherical(E2,ux,uy,theta,phi)；将远场数据转化为球坐标系下的数据。输出结果为维度为（M*N,1）（其中M和N分别为nb和na的投影分辨率）；

（2）示例1：生成E2_far关于立体角theta的图片，方位角phi=0；

代码如下：

m="Monitor1"; # Monitor name
res = 201; # projection resolution
E2 = farfield3d(m,1,res,res);
ux = farfieldux(m,1,res,res);
uy = farfielduy(m,1,res,res);
theta = linspace(-90,90,100); #theta角取值范围[-90,90],等间隔取100个点位
phi = 0;
plot(theta, farfieldspherical(E2,ux,uy,theta,phi) ,"theta", "E^2", "E^2 at phi=0");#以theta为x轴，farfieldspherical为轴，"theta"为x轴名称，“E^2”为y轴名称，“E^2 at phi=0”为图片名称，进行作图
示例2：将厂数据转化为以theta和phi角网格坐标系下。

代码如下：

theta = linspace(-90,90,10);
phi = linspace(0,45,11);
Theta = meshgridx(theta,phi); #将theta数据转化为11*10阶矩阵，其中矩阵中每行数据中元素为原始theta矩阵中的元素，且每行元素均相同
Phi = meshgridy(theta,phi);#将phi数据转化为11*10阶矩阵
E2_angle = farfieldspherical(E2,ux,uy,Theta,Phi);
E2_angle = reshape(E2_angle, [length(theta), length(phi)]);#将E2_angle矩阵11*10转换为10*11阶矩阵
image(theta, phi, E2_angle, "theta","phi","E2");#以theta为x轴，phi为y轴，对z进行作图，“theta”为x轴名称，“phi”为y轴名称，“E2”为图片名称
（3）其它参考

List of commands , farfield3d , farfieldux , farfielduy , Far field projections - Direction unit vector coordinates , meshgridx , meshgridy；
————————————————

版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

原文链接：https://blog.csdn.net/jxxl_1314/article/details/130329297