6 源命令和输出命令

6.1 waveform命令

        功能：在模型中指定用于源的波形。

        命令格式为：#waveform: str1 f1 f2 str2

str1：为波形的类型，它可以是

        gaussian：是一个高斯波形。

        gaussiandot：高斯波形的一阶导数。

        gaussiandotnorm：高斯波形的一阶导数的归一化。

        gaussiandotdot：高斯波形的二阶导数。

        gaussiandotdotnorm：高斯波形的归一化二阶导数。

        ricker：高斯波形的负的、归一化的二阶导数。

        sine：正弦波形的一个单周期。

        contsine ：一个连续的正弦波形。为了避免在计算中引入噪声，对正弦波的第一个周期(斜坡

激励)调制波形的幅度。

f1：为波形最大幅度的缩放(对于#赫兹偶极子，单位为安培，对于#电压源或#传输线，单位为伏

特)。

f2：为波形的频率(赫兹)。

str2：用于将波形分配给源的标识符

        例如，要指定一个振幅为1，中心频率为1.2 GHz的高斯波形。则指令为：

        #waveform: gaussian 1 1.2e9 my_gauss_pulse.

        注：用于创建波形的函数可以在工具部分找到

6.2 excitation_file命令

        功能：指定一个ASCII文件，该文件包含幅度值列，这些幅度值列指定可以与模型中的源一起

使用的自定义波形形状。每列的第一行必须以标识符字符串开头，该字符串将用作每个波形的名

称。如果振幅值少于将要执行的迭代次数，则振幅值序列的末尾将用零值填充，直至迭代次数。如

果指定的幅度值超出需要，则忽略。

        命令格式为：#excitation_file: file1

file1：为与输入文件在同一目录下包含指定波形的文件的名称，或者file1可以是包含指定波形的文

件的完整路径(允许您指定任何位置)。

        例如，要指定包含两个自定义波形形状的文件my_waves.txt。则指令为：

        #excitation_file: my_waves.txt

my_waves.txt文件的内容将采用如下形式:

然后使用my_pulse1自定义波形形状。

        例如，一个z极化赫兹偶极子源。则指令为：

        #hertzian_dipole: z 0.5 0.5 0.5 my_pulse1

       注：使用自定义波形激励时，不需要使用#waveform命令，只需要使用#excitation_file命令和

任何将与自定义波形激励一起使用的源。

6.3 hertzian_dipole命令

        功能：在电场位置指定电流密度项最简单的激励，通常称为加性或软源。

        式中𝐽𝑠为电流密度，𝐼为电流，∆𝑙为无穷小电偶极子的长度，∆x，∆y，∆z为网格的空间分辨

率。根据指定的极化，∆𝑙被设置为∆x，∆y或∆z

        命令格式为：#hertzian_dipole: c1 f1 f2 f3 str1 [f4 f5]

c1：为光源的偏振，可以是x、y或z。

f1、f2、f3：为模型中源的坐标(x、y、z)。

f4、f5：为可选参数。f4是启动源的时间延迟。f5是去除源的时间。

        如果时间窗口长于源去除时间，则源将在源去除时间后停止。如果源移除时间长于时间窗

口，则源将rx间窗口内处于活动状态。如果省略f4 f5，则源将在时间窗口开始时启动，在时间窗口

结束时停止。

str1：为与源一起使用的波形的标识符。

        例如，使用单位振幅的x极化赫兹偶极子和600 MHz中心频率的里克波形。则指令为：

        #waveform: ricker 1 600e6 my_ricker_pulse

        #hertzian_dipole: x 0.05 0.05 0.05 my_ricker_pulse

        当#hertzian_dipole在2d模拟中使用时，它在模拟的无限方向上充当电流的线源。

6.4 magnetic_dipole命令

        功能：模拟一个无穷小的磁偶极子。这通常被称为添加剂或软源。

        命令格式为：#magnetic_dipole: c1 f1 f2 f3 str1 [f4 f5]

c1：为光源的偏振，可以是x、y或z。

f1、f2、f3：为模型中源的坐标(x、y、z)。

f4、f5：为可选参数。f4是启动源的时间延迟。f5是去除源的时间。

        如果时间窗口长于源去除时间，则源将在源去除时间后停止。如果源移除时间长于时间窗

口，则源将rx间窗口内处于活动状态。如果省略f4 f5，则源将在时间窗口开始时启动，在时间窗口

结束时停止。

str1：为与源一起使用的波形的标识符。

6.5 voltage_source命令

        功能：在电场位置引入电压源。如果它的电阻为零，即指定电场分量的时间变化是规定的，

它可以是硬源，或者如果它的电阻不为零，它可以作为电阻电压源。当模型中包含天线的物理特性

时，对天线的激励是有用的。

        命令格式为：#voltage_source: c1 f1 f2 f3 f4 str1 [f5 f6]

c1：为光源的偏振，可以是x、y或z。

f1、f2、f3：为模型中源的坐标(x、y、z)。

f4、f5：为可选参数。f4是启动源的时间延迟。f5是去除源的时间。

        如果时间窗口长于源去除时间，则源将在源去除时间后停止。如果源移除时间长于时间窗

口，则源将rx间窗口内处于活动状态。如果省略f4 f5，则源将在时间窗口开始时启动，在时间窗口

结束时停止。

str1：为与源一起使用的波形的标识符。

        例如，一个y方向，内阻为50欧姆，幅度为5，中心频率为1.2 GHz高斯波形的电压源。则指令为：

        #waveform: gaussian 5 1.2e9my_gauss_pulse

        #voltage_source: y 0.05 0.05 0.05 50 my_gauss_pulse.

6.6 transmission_line:命令

        功能：在电场位置引入一维传输线模型。传输线的指定电阻可以大于零，小于自由空间的阻

抗(376.73欧姆)。当模型中包含天线的物理特性时，对天线的激励是有用的。

        命令格式为：#transmission_line: c1 f1 f2 f3 f4 str1 [f5 f6]

c1：为传输线的极化，可以是x、y或z。

f1、f2、f3：为模型中传输线的坐标(x、y、z)

f4为：传输线源特性电阻，单位为欧姆。可以是大于零且小于自由空间阻抗(376.73欧姆)的值。

f5、f6：为可选参数。f5是启动传输线励磁的时间延迟。f6为消除传输线励磁的时间。如果时间窗

长于输电线拆除时间的激励，则在输电线拆除时间的激励之后，输电线的激励将停止。如果传输线

的激励解除时间长于时间窗，则传输线的激励将在整个时间窗内处于活动状态。如果省略f5 f6，则

传输线的励磁将在时间窗开始时开始，在时间窗结束时停止。

str1：为与源一起使用的波形的标识符。

        传输线中电压和电流值的时间历史记录保存到输出文件中。这些都记录在输出文件部分中。

这些参数对于计算天线的特性有用。

        例如，输入阻抗或s参数。gprMax包含一个Python模块(在工具包中)来帮助您查看使用传输线

馈送的天线模型的输入阻抗和s11参数。

        例如，一个z向传输线源，其电阻为75欧姆，幅度为5，中心频率为1.2 GHz高斯波形。则指令

为：

        #waveform: gaussian 5 1.2e9 my_gauss_pulse

        #transmission_line: z 0.05 0.05 0.05 75 my_gauss_pulse

6.7 rx命令

        功能：在模型中引入输出点。在这些位置上，电场和磁场分量的值在模型迭代次数上将被保

存到文件中。

        命令格式为：#rx: f1 f2 f3 [str1 str2]

f1、f2、f3：为模型中接收机的坐标(x,y,z)。

str1：为接收方的标识符。

str2：为接收器的输出列表。它可以是Ex、Ey、Ez、Hx、Hy、Hz、Ix、Iy或Iz中的任何选项。

        注：当没有可选参数str1和str2时，所有电场和磁场分量将随接收点输出。

6.8 rx_array命令

        功能：在模型中定义多个输出点。

        命令格式为：#rx_array: f1 f2 f3 f4 f5 f6 f7 f8 f9

f1、f2、f3：为输出线/面/体积的左下(x,y,z)坐标，

f4、f5、f6：为输出线/面/体积的右上(x,y,z)坐标。

f7、f8、f9：为(x,y,z)增量，它定义了每个方向上输出点的数量。

        f7、f8、f9可以设置为零，以防止在特定方向上出现任何输出点。否则，f7的最小值为∆x，f8

的最小值为∆y，f9的最小值为∆z。

6.9 src_steps和rx_steps命令

        功能：在模型运行之间移动所有简单源(#src_steps)或所有接收器(#rx_steps)的位置。

        命令的语法为：#src_steps: f1 f2 f3

        命令的语法为：#rx_steps: f1 f2 f3

f1、f2、f3：为增量(x,y,z)来移动所有简单源(#hertzian_dipole或#magnetic_dipole)或所有接收器

(使用#rx或#rx_array命令创建)。

6.10 geometry_view命令

        功能：输出关于模型几何的文件信息。

        命令格式为：#geometry_view: f1 f2 f3 f4 f5 f6 f7 f8 f9 file1 c1

f1、f2、f3：为几何视图体积的左下角(x、y、z)坐标，以米为单位。

f4、f5、f6：为几何视图体积的右上角(x,y,z)坐标，以米为单位。

f7、f8、f9：为几何视图的空间离散，以米为单位。

        通常这些参数与模型的空间离散化相同，也可以改得更小。

file1：为文件的文件名，几何视图将存储在与输入文件相同的目录中。

c1：n (normal)或f (fine)，它指定是以每个单元为基础(n)还是以每个单元边缘为基)输出几何信息。

fine模式应该保留用于查看占用小体积的几何结构的详细部分，因为使用此模式可以生成具有大文

件大小的几何文件。

        注：只检查模型的几何形状时，可以使用可选的命令行参数--geometry-only运行gprMax。构

建模型并生成任何几何视图文件，但不会运行模拟。

6.11 snapshot命令

        功能：在给定的时间瞬间获得关于模型体积内电磁场的信息。

        命令格式为：#snapshot: f1 f2 f3 f4 f5 f6 f7 f8 f9 f10 file1

                             #snapshot: f1 f2 f3 f4 f5 f6 f7 f8 f9 i1 file1

f1、f2、f3：为快照体积的左下角(x,y,z)坐标，单位为米。

f4、f5、f6：为快照体积的右上角(x,y,z)坐标，单位为米。

f7、f8、f9：为快照的空间离散，单位为米。

f10、i1：为以秒为单位的时间(浮点数)或迭代数(整数)，表示快照将被拍摄的时间点。

file1：为要保存快照的文件名。快照文件会自动存储在输入文件名后面加上' _snaps '的目录中。对于多个模型运行，每个模型运行将有自己的目录，即' _snaps1 '， ' snaps2 '等…

        例如，在模拟时间为3纳秒时保存模型中电磁场的快照。则指令为：

        #snapshot: 0 0 0 1 1 1 0.1 0.1 0.1 3e-9 snap1

        可以利用Python脚本创建一系列快照。

        例如，要以0.1ns的间隔从0.1ns开始到3ns创建30个快照，请在输入文件中使用以下代码片

段。相应替换x1 y1 z1 x2 y2 z2 dx dy dz

#python:

        for i in range(1, 31):

        print('#snapshot: x1 y1 z1 x2 y2 z2 dx dy dz {} snapshot{}'.format((i/10)*1e-9,i))

#end_python: