目录
一、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远场投影数据的积分。
参数 | 默认值 | 类型 | 描述 | |
E2 | 必须 | 矩阵 | 来自于farfield3d的电场数据 | |
ux | 必须 | 矢量 | 来自于farfieldux。注意这个结果是个矢量,因此只在一个频率点下执行farfieldux脚本命令就足够了。 | |
uy | 必须 | 矢量 | 来自于farfielduy。注意这个结果是个矢量,因此只在一个频率点下执行farfieldux脚本命令就足够了。 | |
halfangle | 可选 | 90 | 矢量 | 积分锥形区域的半角大小,单位为度,长度必须为L或1(维度相关长度?)。半角应该在0-90°之间。 |
theta0 | 可选 | 0 | 矢量 | 积分锥形区域的中心较theta,单位为度,长度必须为L或1。Theta0必须在0-90°之间。 |
phi0 | 可选 | 0 | 矢量 | 积分的锥形区域的中心角,单位为度,长度必须为L或1。phi0必须在0-360°之间。 |
(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的投影分辨率);
参数 | 数值 | 类型 | 描述 | |
E2 | 必须 | 矩阵 | 来自于farfield3d的电场数据 | |
ux | 必须 | 矢量 | 来自于farfieldux的ux数据。注意这个结果是个矢量,因此只在一个频率点下执行farfieldux脚本命令就足够了。 | |
uy | 必须 | 矢量 | 来自于farfielduy的uy数据。注意这个结果是个矢量,因此只在一个频率点下执行farfielduy脚本命令就足够了。 | |
theta | 必须 | 矢量 | 立体角θ矢量,以度为单位。长度必须为M或者1。 | |
phi | 必须 | 矢量 | 方位角φ矢量,以度为单位。长度必须为N或者1。 |
(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;