FDTD Script命令学习-farfield3dintegrate/farfieldspherical

目录

一、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；