2024年GPRMA2D知识点汇总VIP专享VIP免费

GprMax是爱丁堡大学的AntonisGiannopoulos于1996年推出来的一种基于时域有限差分（FDTD）算法和理想匹配层（PML）边界吸取条件的探地雷达正演数值模拟软件，用于探地雷到达像研究。其中，GprMax2D是二维正演，GprMax3D为三维正演。该软件可以在Windows、Linux和MacOS三个平台上使用。本文重要针对Windows平台进行阐明。一、软件获得该软件为免费软件，可以去GprMax官网下载。也可点此直接下载。二、准备工作软件无需安装，下载后用解压工具解压，找到Windows文献夹，直接双击GprMax2D.exe即可运行。不过为了以便批量模拟，提议把Windows文献夹下的文献（cygwin1.dll、GprMax2D.exe和GprMax3D.exe）所有复制到系统盘系统搜索途径下，例如：C:\Windows。假如只进行二维正演，只复制cygwin1.dll和GprMax2D.exe即可。当然也可放在任意途径下，只要使用时包括所在途径就OK了。笔者习惯放在系统要目录下，即C:\。此外，最佳建立一种输入文献和模拟成果寄存的专用文献夹，并且把tools文献夹下的文献（gprmax.m、gprmax2g.m、gprmax3g.m、gprmaxde.m和gprmaxso.m）所有复制到该专用文献夹下，注意文献夹名最佳使用英文。假如只进行二维正演，只复制gprmax.m和gprmax2g.m即可。笔者习惯使用D:\GPR。正演成果需要用MATLAB进行绘图，因此需要安装有MATLAB软件。软件下载地址和安装措施此处不再阐明，后来使用将假设读者已经成功安装MATLAB软件。三、软件使用一般进行数值模拟时一般都是若干个对比模型进行模拟，因此本文只简介批量模拟的措施。如此一来，笔者将认为读者的软件和笔者同样放在C:\，输入文献和gprmax.m、gprmax2g.m文献已经放在D:\GPR。欲进行正演模拟需要先建立输入文献，输入文献的建立将在背面简介，此处假设在D:\GPR已经有若干输入文献，文献名分别为fname1.in，fname2.in，fname3.in。打开记事本，写入如下两行直线间的内容，另存为*.bat文献，例如：GPR.bat。其中括号里的内容为本条语句的阐明。------------------------------------------------------------------------------------------pathC:\;%path%(把软件所在途径包括进系统搜索途径)cd/dD:\GPR(设置专用文献夹途径为活动途径)gprmax2dD:\GPR\fname1.in（按次序进行批量正演模拟）gprmax2dD:\GPR\fname2.ingprmax2dD:\GPR\fname3.in------------------------------------------------------------------------------------------完毕后保留并关闭，然后双击运行*.bat文献即可进行批量正演模拟四、建立输入文献输入文献是纯文本文献，可以用任何文字处理程序编辑，并且内容所有为英文。打开记事本，写入输入文献内容，另存为*.in即可。输入文献必需包括模拟模型的所有必要信息，每一种有效命令行都必须以符号（#）开始，否则该命令行无效。所有命令的一般语法为：#命令名:参数1参数2参数3……为了以便简介GprMax2D的命令及参数，做了如下约定：·f代表实数，如1.5·i代表整数，如15·c代表单个字母，如y·str代表字条串，如air·file代表文献名，如test.in·所有的长度单位均为m·所有的时间单位均为s·所有的频率单位均为Hz·模型空间的坐标原点（0,0）在左下角输入文献中除了#analysis:和#end_analysis:命令外，先后次序一般是没有规定的。GprMax2D2.0版一共有32个命令：------------------------------------------------------------------------------------------#title:str输出文献的标题，包括在输出文献中。#domain:f1f2指定模型空间大小，单位是米。·f1x方向的大小·f2y方向的大小#dx_dy:f1f2指定网格步长，单位是米。·f1x方向的步长，即Δx·f2y方向的步长，即Δy模型划分的网格数目是由#domain:和#dx_dy:两个命令共同决定的。网格划分越小模型越精确，但模拟的计算量也越在，需要权衡确定。#time_step_stability_factor:f1时间步长的稳定系数，取值范围（0，1］，实际使用时间为f1×Δt。#time_window:f1ori1时窗大小，单位是秒。·f1为时窗，如f1=16e-9，时窗为16ns·i1为时间步长数，时窗为i1×Δt。#messages:c1与否启动屏幕信息。y(es)启动n(o)关闭，默认n。#number_of_media:i1介质数目，i1应...